Editorial

The last year has been a very busy one for the BMB. We have published 140 articles online, which is over twice the publication rate of the previous year. This reflects the increase in the submission rates and Springer's view that papers should be published online as quickly as possible. In addition, the entire archive of the Bulletin of Mathematical Biology is now available on the Springer website for the journal, digitized back to Vol. 1, No. 1, published in 1939

Periodic solutions of a delayed predator-prey model with stage structure for predator

A periodic time-dependent Lotka-Volterra-type predator-prey model with stage structure for the predator and time delays due to negative feedback and gestation is investigated. Sufficient conditions are derived, respectively, for the existence and global stability of positive periodic solutions to the proposed model

Modeling virus pandemics in a globally connected world a challenge towards a mathematics for living systems

This editorial paper presents the papers published in a special issue devoted to the modeling and simulation of mutating virus pandemics in a globally connected world. The presentation is proposed in three parts. First, motivations and objectives are presented according to the idea that mathematical models should go beyond deterministic population dynamics by considering the multiscale, heterogeneous features of the complex system under consideration. Subsequently, the contents of the papers in this issue are presented referring to the aforementioned complexity features. Finally, a critical analysis of the overall contents of the issue is proposed, with the aim of providing a forward look to research perspectives.PostprintPeer reviewe

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

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

The role of spatial variations of abiotic factors in mediating intratumour phenotypic heterogeneity

We present here a space- and phenotype-structured model of selection dynamics between cancer cells within a solid tumour. In the framework of this model, we combine formal analyses with numerical simulations to investigate in silico the role played by the spatial distribution of abiotic components of the tumour microenvironment in mediating phenotypic selection of cancer cells. Numerical simulations are performed both on the 3D geometry of an in silico multicellular tumour spheroid and on the 3D geometry of an in vivo human hepatic tumour, which was imaged using computerised tomography. The results obtained show that inhomogeneities in the spatial distribution of oxygen, currently observed in solid tumours, can promote the creation of distinct local niches and lead to the selection of different phenotypic variants within the same tumour. This process fosters the emergence of stable phenotypic heterogeneity and supports the presence of hypoxic cells resistant to cytotoxic therapy prior to treatment. Our theoretical results demonstrate the importance of integrating spatial data with ecological principles when evaluating the therapeutic response of solid tumours to cytotoxic therapy

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

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

A Mathematical Framework for Modelling the Metastatic Spread of Cancer

Cancer is a complex disease that starts with mutations of key genes in one cell or a small group of cells at a primary site in the body. If these cancer cells continue to grow successfully and, at some later stage, invade the surrounding tissue and acquire a vascular network, they can spread to distant secondary sites in the body. This process, known as metastatic spread, is responsible for around 90% of deaths from cancer and is one of the so-called hallmarks of cancer. To shed light on the metastatic process, we present a mathematical modelling framework that captures for the first time the interconnected processes of invasion and metastatic spread of individual cancer cells in a spatially explicit manner—a multigrid, hybrid, individual-based approach. This framework accounts for the spatiotemporal evolution of mesenchymal- and epithelial-like cancer cells, membrane-type-1 matrix metalloproteinase (MT1-MMP) and the diffusible matrix metalloproteinase-2 (MMP-2), and for their interactions with the extracellular matrix. Using computational simulations, we demonstrate that our model captures all the key steps of the invasion-metastasis cascade, i.e. invasion by both heterogeneous cancer cell clusters and by single mesenchymal-like cancer cells; intravasation of these clusters and single cells both via active mechanisms mediated by matrix-degrading enzymes (MDEs) and via passive shedding; circulation of cancer cell clusters and single cancer cells in the vasculature with the associated risk of cell death and disaggregation of clusters; extravasation of clusters and single cells; and metastatic growth at distant secondary sites in the body. By faithfully reproducing experimental results, our simulations support the evidence-based hypothesis that the membrane-bound MT1-MMP is the main driver of invasive spread rather than diffusible MDEs such as MMP-2