Interaction of Tumor with Its Micro-environment: A Mathematical Model

Mathematical models of intercellular signaling in breast cancer

BACKGROUND AND OBJECTIVES

The development and regulation of healthy and cancerous breast tissue is guided by communication between cells. Diverse signals are exchanged between cancer cells and non-cancerous cells of the tumor microenvironment (TME), influencing all stages of tumor progression. Mathematical models are essential for understanding how this complex network determines cancer progression and the effectiveness of treatment.

METHODOLOGY

We reviewed the current dynamical mathematical models of intercellular signaling in breast cancer, examining models with cancer cells only, fibroblasts, endothelial cells, macrophages and the immune system as whole. We categorized the goals and complexity of these models, to highlight how they can explain many features of cancer emergence and progression.

RESULTS

We found that dynamical models of intercellular signaling can elucidate tissue-level dysregulation in cancer by explaining: i) maintenance of non-heritable intratumor phenotypic heterogeneity, ii) transitions between tumor dormancy and accelerated invasive growth, iii) stromal support of tumor vascularization and growth factor enrichment and iv) suppression of immune infiltration and cancer surveillance. These models also provide a framework to propose novel TME-targeting treatment strategies. However, most models were focused on a highly selected and small set of signaling interactions between a few cell types, and their translational applicability were severely limited by the availability of tumor-specific data for personalized model calibration.

CONCLUSIONS AND IMPLICATIONS

Mathematical models of breast cancer have many challenges and opportunities to incorporate signaling. The four key challenges are: 1) finding ways to treat signaling networks as a context-dependent language that incorporates non-linear and non-additive responses, 2) identifying the key cell phenotypes that signals control and understanding the feedbacks between signals and phenotype that determine the progression of cancer, (3) estimating parameters of specific patient tumors early in treatment, 4) linking models with novel data collection methods that have single cell and spatial resolution. As our approaches advance, it is our hope that dynamical mathematical models of inter-cellular signaling can play a central role in identifying and testing new treatment strategies as well as forecasting impacts of disease treatment.

Investigating the spatial interaction of immune cells in colon cancer

The intricate network of interactions between cells and molecules in the tumor microenvironment creates a heterogeneous ecosystem. The proximity of the cells and molecules to their activators and inhibitors is essential in the progression of tumors. Here, we develop a system of partial differential equations coupled with linear elasticity to investigate the effects of spatial interactions on the tumor microenvironment. We observe interesting cell and cytokine distribution patterns, which are heavily affected by macrophages. We also see that cytotoxic T cells get recruited and suppressed at the site of macrophages. Moreover, we observe that anti-tumor macrophages reorganize the patterns in favor of a more spatially restricted cancer and necrotic core. Furthermore, the adjoint-based sensitivity analysis indicates that the most sensitive model's parameters are directly related to macrophages. The results emphasize the widely acknowledged effect of macrophages in controlling cancer cells population and spatially arranging cells in the tumor microenvironment.

Modeling tumour heterogeneity of PD-L1 expression in tumour progression and adaptive therapy

Although PD-1/PD-L1 inhibitors show potent and durable anti-tumour effects in some refractory tumours, the response rate in overall patients is unsatisfactory, which in part due to the inherent heterogeneity of PD-L1. In order to establish an approach for predicting and estimating the dynamic alternation of PD-L1 heterogeneity during cancer progression and treatment, this study establishes a comprehensive modelling and computational framework based on a mathematical model of cancer cell evolution in the tumour-immune microenvironment, and in combination with epigenetic data and overall survival data of clinical patients from The Cancer Genome Atlas. Through PD-L1 heterogeneous virtual patients obtained by the computational framework, we explore the adaptive therapy of administering anti-PD-L1 according to the dynamic of PD-L1 state among cancer cells. Our results show that in contrast to the continuous maximum tolerated dose treatment, adaptive therapy is more effective for PD-L1 positive patients, in that it prolongs the survival of patients by administration of drugs at lower dosage.

Role of senescent tumor cells in building a cytokine shield in the tumor microenvironment: mathematical modeling

Cellular senescence can induce dual effects (promotion or inhibition) on cancer progression. While immune cells naturally respond and migrate toward various chemotactic sources from the tumor mass, various factors including senescent tumor cells (STCs) in the tumor microenvironment may affect this chemotactic movement. In this work, we investigate the mutual interactions between the tumor cells and the immune cells that either inhibit or facilitate tumor growth by developing a mathematical model that consists of taxis-reaction-diffusion equations and receptor kinetics for the key players in the interaction network. We apply a mathematical model to a transwell Boyden chamber invasion assay used in the experiments to illustrate that STCs can play a pivotal role in negating immune attack through tight regulation of intra- and extra-cellular signaling molecules. In particular, we show that senescent tumor cells in cell cycle arrest can block intratumoral infiltration of CD8+ T cells by secreting a high level of CXCL12, which leads to significant reduction its receptors, CXCR4, on T cells, and thus impaired chemotaxis. The predictions of nonlinear responses to CXCL12 were in good agreement with experimental data. We tested several hypotheses on immune-tumor interactions under various biochemical conditions in the tumor microenvironment and developed new concepts for anti-tumor strategies targeting senescence induced immune impairment.

Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment

The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor's growth.

Multiscale modeling of collective cell migration elucidates the mechanism underlying tumor-stromal interactions in different spatiotemporal scales

Metastasis is the pathogenic spread of cancer cells from a primary tumor to a secondary site which happens at the late stages of cancer. It is caused by a variety of biological, chemical, and physical processes, such as molecular interactions, intercellular communications, and tissue-level activities. Complex interactions of cancer cells with their microenvironment components such as cancer associated fibroblasts (CAFs) and extracellular matrix (ECM) cause them to adopt an invasive phenotype that promotes tumor growth and migration. This paper presents a multiscale model for integrating a wide range of time and space interactions at the molecular, cellular, and tissue levels in a three-dimensional domain. The modeling procedure starts with presenting nonlinear dynamics of cancer cells and CAFs using ordinary differential equations based on TGFβ, CXCL12, and LIF signaling pathways. Unknown kinetic parameters in these models are estimated using hybrid unscented Kalman filter and the models are validated using experimental data. Then, the principal role of CAFs on metastasis is revealed by spatial-temporal modeling of circulating signals throughout the TME. At this stage, the model has evolved into a coupled ODE-PDE system that is capable of determining cancer cells' status in one of the quiescent, proliferating or migratory conditions due to certain metastasis factors and ECM characteristics. At the tissue level, we consider a force-based framework to model the cancer cell proliferation and migration as the final step towards cancer cell metastasis. The ability of the multiscale model to depict cancer cells' behavior in different levels of modeling is confirmed by comparing its outputs with the results of RT PCR and wound scratch assay techniques. Performance evaluation of the model indicates that the proposed multiscale model can pave the way for improving the efficiency of therapeutic methods in metastasis prevention.

