Structured models of cell migration incorporating molecular binding processes

Novel Aspects in Pattern Formation Arise from Coupling Turing Reaction-Diffusion and Chemotaxis

Recent experimental studies on primary hair follicle formation and feather bud morphogenesis indicate a coupling between Turing-type diffusion driven instability and chemotactic patterning. Inspired by these findings we develop and analyse a mathematical model that couples chemotaxis to a reaction-diffusion system exhibiting diffusion-driven (Turing) instability. While both systems, reaction-diffusion systems and chemotaxis, can independently generate spatial patterns, we were interested in how the coupling impacts the stability of the system, parameter region for patterning, pattern geometry, as well as the dynamics of pattern formation. We conduct a classical linear stability analysis for different model structures, and confirm our results by numerical analysis of the system. Our results show that the coupling generally increases the robustness of the patterning process by enlarging the pattern region in the parameter space. Concerning time scale and pattern regularity, we find that an increase in the chemosensitivity can speed up the patterning process for parameters inside and outside of the Turing space, but generally reduces spatial regularity of the pattern. Interestingly, our analysis indicates that pattern formation can also occur when neither the Turing nor the chemotaxis system can independently generate pattern. On the other hand, for some parameter settings, the coupling of the two processes can extinguish the pattern formation, rather than reinforce it. These theoretical findings can be used to corroborate the biological findings on morphogenesis and guide future experimental studies. From a mathematical point of view, this work sheds a light on coupling classical pattern formation systems from the parameter space perspective.

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.

Nonlocal and local models for taxis in cell migration: a rigorous limit procedure

A rigorous limit procedure is presented which links nonlocal models involving adhesion or nonlocal chemotaxis to their local counterparts featuring haptotaxis and classical chemotaxis, respectively. It relies on a novel reformulation of the involved nonlocalities in terms of integral operators applied directly to the gradients of signal-dependent quantities. The proposed approach handles both model types in a unified way and extends the previous mathematical framework to settings that allow for general solution-dependent coefficient functions. The previous forms of nonlocal operators are compared with the new ones introduced in this paper and the advantages of the latter are highlighted by concrete examples. Numerical simulations in 1D provide an illustration of some of the theoretical findings.

Mathematical models for cell migration: a non-local perspective

We provide a review of recent advancements in non-local continuous models for migration, mainly from the perspective of its involvement in embryonal development and cancer invasion. Particular emphasis is placed on spatial non-locality occurring in advection terms, used to characterize a cell's motility bias according to its interactions with other cellular and acellular components in its vicinity (e.g. cell-cell and cell-tissue adhesions, non-local chemotaxis), but we also briefly address spatially non-local source terms. Following a short introduction and description of applications, we give a systematic classification of available PDE models with respect to the type of featured non-localities and review some of the mathematical challenges arising from such models, with a focus on analytical aspects. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo heritable, spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of the evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.

Construction of a competing endogenous RNA network using differentially expressed lncRNAs, miRNAs and mRNAs in non‑small cell lung cancer

The competing endogenous RNA (ceRNA) network is crucial for the development and progression of tumors, including non-small cell lung cancer (NSCLC). However, what type of ceRNA network regulates NSCLC has not been clarified. The present study aimed to elucidate the long non-coding RNA (lncRNA)/microRNA (miRNA)/mRNA ceRNA network in NSCLC, particularly for the significance of lncRNAs in NSCLC. NSCLC-specific differentially expressed lncRNAs, miRNAs and mRNAs in the Cancer Genome Atlas (TCGA) were analyzed and their relationship was analyzed by a ceRNA network. Their potential functions of differentially expressed mRNAs were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Furthermore, the expression levels of four selected lncRNAs in TCGA were determined and their associated survival of patients was examined. In addition, the expression profiles of these four lncRNAs in 48 NSCLC specimens and cell lines, their cellular distribution and associated clinical parameters were examined. We successfully constructed a ceRNA network, including 113 lncRNAs, 12 miRNAs and 36 mRNAs differentially expressed between NSCLC and non-tumor tissues. LINC00525, MED4-AS1, STEAP2-AS1 and SYNPR-AS1 lncRNAs were selected and validated for their association with the survival of NSCLC patients. The expression of these lncRNAs was upregulated in 48 NSCLC tissues and was varying in NSCLC cells. While LINC00525 was mainly expressed in the cytoplasm, MED4-AS1 was in both the nucleus and cytoplasm of A549 cells. In addition, the expression of LINC00525 was significantly associated with smoking history (P<0.05); MED4-AS1 was significantly associated with women, poor differentiation and lymph node metastasis (P<0.05); STEAP2-AS1 was significantly associated with women (P<0.01); and SYNPR-AS1 was significantly associated with women and adenocarcinoma (P<0.05). These lncRNAs may be valuable biomarkers for prognosis of NSCLC and the ceRNA network may provide new insights in the pathogenesis of NSCLC.

Age Structure Can Account for Delayed Logistic Proliferation of Scratch Assays

Scratch assays are in vitro methods for studying cell migration. In these experiments, a scratch is made on a cell monolayer and recolonisation of the scratched region is imaged to quantify cell migration rates. Typically, scratch assays are modelled by reaction diffusion equations depicting cell migration by Fickian diffusion and proliferation by a logistic term. In a recent paper (Jin et al. in Bull Math Biol 79(5):1028-1050, 2017), the authors observed experimentally that during the early stage of the recolonisation process, there is a disturbance phase where proliferation is not logistic, and this is followed by a growth phase where proliferation appears to be logistic. The authors did not identify the precise mechanism that causes the disturbance phase but showed that ignoring it can lead to incorrect parameter estimates. The aim of this work is to show that a nonlinear age-structured population model can account for the two phases of proliferation in scratch assays. The model consists of an age-structured cell cycle model of a cell population, coupled with an ordinary differential equation describing the resource concentration dynamics in the substrate. The model assumes a resource-dependent cell cycle threshold age, above which cells are able to proliferate. By studying the dynamics of the full system in terms of the subpopulations of cells that can proliferate and the ones that can not, we are able to find conditions under which the model captures the two-phase behaviour. Through numerical simulations, we are able to show that the interplay between the resource concentration in the substrate and the cell subpopulations dynamics can explain the biphasic dynamics.

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.

INVESTIGATION OF A STRUCTURED FISHER'S EQUATION WITH APPLICATIONS IN BIOCHEMISTRY

Recent biological research has sought to understand how biochemical signaling pathways, such as the mitogen-activated protein kinase (MAPK) family, influence the migration of a population of cells during wound healing. Fisher's Equation has been used extensively to model experimental wound healing assays due to its simple nature and known traveling wave solutions. This partial differential equation with independent variables of time and space cannot account for the effects of biochemical activity on wound healing, however. To this end, we derive a structured Fisher's Equation with independent variables of time, space, and biochemical pathway activity level and prove the existence of a self-similar traveling wave solution to this equation. We exhibit that these methods also apply to a general structured reaction-diffusion equation and a chemotaxis equation. We also consider a more complicated model with different phenotypes based on MAPK activation and numerically investigate how various temporal patterns of biochemical activity can lead to increased and decreased rates of population migration.

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.