Simulating photodynamic therapy for the treatment of glioblastoma using Monte Carlo radiative transport.

Significance: Glioblastoma (GBM) is a rare but deadly form of brain tumor with a low median survival rate of 14.6 months, due to its resistance to treatment. An independent simulation of the INtraoperative photoDYnamic therapy for GliOblastoma (INDYGO) trial, a clinical trial aiming to treat the GBM resection cavity with photodynamic therapy (PDT) via a laser coupled balloon device, is demonstrated. Aim: To develop a framework providing increased understanding for the PDT treatment, its parameters, and their impact on the clinical outcome. Approach: We use Monte Carlo radiative transport techniques within a computational brain model containing a GBM to simulate light path and PDT effects. Treatment parameters (laser power, photosensitizer concentration, and irradiation time) are considered, as well as PDT's impact on brain tissue temperature. Results: The simulation suggests that 39% of post-resection GBM cells are killed at the end of treatment when using the standard INDYGO trial protocol (light fluence = 200 J/cm2 at balloon wall) and assuming an initial photosensitizer concentration of 5 µM. Increases in treatment time and light power (light fluence = 400 J/cm2 at balloon wall) result in further cell kill but increase brain cell temperature, which potentially affects treatment safety. Increasing the p hotosensitizer concentration produces the most significant increase in cell kill, with 61% of GBM cells killed when doubling concentration to 10 µM and keeping the treatment time and power the same. According to these simulations, the standard trial protocol is reasonably well optimized with improvements in cell kill difficult to achieve without potentially dangerous increases in temperature. To improve treatment outcome, focus should be placed on improving the photosensitizer. Conclusions: With further development and optimization, the simulation could have potential clinical benefit and be used to help plan and optimize intraoperative PDT treatment for GBM.

A computational investigation of COVID-19 transmission inside hospital wards and associated costs.

The COVID-19 pandemic has placed a particular burden on hospitals: from intra-hospital transmission of the infections to reduced admissions of non-COVID-19 patients. There are also high costs associated with the treatment of hospitalised COVID-19 patients, as well as reductions in revenues due to delayed and cancelled treatments. In this study we investigate computationally the transmission of COVID-19 inside a hospital ward that contains multiple-bed bays (with 4 or 6 beds) and multiple single-bed side rooms (that can accommodate the contacts of COVID-19-positive patients). The aim of this study is to investigate the role of 4-bed bays vs. 6-bed bays on the spread of infections and the hospital costs. We show that 4-bed bays are associated with lower infections only when we reduce the discharge time of some patients from 10 days to 5 days. This also leads to lower costs for the treatment of COVID-19 patients. In contrast, 6-bed bays are associated with reduced hospital waiting lists (especially when there are also multiple side rooms available to accommodate the contacts of COVID-19-positive patients identified inside the 6-bed bays).

Nonlocal multiscale modelling of tumour-oncolytic viruses interactions within a heterogeneous fibrous/non-fibrous extracellular matrix.

In this study we investigate computationally tumour-oncolytic virus (OV) interactions that take place within a heterogeneous extracellular matrix (ECM). The ECM is viewed as a mixture of two constitutive phases, namely a fibre phase and a non-fibre phase. The multiscale mathematical model presented here focuses on the nonlocal cell-cell and cell-ECM interactions, and how these interactions might be impacted by the infection of cancer cells with the OV. At macroscale we track the kinetics of cancer cells, virus particles and the ECM. At microscale we track (i) the degradation of ECM by matrix degrading enzymes (MDEs) produced by cancer cells, which further influences the movement of tumour boundary; (ii) the re-arrangement of the microfibres that influences the re-arrangement of macrofibres (i.e., fibres at macroscale). With the help of this new multiscale model, we investigate two questions: (i) whether the infected cancer cell fluxes are the result of local or non-local advection in response to ECM density; and (ii) what is the effect of ECM fibres on the the spatial spread of oncolytic viruses and the outcome of oncolytic virotherapy.

Computational Model of Heterogeneity in Melanoma: Designing Therapies and Predicting Outcomes.