Macrophage polarization as a potential therapeutic target for atherosclerosis: a dynamic stochastic modelling study

We proposed a dynamic stochastic mathematical model to evaluate the role of macrophage polarization in plaque development. The dynamic process of macrophages from proliferation to death was simulated under different lipid microenvironments. The probability of macrophage phenotypic switching was described using a Bernoulli distribution where the stochastic variable was determined by the local lipid level. Moreover, the interactions between macrophages and microenvironmental factors vary with macrophage phenotype. We investigated the distribution of key microenvironmental factors, the dynamics of macrophage polarization and its influence on foam cell formation. M1 macrophages were found to predominate in advanced plaque corresponding to the exacerbated inflammation observed in mice experiments. The imbalance between the deposition of oxidized low-density lipoprotein and phagocytic effects of macrophages governed the formation of foam cells. Furthermore, we simulated targeted therapies by either directly inhibiting the polarization probability to M1 macrophages or indirectly regulating macrophage polarization due to high-density lipoprotein levels. Comparison of simulation results with experimental findings in both therapies indicated that the intervention and regulation of macrophage polarization could influence plaque microenvironment and subsequently induce plaque regression, especially in the early stage. The proposed modelling system can facilitate the evaluation of novel therapies targeting macrophage polarization.

Prediction of the atherosclerotic plaque development in carotid arteries; the effect of T-cells

The carotid artery disease is one of the leading causes of mortality worldwide, as it leads to the progressive arterial stenosis that may result to stroke. To address this issue, the scientific community is attempting not only to enrich our knowledge on the underlying atherosclerotic mechanisms, but also to enable the prediction of the atherosclerotic progression. This study investigates the role of T-cells in the atherosclerotic plaque growth process through the implementation of a computational model in realistic geometries of carotid arteries. T-cells mediate in the inflammatory process by secreting interferon-y that enhances the activation of macrophages. In this analysis, we used 5 realistic human carotid arterial segments as input to the model. In particular, magnetic resonance imaging data, as well as, clinical data were collected from the patients at two time points. Using the baseline data, plaque growth was predicted and correlated to the follow-up arterial geometries. The results exhibited a very good agreement between them, presenting a high coefficient of determination R²=0.64.

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells' location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells' diffusion and immune cells' influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes' diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

A novel nonlocal partial differential equation model of endothelial progenitor cell cluster formation during the early stages of vasculogenesis

The formation of new vascular networks is essential for tissue development and regeneration, in addition to playing a key role in pathological settings such as ischemia and tumour development. Experimental findings in the past two decades have led to the identification of a new mechanism of neovascularisation, known as cluster-based vasculogenesis, during which endothelial progenitor cells (EPCs) mobilised from the bone marrow are capable of bridging distant vascular beds in a variety of hypoxic settings in vivo. This process is characterised by the formation of EPC clusters during its early stages and, while much progress has been made in identifying various mechanisms underlying cluster formation, we are still far from a comprehensive description of such spatio-temporal dynamics. In order to achieve this, we propose a novel mathematical model of the early stages of cluster-based vasculogenesis, comprising of a system of nonlocal partial differential equations including key mechanisms such as endogenous chemotaxis, matrix degradation, cell proliferation and cell-to-cell adhesion. We conduct a linear stability analysis on the system and solve the equations numerically. We then conduct a parametric analysis of the numerical solutions of the one-dimensional problem to investigate the role of underlying dynamics on the speed of cluster formation and the size of clusters, measured via appropriate metrics for the cluster width and compactness. We verify the key results of the parametric analysis with simulations of the two-dimensional problem. Our results, which qualitatively compare with data from in vitro experiments, elucidate the complementary role played by endogenous chemotaxis and matrix degradation in the formation of clusters, suggesting chemotaxis is responsible for the cluster topology while matrix degradation is responsible for the speed of cluster formation. Our results also indicate that the nonlocal cell-to-cell adhesion term in our model, even though it initially causes cells to aggregate, is not sufficient to ensure clusters are stable over long time periods. Consequently, new modelling strategies for cell-to-cell adhesion are required to stabilise in silico clusters. We end the paper with a thorough discussion of promising, fruitful future modelling and experimental research perspectives.

Optimal regulation of tumour-associated neutrophils in cancer progression

In a tumour microenvironment, tumour-associated neutrophils could display two opposing differential phenotypes: anti-tumour (N1) and pro-tumour (N2) effector cells. Converting N2 to N1 neutrophils provides innovative therapies for cancer treatment. In this study, a mathematical model for N1-N2 dynamics describing the cancer survival and immune inhibition in response to TGF-β and IFN-β is considered. The effects of exogenous intervention of TGF-β inhibitor and IFN-β are examined in order to enhance N1 recruitment to combat tumour progression. Our approach employs optimal control theory to determine drug infusion protocols that could minimize tumour volume with least administration cost possible. Four optimal control scenarios corresponding to different therapeutic strategies are explored, namely, TGF-β inhibitor control only, IFN-β control only, concomitant TGF-β inhibitor and IFN-β controls, and alternating TGF-β inhibitor and IFN-β controls. For each scheme, different initial conditions are varied to depict different pathophysiological condition of a cancer patient, leading to adaptive treatment schedule. TGF-β inhibitor and IFN-β drug dosages, total drug amount, infusion times and relative cost of drug administrations are obtained under various circumstances. The control strategies achieved could guide in designing individualized therapeutic protocols.

Quantitative predictive approaches for Dupuytren disease: a brief review and future perspectives

