Growth:Visualisation of predictive mathematical models using 3D computer graphics and animation

‘Growth’ is a short film created as a result of a multidisciplinary collaboration between artists and scientists from the University of Dundee. The author worked with mathematicians to investigate how techniques and technologies used in 3D computer animation and visual effects industries might support the enhanced visualisation of predictive mathematical modelling of solid tumour growth. This paper will discuss the practical and artistic processes of this visualisation research, including the technical innovations required in undertaking this work – such as the creation of custom tools for ‘reading’ data, or the addition of stereoscopic output. By transforming numerical data into three-dimensional ‘objects’, artists can provide new ways of ‘seeing’ information and identifying patterns or results. Developing visualisation techniques can be used to improve the communication of cancer growth to patients (by increasing patient understanding and relieving their fears) through the use of new and innovative visual material.<br/

Initial/boundary-value problems of tumor growth within a host tissue

This paper concerns multiphase models of tumor growth in interaction with a surrounding tissue, taking into account also the interplay with diffusible nutrients feeding the cells. Models specialize in nonlinear systems of possibly degenerate parabolic equations, which include phenomenological terms related to specific cell functions. The paper discusses general modeling guidelines for such terms, as well as for initial and boundary conditions, aiming at both biological consistency and mathematical robustness of the resulting problems. Particularly, it addresses some qualitative properties such as a priori nonnegativity, boundedness, and uniqueness of the solutions. Existence of the solutions is studied in the one-dimensional time-independent case.Comment: 30 pages, 5 figure

Spatio-temporal Models of Lymphangiogenesis in Wound Healing

