Traveling Wave Model to Interpret a Wound-Healing Cell Migration Assay for Human Peritoneal Mesothelial Cells

Osteoarthritis as a Systemic Disease Promoted Prostate Cancer In Vivo and In Vitro

Osteoarthritis (OA) is increasing worldwide, and previous work found that OA increases systemic cartilage oligomeric matrix protein (COMP), which has also been implicated in prostate cancer (PCa). As such, we sought to investigate whether OA augments PCa progression. Cellular proliferation and migration of RM1 murine PCa cells treated with interleukin (IL)-1α, COMP, IL-1α + COMP, or conditioned media from cartilage explants treated with IL-1α (representing OA media) and with inhibitors of COMP were assessed. A validated murine model was used for tumor growth and marker expression analysis. Both proliferation and migration were greater in PCa cells treated with OA media compared to controls (p < 0.001), which was not seen with direct application of the stimulants. Migration and proliferation were not negatively affected when OA media was mixed with downstream and COMP inhibitors compared to controls (p > 0.05 for all). Mice with OA developed tumors 100% of the time, whereas mice without OA only 83.4% (p = 0.478). Tumor weight correlated with OA severity (Pearson correlation = 0.813, p = 0.002). Moreover, tumors from mice with OA demonstrated increased Ki-67 expression compared to controls (mean 24.56% vs. 6.91%, p = 0.004) but no difference in CD31, PSMA, or COMP expression (p > 0.05). OA appears to promote prostate cancer in vitro and in vivo.

Making Predictions Using Poorly Identified Mathematical Models

Many commonly used mathematical models in the field of mathematical biology involve challenges of parameter non-identifiability. Practical non-identifiability, where the quality and quantity of data does not provide sufficiently precise parameter estimates is often encountered, even with relatively simple models. In particular, the situation where some parameters are identifiable and others are not is often encountered. In this work we apply a recent likelihood-based workflow, called Profile-Wise Analysis (PWA), to non-identifiable models for the first time. The PWA workflow addresses identifiability, parameter estimation, and prediction in a unified framework that is simple to implement and interpret. Previous implementations of the workflow have dealt with idealised identifiable problems only. In this study we illustrate how the PWA workflow can be applied to both structurally non-identifiable and practically non-identifiable models in the context of simple population growth models. Dealing with simple mathematical models allows us to present the PWA workflow in a didactic, self-contained document that can be studied together with relatively straightforward Julia code provided on GitHub . Working with simple mathematical models allows the PWA workflow prediction intervals to be compared with gold standard full likelihood prediction intervals. Together, our examples illustrate how the PWA workflow provides us with a systematic way of dealing with non-identifiability, especially compared to other approaches, such as seeking ad hoc parameter combinations, or simply setting parameter values to some arbitrary default value. Importantly, we show that the PWA workflow provides insight into the commonly-encountered situation where some parameters are identifiable and others are not, allowing us to explore how uncertainty in some parameters, and combinations of parameters, regardless of their identifiability status, influences model predictions in a way that is insightful and interpretable.

Discrete and continuous mathematical models of sharp-fronted collective cell migration and invasion

Mathematical models describing the spatial spreading and invasion of populations of biological cells are often developed in a continuum modelling framework using reaction-diffusion equations. While continuum models based on linear diffusion are routinely employed and known to capture key experimental observations, linear diffusion fails to predict well-defined sharp fronts that are often observed experimentally. This observation has motivated the use of nonlinear degenerate diffusion; however, these nonlinear models and the associated parameters lack a clear biological motivation and interpretation. Here, we take a different approach by developing a stochastic discrete lattice-based model incorporating biologically inspired mechanisms and then deriving the reaction-diffusion continuum limit. Inspired by experimental observations, agents in the simulation deposit extracellular material, which we call a substrate, locally onto the lattice, and the motility of agents is taken to be proportional to the substrate density. Discrete simulations that mimic a two-dimensional circular barrier assay illustrate how the discrete model supports both smooth and sharp-fronted density profiles depending on the rate of substrate deposition. Coarse-graining the discrete model leads to a novel partial differential equation (PDE) model whose solution accurately approximates averaged data from the discrete model. The new discrete model and PDE approximation provide a simple, biologically motivated framework for modelling the spreading, growth and invasion of cell populations with well-defined sharp fronts. Open-source Julia code to replicate all results in this work is available on GitHub.

Bayesian inverse problem for a fractional diffusion model of cell migration

In the present work, both direct and inverse problems are considered for a Fisher-type fractional diffusion equation, which is proposed to describe the phenomenon of cell migration. For the direct problem, a solution is given via the Fourier method and the Laplace transform. On the other hand, we solved the inverse problem from a Bayesian statistical framework using a set of data that are the result of a cell migration experiment on a wound closure assay. We estimated the parameters of the mathematical model via Markov Chain Monte Carlo methods.

Modelling count data with partial differential equation models in biology

Partial differential equation (PDE) models are often used to study biological phenomena involving movement-birth-death processes, including ecological population dynamics and the invasion of populations of biological cells. Count data, by definition, is non-negative, and count data relating to biological populations is often bounded above by some carrying capacity that arises through biological competition for space or nutrients. Parameter estimation, parameter identifiability, and making model predictions usually involves working with a measurement error model that explicitly relating experimental measurements with the solution of a mathematical model. In many biological applications, a typical approach is to assume the data are normally distributed about the solution of the mathematical model. Despite the widespread use of the standard additive Gaussian measurement error model, the assumptions inherent in this approach are rarely explicitly considered or compared with other options. Here, we interpret scratch assay data, involving migration, proliferation and delays in a population of cancer cells using a reaction-diffusion PDE model. We consider relating experimental measurements to the PDE solution using a standard additive Gaussian measurement error model alongside a comparison to a more biologically realistic binomial measurement error model. While estimates of model parameters are relatively insensitive to the choice of measurement error model, model predictions for data realisations are very sensitive. The standard additive Gaussian measurement error model leads to biologically inconsistent predictions, such as negative counts and counts that exceed the carrying capacity across a relatively large spatial region within the experiment. Furthermore, the standard additive Gaussian measurement error model requires estimating an additional parameter compared to the binomial measurement error model. In contrast, the binomial measurement error model leads to biologically plausible predictions and is simpler to implement. We provide open source Julia software on GitHub to replicate all calculations in this work, and we explain how to generalise our approach to deal with coupled PDE models with several dependent variables through a multinomial measurement error model, as well as pointing out other potential generalisations by linking our work with established practices in the field of generalised linear models.

Modeling the extracellular matrix in cell migration and morphogenesis: a guide for the curious biologist