In this study we review the current state of the art for Dupuytren's disease (DD), while emphasising the need for a better integration of clinical, experimental and quantitative predictive approaches to understand the evolution of the disease and improve current treatments. We start with a brief review of the biology of this disease and current treatment approaches. Then, since certain aspects in the pathogenesis of this disorder have been compared to various biological aspects of wound healing and malignant processes, next we review some in silico (mathematical modelling and simulations) predictive approaches for complex multi-scale biological interactions occurring in wound healing and cancer. We also review the very few in silico approaches for DD, and emphasise the applicability of these approaches to address more biological questions related to this disease. We conclude by proposing new mathematical modelling and computational approaches for DD, which could be used in the absence of animal models to make qualitative and quantitative predictions about the evolution of this disease that could be further tested in vitro.

Atorvastatin-mediated rescue of cancer-related cognitive changes in combined anticancer therapies

Acute administration of trastuzumab (TZB) may induce various forms of cognitive impairment. These cancer-related cognitive changes (CRCC) are regulated by an adverse biological process involving cancer stem cells (CSCs) and IL-6. Recent studies have reported that atorvastatin (ATV) may change the dynamic of cognitive impairment in a combination (TZB+ATV) therapy. In this study, we investigate the mutual interactions between cancer stem cells and the tumor cells that facilitate cognitive impairment during long term TZB therapy by developing a mathematical model that involves IL-6 and the key apoptotic regulation. These include the densities of tumor cells and CSCs, and the concentrations of intracellular signaling molecules (NFκB, Bcl-2, BAX). We apply the mathematical model to a single or combination (ATV+TZB) therapy used in the experiments to demonstrate that the CSCs can enhance CRCC by secreting IL-6 and ATV may interfere the whole regulation. We show that the model can both reproduce the major experimental observation on onset and prevention of CRCC, and suggest several important predictions to guide future experiments with the goal of the development of new anti-tumor and anti-CRCC strategies. Moreover, using this model, we investigate the fundamental mechanism of onset of cognitive impairment in TZB-treated patients and the impact of alternating therapies on the anti-tumor efficacy and intracellular response to different treatment schedules.

A conventional drug, trastuzumab (TZB), was shown to be an effective weapon in killing cancer cells in brain. However, long term treatment of TZB increases the proportion of cancer stem cells (CSCs) in the tumour microenvironment (TME) and induces up-regulation of pro-tumoral molecules such as IL-6 in TME. These cancer cells then become more resistant to this chemotherapy through the IL-mediated up-regulation of NFκB and CSCs. More importantly, these changes in TME result in a serious side effect, cognitive impairment called cancer-related cognitive changes (CRCC). The detailed mechanism of CRCC is still poorly understood. However, cancer patients with chemotherapy-induced cognitive impairment can have long-term or delayed mental changes. In this study, we investigated the fundamental mechanism of CRCC in cancer patients based on experiments and a mathematical model that describes how tumor cells interact with CSCs in response to chemo drugs. In particular, we investigate how TZB-induced CSCs with modified IL-6 landscapes shape the cognitive functions in cancer patients. We showed that the combination treatment with another drug, atorvastatin (ATV), can abrogate the TZB-induced CRCC and enhance the survival probability of cancer patients by synergistic anti-tumor effect. We demonstrate that the cognitive functions and survival rates in cancer patients depend on the apoptotic signaling pathways via the critical communication and IL-6 landscapes of stimulated CTCs.

A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration

Breast cancer is the most prominent type of cancer among women. Understanding the microenvironment of breast cancer and the interactions between cells and cytokines will lead to better treatment approaches for patients. In this study, we developed a data-driven mathematical model to investigate the dynamics of key cells and cytokines involved in breast cancer development. We used gene expression profiles of tumors to estimate the relative abundance of each immune cell and group patients based on their immune patterns. Dynamical results show the complex interplay between cells and molecules, and sensitivity analysis emphasizes the direct effects of macrophages and adipocytes on cancer cell growth. In addition, we observed the dual effect of IFN-γ on cancer proliferation, either through direct inhibition of cancer cells or by increasing the cytotoxicity of CD8+ T-cells.

Mathematical model of STAT signalling pathways in cancer development and optimal control approaches

In various diseases, the STAT family display various cellular controls over various challenges faced by the immune system and cell death programs. In this study, we investigate how an intracellular signalling network (STAT1, STAT3, Bcl-2 and BAX) regulates important cellular states, either anti-apoptosis or apoptosis of cancer cells. We adapt a mathematical framework to illustrate how the signalling network can generate a bi-stability condition so that it will induce either apoptosis or anti-apoptosis status of tumour cells. Then, we use this model to develop several anti-tumour strategies including IFN-β infusion. The roles of JAK-STATs signalling in regulation of the cell death program in cancer cells and tumour growth are poorly understood. The mathematical model unveils the structure and functions of the intracellular signalling and cellular outcomes of the anti-tumour drugs in the presence of IFN-β and JAK stimuli. We identify the best injection order of IFN-β and DDP among many possible combinations, which may suggest better infusion strategies of multiple anti-cancer agents at clinics. We finally use an optimal control theory in order to maximize anti-tumour efficacy and minimize administrative costs. In particular, we minimize tumour volume and maximize the apoptotic potential by minimizing the Bcl-2 concentration and maximizing the BAX level while minimizing total injection amount of both IFN-β and JAK2 inhibitors (DDP).

Data Driven Mathematical Model of FOLFIRI Treatment for Colon Cancer

Many colon cancer patients show resistance to their treatments. Therefore, it is important to consider unique characteristic of each tumor to find the best treatment options for each patient. In this study, we develop a data driven mathematical model for interaction between the tumor microenvironment and FOLFIRI drug agents in colon cancer. Patients are divided into five distinct clusters based on their estimated immune cell fractions obtained from their primary tumors' gene expression data. We then analyze the effects of drugs on cancer cells and immune cells in each group, and we observe different responses to the FOLFIRI drugs between patients in different immune groups. For instance, patients in cluster 3 with the highest T-reg/T-helper ratio respond better to the FOLFIRI treatment, while patients in cluster 2 with the lowest T-reg/T-helper ratio resist the treatment. Moreover, we use ROC curve to validate the model using the tumor status of the patients at their follow up, and the model predicts well for the earlier follow up days.

Role of neutrophil extracellular traps in regulation of lung cancer invasion and metastasis: Structural insights from a computational model