Several studies suggest that one possible cause of impaired wound healing is failed or insufficient lymphangiogenesis, that is the formation of new lymphatic capillaries. Although many mathematical models have been developed to describe the formation of blood capillaries (angiogenesis), very few have been proposed for the regeneration of the lymphatic network. Lymphangiogenesis is a markedly different process from angiogenesis, occurring at different times and in response to different chemical stimuli. Two main hypotheses have been proposed: 1) lymphatic capillaries sprout from existing interrupted ones at the edge of the wound in analogy to the blood angiogenesis case; 2) lymphatic endothelial cells first pool in the wound region following the lymph flow and then, once sufficiently populated, start to form a network. Here we present two PDE models describing lymphangiogenesis according to these two different hypotheses. Further, we include the effect of advection due to interstitial flow and lymph flow coming from open capillaries. The variables represent different cell densities and growth factor concentrations, and where possible the parameters are estimated from biological data. The models are then solved numerically and the results are compared with the available biological literature.Comment: 29 pages, 9 Figures, 6 Tables (39 figure files in total

Modeling the emergence of phenotypic heterogeneity in vascularized tumors

We present a mathematical study of the emergence of phenotypic heterogeneity in vascularized tumors. Our study is based on formal asymptotic analysis and numerical simulations of a system of nonlocal parabolic equations that describes the phenotypic evolution of tumor cells and their nonlinear dynamic interactions with the oxygen, which is released from the intratumoral vascular network. Numerical simulations are carried out both in the case of arbitrary distributions of intratumor blood vessels and in the case where the intratumoral vascular network is reconstructed from clinical images obtained using dynamic optical coherence tomography. The results obtained support a more in-depth theoretical understanding of the eco-evolutionary process which underpins the emergence of phenotypic heterogeneity in vascularized tumors. In particular, our results offer a theoretical basis for empirical evidence indicating that the phenotypic properties of cancer cells in vascularized tumors vary with the distance from the blood vessels, and establish a relation between the degree of tumor tissue vascularization and the level of intratumor phenotypic heterogeneity

Value-Driven Analysis of New Paradigms in Space Architectures: An Ilities-Based Approach

Current commercial, civil, and military space architecture designs perform exquisitely and reliably. However, today’s architecture paradigms are also characterized by expensive launches, large and expensive high-performance spacecraft, long development cycles, and wide variations in ground architectures. While current assets provide high-quality services, and future assets are slated to improve performance within the same design frameworks, proposed future architectures may not be capitalizing on technology improvements, system innovations, or policy alternatives explored during the last two decades. This paper identifies five “trends” along which space architectures may develop, aimed at granting systems several “ilities,” such as resiliency, robustness, flexibility, scalability, and affordability. The trends examined include: commercialization of space, significant reductions in launch costs and the development of hybrid or reusable launch systems, development of on-orbit infrastructure and servicing, aggregation or disaggregation of orbital assets, and the automation and standardization of ground architectures. Further refinement of these key technological and system trends could result in major paradigm shifts in the development and fielding of space operations as well as lead to space architecture designs in the future that are radically different from those today. Within the framework of systems engineering ilities and risk management, this paper reviews current literature surrounding these new change trends and justifies their potential to cause significant paradigm shifts. By examining the work and research conducted so far through an ilities-based approach, systems engineers can more fully appreciate the value being offered by these trends

Evolutionary Dynamics in Vascularised Tumours under Chemotherapy: Mathematical Modelling, Asymptotic Analysis and Numerical Simulations

We consider a mathematical model for the evolutionary dynamics of tumour cells in vascularised tumours under chemotherapy. The model comprises a system of coupled partial integro-differential equations for the phenotypic distribution of tumour cells, the concentration of oxygen and the concentration of a chemotherapeutic agent. In order to disentangle the impact of different evolutionary parameters on the emergence of intra-tumour phenotypic heterogeneity and the development of resistance to chemotherapy, we construct explicit solutions to the equation for the phenotypic distribution of tumour cells and provide a detailed quantitative characterisation of the long-time asymptotic behaviour of such solutions. Analytical results are integrated with numerical simulations of a calibrated version of the model based on biologically consistent parameter values. The results obtained provide a theoretical explanation for the observation that the phenotypic properties of tumour cells in vascularised tumours vary with the distance from the blood vessels. Moreover, we demonstrate that lower oxygen levels may correlate with higher levels of phenotypic variability, which suggests that the presence of hypoxic regions supports intra-tumour phenotypic heterogeneity. Finally, the results of our analysis put on a rigorous mathematical basis the idea, previously suggested by formal asymptotic results and numerical simulations, that hypoxia favours the selection for chemoresistant phenotypic variants prior to treatment. Consequently, this facilitates the development of resistance following chemotherapy

A stochastic individual-based model to explore the role of spatial interactions and antigen recognition in the immune response against solid tumours

Spatial interactions between cancer and immune cells, as well as the recognition of tumour antigens by cells of the immune system, play a key role in the immune response against solid tumours. The existing mathematical models generally focus only on one of these key aspects. We present here a spatial stochastic individual-based model that explicitly captures antigen expression and recognition. In our model, each cancer cell is characterised by an antigen profile which can change over time due to either epimutations or mutations. The immune response against the cancer cells is initiated by the dendritic cells that recognise the tumour antigens and present them to the cytotoxic T cells. Consequently, T cells become activated against the tumour cells expressing such antigens. Moreover, the differences in movement between inactive and active immune cells are explicitly taken into account by the model. Computational simulations of our model clarify the conditions for the emergence of tumour clearance, dormancy or escape, and allow us to assess the impact of antigenic heterogeneity of cancer cells on the efficacy of immune action. Ultimately, our results highlight the complex interplay between spatial interactions and adaptive mechanisms that underpins the immune response against solid tumours, and suggest how this may be exploited to further develop cancer immunotherapies

Bridging the gap between individual-based and continuum models of growing cell populations

Continuum models for the spatial dynamics of growing cell populations have been widely used to investigate the mechanisms underpinning tissue development and tumour invasion. These models consist of nonlinear partial differential equations that describe the evolution of cellular densities in response to pressure gradients generated by population growth. Little prior work has explored the relation between such continuum models and related single-cell-based models. We present here a simple stochastic individual-based model for the spatial dynamics of multicellular systems whereby cells undergo pressure-driven movement and pressure-dependent proliferation. We show that nonlinear partial differential equations commonly used to model the spatial dynamics of growing cell populations can be formally derived from the branching random walk that underlies our discrete model. Moreover, we carry out a systematic comparison between the individual-based model and its continuum counterparts, both in the case of one single cell population and in the case of multiple cell populations with different biophysical properties. The outcomes of our comparative study demonstrate that the results of computational simulations of the individual-based model faithfully mirror the qualitative and quantitative properties of the solutions to the corresponding nonlinear partial differential equations. Ultimately, these results illustrate how the simple rules governing the dynamics of single cells in our individual-based model can lead to the emergence of complex spatial patterns of population growth observed in continuum models

A hybrid discrete-continuum approach to model Turing pattern formation

Since its introduction in 1952, with a further refinement in 1972 by Gierer and Meinhardt, Turing's (pre-)pattern theory (the chemical basis of morphogenesis) has been widely applied to a number of areas in developmental biology, where evolving cell and tissue structures are naturally observed. The related pattern formation models normally comprise a system of reaction-diffusion equations for interacting chemical species (morphogens), whose heterogeneous distribution in some spatial domain acts as a template for cells to form some kind of pattern or structure through, for example, differentiation or proliferation induced by the chemical pre-pattern. Here we develop a hybrid discrete-continuum modelling framework for the formation of cellular patterns via the Turing mechanism. In this framework, a stochastic individual-based model of cell movement and proliferation is combined with a reaction-diffusion system for the concentrations of some morphogens. As an illustrative example, we focus on a model in which the dynamics of the morphogens are governed by an activator-inhibitor system that gives rise to Turing pre-patterns. The cells then interact with the morphogens in their local area through either of two forms of chemically-dependent cell action: Chemotaxis and chemically-controlled proliferation. We begin by considering such a hybrid model posed on static spatial domains, and then turn to the case of growing domains. In both cases, we formally derive the corresponding deterministic continuum limit and show that that there is an excellent quantitative match between the spatial patterns produced by the stochastic individual-based model and its deterministic continuum counterpart, when sufficiently large numbers of cells are considered. This paper is intended to present a proof of concept for the ideas underlying the modelling framework, with the aim to then apply the related methods to the study of specific patterning and morphogenetic processes in the future

Modelling the Immune Response to Cancer: An Individual-Based Approach Accounting for the Difference in Movement Between Inactive and Activated T Cells

A growing body of experimental evidence indicates that immune cells move in an unrestricted search pattern if they are in the pre-activated state, whilst they tend to stay within a more restricted area upon activation induced by the presence of tumour antigens. This change in movement is not often considered in the existing mathematical models of the interactions between immune cells and cancer cells. With the aim to fill such a gap in the existing literature, in this work we present a spatially structured individual-based model of tumour–immune competition that takes explicitly into account the difference in movement between inactive and activated immune cells. In our model, a Lévy walk is used to capture the movement of inactive immune cells, whereas Brownian motion is used to describe the movement of antigen-activated immune cells. The effects of activation of immune cells, the proliferation of cancer cells and the immune destruction of cancer cells are also modelled. We illustrate the ability of our model to reproduce qualitatively the spatial trajectories of immune cells observed in experimental data of single-cell tracking. Computational simulations of our model further clarify the conditions for the onset of a successful immune action against cancer cells and may suggest possible targets to improve the efficacy of cancer immunotherapy. Overall, our theoretical work highlights the importance of taking into account spatial interactions when modelling the immune response to cancer cells