The extracellular matrix (ECM) is a highly complex structure through which biochemical and mechanical signals are transmitted. In processes of cell migration, the ECM also acts as a scaffold, providing structural support to cells as well as points of potential attachment. Although the ECM is a well-studied structure, its role in many biological processes remains difficult to investigate comprehensively due to its complexity and structural variation within an organism. In tandem with experiments, mathematical models are helpful in refining and testing hypotheses, generating predictions, and exploring conditions outside the scope of experiments. Such models can be combined and calibrated with in vivo and in vitro data to identify critical cell-ECM interactions that drive developmental and homeostatic processes, or the progression of diseases. In this review, we focus on mathematical and computational models of the ECM in processes such as cell migration including cancer metastasis, and in tissue structure and morphogenesis. By highlighting the predictive power of these models, we aim to help bridge the gap between experimental and computational approaches to studying the ECM and to provide guidance on selecting an appropriate model framework to complement corresponding experimental studies.

Understanding the ideal wound healing mechanistic behavior using in silico modelling perspectives: A review

Complexity of the entire body precludes an accurate assessment of the specific contributions of tissues or cells during the healing process, which might be expensive and time consuming. Because of this, controlling the wound's size, depth, and dimensions may be challenging, and there is not yet an efficient and reliable chronic wound model representation. Furthermore, given the inherent challenges associated with conducting non-invasive in vivo investigations, it becomes peremptory to explore alternative methodologies for studying wound healing. In this context, biologically-realistic mathematical and computational models emerge as a valuable framework that can effectively address this need. Therefore, it might improve our approach to understanding the process at its core. This article will examines all facets of wound healing, including the kinds, pathways, and most current developments in wound treatment worldwide, particularly in silico modelling utilizing both mathematical and structure-based modelling techniques. It may be helpful to identify the crucial traits through the feedback loop of computer models and experimental investigations in order to build innovative therapies to cure wounds. Hence the effectiveness of personalised medicine and more targeted therapy in the healing of wounds may be enhanced by this interdisciplinary expertise.

Modelling the Effect of Matrix Metalloproteinases in Dermal Wound Healing

With over 2 million people in the UK suffering from chronic wounds, understanding the biochemistry and pharmacology that underpins these wounds and wound healing is of high importance. Chronic wounds are characterised by high levels of matrix metalloproteinases (MMPs), which are necessary for the modification of healthy tissue in the healing process. Overexposure of MMPs, however, adversely affects healing of the wound by causing further destruction of the surrounding extracellular matrix. In this work, we propose a mathematical model that focuses on the interaction of MMPs with dermal cells using a system of partial differential equations. Using biologically realistic parameter values, this model gives rise to travelling waves corresponding to a front of healthy cells invading a wound. From the arising travelling wave analysis, we observe that deregulated apoptosis results in the emergence of chronic wounds, characterised by elevated MMP concentrations. We also observe hysteresis effects when both the apoptotic rate and MMP production rate are varied, providing further insight into the management (and potential reversal) of chronic wounds.

Quantifying tissue growth, shape and collision via continuum models and Bayesian inference

Although tissues are usually studied in isolation, this situation rarely occurs in biology, as cells, tissues and organs coexist and interact across scales to determine both shape and function. Here, we take a quantitative approach combining data from recent experiments, mathematical modelling and Bayesian parameter inference, to describe the self-assembly of multiple epithelial sheets by growth and collision. We use two simple and well-studied continuum models, where cells move either randomly or following population pressure gradients. After suitable calibration, both models prove to be practically identifiable, and can reproduce the main features of single tissue expansions. However, our findings reveal that whenever tissue-tissue interactions become relevant, the random motion assumption can lead to unrealistic behaviour. Under this setting, a model accounting for population pressure from different cell populations is more appropriate and shows a better agreement with experimental measurements. Finally, we discuss how tissue shape and pressure affect multi-tissue collisions. Our work thus provides a systematic approach to quantify and predict complex tissue configurations with applications in the design of tissue composites and more generally in tissue engineering.

A mathematical model of wound healing in bovine corneal endothelium

We developed a mathematical model to describe healing processes in bovine corneal endothelial (BCE) cells in culture, triggered by mechanical wounds with parallel edges. Previous findings from our laboratory show that, in these cases, BCE monolayers exhibit an approximately constant healing velocity. Also, that caspase-dependent apoptosis occurs, with the fraction of apoptotic cells increasing with the distance traveled by the healing edge. In addition, in this study we report the novel findings that, for wound scratch assays performed preserving the basal extracellular matrix: i) the healing cells increase their en face surface area in a characteristic fashion, and ii) the average length of the segments of the cell columns actively participating in the healing process increases linearly with time. These latter observations preclude the utilization of standard traveling wave formalisms to model wound healing in BCE cells. Instead, we developed and studied a simple phenomenological model based on a plausible formula for the spreading dynamics of the individual healing cells, that incorporates original evidence about the process in BCE cells. The model can be simulated to: i) obtain an approximately constant healing velocity; ii) reproduce the profile of the healing cell areas, and iii) obtain approximately linear time dependences of the mean cell area and average length of the front active segments per column. In view of its accuracy to account for the experimental observations, the model can also be acceptably employed to quantify the appearance of apoptotic cells during BCE wound healing. The strategy utilized here could offer a novel formal framework to represent modifications undergone by some epithelial cell lines during wound healing.

Geometric analysis enables biological insight from complex non-identifiable models using simple surrogates

An enduring challenge in computational biology is to balance data quality and quantity with model complexity. Tools such as identifiability analysis and information criterion have been developed to harmonise this juxtaposition, yet cannot always resolve the mismatch between available data and the granularity required in mathematical models to answer important biological questions. Often, it is only simple phenomenological models, such as the logistic and Gompertz growth models, that are identifiable from standard experimental measurements. To draw insights from complex, non-identifiable models that incorporate key biological mechanisms of interest, we study the geometry of a map in parameter space from the complex model to a simple, identifiable, surrogate model. By studying how non-identifiable parameters in the complex model quantitatively relate to identifiable parameters in surrogate, we introduce and exploit a layer of interpretation between the set of non-identifiable parameters and the goodness-of-fit metric or likelihood studied in typical identifiability analysis. We demonstrate our approach by analysing a hierarchy of mathematical models for multicellular tumour spheroid growth experiments. Typical data from tumour spheroid experiments are limited and noisy, and corresponding mathematical models are very often made arbitrarily complex. Our geometric approach is able to predict non-identifiabilities, classify non-identifiable parameter spaces into identifiable parameter combinations that relate to features in the data characterised by parameters in a surrogate model, and overall provide additional biological insight from complex non-identifiable models.

Proper Orthogonal Decomposition Analysis Reveals Cell Migration Directionality During Wound Healing

A proper orthogonal decomposition (POD) order reduction method was implemented to reduce the full high dimensional dynamical system associated with a wound healing cell migration assay to a low-dimensional approximation that identified the prevailing cell trajectories. The POD analysis generated POD modes that were representative of the prevalent cell trajectories. The shapes of the POD modes depended on the location of the cells with respect to the wound and exposure to disturbed (DF) or undisturbed (UF) fluid flow where the net flow was in the antegrade direction with a retrograde component or fully antegrade, respectively. For DF and UF, the POD modes of the downstream cells identified trajectories that moved upstream against the flow, while upstream POD modes exhibited sideways cell migrations. In the absence of flow, no major shape differences were observed in the POD modes on either side of the wound. The POD modes also served to reconstruct the cell migration of individual cells. With as few as three modes, the predominant cell migration trajectories were reconstructed, while the level of accuracy increased with the inclusion of more POD modes. The POD order reduction method successfully identified the predominant cell migratory trajectories under static and varying pulsatile fluid flow conditions serving as a first step in the development of artificial intelligence models of cell migration in disease and development.