Lung cancer is one of the leading causes of cancer-related deaths worldwide and is characterized by hijacking immune system for active growth and aggressive metastasis. Neutrophils, which in their original form should establish immune activities to the tumor as a first line of defense, are undermined by tumor cells to promote tumor invasion in several ways. In this study, we investigate the mutual interactions between the tumor cells and the neutrophils that facilitate tumor invasion by developing a mathematical model that involves taxis-reaction-diffusion equations for the critical components in the interaction. These include the densities of tumor and neutrophils, and the concentrations of signaling molecules and structure such as neutrophil extracellular traps (NETs). We apply the mathematical model to a Boyden invasion assay used in the experiments to demonstrate that the tumor-associated neutrophils can enhance tumor cell invasion by secreting the neutrophil elastase. We show that the model can both reproduce the major experimental observation on NET-mediated cancer invasion and make several important predictions to guide future experiments with the goal of the development of new anti-tumor strategies. Moreover, using this model, we investigate the fundamental mechanism of NET-mediated invasion of cancer cells and the impact of internal and external heterogeneity on the migration patterning of tumour cells and their response to different treatment schedules.

Data Driven Mathematical Model of Colon Cancer Progression

Every colon cancer has its own unique characteristics, and therefore may respond differently to identical treatments. Here, we develop a data driven mathematical model for the interaction network of key components of immune microenvironment in colon cancer. We estimate the relative abundance of each immune cell from gene expression profiles of tumors, and group patients based on their immune patterns. Then we compare the tumor sensitivity and progression in each of these groups of patients, and observe differences in the patterns of tumor growth between the groups. For instance, in tumors with a smaller density of naive macrophages than activated macrophages, a higher activation rate of macrophages leads to an increase in cancer cell density, demonstrating a negative effect of macrophages. Other tumors however, exhibit an opposite trend, showing a positive effect of macrophages in controlling tumor size. Although the results indicate that for all patients the size of the tumor is sensitive to the parameters related to macrophages, such as their activation and death rate, this research demonstrates that no single biomarker could predict the dynamics of tumors.

Synthesis and secretome release by human bone marrow mesenchymal stem cell spheroids within three-dimensional collagen hydrogels: Integrating experiments and modelling

Myocardial infarction results in loss of cardiac cell types, inflammation, extracellular matrix (ECM) degradation, and fibrotic scar. Transplantation of bone marrow-derived mesenchymal stem cells (BM-MSCs) is being explored as they could differentiate into cardiomyocyte-like cells, integrate into host tissue, and enhance resident cell activity. The ability of these cells to restore lost ECM, remodel the inflammatory scar tissue, and repair the injured myocardium remains unexplored. We here elucidated the synthesis and deposition of ECM (e.g., elastin, sulfated glycosaminoglycans, hyaluronan, collagen type III, laminin, fibrillin, lysyl oxidase, and nitric oxide synthases), matrix metalloproteinases (MMPs) and their inhibitors (TIMPs), and other secretome (cytokines, chemokines, and growth factors) in adult human BM-MSC spheroid cultures within three-dimensional collagen gels. The roles of species-specific type I collagen and 5-azacytadine were assessed over a 28-day period. Results revealed that human collagen (but not rat-derived) suppressed MSC proliferation and survival, and MSCs synthesized and released a variety of ECM proteins and secretome over the 28 days. Matrix deposition is at least an order of magnitude lower than their release levels at every time point, most possibly due to elevated MMP levels and interleukins with a concomitant decrease in TIMPs. Matrix synthesis over the 28-day period was fitted to a competitive inhibition form of Michaelis-Menten kinetics, and the production and decay rates of ECM, MMPs, and TIMPs, along with the kinetic model parameters quantified. Such an integrated experimental and modelling approach would help elucidate the critical roles of various parameters (e.g., cell encapsulation and delivery vehicles) in stem cell-based transplantation therapies.

Synergistic Effects of Bortezomib-OV Therapy and Anti-Invasive Strategies in Glioblastoma: A Mathematical Model

It is well-known that the tumor microenvironment (TME) plays an important role in the regulation of tumor growth and the efficacy of anti-tumor therapies. Recent studies have demonstrated the potential of combination therapies, using oncolytic viruses (OVs) in conjunction with proteosome inhibitors for the treatment of glioblastoma, but the role of the TME in such therapies has not been studied. In this paper, we develop a mathematical model for combination therapies based on the proteosome inhibitor bortezomib and the oncolytic herpes simplex virus (oHSV), with the goal of understanding their roles in bortezomib-induced endoplasmic reticulum (ER) stress, and how the balance between apoptosis and necroptosis is affected by the treatment protocol. We show that the TME plays a significant role in anti-tumor efficacy in OV combination therapy, and illustrate the effect of different spatial patterns of OV injection. The results illustrate a possible phenotypic switch within tumor populations in a given microenvironment, and suggest new anti-invasion therapies.

Role of tumor-associated neutrophils in regulation of tumor growth in lung cancer development: A mathematical model

Neutrophils display rapid and potent innate immune responses in various diseases. Tumor-associated neutrophils (TANs) however either induce or overcome immunosuppressive functions of the tumor microenvironment through complex tumor-stroma crosstalk. We developed a mathematical model to address the question of how phenotypic alterations between tumor suppressive N1 TANS, and tumor promoting N2 TANs affect nonlinear tumor growth in a complex tumor microenvironment. The model provides a visual display of the complex behavior of populations of TANs and tumors in response to various TGF-β and IFN-β stimuli. In addition, the effect of anti-tumor drug administration is incorporated in the model in an effort to achieve optimal anti-tumor efficacy. The simulation results from the mathematical model were in good agreement with experimental data. We found that the N2-to-N1 ratio (N21R) index is positively correlated with aggressive tumor growth, suggesting that this may be a good prognostic factor. We also found that the antitumor efficacy increases when the relative ratio (Dap) of delayed apoptotic cell death of N1 and N2 TANs is either very small or relatively large, providing a basis for therapeutically targeting prometastatic N2 TANs.

Modeling three-dimensional invasive solid tumor growth in heterogeneous microenvironment under chemotherapy

A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations indicate that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.

Modeling triple-negative breast cancer heterogeneity: Effects of stromal macrophages, fibroblasts and tumor vasculature

