An individual-based model to explore the impact of psychological stress on immune infiltration into tumour spheroids.

In recentin vitroexperiments on co-culture between breast tumour spheroids and activated immune cells, it was observed that the introduction of the stress hormone cortisol resulted in a decreased immune cell infiltration into the spheroids. Moreover, the presence of cortisol deregulated the normal levels of the pro- and anti-inflammatory cytokines IFN-γand IL-10. We present an individual-based model to explore the interaction dynamics between tumour and immune cells under psychological stress conditions. With our model, we explore the processes underlying the emergence of different levels of immune infiltration, with particular focus on the biological mechanisms regulated by IFN-γand IL-10. The set-up of numerical simulations is defined to mimic the scenarios considered in the experimental study. Similarly to the experimental quantitative analysis, we compute a score that quantifies the level of immune cell infiltration into the tumour. The results of numerical simulations indicate that the motility of immune cells, their capability to infiltrate through tumour cells, their growth rate and the interplay between these cell parameters can affect the level of immune cell infiltration in different ways. Ultimately, numerical simulations of this model support a deeper understanding of the impact of biological stress-induced mechanisms on immune infiltration.

Force Estimation during Cell Migration Using Mathematical Modelling.

Cell migration is essential for physiological, pathological and biomedical processes such as, in embryogenesis, wound healing, immune response, cancer metastasis, tumour invasion and inflammation. In light of this, quantifying mechanical properties during the process of cell migration is of great interest in experimental sciences, yet few theoretical approaches in this direction have been studied. In this work, we propose a theoretical and computational approach based on the optimal control of geometric partial differential equations to estimate cell membrane forces associated with cell polarisation during migration. Specifically, cell membrane forces are inferred or estimated by fitting a mathematical model to a sequence of images, allowing us to capture dynamics of the cell migration. Our approach offers a robust and accurate framework to compute geometric mechanical membrane forces associated with cell polarisation during migration and also yields geometric information of independent interest, we illustrate one such example that involves quantifying cell proliferation levels which are associated with cell division, cell fusion or cell death.

An integrated framework for quantifying immune-tumour interactions in a 3D co-culture model.

Investigational in vitro models that reflect the complexity of the interaction between the immune system and tumours are limited and difficult to establish. Herein, we present a platform to study the tumour-immune interaction using a co-culture between cancer spheroids and activated immune cells. An algorithm was developed for analysis of confocal images of the co-culture to evaluate the following quantitatively; immune cell infiltration, spheroid roundness and spheroid growth. As a proof of concept, the effect of the glucocorticoid stress hormone, cortisol was tested on 66CL4 co-culture model. Results were comparable to 66CL4 syngeneic in vivo mouse model undergoing psychological stress. Furthermore, administration of glucocorticoid receptor antagonists demonstrated the use of this model to determine the effect of treatments on the immune-tumour interplay. In conclusion, we provide a method of quantifying the interaction between the immune system and cancer, which can become a screening tool in immunotherapy design.

A sensor kinase controls turgor-driven plant infection by the rice blast fungus.

The blast fungus Magnaporthe oryzae gains entry to its host plant by means of a specialized pressure-generating infection cell called an appressorium, which physically ruptures the leaf cuticle1,2. Turgor is applied as an enormous invasive force by septin-mediated reorganization of the cytoskeleton and actin-dependent protrusion of a rigid penetration hypha3. However, the molecular mechanisms that regulate the generation of turgor pressure during appressorium-mediated infection of plants remain poorly understood. Here we show that a turgor-sensing histidine-aspartate kinase, Sln1, enables the appressorium to sense when a critical turgor threshold has been reached and thereby facilitates host penetration. We found that the Sln1 sensor localizes to the appressorium pore in a pressure-dependent manner, which is consistent with the predictions of a mathematical model for plant infection. A Δsln1 mutant generates excess intracellular appressorium turgor, produces hyper-melanized non-functional appressoria and does not organize the septins and polarity determinants that are required for leaf infection. Sln1 acts in parallel with the protein kinase C cell-integrity pathway as a regulator of cAMP-dependent signalling by protein kinase A. Pkc1 phosphorylates the NADPH oxidase regulator NoxR and, collectively, these signalling pathways modulate appressorium turgor and trigger the generation of invasive force to cause blast disease.

Cell migration through three-dimensional confining pores: speed accelerations by deformation and recoil of the nucleus.