Revisiting Fisher-KPP model to interpret the spatial spreading of invasive cell population in biology

In this paper, the homotopy analysis method, a powerful analytical technique, is applied to obtain analytical solutions to the Fisher-KPP equation in studying the spatial spreading of invasive species in ecology and to extract the nature of the spatial spreading of invasive cell populations in biology. The effect of the proliferation rate of the model of interest on the entire population is studied. It is observed that the invasive cell or the invasive population is decreased within a short time with the minimum proliferation rate. The homotopy analysis method is found to be superior to other analytical methods, namely the Adomian decomposition method, the homotopy perturbation method, etc. because of containing an auxiliary parameter, which provides us with a convenient way to adjust and control the region of convergence of the series solution. Graphical representation of the approximate series solutions obtained by the homotopy analysis method, the Adomian decomposition method, and the Homotopy perturbation method is illustrated, which shows the superiority of the homotopy analysis method. The method is examined on several examples, which reveal the ingenuousness and the effectiveness of the method of interest.

•

Closed-form solutions are obtained for the Fisher-KPP equation through the Homotopy analysis method.

•

The effect of the proliferation rate of the model of interest on the entire population is studied.

•

The invasive cell or the invasive population decreases in short time with the minimum proliferation rate.

•

The Homotopy analysis method is found superior over other analytical methods.

Non-vanishing sharp-fronted travelling wave solutions of the Fisher-Kolmogorov model

The Fisher-Kolmogorov-Petrovsky-Piskunov (KPP) model, and generalizations thereof, involves simple reaction-diffusion equations for biological invasion that assume individuals in the population undergo linear diffusion with diffusivity $D$, and logistic proliferation with rate $\lambda $. For the Fisher-KPP model, biologically relevant initial conditions lead to long-time travelling wave solutions that move with speed $c=2\sqrt {\lambda D}$. Despite these attractive features, there are several biological limitations of travelling wave solutions of the Fisher-KPP model. First, these travelling wave solutions do not predict a well-defined invasion front. Second, biologically relevant initial conditions lead to travelling waves that move with speed $c=2\sqrt {\lambda D}> 0$. This means that, for biologically relevant initial data, the Fisher-KPP model cannot be used to study invasion with $c \ne 2\sqrt {\lambda D}$, or retreating travelling waves with $c < 0$. Here, we reformulate the Fisher-KPP model as a moving boundary problem and show that this reformulated model alleviates the key limitations of the Fisher-KPP model. Travelling wave solutions of the moving boundary problem predict a well-defined front that can propagate with any wave speed, $-\infty < c < \infty $. Here, we establish these results using a combination of high-accuracy numerical simulations of the time-dependent partial differential equation, phase plane analysis and perturbation methods. All software required to replicate this work is available on GitHub.

A Model for Cell Proliferation in a Developing Organism

In mathematical biology, there is a great deal of interest in producing continuum models by scaling discrete agent-based models governed by local stochastic rules. We discuss a particular example of this approach: a model for the proliferation of neural crest cells that can help us understand the development of Hirschprung's disease, a potentially-fatal condition in which the enteric nervous system of a new-born child does not extend all the way through the intestine and colon. Our starting point is a discrete-state, continuous-time Markov chain model proposed by Hywood et al. (2013a) for the location of the neural crest cells that make up the enteric nervous system. Hywood et al. (2013a) scaled their model to derive an approximate second order partial differential equation describing how the limiting expected number of neural crest cells evolve in space and time. In contrast, we exploit the relationship between the above-mentioned Markov chain model and the well-known Yule-Furry process to derive the exact form of the scaled version of the process. Furthermore, we provide expressions for other features of the domain agent occupancy process, such as the variance of the marginal occupancy at a particular site, the distribution of the number of agents that are yet to reach a given site and a stochastic description of the process itself.

Reliable and efficient parameter estimation using approximate continuum limit descriptions of stochastic models

Stochastic individual-based mathematical models are attractive for modelling biological phenomena because they naturally capture the stochasticity and variability that is often evident in biological data. Such models also allow us to track the motion of individuals within the population of interest. Unfortunately, capturing this microscopic detail means that simulation and parameter inference can become computationally expensive. One approach for overcoming this computational limitation is to coarse-grain the stochastic model to provide an approximate continuum model that can be solved using far less computational effort. However, coarse-grained continuum models can be biased or inaccurate, particularly for certain parameter regimes. In this work, we combine stochastic and continuum mathematical models in the context of lattice-based models of two-dimensional cell biology experiments by demonstrating how to simulate two commonly used experiments: cell proliferation assays and barrier assays. Our approach involves building a simple statistical model of the discrepancy between the expensive stochastic model and the associated computationally inexpensive coarse-grained continuum model. We form this statistical model based on a limited number of expensive stochastic model evaluations at design points sampled from a user-chosen distribution, corresponding to a computer experiment design problem. With straightforward design point selection schemes, we show that using the statistical model of the discrepancy in tandem with the computationally inexpensive continuum model allows us to carry out prediction and inference while correcting for biases and inaccuracies due to the continuum approximation. We demonstrate this approach by simulating a proliferation assay, where the continuum limit model is the well-known logistic ordinary differential equation, as well as a barrier assay where the continuum limit model is closely related to the well-known Fisher-KPP partial differential equation. We construct an approximate likelihood function for parameter inference, both with and without discrepancy correction terms. Using maximum likelihood estimation, we provide point estimates of the unknown parameters, and use the profile likelihood to characterise the uncertainty in these estimates and form approximate confidence intervals. For the range of inference problems considered, working with the continuum limit model alone leads to biased parameter estimation and confidence intervals with poor coverage. In contrast, incorporating correction terms arising from the statistical model of the model discrepancy allows us to recover the parameters accurately with minimal computational overhead. The main tradeoff is that the associated confidence intervals are typically broader, reflecting the additional uncertainty introduced by the approximation process. All algorithms required to replicate the results in this work are written in the open source Julia language and are available at GitHub.

Interpreting how nonlinear diffusion affects the fate of bistable populations using a discrete modelling framework

Understanding whether a population will survive or become extinct is a central question in population biology. One way of exploring this question is to study population dynamics using reaction–diffusion equations, where migration is usually represented as a linear diffusion term, and birth–death is represented with a nonlinear source term. While linear diffusion is most commonly employed to study migration, there are several limitations of this approach, such as the inability of linear diffusion-based models to predict a well-defined population front. One way to overcome this is to generalize the constant diffusivity, D, to a nonlinear diffusivity function D(C), where C>0 is the population density. While the choice of D(C) affects long-term survival or extinction of a bistable population, working solely in a continuum framework makes it difficult to understand how the choice of D(C) affects survival or extinction. We address this question by working with a discrete simulation model that is easy to interpret. This approach provides clear insight into how the choice of D(C) either encourages or suppresses population extinction relative to the classical linear diffusion model.