A hallmark of breast tumors is its spatial heterogeneity that includes its distribution of cancer stem cells and progenitor cells, but also heterogeneity in the tumor microenvironment. In this study we focus on the contributions of stromal cells, specifically macrophages, fibroblasts, and endothelial cells on tumor progression. We develop a computational model of triple-negative breast cancer based on our previous work and expand it to include macrophage infiltration, fibroblasts, and angiogenesis. In vitro studies have shown that the secretomes of tumor-educated macrophages and fibroblasts increase both the migration and proliferation rates of triple-negative breast cancer cells. In vivo studies also demonstrated that blocking signaling of selected secreted factors inhibits tumor growth and metastasis in mouse xenograft models. We investigate the influences of increased migration and proliferation rates on tumor growth, the effect of the presence on fibroblasts or macrophages on growth and morphology, and the contributions of macrophage infiltration on tumor growth. We find that while the presence of macrophages increases overall tumor growth, the increase in macrophage infiltration does not substantially increase tumor growth and can even stifle tumor growth at excessive rates.

A mathematical model of aortic aneurysm formation

Abdominal aortic aneurysm (AAA) is a localized enlargement of the abdominal aorta, such that the diameter exceeds 3 cm. The natural history of AAA is progressive growth leading to rupture, an event that carries up to 90% risk of mortality. Hence there is a need to predict the growth of the diameter of the aorta based on the diameter of a patient’s aneurysm at initial screening and aided by non-invasive biomarkers. IL-6 is overexpressed in AAA and was suggested as a prognostic marker for the risk in AAA. The present paper develops a mathematical model which relates the growth of the abdominal aorta to the serum concentration of IL-6. Given the initial diameter of the aorta and the serum concentration of IL-6, the model predicts the growth of the diameter at subsequent times. Such a prediction can provide guidance to how closely the patient’s abdominal aorta should be monitored. The mathematical model is represented by a system of partial differential equations taking place in the aortic wall, where the media is assumed to have the constituency of an hyperelastic material.

The role of myosin II in glioma invasion: A mathematical model

Gliomas are malignant tumors that are commonly observed in primary brain cancer. Glioma cells migrate through a dense network of normal cells in microenvironment and spread long distances within brain. In this paper we present a two-dimensional multiscale model in which a glioma cell is surrounded by normal cells and its migration is controlled by cell-mechanical components in the microenvironment via the regulation of myosin II in response to chemoattractants. Our simulation results show that the myosin II plays a key role in the deformation of the cell nucleus as the glioma cell passes through the narrow intercellular space smaller than its nuclear diameter. We also demonstrate that the coordination of biochemical and mechanical components within the cell enables a glioma cell to take the mode of amoeboid migration. This study sheds lights on the understanding of glioma infiltration through the narrow intercellular spaces and may provide a potential approach for the development of anti-invasion strategies via the injection of chemoattractants for localization.

A Macroscopic Mathematical Model for Cell Migration Assays Using a Real-Time Cell Analysis

Experiments of cell migration and chemotaxis assays have been classically performed in the so-called Boyden Chambers. A recent technology, xCELLigence Real Time Cell Analysis, is now allowing to monitor the cell migration in real time. This technology measures impedance changes caused by the gradual increase of electrode surface occupation by cells during the course of time and provide a Cell Index which is proportional to cellular morphology, spreading, ruffling and adhesion quality as well as cell number. In this paper we propose a macroscopic mathematical model, based on advection-reaction-diffusion partial differential equations, describing the cell migration assay using the real-time technology. We carried out numerical simulations to compare simulated model dynamics with data of observed biological experiments on three different cell lines and in two experimental settings: absence of chemotactic signals (basal migration) and presence of a chemoattractant. Overall we conclude that our minimal mathematical model is able to describe the phenomenon in the real time scale and numerical results show a good agreement with the experimental evidences.

Dynamics and control of a mathematical model for metronomic chemotherapy

A 3-compartment model for metronomic chemotherapy that takes into account cancerous cells, the tumor vasculature and tumor immune-system interactions is considered as an optimal control problem. Metronomic chemo-therapy is the regular, almost continuous administration of chemotherapeutic agents at low dose, possibly with small interruptions to increase the efficacy of the drugs. There exists medical evidence that such administrations of specific cytotoxic agents (e.g., cyclophosphamide) have both antiangiogenic and immune stimulatory effects. A mathematical model for angiogenic signaling formulated by Hahnfeldt et al. is combined with the classical equations for tumor immune system interactions by Stepanova to form a minimally parameterized model to capture these effects of low dose chemotherapy. The model exhibits bistable behavior with the existence of both benign and malignant locally asymptotically stable equilibrium points. In this paper, the transfer of states from the malignant into the benign regions is used as a motivation for the construction of an objective functional that induces this process and the analysis of the corresponding optimal control problem is initiated.

Strategies of eradicating glioma cells: a multi-scale mathematical model with MiR-451-AMPK-mTOR control

The cellular dispersion and therapeutic control of glioblastoma, the most aggressive type of primary brain cancer, depends critically on the migration patterns after surgery and intracellular responses of the individual cancer cells in response to external biochemical and biomechanical cues in the microenvironment. Recent studies have shown that a particular microRNA, miR-451, regulates downstream molecules including AMPK and mTOR to determine the balance between rapid proliferation and invasion in response to metabolic stress in the harsh tumor microenvironment. Surgical removal of main tumor is inevitably followed by recurrence of the tumor due to inaccessibility of dispersed tumor cells in normal brain tissue. In order to address this multi-scale nature of glioblastoma proliferation and invasion and its response to conventional treatment, we propose a hybrid model of glioblastoma that analyses spatio-temporal dynamics at the cellular level, linking individual tumor cells with the macroscopic behaviour of cell organization and the microenvironment, and with the intracellular dynamics of miR-451-AMPK-mTOR signaling within a tumour cell. The model identifies a key mechanism underlying the molecular switches between proliferative phase and migratory phase in response to metabolic stress and biophysical interaction between cells in response to fluctuating glucose levels in the presence of blood vessels (BVs). The model predicts that cell migration, therefore efficacy of the treatment, not only depends on oxygen and glucose availability but also on the relative balance between random motility and strength of chemoattractants. Effective control of growing cells near BV sites in addition to relocalization of invisible migratory cells back to the resection site was suggested as a way of eradicating these migratory cells.

Mathematical modeling of Interleukin-35 promoting tumor growth and angiogenesis