Directional cell migration in dense three-dimensional (3D) environments critically depends upon shape adaptation and is impeded depending on the size and rigidity of the nucleus. Accordingly, the nucleus is primarily understood as a physical obstacle; however, its pro-migratory functions by stepwise deformation and reshaping remain unclear. Using atomic force spectroscopy, time-lapse fluorescence microscopy and shape change analysis tools, we determined the nuclear size, deformability, morphology and shape change of HT1080 fibrosarcoma cells expressing the Fucci cell cycle indicator or being pre-treated with chromatin-decondensating agent TSA. We show oscillating peak accelerations during migration through 3D collagen matrices and microdevices that occur during shape reversion of deformed nuclei (recoil), and increase with confinement. During G1 cell-cycle phase, nucleus stiffness was increased and yielded further increased speed fluctuations together with sustained cell migration rates in confinement when compared to interphase populations or to periods of intrinsic nuclear softening in the S/G2 cell-cycle phase. Likewise, nuclear softening by pharmacological chromatin decondensation or after lamin A/C depletion reduced peak oscillations in confinement. In conclusion, deformation and recoil of the stiff nucleus contributes to saltatory locomotion in dense tissues. This article is part of a discussion meeting issue 'Forces in cancer: interdisciplinary approaches in tumour mechanobiology'.

Bayesian Parameter Identification for Turing Systems on Stationary and Evolving Domains.

In this study, we apply the Bayesian paradigm for parameter identification to a well-studied semi-linear reaction-diffusion system with activator-depleted reaction kinetics, posed on stationary as well as evolving domains. We provide a mathematically rigorous framework to study the inverse problem of finding the parameters of a reaction-diffusion system given a final spatial pattern. On the stationary domain the parameters are finite-dimensional, but on the evolving domain we consider the problem of identifying the evolution of the domain, i.e. a time-dependent function. Whilst others have considered these inverse problems using optimisation techniques, the Bayesian approach provides a rigorous mathematical framework for incorporating the prior knowledge on uncertainty in the observation and in the parameters themselves, resulting in an approximation of the full probability distribution for the parameters, given the data. Furthermore, using previously established results, we can prove well-posedness results for the inverse problem, using the well-posedness of the forward problem. Although the numerical approximation of the full probability is computationally expensive, parallelised algorithms make the problem solvable using high-performance computing.

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.

A computational framework for particle and whole cell tracking applied to a real biological dataset.

Cell tracking is becoming increasingly important in cell biology as it provides a valuable tool for analysing experimental data and hence furthering our understanding of dynamic cellular phenomena. The advent of high-throughput, high-resolution microscopy and imaging techniques means that a wealth of large data is routinely generated in many laboratories. Due to the sheer magnitude of the data involved manual tracking is often cumbersome and the development of computer algorithms for automated cell tracking is thus highly desirable. In this work, we describe two approaches for automated cell tracking. Firstly, we consider particle tracking. We propose a few segmentation techniques for the detection of cells migrating in a non-uniform background, centroids of the segmented cells are then calculated and linked from frame to frame via a nearest-neighbour approach. Secondly, we consider the problem of whole cell tracking in which one wishes to reconstruct in time whole cell morphologies. Our approach is based on fitting a mathematical model to the experimental imaging data with the goal being that the physics encoded in the model is reflected in the reconstructed data. The resulting mathematical problem involves the optimal control of a phase-field formulation of a geometric evolution law. Efficient approximation of this challenging optimal control problem is achieved via advanced numerical methods for the solution of semilinear parabolic partial differential equations (PDEs) coupled with parallelisation and adaptive resolution techniques. Along with a detailed description of our algorithms, a number of simulation results are reported on. We focus on illustrating the effectivity of our approaches by applying the algorithms to the tracking of migrating cells in a dataset which reflects many of the challenges typically encountered in microscopy data.

A Model for Selection of Eyespots on Butterfly Wings.

UNSOLVED PROBLEM: The development of eyespots on the wing surface of butterflies of the family Nympalidae is one of the most studied examples of biological pattern formation.However, little is known about the mechanism that determines the number and precise locations of eyespots on the wing. Eyespots develop around signaling centers, called foci, that are located equidistant from wing veins along the midline of a wing cell (an area bounded by veins). A fundamental question that remains unsolved is, why a certain wing cell develops an eyespot, while other wing cells do not. KEY IDEA AND MODEL: We illustrate that the key to understanding focus point selection may be in the venation system of the wing disc. Our main hypothesis is that changes in morphogen concentration along the proximal boundary veins of wing cells govern focus point selection. Based on previous studies, we focus on a spatially two-dimensional reaction-diffusion system model posed in the interior of each wing cell that describes the formation of focus points. Using finite element based numerical simulations, we demonstrate that variation in the proximal boundary condition is sufficient to robustly select whether an eyespot focus point forms in otherwise identical wing cells. We also illustrate that this behavior is robust to small perturbations in the parameters and geometry and moderate levels of noise. Hence, we suggest that an anterior-posterior pattern of morphogen concentration along the proximal vein may be the main determinant of the distribution of focus points on the wing surface. In order to complete our model, we propose a two stage reaction-diffusion system model, in which an one-dimensional surface reaction-diffusion system, posed on the proximal vein, generates the morphogen concentrations that act as non-homogeneous Dirichlet (i.e., fixed) boundary conditions for the two-dimensional reaction-diffusion model posed in the wing cells. The two-stage model appears capable of generating focus point distributions observed in nature. RESULT: We therefore conclude that changes in the proximal boundary conditions are sufficient to explain the empirically observed distribution of eyespot focus points on the entire wing surface. The model predicts, subject to experimental verification, that the source strength of the activator at the proximal boundary should be lower in wing cells in which focus points form than in those that lack focus points. The model suggests that the number and locations of eyespot foci on the wing disc could be largely controlled by two kinds of gradients along two different directions, that is, the first one is the gradient in spatially varying parameters such as the reaction rate along the anterior-posterior direction on the proximal boundary of the wing cells, and the second one is the gradient in source values of the activator along the veins in the proximal-distal direction of the wing cell.