A Continuum Mathematical Model of Substrate-Mediated Tissue Growth

We consider a continuum mathematical model of biological tissue formation inspired by recent experiments describing thin tissue growth in 3D-printed bioscaffolds. The continuum model, which we call the substrate model, involves a partial differential equation describing the density of tissue, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{u}}(\hat{{\mathbf {x}}},{\hat{t}})$$\end{document}u^(x^,t^) that is coupled to the concentration of an immobile extracellular substrate, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}(\hat{{\mathbf {x}}},{\hat{t}})$$\end{document}s^(x^,t^). Cell migration is modelled with a nonlinear diffusion term, where the diffusive flux is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}$$\end{document}s^, while a logistic growth term models cell proliferation. The extracellular substrate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}$$\end{document}s^ is produced by cells and undergoes linear decay. Preliminary numerical simulations show that this mathematical model is able to recapitulate key features of recent tissue growth experiments, including the formation of sharp fronts. To provide a deeper understanding of the model we analyse travelling wave solutions of the substrate model, showing that the model supports both sharp-fronted travelling wave solutions that move with a minimum wave speed, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c = c_{\mathrm{min}}$$\end{document}c=cmin, as well as smooth-fronted travelling wave solutions that move with a faster travelling wave speed, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c > c_{\mathrm{min}}$$\end{document}c>cmin. We provide a geometric interpretation that explains the difference between smooth and sharp-fronted travelling wave solutions that is based on a slow manifold reduction of the desingularised three-dimensional phase space. In addition, we also develop and test a series of useful approximations that describe the shape of the travelling wave solutions in various limits. These approximations apply to both the sharp-fronted and smooth-fronted travelling wave solutions. Software to implement all calculations is available at GitHub.

Parameter identifiability and model selection for sigmoid population growth models

Sigmoid growth models, such as the logistic, Gompertz and Richards' models, are widely used to study population dynamics ranging from microscopic populations of cancer cells, to continental-scale human populations. Fundamental questions about model selection and parameter estimation are critical if these models are to be used to make practical inferences. However, the question of parameter identifiability - whether a data set contains sufficient information to give unique or sufficiently precise parameter estimates - is often overlooked. We use a profile-likelihood approach to explore practical parameter identifiability using data describing the re-growth of hard coral. With this approach, we explore the relationship between parameter identifiability and model misspecification, finding that the logistic growth model does not suffer identifiability issues for the type of data we consider whereas the Gompertz and Richards' models encounter practical non-identifiability issues. This analysis of parameter identifiability and model selection is important because different growth models are in biological modelling without necessarily considering whether parameters are identifiable. Standard practices that do not consider parameter identifiability can lead to unreliable or imprecise parameter estimates and potentially misleading mechanistic interpretations. For example, using the Gompertz model, the estimate of the time scale of coral re-growth is 625 days when we estimate the initial density from the data, whereas it is 1429 days using a more standard approach where variability in the initial density is ignored. While tools developed here focus on three standard sigmoid growth models only, our theoretical developments are applicable to any sigmoid growth model and any appropriate data set. MATLAB implementations of all software are available on GitHub.

Extinction of Bistable Populations is Affected by the Shape of their Initial Spatial Distribution

The question of whether biological populations survive or are eventually driven to extinction has long been examined using mathematical models. In this work, we study population survival or extinction using a stochastic, discrete lattice-based random walk model where individuals undergo movement, birth and death events. The discrete model is defined on a two-dimensional hexagonal lattice with periodic boundary conditions. A key feature of the discrete model is that crowding effects are introduced by specifying two different crowding functions that govern how local agent density influences movement events and birth/death events. The continuum limit description of the discrete model is a nonlinear reaction-diffusion equation, and we focus on crowding functions that lead to linear diffusion and a bistable source term that is often associated with the strong Allee effect. Using both the discrete and continuum modelling tools, we explore the complicated relationship between the long-term survival or extinction of the population and the initial spatial arrangement of the population. In particular, we study different spatial arrangements of initial distributions: (i) a well-mixed initial distribution where the initial density is independent of position in the domain; (ii) a vertical strip initial distribution where the initial density is independent of vertical position in the domain; and, (iii) several forms of two-dimensional initial distributions where the initial population is distributed in regions with different shapes. Our results indicate that the shape of the initial spatial distribution of the population affects extinction of bistable populations. All software required to solve the discrete and continuum models used in this work are available on GitHub .

Physics-informed learning of governing equations from scarce data

Harnessing data to discover the underlying governing laws or equations that describe the behavior of complex physical systems can significantly advance our modeling, simulation and understanding of such systems in various science and engineering disciplines. This work introduces a novel approach called physics-informed neural network with sparse regression to discover governing partial differential equations from scarce and noisy data for nonlinear spatiotemporal systems. In particular, this discovery approach seamlessly integrates the strengths of deep neural networks for rich representation learning, physics embedding, automatic differentiation and sparse regression to approximate the solution of system variables, compute essential derivatives, as well as identify the key derivative terms and parameters that form the structure and explicit expression of the equations. The efficacy and robustness of this method are demonstrated, both numerically and experimentally, on discovering a variety of partial differential equation systems with different levels of data scarcity and noise accounting for different initial/boundary conditions. The resulting computational framework shows the potential for closed-form model discovery in practical applications where large and accurate datasets are intractable to capture.

Combining experiments and in silico modeling to infer the role of adhesion and proliferation on the collective dynamics of cells

The collective dynamics of cells on surfaces and interfaces poses technological and theoretical challenges in the study of morphogenesis, tissue engineering, and cancer. Different mechanisms are at play, including, cell–cell adhesion, cell motility, and proliferation. However, the relative importance of each one is elusive. Here, experiments with a culture of glioblastoma multiforme cells on a substrate are combined with in silico modeling to infer the rate of each mechanism. By parametrizing these rates, the time-dependence of the spatial correlation observed experimentally is reproduced. The obtained results suggest a reduction in cell–cell adhesion with the density of cells. The reason for such reduction and possible implications for the collective dynamics of cancer cells are discussed.

Estimating parameters of a stochastic cell invasion model with fluorescent cell cycle labelling using approximate Bayesian computation

We develop a parameter estimation method based on approximate Bayesian computation (ABC) for a stochastic cell invasion model using fluorescent cell cycle labelling with proliferation, migration and crowding effects. Previously, inference has been performed on a deterministic version of the model fitted to cell density data, and not all parameters were identifiable. Considering the stochastic model allows us to harness more features of experimental data, including cell trajectories and cell count data, which we show overcomes the parameter identifiability problem. We demonstrate that, while difficult to collect, cell trajectory data can provide more information about the parameters of the cell invasion model. To handle the intractability of the likelihood function of the stochastic model, we use an efficient ABC algorithm based on sequential Monte Carlo. Rcpp and MATLAB implementations of the simulation model and ABC algorithm used in this study are available at https://github.com/michaelcarr-stats/FUCCI.