Interleukin-35 (IL-35), a cytokine from the Interleukin-12 cytokine family, has been considered as an anti-inflammatory cytokine which promotes tumor progression and tumor immune evasion. It has also been demonstrated that IL-35 is secreted by regulatory T cells. Recent mouse experiments have shown that IL-35 produced by cancer cells promotes tumor growth via enhancing myeloid cell accumulation and angiogenesis, and reducing the infiltration of activated CD8[Formula: see text] T cells into tumor microenvironment. In the present paper we develop a mathematical model based on these experimental results. We include in the model an anti-IL-35 drug as treatment. The extended model (with drug) is used to design protocols of anti-IL-35 injections for treatment of cancer. We find that with a fixed total amount of drug, continuous injection has better efficacy than intermittent injections in reducing the tumor load while the treatment is ongoing. We also find that the percentage of tumor reduction under anti-IL-35 treatment improves when the production of IL-35 by cancer is increased.

Choindroitinase ABC I-mediated enhancement of oncolytic virus spread and anti tumor efficacy: a mathematical model

Oncolytic viruses are genetically engineered viruses that are designed to kill cancer cells while doing minimal damage to normal healthy tissue. After being injected into a tumor, they infect cancer cells, multiply inside them, and when a cancer cell is killed they move on to spread and infect other cancer cells. Chondroitinase ABC (Chase-ABC) is a bacterial enzyme that can remove a major glioma ECM component, chondroitin sulfate glycosoamino glycans from proteoglycans without any deleterious effects in vivo. It has been shown that Chase-ABC treatment is able to promote the spread of the viruses, increasing the efficacy of the viral treatment. In this paper we develop a mathematical model to investigate the effect of the Chase-ABC on the treatment of glioma by oncolytic viruses (OV). We show that the model's predictions agree with experimental results for a spherical glioma. We then use the model to test various treatment options in the heterogeneous microenvironment of the brain. The model predicts that separate injections of OV, one into the center of the tumor and another outside the tumor will result in better outcome than if the total injection is outside the tumor. In particular, the injection of the ECM-degrading enzyme (Chase-ABC) on the periphery of the main tumor core need to be administered in an optimal strategy in order to infect and eradicate the infiltrating glioma cells outside the tumor core in addition to proliferative cells in the bulk of tumor core. The model also predicts that the size of tumor satellites and distance between the primary tumor and multifocal/satellite lesions may be an important factor for the efficacy of the viral therapy with Chase treatment.

Mathematical model of macrophage-facilitated breast cancer cells invasion

Mortality from breast cancer stems from its tendency to invade into surrounding tissues and organs. Experiments have shown that this metastatic process is facilitated by macrophages in a short-ranged chemical signalling loop. Macrophages secrete epidermal growth factor, EGF, and respond to the colony stimulating factor 1, CSF-1. Tumor cells secrete CSF-1 and respond to EGF. In this way, the cells coordinate aggregation and cooperative migration. Here we investigate this process in a model for in vitro interactions using two distinct but related mathematical approaches. In the first, we analyze and simulate a set of partial differential equations to determine conditions for aggregation. In the second, we use a cell-based discrete 3D simulation to follow the fates and motion of individual cells during aggregation. Linear stability analysis of the PDE model reveals that decreasing the chemical secretion, chemotaxis coefficients or density of cells or increasing the chemical degradation in the model could eliminate the spontaneous aggregation of cells. Simulations with the discrete model show that the ratio between tumor cells and macrophages in aggregates increases when the EGF secretion parameter is increased. The results also show how CSF-1/CSF-1R autocrine signalling in tumor cells affects the ratio between the two cell types. Comparing the continuum results with simulations of a discrete cell-based model, we find good qualitative agreement.

The LDL-HDL profile determines the risk of atherosclerosis: a mathematical model

Atherosclerosis, the leading death in the United State, is a disease in which a plaque builds up inside the arteries. As the plaque continues to grow, the shear force of the blood flow through the decreasing cross section of the lumen increases. This force may eventually cause rupture of the plaque, resulting in the formation of thrombus, and possibly heart attack. It has long been recognized that the formation of a plaque relates to the cholesterol concentration in the blood. For example, individuals with LDL above 190 mg/dL and HDL below 40 mg/dL are at high risk, while individuals with LDL below 100 mg/dL and HDL above 50 mg/dL are at no risk. In this paper, we developed a mathematical model of the formation of a plaque, which includes the following key variables: LDL and HDL, free radicals and oxidized LDL, MMP and TIMP, cytockines: MCP-1, IFN-γ, IL-12 and PDGF, and cells: macrophages, foam cells, T cells and smooth muscle cells. The model is given by a system of partial differential equations with in evolving plaque. Simulations of the model show how the combination of the concentrations of LDL and HDL in the blood determine whether a plaque will grow or disappear. More precisely, we create a map, showing the risk of plaque development for any pair of values (LDL,HDL).

Modeling the role of TGF-β in regulation of the Th17 phenotype in the LPS-driven immune system

Airway exposure levels of lipopolysaccharide (LPS) are known to determine type I versus type II helper T cell induced experimental asthma. While low doses of LPS derive Th2 inflammatory responses, high (and/or intermediate) LPS levels induce Th1- or Th17-dominant responses. The present paper develops a mathematical model of the phenotypic switches among three Th phenotypes (Th1, Th2, and Th17) in response to various LPS levels. In the present work, we simplify the complex network of the interactions between cells and regulatory molecules. The model describes the nonlinear cross-talks between the IL-4/Th2 activities and a key regulatory molecule, transforming growth factor β (TGF-β), in response to high, intermediate, and low levels of LPS. The model characterizes development of three phenotypes (Th1, Th2, and Th17) and predicts the onset of a new phenotype, Th17, under the tight control of TGF-β. Analysis of the model illustrates the mono-, bi-, and oneway-switches in the key regulatory parameter sets in the absence or presence of time delays. The model also predicts coexistence of those phenotypes and Th1- or Th2-dominant immune responses in a spatial domain under various biochemical and bio-mechanical conditions in the microenvironment.

The Need for Integrative Computational Oncology: An Illustrated Example through MMP-Mediated Tissue Degradation