Cutaneous melanoma is a highly invasive tumor and, despite the development of recent therapies, most patients with advanced metastatic melanoma have a poor clinical outcome. The most frequent mutations in melanoma affect the BRAF oncogene, a protein kinase of the MAPK signaling pathway. Therapies targeting both BRAF and MEK are effective for only 50% of patients and, almost systematically, generate drug resistance. Genetic and non-genetic mechanisms associated with the strong heterogeneity and plasticity of melanoma cells have been suggested to favor drug resistance but are still poorly understood. Recently, we have introduced a novel mathematical formalism allowing the representation of the relation between tumor heterogeneity and drug resistance and proposed several models for the development of resistance of melanoma treated with BRAF/MEK inhibitors. In this paper, we further investigate this relationship by using a new computational model that copes with multiple cell states identified by single cell mRNA sequencing data in melanoma treated with BRAF/MEK inhibitors. We use this model to predict the outcome of different therapeutic strategies. The reference therapy, referred to as "continuous" consists in applying one or several drugs without disruption. In "combination therapy", several drugs are used sequentially. In "adaptive therapy" drug application is interrupted when the tumor size is below a lower threshold and resumed when the size goes over an upper threshold. We show that, counter-intuitively, the optimal protocol in combination therapy of BRAF/MEK inhibitors with a hypothetical drug targeting cell states that develop later during the tumor response to kinase inhibitors, is to treat first with this hypothetical drug. Also, even though there is little difference in the timing of emergence of the resistance between continuous and adaptive therapies, the spatial distribution of the different melanoma subpopulations is more zonated in the case of adaptive therapy.

Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics.

Cancer cell mutations occur when cells undergo multiple cell divisions, and these mutations can be spontaneous or environmentally-induced. The mechanisms that promote and sustain these mutations are still not fully understood. This study deals with the identification (or reconstruction) of the usually unknown cancer cell mutation law, which lead to the transformation of a primary tumour cell population into a secondary, more aggressive cell population. We focus on local and nonlocal mathematical models for cell dynamics and movement, and identify these mutation laws from macroscopic tumour snapshot data collected at some later stage in the tumour evolution. In a local cancer invasion model, we first reconstruct the mutation law when we assume that the mutations depend only on the surrounding cancer cells (i.e., the ECM plays no role in mutations). Second, we assume that the mutations depend on the ECM only, and we reconstruct the mutation law in this case. Third, we reconstruct the mutation when we assume that there is no prior knowledge about the mutations. Finally, for the nonlocal cancer invasion model, we reconstruct the mutation law that depends on the cancer cells and on the ECM. For these numerical reconstructions, our approximations are based on the finite difference method combined with the finite elements method. As the inverse problem is ill-posed, we use the Tikhonov regularisation technique in order to regularise the solution. Stability of the solution is examined by adding additive noise into the measurements.

Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions.

We propose and study computationally a novel non-local multiscale moving boundary mathematical model for tumour and oncolytic virus (OV) interactions when we consider the go or grow hypothesis for cancer dynamics. This spatio-temporal model focuses on two cancer cell phenotypes that can be infected with the OV or remain uninfected, and which can either move in response to the extracellular-matrix (ECM) density or proliferate. The interactions between cancer cells, those among cancer cells and ECM, and those among cells and OV occur at the macroscale. At the micro-scale, we focus on the interactions between cells and matrix degrading enzymes (MDEs) that impact the movement of tumour boundary. With the help of this multiscale model we explore the impact on tumour invasion patterns of two different assumptions that we consider in regard to cell-cell and cell-matrix interactions. In particular we investigate model dynamics when we assume that cancer cell fluxes are the result of local advection in response to the density of extracellular matrix (ECM), or of non-local advection in response to cell-ECM adhesion. We also investigate the role of the transition rates between mainly-moving and mainly-growing cancer cell sub-populations, as well as the role of virus infection rate and virus replication rate on the overall tumour dynamics.

Collective Cell Migration in a Fibrous Environment: A Hybrid Multiscale Modelling Approach.

The specific structure of the extracellular matrix (ECM), and in particular the density and orientation of collagen fibres, plays an important role in the evolution of solid cancers. While many experimental studies discussed the role of ECM in individual and collective cell migration, there are still unanswered questions about the impact of nonlocal cell sensing of other cells on the overall shape of tumour aggregation and its migration type. There are also unanswered questions about the migration and spread of tumour that arises at the boundary between different tissues with different collagen fibre orientations. To address these questions, in this study we develop a hybrid multi-scale model that considers the cells as individual entities and ECM as a continuous field. The numerical simulations obtained through this model match experimental observations, confirming that tumour aggregations are not moving if the ECM fibres are distributed randomly, and they only move when the ECM fibres are highly aligned. Moreover, the stationary tumour aggregations can have circular shapes or irregular shapes (with finger-like protrusions), while the moving tumour aggregations have elongate shapes (resembling to clusters, strands or files). We also show that the cell sensing radius impacts tumour shape only when there is a low ratio of fibre to non-fibre ECM components. Finally, we investigate the impact of different ECM fibre orientations corresponding to different tissues, on the overall tumour invasion of these neighbouring tissues.

