Choosing an Appropriate Modelling Framework for Analysing Multispecies Co-culture Cell Biology Experiments

In vitro cell biology assays play a crucial role in informing our understanding of the migratory, proliferative and invasive properties of many cell types in different biological contexts. While mono-culture assays involve the study of a population of cells composed of a single cell type, co-culture assays study a population of cells composed of multiple cell types (or subpopulations of cells). Such co-culture assays can provide more realistic insights into many biological processes including tissue repair, tissue regeneration and malignant spreading. Typically, system parameters, such as motility and proliferation rates, are estimated by calibrating a mathematical or computational model to the observed experimental data. However, parameter estimates can be highly sensitive to the choice of model and modelling framework. This observation motivates us to consider the fundamental question of how we can best choose a model to facilitate accurate parameter estimation for a particular assay. In this work we describe three mathematical models of mono-culture and co-culture assays that include different levels of spatial detail. We study various spatial summary statistics to explore if they can be used to distinguish between the suitability of each model over a range of parameter space. Our results for mono-culture experiments are promising, in that we suggest two spatial statistics that can be used to direct model choice. However, co-culture experiments are far more challenging: we show that these same spatial statistics which provide useful insight into mono-culture systems are insufficient for co-culture systems. Therefore, we conclude that great care ought to be exercised when estimating the parameters of co-culture assays.

[1]  Ying Wang,et al.  Enhanced Keratinocyte Proliferation and Migration in Co-culture with Fibroblasts , 2012, PloS one.

[2]  Zee Upton,et al.  Development of a three-dimensional human skin equivalent wound model for investigating novel wound healing therapies. , 2010, Tissue engineering. Part C, Methods.

[3]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[4]  Matthew J Simpson,et al.  Corrected mean-field models for spatially dependent advection-diffusion-reaction phenomena. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[5]  Peter J. Diggle,et al.  Statistical Analysis of Spatial and Spatio-Temporal Point Patterns , 2013 .

[6]  C. Please,et al.  Experimental characterization and computational modelling of two-dimensional cell spreading for skeletal regeneration , 2007, Journal of The Royal Society Interface.

[7]  H. Green,et al.  QUANTITATIVE STUDIES OF THE GROWTH OF MOUSE EMBRYO CELLS IN CULTURE AND THEIR DEVELOPMENT INTO ESTABLISHED LINES , 1963, The Journal of cell biology.

[8]  D. McElwain,et al.  Do Pioneer Cells Exist? , 2014, PloS one.

[9]  A. Singer,et al.  Maximum entropy formulation of the Kirkwood superposition approximation. , 2004, The Journal of chemical physics.

[10]  B. Bolker,et al.  Spatial Moment Equations for Plant Competition: Understanding Spatial Strategies and the Advantages of Short Dispersal , 1999, The American Naturalist.

[11]  C. Zheng,et al.  Applied contaminant transport modeling , 2002 .

[12]  Christophe Deroulers,et al.  Modeling tumor cell migration: From microscopic to macroscopic models. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Charles L. Tucker,et al.  Lagrangian particle calculations of distributive mixing: Limitations and applications , 2006 .

[14]  Ulf Dieckmann,et al.  Relaxation Projections and the Method of Moments , 1999 .

[15]  Kerry A Landman,et al.  Multi-scale modeling of a wound-healing cell migration assay. , 2007, Journal of theoretical biology.

[16]  J. Schalkwijk,et al.  Keratinocytes drive melanoma invasion in a reconstructed skin model , 2010, Melanoma research.

[17]  Matthew J. Simpson,et al.  Multi-species simple exclusion processes , 2009 .

[18]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[19]  William H. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[20]  P. M. Lee,et al.  Random Walks and Random Environments: Volume 1: Random Walks , 1995 .

[21]  J. Murray,et al.  Virtual and real brain tumors: using mathematical modeling to quantify glioma growth and invasion , 2003, Journal of the Neurological Sciences.

[22]  Geoffrey W. Stevens,et al.  Cell migration and proliferation during monolayer formation and wound healing , 2009 .

[23]  Matthew J Simpson,et al.  Correcting mean-field approximations for birth-death-movement processes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Matthew J. Simpson,et al.  Distinguishing between mean-field, moment dynamics and stochastic descriptions of birth death movement processes , 2014 .

[25]  Ulf Dieckmann,et al.  A multiscale maximum entropy moment closure for locally regulated space–time point process models of population dynamics , 2011, Journal of mathematical biology.

[26]  Deborah C Markham,et al.  Simplified method for including spatial correlations in mean-field approximations. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  Kieran J. Sharkey,et al.  Pair-level approximations to the spatio-temporal dynamics of epidemics on asymmetric contact networks , 2006, Journal of mathematical biology.

[28]  J. Filipe,et al.  Hybrid closure-approximation to epidemic models , 1999 .

[29]  Ulf Dieckmann,et al.  On moment closures for population dynamics in continuous space. , 2004, Journal of theoretical biology.

[30]  Matthew J Simpson,et al.  Experimental and Modelling Investigation of Monolayer Development with Clustering , 2013, Bulletin of mathematical biology.

[31]  Edward A. Codling,et al.  Random walk models in biology , 2008, Journal of The Royal Society Interface.

[32]  Kerry A Landman,et al.  Quantifying evenly distributed states in exclusion and nonexclusion processes. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[33]  J A N Filipe,et al.  Analytical methods for predicting the behaviour of population models with general spatial interactions. , 2003, Mathematical biosciences.

[34]  B. Bolker,et al.  Using Moment Equations to Understand Stochastically Driven Spatial Pattern Formation in Ecological Systems , 1997, Theoretical population biology.

[35]  Jonathan A. Sherratt,et al.  Models of epidermal wound healing , 1990, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[36]  Christine Unger,et al.  In vitro cell migration and invasion assays. , 2013, Mutation research.

[37]  T. Liggett,et al.  Stochastic Interacting Systems: Contact, Voter and Exclusion Processes , 1999 .

[38]  U. Dieckmann,et al.  POPULATION GROWTH IN SPACE AND TIME: SPATIAL LOGISTIC EQUATIONS , 2003 .

[39]  Deborah C Markham,et al.  Incorporating spatial correlations into multispecies mean-field models. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  Steven C. Chapra,et al.  Numerical Methods for Engineers , 1986 .

[41]  Ruth E. Baker,et al.  Assessing the role of spatial correlations during collective cell spreading , 2014, Scientific Reports.