Physical oncology is a growing force in cancer research, and it is enhanced by integrative computational oncology: the fusion of novel experiments with mathematical and computational modeling. Computational models must be assessed with accurate numerical methods on correctly scaled tissues to avoid numerical artifacts that can cloud analysis. Simulation-driven analyses can only be validated by careful experiments. In this perspectives piece, we evaluate a current, widespread model of matrix metalloproteinase-driven tissue degradation during cancer invasion to illustrate that integrative computational oncology may not realize its fullest potential if either of these critical steps is neglected.

Evolution and morphology of microenvironment-enhanced malignancy of three-dimensional invasive solid tumors

The emergence of invasive and metastatic behavior in malignant tumors can often lead to fatal outcomes for patients. The collective malignant tumor behavior resulting from the complex tumor-host interactions and the interactions between the tumor cells is currently poorly understood. In this paper, we employ a cellular automaton (CA) model to investigate microenvironment-enhanced malignant behaviors and morphologies of in vitro avascular invasive solid tumors in three dimensions. Our CA model incorporates a variety of microscopic-scale tumor-host interactions, including the degradation of the extracellular matrix by the malignant cells, nutrient-driven cell migration, pressure buildup due to the deformation of the microenvironment by the growing tumor, and its effect on the local tumor-host interface stability. Moreover, the effects of cell-cell adhesion on tumor growth are explicitly taken into account. Specifically, we find that while strong cell-cell adhesion can suppress the invasive behavior of the tumors growing in soft microenvironments, cancer malignancy can be significantly enhanced by harsh microenvironmental conditions, such as exposure to high pressure levels. We infer from the simulation results a qualitative phase diagram that characterizes the expected malignant behavior of invasive solid tumors in terms of two competing malignancy effects: the rigidity of the microenvironment and cell-cell adhesion. This diagram exhibits phase transitions between noninvasive and invasive behaviors. We also discuss the implications of our results for the diagnosis, prognosis, and treatment of malignant tumors.

Inter-cellular signaling network reveals a mechanistic transition in tumor microenvironment

We conducted inter-cellular cytokine correlation and network analysis based upon a stochastic population dynamics model that comprises five cell types and fifteen signaling molecules inter-connected through a large number of cell-cell communication pathways. We observed that the signaling molecules are tightly correlated even at very early stages (e.g. the first month) of human glioma, but such correlation rapidly diminishes when tumor grows to a size that can be clinically detected. Further analysis suggests that paracrine is shown to be the dominant force during tumor initiation and priming, while autocrine supersedes it and supports a robust tumor expansion. In correspondence, the cytokine correlation network evolves through an increasing to decreasing complexity. This study indicates a possible mechanistic transition from the microenvironment-controlled, paracrine-based regulatory mechanism to self-sustained rapid progression to fetal malignancy. It also reveals key nodes that are responsible for such transition and can be potentially harnessed for the design of new anti-cancer therapies.

A hybrid model of tumor-stromal interactions in breast cancer

Ductal carcinoma in situ (DCIS) is an early stage noninvasive breast cancer that originates in the epithelial lining of the milk ducts, but it can evolve into comedo DCIS and ultimately, into the most common type of breast cancer, invasive ductal carcinoma. Understanding the progression and how to effectively intervene in it presents a major scientific challenge. The extracellular matrix (ECM) surrounding a duct contains several types of cells and several types of growth factors that are known to individually affect tumor growth, but at present the complex biochemical and mechanical interactions of these stromal cells and growth factors with tumor cells is poorly understood. Here we develop a mathematical model that incorporates the cross-talk between stromal and tumor cells, which can predict how perturbations of the local biochemical and mechanical state influence tumor evolution. We focus on the EGF and TGF-β signaling pathways and show how up- or down-regulation of components in these pathways affects cell growth and proliferation. We then study a hybrid model for the interaction of cells with the tumor microenvironment (TME), in which epithelial cells (ECs) are modeled individually while the ECM is treated as a continuum, and show how these interactions affect the early development of tumors. Finally, we incorporate breakdown of the epithelium into the model and predict the early stages of tumor invasion into the stroma. Our results shed light on the interactions between growth factors, mechanical properties of the ECM, and feedback signaling loops between stromal and tumor cells, and suggest how epigenetic changes in transformed cells affect tumor progression.

Mathematical model of colitis-associated colon cancer

As a result of chronic inflammation of their colon, patients with ulcerative colitis or Crohn's disease are at risk of developing colon cancer. In this paper, we consider the progression of colitis-associated colon cancer. Unlike normal colon mucosa, the inflammed colon mucosa undergoes genetic mutations, affecting, in particular, tumor suppressors TP53 and adenomatous polyposis coli (APC) gene. We develop a mathematical model that involves these genes, under chronic inflammation, as well as NF-κB, β-catenin, MUC1 and MUC2. The model demonstrates that increased level of cells with TP53 mutations results in abnormal growth and proliferation of the epithelium; further increase in the epithelium proliferation results from additional APC mutations. The model may serve as a conceptual framework for further data-based study of the early stage of colon cancer.

Study of tumor growth under hyperthermia condition

The new concept of keeping primary tumor under control in situ to suppress distant foci sheds light on the treatment of metastatic tumor. Hyperthermia is considered as one of the means for controlling tumor growth. To simulate the tumor growth, a continuum mathematical model has been introduced. The newest understanding of the Warburg effect on the cellular metabolism and diffusion of the nutrients in the tissue has been taken into consideration. The numerical results are compared with the in vivo experimental data by fitting the tumor cell doubling time/tumor cell growth rate under different thermal conditions. Both the tumor growth curve and corresponding average glucose concentration have been predicted. The numerical results have quantitatively illustrated the controlling effect on tumor growth under hyperthermia condition in the initial stage.

Diversity of dynamics and morphologies of invasive solid tumors

