Abstract: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.

Optimal quantification of contact inhibition in cell populations

The impact of contact inhibition on collective cell migration and proliferation

Minimal cellular automaton model with heterogeneous cell sizes predicts epithelial colony growth

Quantifying assays: inhibition of signalling pathways of cancer

Local contact inhibition leads to universal principles of cell population growth

Quantifying cell cycle regulation by tissue crowding

Quantifying assays: A Modeling tale of variability in cancer therapeutics assessed on cancer cells

Control of cell colony growth by contact inhibition

Influence of cell-cell interactions on the population growth rate in a tumor

Limits on the accuracy of contact inhibition of locomotion

Effect of Contact Inhibition Locomotion on confined cellular organization

Modeling tumor cell migration: From microscopic to macroscopic models

Discrete and continuum approximations for collective cell migration in a scratch assay with cell size dynamics

Inference of weak-form partial differential equations describing migration and proliferation mechanisms in wound healing experiments on cancer cells

Using Experimental Data and Information Criteria to Guide Model Selection for Reaction–Diffusion Problems in Mathematical Biology

Modelling count data with partial differential equation models in biology

A Reaction-Diffusion Model of the Cadherin-Catenin System: A Possible Mechanism for Contact Inhibition and Implications for Tumorigenesis

Partial differential equation-based inference of migration and proliferation mechanisms in cancer cell populations

Collective motion of cells crawling on a substrate: roles of cell shape and contact inhibition

A Mathematical Tumor Model with Immune Resistance and Drug Therapy: An Optimal Control Approach