Stability analysis and simulations of coupled bulk-surface reaction-diffusion systems.

In this article, we formulate new models for coupled systems of bulk-surface reaction-diffusion equations on stationary volumes. The bulk reaction-diffusion equations are coupled to the surface reaction-diffusion equations through linear Robin-type boundary conditions. We then state and prove the necessary conditions for diffusion-driven instability for the coupled system. Owing to the nature of the coupling between bulk and surface dynamics, we are able to decouple the stability analysis of the bulk and surface dynamics. Under a suitable choice of model parameter values, the bulk reaction-diffusion system can induce patterning on the surface independent of whether the surface reaction-diffusion system produces or not, patterning. On the other hand, the surface reaction-diffusion system cannot generate patterns everywhere in the bulk in the absence of patterning from the bulk reaction-diffusion system. For this case, patterns can be induced only in regions close to the surface membrane. Various numerical experiments are presented to support our theoretical findings. Our most revealing numerical result is that, Robin-type boundary conditions seem to introduce a boundary layer coupling the bulk and surface dynamics.

Parameter identification problems in the modelling of cell motility.

We present a novel parameter identification algorithm for the estimation of parameters in models of cell motility using imaging data of migrating cells. Two alternative formulations of the objective functional that measures the difference between the computed and observed data are proposed and the parameter identification problem is formulated as a minimisation problem of nonlinear least squares type. A Levenberg-Marquardt based optimisation method is applied to the solution of the minimisation problem and the details of the implementation are discussed. A number of numerical experiments are presented which illustrate the robustness of the algorithm to parameter identification in the presence of large deformations and noisy data and parameter identification in three dimensional models of cell motility. An application to experimental data is also presented in which we seek to identify parameters in a model for the monopolar growth of fission yeast cells using experimental imaging data. Our numerical tests allow us to compare the method with the two different formulations of the objective functional and we conclude that the results with both objective functionals seem to agree.

Modelling cell motility and chemotaxis with evolving surface finite elements.

We present a mathematical and a computational framework for the modelling of cell motility. The cell membrane is represented by an evolving surface, with the movement of the cell determined by the interaction of various forces that act normal to the surface. We consider external forces such as those that may arise owing to inhomogeneities in the medium and a pressure that constrains the enclosed volume, as well as internal forces that arise from the reaction of the cells' surface to stretching and bending. We also consider a protrusive force associated with a reaction-diffusion system (RDS) posed on the cell membrane, with cell polarization modelled by this surface RDS. The computational method is based on an evolving surface finite-element method. The general method can account for the large deformations that arise in cell motility and allows the simulation of cell migration in three dimensions. We illustrate applications of the proposed modelling framework and numerical method by reporting on numerical simulations of a model for eukaryotic chemotaxis and a model for the persistent movement of keratocytes in two and three space dimensions. Movies of the simulated cells can be obtained from http://homepages.warwick.ac.uk/∼maskae/CV_Warwick/Chemotaxis.html.

Global existence for semilinear reaction-diffusion systems on evolving domains.

We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly arise in the theory of pattern formation. We present numerical results illustrating our theoretical findings.

Modeling parr-mark pattern formation during the early development of Amago trout.

This paper studies the formation of the large dark patterns, known as parr marks, that form on the Amago trout as it grows from the early larval stages to adulthood. The Amago trout, known as Oncorhynchus masou ishikawa, exhibits stripes during the early stages of development that in turn evolve (through reorientation and peak insertion) to form zigzag spot patterns as the fish grows to adulthood. By considering a standard representation of the Turing model for biological self-organization via interacting and diffusing morphogens, we illustrate that a diffusively driven instability can generate transient patterns consistent with those experimentally observed during the process of parr-mark formation in the early development of the Amago trout. Surface evolution is modeled through an experimentally driven growth function. Our studies conclude that the surface evolution profile, the surface geometry, and the curvature are key factors that play a pivotal role in reaction-diffusion systems in a study motivated by observations of Amago trout parr-mark pattern formation.