Complex tumor-host interactions can significantly affect the growth dynamics and morphologies of progressing neoplasms. The growth of a confined solid tumor induces mechanical pressure and deformation of the surrounding microenvironment, which in turn influences tumor growth. In this paper, we generalize a recently developed cellular automaton model for invasive tumor growth in heterogeneous microenvironments [Y. Jiao and S. Torquato, PLoS Comput. Biol.7, e1002314 (2011)] by incorporating the effects of pressure. Specifically, we explicitly model the pressure exerted on the growing tumor due to the deformation of the microenvironment and its effect on the local tumor-host interface instability. Both noninvasive-proliferative growth and invasive growth with individual cells that detach themselves from the primary tumor and migrate into the surrounding microenvironment are investigated. We find that while noninvasive tumors growing in "soft" homogeneous microenvironments develop almost isotropic shapes, both high pressure and host heterogeneity can strongly enhance malignant behavior, leading to finger-like protrusions of the tumor surface. Moreover, we show that individual invasive cells of an invasive tumor degrade the local extracellular matrix at the tumor-host interface, which diminishes the fingering growth of the primary tumor. The implications of our results for cancer diagnosis, prognosis and therapy are discussed.

Emergent behaviors from a cellular automaton model for invasive tumor growth in heterogeneous microenvironments

Understanding tumor invasion and metastasis is of crucial importance for both fundamental cancer research and clinical practice. In vitro experiments have established that the invasive growth of malignant tumors is characterized by the dendritic invasive branches composed of chains of tumor cells emanating from the primary tumor mass. The preponderance of previous tumor simulations focused on non-invasive (or proliferative) growth. The formation of the invasive cell chains and their interactions with the primary tumor mass and host microenvironment are not well understood. Here, we present a novel cellular automaton (CA) model that enables one to efficiently simulate invasive tumor growth in a heterogeneous host microenvironment. By taking into account a variety of microscopic-scale tumor-host interactions, including the short-range mechanical interactions between tumor cells and tumor stroma, degradation of the extracellular matrix by the invasive cells and oxygen/nutrient gradient driven cell motions, our CA model predicts a rich spectrum of growth dynamics and emergent behaviors of invasive tumors. Besides robustly reproducing the salient features of dendritic invasive growth, such as least-resistance paths of cells and intrabranch homotype attraction, we also predict nontrivial coupling between the growth dynamics of the primary tumor mass and the invasive cells. In addition, we show that the properties of the host microenvironment can significantly affect tumor morphology and growth dynamics, emphasizing the importance of understanding the tumor-host interaction. The capability of our CA model suggests that sophisticated in silico tools could eventually be utilized in clinical situations to predict neoplastic progression and propose individualized optimal treatment strategies.

The goal of the present work is to develop an efficient single-cell based cellular automaton (CA) model that enables one to investigate the growth dynamics and morphology of invasive solid tumors. Recent experiments have shown that highly malignant tumors develop dendritic branches composed of tumor cells that follow each other, which massively invade into the host microenvironment and ultimately lead to cancer metastasis. Previous theoretical/computational cancer modeling neither addressed the question of how such chain-like invasive branches form nor how they interact with the host microenvironment and the primary tumor. Our CA model, which incorporates a variety of microscopic-scale tumor-host interactions (e.g., the mechanical interactions between tumor cells and tumor stroma, degradation of the extracellular matrix by the tumor cells and oxygen/nutrient gradient driven cell motions), can robustly reproduce experimentally observed invasive tumor evolution and predict a wide spectrum of invasive tumor growth dynamics and emergent behaviors in various different heterogeneous environments. Further refinement of our CA model could eventually lead to the development of a powerful simulation tool for clinical purposes capable of predicting neoplastic progression and suggesting individualized optimal treatment strategies.

Mechanistic modeling of the effects of myoferlin on tumor cell invasion

Myoferlin (MYOF) is a member of the evolutionarily conserved ferlin family of proteins, noted for their role in a variety of membrane processes, including endocytosis, repair, and vesicular transport. Notably, ferlins are implicated in Caenorhabditis elegans sperm motility (Fer-1), mammalian skeletal muscle development and repair (MYOF and dysferlin), and presynaptic transmission in the auditory system (otoferlin). In this paper, we demonstrate that MYOF plays a previously unrecognized role in cancer cell invasion, using a combination of mathematical modeling and in vitro experiments. Using a real-time impedance-based invasion assay (xCELLigence), we have shown that lentiviral-based knockdown of MYOF significantly reduced invasion of MDA-MB-231 breast cancer cells in Matrigel bioassays. Based on these experimental data, we developed a partial differential equation model of MYOF effects on cancer cell invasion, which we used to generate mechanistic hypotheses. The mathematical model predictions revealed that matrix metalloproteinases (MMPs) may play a key role in modulating this invasive property, which was supported by experimental data using quantitative RT-PCR screens. These results suggest that MYOF may be a promising target for biomarkers or drug target for metastatic cancer diagnosis and therapy, perhaps mediated through MMPs.

The role of the microenvironment in tumor growth and invasion

Mathematical modeling and computational analysis are essential for understanding the dynamics of the complex gene networks that control normal development and homeostasis, and can help to understand how circumvention of that control leads to abnormal outcomes such as cancer. Our objectives here are to discuss the different mechanisms by which the local biochemical and mechanical microenvironment, which is comprised of various signaling molecules, cell types and the extracellular matrix (ECM), affects the progression of potentially-cancerous cells, and to present new results on two aspects of these effects. We first deal with the major processes involved in the progression from a normal cell to a cancerous cell at a level accessible to a general scientific readership, and we then outline a number of mathematical and computational issues that arise in cancer modeling. In Section 2 we present results from a model that deals with the effects of the mechanical properties of the environment on tumor growth, and in Section 3 we report results from a model of the signaling pathways and the tumor microenvironment (TME), and how their interactions affect the development of breast cancer. The results emphasize anew the complexities of the interactions within the TME and their effect on tumor growth, and show that tumor progression is not solely determined by the presence of a clone of mutated immortal cells, but rather that it can be 'community-controlled'.

Tumor cells proliferation and migration under the influence of their microenvironment

It is well known that tumor microenvironment affects tumor growth and metastasis: Tumor cells may proliferate at different rates and migrate in different patterns depending on the microenvironment in which they are embedded. There is a huge literature that deals with mathematical models of tumor growth and proliferation, in both the avascular and vascular phases. In particular, a review of the literature of avascular tumor growth (up to 2006) can be found in Lolas (G. Lolas, Lecture Notes in Mathematics, Springer Berlin / Heidelberg, 1872, 77 (2006)). In this article we report on some of our recent work. We consider two aspects, proliferation and of migration, and describe mathematical models based on in vitro experiments. Simulations of the models are in agreement with experimental results. The models can be used to generate hypotheses regarding the development of drugs which will confine tumor growth.