Model-based data analysis of tissue growth in thin 3D printed scaffolds

Tissue growth in three-dimensional (3D) printed scaffolds enables exploration and control of cell behaviour in more biologically realistic geometries than that allowed by traditional 2D cell culture. Cell proliferation and migration in these experiments have yet to be explicitly characterised, limiting the ability of experimentalists to determine the effects of various experimental conditions, such as scaffold geometry, on cell behaviour. We consider tissue growth by osteoblastic cells in melt electro-written scaffolds that comprise thin square pores with sizes that were deliberately increased between experiments. We collect highly detailed temporal measurements of the average cell density, tissue coverage, and tissue geometry. To quantify tissue growth in terms of the underlying cell proliferation and migration processes, we introduce and calibrate a mechanistic mathematical model based on the Porous-Fisher reaction-diffusion equation. Parameter estimates and uncertainty quantification through profile likelihood analysis reveal consistency in the rate of cell proliferation and steady-state cell density between pore sizes. This analysis also serves as an important model verification tool: while the use of reaction-diffusion models in biology is widespread, the appropriateness of these models to describe tissue growth in 3D scaffolds has yet to be explored. We find that the Porous-Fisher model is able to capture features relating to the cell density and tissue coverage, but is not able to capture geometric features relating to the circularity of the tissue interface. Our analysis identifies two distinct stages of tissue growth, suggests several areas for model refinement, and provides guidance for future experimental work that explores tissue growth in 3D printed scaffolds.

Mathematical Modeling Can Advance Wound Healing Research

Significance: For over 30 years, there has been sustained interest in the development of mathematical models for investigating the complex mechanisms underlying each stage of the wound healing process. Despite the immense associated challenges, such models have helped usher in a paradigm shift in wound healing research. Recent Advances: In this article, we review contributions in the field that span epidermal, dermal, and corneal wound healing, and treatments of nonhealing wounds. The recent influence of mathematical models on biological experiments is detailed, with a focus on wound healing assays and fibroblast-populated collagen lattices. Critical Issues: We provide an overview of the field of mathematical modeling of wound healing, highlighting key advances made in recent decades, and discuss how such models have contributed to the development of improved treatment strategies and/or an enhanced understanding of the tightly regulated steps that comprise the healing process. Future Directions: We detail some of the open problems in the field that could be addressed through a combination of theoretical and/or experimental approaches. To move the field forward, we need to have a common language between scientists to facilitate cross-collaboration, which we hope this review can support by highlighting progress to date.

The role of mechanical interactions in EMT

The detachment of cells from the boundary of an epithelial tissue and the subsequent invasion of these cells into surrounding tissues is important for cancer development and wound healing, and is strongly associated with the epithelial-mesenchymal transition (EMT). Chemical signals, such as TGF-β, produced by surrounding tissue can be uptaken by cells and induce EMT. In this work, we present a novel cell-based discrete mathematical model of mechanical cellular relaxation, cell proliferation, and cell detachment driven by chemically-dependent EMT in an epithelial tissue. A continuum description of the model is then derived in the form of a novel nonlinear free boundary problem. Using the discrete and continuum models we explore how the coupling of chemical transport and mechanical interactions influences EMT, and postulate how this could be used to help control EMT in pathological situations.

Mathematical Model of Tumour Spheroid Experiments with Real-Time Cell Cycle Imaging