Directionality of Macrophages Movement in Tumour Invasion: A Multiscale Moving-Boundary Approach.

Invasion of the surrounding tissue is one of the recognised hallmarks of cancer (Hanahan and Weinberg in Cell 100: 57-70, 2000. https://doi.org/10.1016/S0092-8674(00)81683-9 ), which is accomplished through a complex heterotypic multiscale dynamics involving tissue-scale random and directed movement of the population of both cancer cells and other accompanying cells (including here, the family of tumour-associated macrophages) as well as the emerging cell-scale activity of both the matrix-degrading enzymes and the rearrangement of the cell-scale constituents of the extracellular matrix (ECM) fibres. The involved processes include not only the presence of cell proliferation and cell adhesion (to other cells and to the extracellular matrix), but also the secretion of matrix-degrading enzymes. This is as a result of cancer cells as well as macrophages, which are one of the most abundant types of immune cells in the tumour micro-environment. In large tumours, these tumour-associated macrophages (TAMs) have a tumour-promoting phenotype, contributing to tumour proliferation and spread. In this paper, we extend a previous multiscale moving-boundary mathematical model for cancer invasion, by considering also the multiscale effects of TAMs, with special focus on the influence that their directional movement exerts on the overall tumour progression. Numerical investigation of this new model shows the importance of the interactions between pro-tumour TAMs and the fibrous ECM, highlighting the impact of the fibres on the spatial structure of solid tumour.

Cell-Scale Degradation of Peritumoural Extracellular Matrix Fibre Network and Its Role Within Tissue-Scale Cancer Invasion.

Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. During the complex interaction between cancer cell population and the extracellular matrix (ECM), of key importance is the role played by both bulk two-scale dynamics of ECM fibres within collective movement of the tumour cells and the multiscale leading edge dynamics driven by proteolytic activity of the matrix-degrading enzymes (MDEs) that are secreted by the cancer cells. As these two multiscale subsystems share and contribute to the same tumour macro-dynamics, in this work we develop further the model introduced in Shuttleworth and Trucu (Bull Math Biol 81:2176-2219, 2019. https://doi.org/10.1007/s11538-019-00598-w) by exploring a new aspect of their interaction that occurs at the cell scale. Specifically, here we will focus on understanding the cell-scale cross talk between the micro-scale parts of these two multiscale subsystems which get to interact directly in the peritumoural region, with immediate consequences both for MDE micro-dynamics occurring at the leading edge of the tumour and for the cell-scale rearrangement of the naturally oriented ECM fibres in the peritumoural region, ultimately influencing the way tumour progresses in the surrounding tissue. To that end, we will propose a new modelling that captures the ECM fibres degradation not only at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour.

Oncolytic viral therapies and the delicate balance between virus-macrophage-tumour interactions: A mathematical approach.

The success of oncolytic virotherapies depends on the tumour microenvironment, which contains a large number of infiltrating immune cells. In this theoretical study, we derive an ODE model to investigate the interactions between breast cancer tumour cells, an oncolytic virus (Vesicular Stomatitis Virus), and tumour-infiltrating macrophages with different phenotypes which can impact the dynamics of oncolytic viruses. The complexity of the model requires a combined analytical-numerical approach to understand the transient and asymptotic dynamics of this model. We use this model to propose new biological hypotheses regarding the impact on tumour elimination/relapse/persistence of: (i) different macrophage polarisation/re-polarisation rates; (ii) different infection rates of macrophages and tumour cells with the oncolytic virus; (iii) different viral burst sizes for macrophages and tumour cells. We show that increasing the rate at which the oncolytic virus infects the tumour cells can delay tumour relapse and even eliminate tumour. Increasing the rate at which the oncolytic virus particles infect the macrophages can trigger transitions between steady-state dynamics and oscillatory dynamics, but it does not lead to tumour elimination unless the tumour infection rate is also very large. Moreover, we confirm numerically that a large tumour-induced M1→M2 polarisation leads to fast tumour growth and fast relapse (if the tumour was reduced before by a strong anti-tumour immune and viral response). The increase in viral-induced M2→M1 re-polarisation reduces temporarily the tumour size, but does not lead to tumour elimination. Finally, we show numerically that the tumour size is more sensitive to the production of viruses by the infected macrophages.

Multiscale dynamics of a heterotypic cancer cell population within a fibrous extracellular matrix.

Local cancer cell invasion is a complex process involving many cellular and tissue interactions and is an important prerequisite for metastatic spread, the main cause of cancer related deaths. As a tumour increases in malignancy, the cancer cells adopt the ability to mutate into secondary cell subpopulations giving rise to a heterogeneous tumour. This new cell subpopulation often carries higher invasive abilities and permits a quicker spread of the tumour. Building upon the recent multiscale modelling framework for cancer invasion within a fibrous ECM introduced in Shuttleworth and Trucu, (2019), in this paper we consider the process of local invasion by a heterotypic tumour consisting of two cancer cell populations mixed with a two-phase ECM. To that end, we address the double feedback link between the tissue-scale cancer dynamics and the cell-scale molecular processes through the development of a two-part modelling framework that crucially incorporates the multiscale dynamic redistribution of oriented fibres occurring within a two-phase extra-cellular matrix and combines this with the multiscale leading edge dynamics exploring key matrix-degrading enzymes molecular processes along the tumour interface that drive the movement of the cancer boundary. The modelling framework will be accompanied by computational results that explore the effects of the underlying fibre network on the overall pattern of cancer invasion.

Multiscale moving boundary modelling of cancer interactions with a fusogenic oncolytic virus: The impact of syncytia dynamics.

Oncolytic viral therapies is one of the new promising strategies against cancer, due to the ability of oncolytic viruses to specifically replicate inside cancer cells and kill them. There is increasing evidence that a sub-class of viruses that contain fusion proteins (triggering the formation of syncytia) can lead to better oncolytic results. Since the details of the tumour dynamics following syncytia formation are not fully understood, in this study we consider a modelling and computational approach to describe the effect of a fusogenic oncolytic virus on the multiscale dynamics of a spreading tumour. We show that for the baseline parameter values considered here, small syncytia diffusion coefficient leads to tumour reduction. Further tumour reduction can be obtained when we increase the probability of syncytia formation, in the context of different viral burst rates and death rates for individually-infected tumour cells and syncytia structures. Finally, we show that the type of syncytia diffusion coefficient (i.e., constant or density dependent) also impacts the outcome of the oncolytic viral therapy.

Multiscale Modelling of Fibres Dynamics and Cell Adhesion within Moving Boundary Cancer Invasion.

Recognised as one of the hallmarks of cancer, local cancer cell invasion is a complex multiscale process that combines the secretion of matrix-degrading enzymes with a series of altered key cell processes (such as abnormal cell proliferation and changes in cell-cell and cell-matrix adhesion leading to enhanced migration) to degrade important components of the surrounding extracellular matrix (ECM) and this way spread further in the human tissue. In order to gain a deeper understanding of the invasion process, we pay special attention to the interacting dynamics between the cancer cell population and various constituents of the surrounding tumour microenvironment. To that end, we consider the key role that ECM plays within the human body tissue, and in particular we focus on the special contribution of its fibrous proteins components, such as collagen and fibronectin, which play an important part in cell proliferation and migration. In this work, we consider the two-scale dynamic cross-talk between cancer cells and a two-component ECM (consisting of both a fibre and a non-fibre phase). To that end, we incorporate the interlinked two-scale dynamics of cell-ECM interactions within the tumour support that contributes simultaneously both to cell adhesion and to the dynamic rearrangement and restructuring of the ECM fibres. Furthermore, this is embedded within a multiscale moving boundary approach for the invading cancer cell population, in the presence of cell adhesion at the tissue scale and cell-scale fibre redistribution activity and leading edge matrix-degrading enzyme molecular proteolytic processes. The overall modelling framework will be accompanied by computational results that will explore the impact on cancer invasion patterns of different levels of cell adhesion in conjunction with the continuous ECM fibres rearrangement.

Multiscale modelling of cancer response to oncolytic viral therapy.

Oncolytic viruses (OV) are viruses that can replicate selectively within cancer cells and destroy them. While the past few decades have seen significant progress related to the use of these viruses in clinical contexts, the success of oncolytic therapies is dampened by the complex spatial tumour-OV interactions. In this work, we present a novel multiscale moving boundary modelling for the tumour-OV interactions, which is based on coupled systems of partial differential equations both at macro-scale (tissue-scale) and at micro-scale (cell-scale) that are connected through a double feedback link. At the macro-scale, we account for the coupled dynamics of uninfected cancer cells, OV-infected cancer cells, extracellular matrix (ECM) and oncolytic viruses. At the same time, at the micro scale, we focus on essential dynamics of urokinase plasminogen activator (uPA) system which is one of the important proteolytic systems responsible for the degradation of the ECM, with notable influence in cancer invasion. While sourced by the cancer cells that arrive during their macro-dynamics within the outer proliferating rim of the tumour, the uPA micro-dynamics is crucial in determining the movement of the macro-scale tumour boundary (both in terms of direction and displacement magnitude). In this investigation, we consider three scenarios for the macro-scale tumour-OV interactions. While assuming the usual context of reaction-diffusion-taxis coupled PDEs, the three macro-dynamics scenarios gradually explore the influence of the ECM taxis over the tumour - OV interaction, in the form of haptotaxis of both uninfected and infected cells populations as well as the indirect ECM taxis for the oncolytic virus. Finally, the complex tumour-OV interactions is investigated numerically through the development a new multiscale moving boundary computational framework. While further investigation is needed to validate the findings of our modelling, for the parameter regimes that we considered, our numerical simulations indicate that the viral therapy leads to control and decrease of the overall cancer expansion and in certain cases this can result even in the elimination of the tumour.

Spatio-Genetic and phenotypic modelling elucidates resistance and re-sensitisation to treatment in heterogeneous melanoma.

Although novel targeted therapies have significantly improved the overall survival of patients with advanced melanoma, understanding and combatting drug resistance remains a major clinical challenge. Using partial differential equations, we describe the evolution of a cellular population through time, space, and phenotype dimensions, in the presence of various drug species. We then use this framework to explore models in which resistance is attained by either mutations (irreversible) or plasticity (reversible). Numerical results suggest that punctuated evolutionary assumptions are more consistent with results obtained from murine melanoma models than gradual evolution. Furthermore, in the context of an evolving tumour cell population, sequencing the treatment, for instance applying immunotherapy before BRAF inhibitors, can increase treatment effectiveness. However, drug strategies which showed success within a spatially homogeneous tumour environment were unsuccessful under heterogeneous conditions, suggesting that spatio-environmental heterogeneity may be the greatest challenge to tumour therapies. Plastic metabolic models are additionally capable of reproducing the characteristic resistant tumour volume curves and predicting re-sensitisation to secondary waves of treatment observed in patient derived xenograft (PDX) melanomas treated with MEK and BRAF inhibitors. Nevertheless, secondary relapse due to a pre-adapted subpopulation, remaining after the first wave of treatment, results in a more rapid development of resistance. Our model provides a framework through which tumour resistance can be understood and would suggest that carefully phased treatments may be able to overcome the development of long-term resistance in melanoma.

Role of extracellular matrix and microenvironment in regulation of tumor growth and LAR-mediated invasion in glioblastoma.

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 cues in the microenvironment. Recent studies have shown that miR-451 regulates downstream molecules including AMPK/CAB39/MARK and mTOR to determine the balance between rapid proliferation and invasion in response to metabolic stress in the harsh tumor microenvironment. Surgical removal of the 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 complex process of cell proliferation and invasion and its response to conventional treatment, we propose a mathematical model that analyzes the intracellular dynamics of the miR-451-AMPK- mTOR-cell cycle signaling pathway within a cell. The model identifies a key mechanism underlying the molecular switches between proliferative phase and migratory phase in response to metabolic stress in response to fluctuating glucose levels. We show how up- or down-regulation of components in these pathways affects the key cellular decision to infiltrate or proliferate in a complex microenvironment in the absence and presence of time delays and stochastic noise. Glycosylated chondroitin sulfate proteoglycans (CSPGs), a major component of the extracellular matrix (ECM) in the brain, contribute to the physical structure of the local brain microenvironment but also induce or inhibit glioma invasion by regulating the dynamics of the CSPG receptor LAR as well as the spatiotemporal activation status of resident astrocytes and tumor-associated microglia. Using a multi-scale mathematical model, we investigate a CSPG-induced switch between invasive and non-invasive tumors through the coordination of ECM-cell adhesion and dynamic changes in stromal cells. We show that the CSPG-rich microenvironment is associated with non-invasive tumor lesions through LAR-CSGAG binding while the absence of glycosylated CSPGs induce the critical glioma invasion. We illustrate how high molecular weight CSPGs can regulate the exodus of local reactive astrocytes from the main tumor lesion, leading to encapsulation of non-invasive tumor and inhibition of tumor invasion. These different CSPG conditions also change the spatial profiles of ramified and activated microglia. The complex distribution of CSPGs in the tumor microenvironment can determine the nonlinear invasion behaviors of glioma cells, which suggests the need for careful therapeutic strategies.

Signal Propagation in Sensing and Reciprocating Cellular Systems with Spatial and Structural Heterogeneity.

Sensing and reciprocating cellular systems (SARs) are important for the operation of many biological systems. Production in interferon (IFN) SARs is achieved through activation of the Jak-Stat pathway, and downstream upregulation of IFN regulatory factor (IRF)-7 and IFN transcription, but the role that high- and low-affinity IFNs play in this process remains unclear. We present a comparative between a minimal spatio-temporal partial differential equation model and a novel spatio-structural-temporal (SST) model for the consideration of receptor, binding, and metabolic aspects of SAR behaviour. Using the SST framework, we simulate single- and multi-cluster paradigms of IFN communication. Simulations reveal a cyclic process between the binding of IFN to the receptor, and the consequent increase in metabolism, decreasing the propensity for binding due to the internal feedback mechanism. One observes the effect of heterogeneity between cellular clusters, allowing them to individualise and increase local production, and within clusters, where we observe 'subpopular quiescence'; a process whereby intra-cluster subpopulations reduce their binding and metabolism such that other such subpopulations may augment their production. Finally, we observe the ability for low-affinity IFN to communicate a long range signal, where high affinity cannot, and the breakdown of this relationship through the introduction of cell motility. Biological systems may utilise cell motility where environments are unrestrictive and may use fixed system, with low-affinity communication, where a localised response is desirable.

Computational Approaches and Analysis for a Spatio-Structural-Temporal Invasive Carcinoma Model.

Spatio-temporal models have long been used to describe biological systems of cancer, but it has not been until very recently that increased attention has been paid to structural dynamics of the interaction between cancer populations and the molecular mechanisms associated with local invasion. One system that is of particular interest is that of the urokinase plasminogen activator (uPA) wherein uPA binds uPA receptors on the cancer cell surface, allowing plasminogen to be cleaved into plasmin, which degrades the extracellular matrix and this way leads to enhanced cancer cell migration. In this paper, we develop a novel numerical approach and associated analysis for spatio-structuro-temporal modelling of the uPA system for up to two-spatial and two-structural dimensions. This is accompanied by analytical exploration of the numerical techniques used in simulating this system, with special consideration being given to the proof of stability within numerical regimes encapsulating a central differences approach to approximating numerical gradients. The stability analysis performed here reveals instabilities induced by the coupling of the structural binding and proliferative processes. The numerical results expound how the uPA system aids the tumour in invading the local stroma, whilst the inhibitor to this system may impede this behaviour and encourage a more sporadic pattern of invasion.

Aggregation and travelling wave dynamics in a two-population model of cancer cell growth and invasion.

Cells adhere to each other and to the extracellular matrix (ECM) through protein molecules on the surface of the cells. The breaking and forming of adhesive bonds, a process critical in cancer invasion and metastasis, can be influenced by the mutation of cancer cells. In this paper, we develop a nonlocal mathematical model describing cancer cell invasion and movement as a result of integrin-controlled cell-cell adhesion and cell-matrix adhesion, for two cancer cell populations with different levels of mutation. The partial differential equations for cell dynamics are coupled with ordinary differential equations describing the ECM degradation and the production and decay of integrins. We use this model to investigate the role of cancer mutation on the possibility of cancer clonal competition with alternating dominance, or even competitive exclusion (phenomena observed experimentally). We discuss different possible cell aggregation patterns, as well as travelling wave patterns. In regard to the travelling waves, we investigate the effect of cancer mutation rate on the speed of cancer invasion.

Structured models of cell migration incorporating molecular binding processes.

The dynamic interplay between collective cell movement and the various molecules involved in the accompanying cell signalling mechanisms plays a crucial role in many biological processes including normal tissue development and pathological scenarios such as wound healing and cancer. Information about the various structures embedded within these processes allows a detailed exploration of the binding of molecular species to cell-surface receptors within the evolving cell population. In this paper we establish a general spatio-temporal-structural framework that enables the description of molecular binding to cell membranes coupled with the cell population dynamics. We first provide a general theoretical description for this approach and then illustrate it with three examples arising from cancer invasion.