Three-dimensional (3D) in vitro tumour spheroid experiments are an important tool for studying cancer progression and potential cancer drug therapies. Standard experiments involve growing and imaging spheroids to explore how different conditions lead to different rates of spheroid growth. These kinds of experiments, however, do not reveal any information about the spatial distribution of the cell cycle within the expanding spheroid. Since 2008, a new experimental technology called fluorescent ubiquitination-based cell cycle indicator (FUCCI) has enabled real-time in situ visualisation of the cell cycle progression. Observations of 3D tumour spheroids with FUCCI labelling reveal significant intratumoural structure, as the cell cycle status can vary with location. Although many mathematical models of tumour spheroid growth have been developed, none of the existing mathematical models are designed to interpret experimental observations with FUCCI labelling. In this work, we adapt the mathematical framework originally proposed by Ward and King (Math Med Biol 14:39-69, 1997. https://doi.org/10.1093/imammb/14.1.39 ) to produce a new mathematical model of FUCCI-labelled tumour spheroid growth. The mathematical model treats the spheroid as being composed of three subpopulations: (i) living cells in G1 phase that fluoresce red; (ii) living cells in S/G2/M phase that fluoresce green; and (iii) dead cells that are not fluorescent. We assume that the rates at which cells pass through different phases of the cell cycle, and the rate of cell death, depend upon the local oxygen concentration. Parameterising the new mathematical model using experimental measurements of cell cycle transition times, we show that the model can qualitatively capture important experimental observations that cannot be addressed using previous mathematical models. Further, we show that the mathematical model can be used to qualitatively mimic the action of anti-mitotic drugs applied to the spheroid. All software programs required to solve the nonlinear moving boundary problem associated with the new mathematical model are available on GitHub. at https://github.com/wang-jin-mathbio/Jin2021.

A novel mathematical model of heterogeneous cell proliferation

We present a novel mathematical model of heterogeneous cell proliferation where the total population consists of a subpopulation of slow-proliferating cells and a subpopulation of fast-proliferating cells. The model incorporates two cellular processes, asymmetric cell division and induced switching between proliferative states, which are important determinants for the heterogeneity of a cell population. As motivation for our model we provide experimental data that illustrate the induced-switching process. Our model consists of a system of two coupled delay differential equations with distributed time delays and the cell densities as functions of time. The distributed delays are bounded and allow for the choice of delay kernel. We analyse the model and prove the nonnegativity and boundedness of solutions, the existence and uniqueness of solutions, and the local stability characteristics of the equilibrium points. We find that the parameters for induced switching are bifurcation parameters and therefore determine the long-term behaviour of the model. Numerical simulations illustrate and support the theoretical findings, and demonstrate the primary importance of transient dynamics for understanding the evolution of many experimental cell populations.

Invading and Receding Sharp-Fronted Travelling Waves

Biological invasion, whereby populations of motile and proliferative individuals lead to moving fronts that invade vacant regions, is routinely studied using partial differential equation models based upon the classical Fisher-KPP equation. While the Fisher-KPP model and extensions have been successfully used to model a range of invasive phenomena, including ecological and cellular invasion, an often-overlooked limitation of the Fisher-KPP model is that it cannot be used to model biological recession where the spatial extent of the population decreases with time. In this work, we study the Fisher-Stefan model, which is a generalisation of the Fisher-KPP model obtained by reformulating the Fisher-KPP model as a moving boundary problem. The nondimensional Fisher-Stefan model involves just one parameter, [Formula: see text], which relates the shape of the density front at the moving boundary to the speed of the associated travelling wave, c. Using numerical simulation, phase plane and perturbation analysis, we construct approximate solutions of the Fisher-Stefan model for both slowly invading and receding travelling waves, as well as for rapidly receding travelling waves. These approximations allow us to determine the relationship between c and [Formula: see text] so that commonly reported experimental estimates of c can be used to provide estimates of the unknown parameter [Formula: see text]. Interestingly, when we reinterpret the Fisher-KPP model as a moving boundary problem, many overlooked features of the classical Fisher-KPP phase plane take on a new interpretation since travelling waves solutions with [Formula: see text] are normally disregarded. This means that our analysis of the Fisher-Stefan model has both practical value and an inherent mathematical value.

Vibration enhanced cell growth induced by surface acoustic waves as in vitro wound-healing model

We report on in vitro wound-healing and cell-growth studies under the influence of radio-frequency (rf) cell stimuli. These stimuli are supplied either by piezoactive surface acoustic waves (SAWs) or by microelectrode-generated electric fields, both at frequencies around 100 MHz. Employing live-cell imaging, we studied the time- and power-dependent healing of artificial wounds on a piezoelectric chip for different cell lines. If the cell stimulation is mediated by piezomechanical SAWs, we observe a pronounced, significant maximum of the cell-growth rate at a specific SAW amplitude, resulting in an increase of the wound-healing speed of up to 135 ± 85% as compared to an internal reference. In contrast, cells being stimulated only by electrical fields of the same magnitude as the ones exposed to SAWs exhibit no significant effect. In this study, we investigate this effect for different wavelengths, amplitude modulation of the applied electrical rf signal, and different wave modes. Furthermore, to obtain insight into the biological response to the stimulus, we also determined both the cell-proliferation rate and the cellular stress levels. While the proliferation rate is significantly increased for a wide power range, cell stress remains low and within the normal range. Our findings demonstrate that SAW-based vibrational cell stimulation bears the potential for an alternative method to conventional ultrasound treatment, overcoming some of its limitations.

Guidelines for establishing a 3-D printing biofabrication laboratory

Advanced manufacturing and 3D printing are transformative technologies currently undergoing rapid adoption in healthcare, a traditionally non-manufacturing sector. Recent development in this field, largely enabled by merging different disciplines, has led to important clinical applications from anatomical models to regenerative bioscaffolding and devices. Although much research to-date has focussed on materials, designs, processes, and products, little attention has been given to the design and requirements of facilities for enabling clinically relevant biofabrication solutions. These facilities are critical to overcoming the major hurdles to clinical translation, including solving important issues such as reproducibility, quality control, regulations, and commercialization. To improve process uniformity and ensure consistent development and production, large-scale manufacturing of engineered tissues and organs will require standardized facilities, equipment, qualification processes, automation, and information systems. This review presents current and forward-thinking guidelines to help design biofabrication laboratories engaged in engineering model and tissue constructs for therapeutic and non-therapeutic applications.

Cell proliferation and migration explain pore bridging dynamics in 3D printed scaffolds of different pore size

Tissue growth in bioscaffolds is influenced significantly by pore geometry, but how this geometric dependence emerges from dynamic cellular processes such as cell proliferation and cell migration remains poorly understood. Here we investigate the influence of pore size on the time required to bridge pores in thin 3D-printed scaffolds. Experimentally, new tissue infills the pores continually from their perimeter under strong curvature control, which leads the tissue front to round off with time. Despite the varied shapes assumed by the tissue during this evolution, we find that time to bridge a pore simply increases linearly with the overall pore size. To disentangle the biological influence of cell behaviour and the mechanistic influence of geometry in this experimental observation, we propose a simple reaction-diffusion model of tissue growth based on Porous-Fisher invasion of cells into the pores. First, this model provides a good qualitative representation of the evolution of the tissue; new tissue in the model grows at an effective rate that depends on the local curvature of the tissue substrate. Second, the model suggests that a linear dependence of bridging time with pore size arises due to geometric reasons alone, not to differences in cell behaviours across pores of different sizes. Our analysis suggests that tissue growth dynamics in these experimental constructs is dominated by mechanistic crowding effects that influence collective cell proliferation and migration processes, and that can be predicted by simple reaction-diffusion models of cells that have robust, consistent behaviours.

Identifying density-dependent interactions in collective cell behaviour

Scratch assays are routinely used to study collective cell behaviour in vitro. Typical experimental protocols do not vary the initial density of cells, and typical mathematical modelling approaches describe cell motility and proliferation based on assumptions of linear diffusion and logistic growth. Jin et al. (Jin et al. 2016 J. Theor. Biol. 390, 136-145 (doi:10.1016/j.jtbi.2015.10.040)) find that the behaviour of cells in scratch assays is density-dependent, and show that standard modelling approaches cannot simultaneously describe data initiated across a range of initial densities. To address this limitation, we calibrate an individual-based model to scratch assay data across a large range of initial densities. Our model allows proliferation, motility, and a direction bias to depend on interactions between neighbouring cells. By considering a hierarchy of models where we systematically and sequentially remove interactions, we perform model selection analysis to identify the minimum interactions required for the model to simultaneously describe data across all initial densities. The calibrated model is able to match the experimental data across all densities using a single parameter distribution, and captures details about the spatial structure of cells. Our results provide strong evidence to suggest that motility is density-dependent in these experiments. On the other hand, we do not see the effect of crowding on proliferation in these experiments. These results are significant as they are precisely the opposite of the assumptions in standard continuum models, such as the Fisher-Kolmogorov equation and its generalizations.

In vitro cell migration quantification method for scratch assays

The scratch assay is an in vitro technique used to assess the contribution of molecular and cellular mechanisms to cell migration. The assay can also be used to evaluate therapeutic compounds before clinical use. Current quantification methods of scratch assays deal poorly with irregular cell-free areas and crooked leading edges which are features typically present in the experimental data. We introduce a new migration quantification method, called 'monolayer edge velocimetry', that permits analysis of low-quality experimental data and better statistical classification of migration rates than standard quantification methods. The new method relies on quantifying the horizontal component of the cell monolayer velocity across the leading edge. By performing a classification test on in silico data, we show that the method exhibits significantly lower statistical errors than standard methods. When applied to in vitro data, our method outperforms standard methods by detecting differences in the migration rates between different cell groups that the other methods could not detect. Application of this new method will enable quantification of migration rates from in vitro scratch assay data that cannot be analysed using existing methods.

Practical parameter identifiability for spatio-temporal models of cell invasion

We examine the practical identifiability of parameters in a spatio-temporal reaction-diffusion model of a scratch assay. Experimental data involve fluorescent cell cycle labels, providing spatial information about cell position and temporal information about the cell cycle phase. Cell cycle labelling is incorporated into the reaction-diffusion model by treating the total population as two interacting subpopulations. Practical identifiability is examined using a Bayesian Markov chain Monte Carlo (MCMC) framework, confirming that the parameters are identifiable when we assume the diffusivities of the subpopulations are identical, but that the parameters are practically non-identifiable when we allow the diffusivities to be distinct. We also assess practical identifiability using a profile likelihood approach, providing similar results to MCMC with the advantage of being an order of magnitude faster to compute. Therefore, we suggest that the profile likelihood ought to be adopted as a screening tool to assess practical identifiability before MCMC computations are performed.

Neuron dynamics on directional surfaces

Geometrical features play a very important role in neuronal growth and the formation of functional connections between neuronal cells. Here, we analyze the dynamics of axonal growth for neuronal cells cultured on micro-patterned polydimethylsiloxane surfaces. We utilize fluorescence microscopy to image axons, quantify their dynamics, and demonstrate that periodic geometrical patterns impart strong directional bias to neuronal growth. We quantify axonal alignment and present a general stochastic approach that quantitatively describes the dynamics of the growth cones. Neuronal growth is described by a general phenomenological model, based on a simple automatic controller with a closed-loop feedback system. We demonstrate that axonal alignment on these substrates is determined by the surface geometry, and it is quantified by the deterministic part of the stochastic (Langevin and Fokker-Planck) equations. We also show that the axonal alignment with the surface patterns is greatly suppressed by the neuron treatment with Blebbistatin, a chemical compound that inhibits the activity of myosin II. These results give new insight into the role played by the molecular motors and external geometrical cues in guiding axonal growth, and could lead to novel approaches for bioengineering neuronal regeneration platforms.

Connecting gene expression to cellular movement: A transport model for cell migration

The adhesion properties and the mobility of biological cells play key roles in the propagation of cancer. These properties are expected to depend on intracellular processes and on the concentrations of chemicals inside the cell. While most existing reaction-diffusion models for cell migration consider that cell mobility and proliferation rate are constant or depend on an external diffusing species, they do not include the gene expression dynamics taking place in moving cells that affect cellular transport. In this work, we propose a multiscale model where mobility and proliferation depend explicitly on the cell's internal state. We focus more specifically on the case of cellular mobility in epithelial tissues. Wound-healing experiments have demonstrated that the loss of a key protein, E-cadherin, results in a significant increase in both mobility and invasiveness of epithelial cells, with dramatic consequences on cancer progression. We can reproduce the results of these experiments under various genetic conditions with a single set of parameters.

Accurate and efficient discretizations for stochastic models providing near agent-based spatial resolution at low computational cost

Understanding how cells proliferate, migrate and die in various environments is essential in determining how organisms develop and repair themselves. Continuum mathematical models, such as the logistic equation and the Fisher-Kolmogorov equation, can describe the global characteristics observed in commonly used cell biology assays, such as proliferation and scratch assays. However, these continuum models do not account for single-cell-level mechanics observed in high-throughput experiments. Mathematical modelling frameworks that represent individual cells, often called agent-based models, can successfully describe key single-cell-level features of these assays but are computationally infeasible when dealing with large populations. In this work, we propose an agent-based model with crowding effects that is computationally efficient and matches the logistic and Fisher-Kolmogorov equations in parameter regimes relevant to proliferation and scratch assays, respectively. This stochastic agent-based model allows multiple agents to be contained within compartments on an underlying lattice, thereby reducing the computational storage compared to existing agent-based models that allow one agent per site only. We propose a systematic method to determine a suitable compartment size. Implementing this compartment-based model with this compartment size provides a balance between computational storage, local resolution of agent behaviour and agreement with classical continuum descriptions.

Revisiting the Fisher-Kolmogorov-Petrovsky-Piskunov equation to interpret the spreading-extinction dichotomy

The Fisher-Kolmogorov-Petrovsky-Piskunov model, also known as the Fisher-KPP model, supports travelling wave solutions that are successfully used to model numerous invasive phenomena with applications in biology, ecology and combustion theory. However, there are certain phenomena that the Fisher-KPP model cannot replicate, such as the extinction of invasive populations. The Fisher-Stefan model is an adaptation of the Fisher-KPP model to include a moving boundary whose evolution is governed by a Stefan condition. The Fisher-Stefan model also supports travelling wave solutions; however, a key additional feature of the Fisher-Stefan model is that it is able to simulate population extinction, giving rise to a spreading-extinction dichotomy. In this work, we revisit travelling wave solutions of the Fisher-KPP model and show that these results provide new insight into travelling wave solutions of the Fisher-Stefan model and the spreading-extinction dichotomy. Using a combination of phase plane analysis, perturbation analysis and linearization, we establish a concrete relationship between travelling wave solutions of the Fisher-Stefan model and often-neglected travelling wave solutions of the Fisher-KPP model. Furthermore, we give closed-form approximate expressions for the shape of the travelling wave solutions of the Fisher-Stefan model in the limit of slow travelling wave speeds, c≪1.

Continuum descriptions of spatial spreading for heterogeneous cell populations: Theory and experiment

Variability in cell populations is frequently observed in both in vitro and in vivo settings. Intrinsic differences within populations of cells, such as differences in cell sizes or differences in rates of cell motility, can be present even within a population of cells from the same cell line. We refer to this variability as cell heterogeneity. Mathematical models of cell migration, for example, in the context of tumour growth and metastatic invasion, often account for both undirected (random) migration and directed migration that is mediated by cell-to-cell contacts and cell-to-cell adhesion. A key feature of standard models is that they often assume that the population is composed of identical cells with constant properties. This leads to relatively simple single-species homogeneous models that neglect the role of heterogeneity. In this work, we use a continuum modelling approach to explore the role of heterogeneity in spatial spreading of cell populations. We employ a three-species heterogeneous model of cell motility that explicitly incorporates different types of experimentally-motivated heterogeneity in cell sizes: (i) monotonically decreasing; (ii) uniform; (iii) non-monotonic; and (iv) monotonically increasing distributions of cell size. Comparing the density profiles generated by the three-species heterogeneous model with density profiles predicted by a more standard single-species homogeneous model reveals that when we are dealing with monotonically decreasing and uniform distributions a simple and computationally efficient single-species homogeneous model can be remarkably accurate in describing the evolution of a heterogeneous cell population. In contrast, we find that the simpler single-species homogeneous model performs relatively poorly when applied to non-monotonic and monotonically increasing distributions of cell sizes. Additional results for heterogeneity in parameters describing both undirected and directed cell migration are also considered, and we find that similar results apply.

Polarization wave at the onset of collective cell migration

Collective cell migration underlies morphogenesis, tissue regeneration, and cancer progression. How the biomechanical coupling between epithelial cells triggers and coordinates the collective migration is an open question. Here, we develop a one-dimensional model for an epithelial monolayer which predicts that after the onset of migration at an open boundary, cells in the bulk of the epithelium are gradually recruited into outward-directed motility, exhibiting traveling-wave-like behavior. We find an exact formula for the speed of this motility wave proportional to the square root of the cells' contractility, which accounts for cortex tension and adhesion between adjacent cells.

Mathematical models incorporating a multi-stage cell cycle replicate normally-hidden inherent synchronization in cell proliferation

We present a suite of experimental data showing that cell proliferation assays, prepared using standard methods thought to produce asynchronous cell populations, persistently exhibit inherent synchronization. Our experiments use fluorescent cell cycle indicators to reveal the normally hidden cell synchronization, by highlighting oscillatory subpopulations within the total cell population. These oscillatory subpopulations would never be observed without these cell cycle indicators. On the other hand, our experimental data show that the total cell population appears to grow exponentially, as in an asynchronous population. We reconcile these seemingly inconsistent observations by employing a multi-stage mathematical model of cell proliferation that can replicate the oscillatory subpopulations. Our study has important implications for understanding and improving experimental reproducibility. In particular, inherent synchronization may affect the experimental reproducibility of studies aiming to investigate cell cycle-dependent mechanisms, including changes in migration and drug response.

Using a continuum model to decipher the mechanics of embryonic tissue spreading from time-lapse image sequences: An approximate Bayesian computation approach

Advanced imaging techniques generate large datasets capable of describing the structure and kinematics of tissue spreading in embryonic development, wound healing, and the progression of many diseases. These datasets can be integrated with mathematical models to infer biomechanical properties of the system, typically identifying an optimal set of parameters for an individual experiment. However, these methods offer little information on the robustness of the fit and are generally ill-suited for statistical tests of multiple experiments. To overcome this limitation and enable efficient use of large datasets in a rigorous experimental design, we use the approximate Bayesian computation rejection algorithm to construct probability density distributions that estimate model parameters for a defined theoretical model and set of experimental data. Here, we demonstrate this method with a 2D Eulerian continuum mechanical model of spreading embryonic tissue. The model is tightly integrated with quantitative image analysis of different sized embryonic tissue explants spreading on extracellular matrix (ECM) and is regulated by a small set of parameters including forces on the free edge, tissue stiffness, strength of cell-ECM adhesions, and active cell shape changes. We find statistically significant trends in key parameters that vary with initial size of the explant, e.g., for larger explants cell-ECM adhesion forces are weaker and free edge forces are stronger. Furthermore, we demonstrate that estimated parameters for one explant can be used to predict the behavior of other similarly sized explants. These predictive methods can be used to guide further experiments to better understand how collective cell migration is regulated during development.

Colloidal Silver Induces Cytoskeleton Reorganization and E-Cadherin Recruitment at Cell-Cell Contacts in HaCaT Cells

Up until the first half of the 20th century, silver found significant employment in medical applications, particularly in the healing of open wounds, thanks to its antibacterial and antifungal properties. Wound repair is a complex and dynamic biological process regulated by several pathways that cooperate to restore tissue integrity and homeostasis. To facilitate healing, injuries need to be promptly treated. Recently, the interest in alternatives to antibiotics has been raised given the widespread phenomenon of antibiotic resistance. Among these alternatives, the use of silver appears to be a valid option, so a resurgence in its use has been recently observed. In particular, in contrast to ionic silver, colloidal silver, a suspension of metallic silver particles, shows antibacterial activity displaying less or no toxicity. However, the human health risks associated with exposure to silver nanoparticles (NP) appear to be conflicted, and some studies have suggested that it could be toxic in different cellular contexts. These potentially harmful effects of silver NP depend on various parameters including NP size, which commonly range from 1 to 100 nm. In this study, we analyzed the effect of a colloidal silver preparation composed of very small and homogeneous nanoparticles of 0.62 nm size, smaller than those previously tested. We found no adverse effect on the cell proliferation of HaCaT cells, even at high NP concentration. Time-lapse microscopy and indirect immunofluorescence experiments demonstrated that this preparation of colloidal silver strongly increased cell migration, re-modeled the cytoskeleton, and caused recruitment of E-cadherin at cell-cell junctions of human cultured keratinocytes.

Anomalous diffusion for neuronal growth on surfaces with controlled geometries

Geometrical cues are known to play a very important role in neuronal growth and the formation of neuronal networks. Here, we present a detailed analysis of axonal growth and dynamics for neuronal cells cultured on patterned polydimethylsiloxane surfaces. We use fluorescence microscopy to image neurons, quantify their dynamics, and demonstrate that the substrate geometrical patterns cause strong directional alignment of axons. We quantify axonal growth and report a general stochastic approach that quantitatively describes the motion of growth cones. The growth cone dynamics is described by Langevin and Fokker-Planck equations with both deterministic and stochastic contributions. We show that the deterministic terms contain both the angular and speed dependence of axonal growth, and that these two contributions can be separated. Growth alignment is determined by surface geometry, and it is quantified by the deterministic part of the Langevin equation. We combine experimental data with theoretical analysis to measure the key parameters of the growth cone motion: speed and angular distributions, correlation functions, diffusion coefficients, characteristics speeds and damping coefficients. We demonstrate that axonal dynamics displays a cross-over from Brownian motion (Ornstein-Uhlenbeck process) at earlier times to anomalous dynamics (superdiffusion) at later times. The superdiffusive regime is characterized by non-Gaussian speed distributions and power law dependence of the axonal mean square length and the velocity correlation functions. These results demonstrate the importance of geometrical cues in guiding axonal growth, and could lead to new methods for bioengineering novel substrates for controlling neuronal growth and regeneration.

Role of geometrical cues in neuronal growth

Geometrical cues play an essential role in neuronal growth. Here, we quantify axonal growth on surfaces with controlled geometries and report a general stochastic approach that quantitatively describes the motion of growth cones. We show that axons display a strong directional alignment on micropatterned surfaces when the periodicity of the patterns matches the dimension of the growth cone. The growth cone dynamics on surfaces with uniform geometry is described by a linear Langevin equation with both deterministic and stochastic contributions. In contrast, axonal growth on surfaces with periodic patterns is characterized by a system of two generalized Langevin equations with both linear and quadratic velocity dependence and stochastic noise. We combine experimental data with theoretical analysis to measure the key parameters of the growth cone motion: angular distributions, correlation functions, diffusion coefficients, characteristics speeds, and damping coefficients. We demonstrate that axonal dynamics displays a crossover from an Ornstein-Uhlenbeck process to a nonlinear stochastic regime when the geometrical periodicity of the pattern approaches the linear dimension of the growth cone. Growth alignment is determined by surface geometry, which is fully quantified by the deterministic part of the Langevin equation. These results provide insight into the role of curvature sensing proteins and their interactions with geometrical cues.

A free boundary model of epithelial dynamics

In this work we analyse a one-dimensional, cell-based model of an epithelial sheet. In the model, cells interact with their nearest neighbouring cells and move deterministically. Cells also proliferate stochastically, with the rate of proliferation specified as a function of the cell length. This mechanical model of cell dynamics gives rise to a free boundary problem. We construct a corresponding continuum-limit description where the variables in the continuum limit description are expanded in powers of the small parameter 1/N, where N is the number of cells in the population. By carefully constructing the continuum limit description we obtain a free boundary partial differential equation description governing the density of the cells within the evolving domain, as well as a free boundary condition that governs the evolution of the domain. We show that care must be taken to arrive at a free boundary condition that conserves mass. By comparing averaged realisations of the cell-based model with the numerical solution of the free boundary partial differential equation, we show that the new mass-conserving boundary condition enables the coarse-grained partial differential equation model to provide very accurate predictions of the behaviour of the cell-based model, including both evolution of the cell density, and the position of the free boundary, across a range of interaction potentials and proliferation functions in the cell based model.