Inferring parameters for a lattice-free model of cell migration and proliferation using experimental data

Collective cell spreading takes place in spatially continuous environments, yet it is often modelled using discrete lattice-based approaches. Here, we use data from a series of cell proliferation assays, with a prostate cancer cell line, to calibrate a spatially continuous individual based model (IBM) of collective cell migration and proliferation. The IBM explicitly accounts for crowding effects by modifying the rate of movement, direction of movement, and the rate of proliferation by accounting for pair-wise interactions. Taking a Bayesian approach we estimate the free parameters in the IBM using rejection sampling on three separate, independent experimental data sets. Since the posterior distributions for each experiment are similar, we perform simulations with parameters sampled from a new posterior distribution generated by combining the three data sets. To explore the predictive power of the calibrated IBM, we forecast the evolution of a fourth experimental data set. Overall, we show how to calibrate a lattice-free IBM to experimental data, and our work highlights the importance of interactions between individuals. Despite great care taken to distribute cells as uniformly as possible experimentally, we find evidence of significant spatial clustering over short distances, suggesting that standard mean-field models could be inappropriate.

[1]  Matthew J Simpson,et al.  Models of collective cell behaviour with crowding effects: comparing lattice-based and lattice-free approaches , 2012, Journal of The Royal Society Interface.

[2]  Matthew J Simpson,et al.  Looking inside an invasion wave of cells using continuum models: proliferation is the key. , 2006, Journal of theoretical biology.

[3]  D. L. Sean McElwain,et al.  Interpreting scratch assays using pair density dynamics and approximate Bayesian computation , 2014, Open Biology.

[4]  Ramon Grima,et al.  Approximation and inference methods for stochastic biochemical kinetics—a tutorial review , 2016, 1608.06582.

[5]  Joe Pitt-Francis,et al.  Bayesian Calibration, Validation and Uncertainty Quantification for Predictive Modelling of Tumour Growth: A Tutorial , 2017, Bulletin of Mathematical Biology.

[6]  M. Evans Statistical Distributions , 2000 .

[7]  Philip K Maini,et al.  Traveling wave model to interpret a wound-healing cell migration assay for human peritoneal mesothelial cells. , 2004, Tissue engineering.

[8]  D. L. Sean McElwain,et al.  Estimating cell diffusivity and cell proliferation rate by interpreting IncuCyte ZOOM™ assay data using the Fisher-Kolmogorov model , 2015, BMC Systems Biology.

[9]  Fabian J. Theis,et al.  Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion , 2016, PLoS Comput. Biol..

[10]  L. G. Pillis,et al.  A Comparison and Catalog of Intrinsic Tumor Growth Models , 2013, Bulletin of mathematical biology.

[11]  C. Bertolotto,et al.  FANCD2 functions as a critical factor downstream of MiTF to maintain the proliferation and survival of melanoma cells , 2016, Scientific Reports.

[12]  T. Skalak,et al.  The FASEB Journal express article 10.1096/fj.03-0933fje. Published online February 6, 2004. Multicellular simulation predicts microvascular patterning and in silico tissue assembly , 2022 .

[13]  Matthew J. Simpson,et al.  Spatial structure arising from neighbour-dependent bias in collective cell movement , 2016, PeerJ.

[14]  Matthew J Simpson,et al.  Quantifying the effect of experimental design choices for in vitro scratch assays. , 2016, Journal of theoretical biology.

[15]  Esha T. Shah,et al.  Logistic Proliferation of Cells in Scratch Assays is Delayed , 2017, Bulletin of Mathematical Biology.

[16]  Matthew J Simpson,et al.  Optimal Quantification of Contact Inhibition in Cell Populations. , 2017, Biophysical journal.

[17]  Xiaobo Zhou,et al.  Computational Modeling of 3D Tumor Growth and Angiogenesis for Chemotherapy Evaluation , 2014, PloS one.

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

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

[20]  Philip K Maini,et al.  Models, measurement and inference in epithelial tissue dynamics. , 2015, Mathematical biosciences and engineering : MBE.

[21]  A. Kabla,et al.  Collective cell migration: leadership, invasion and segregation , 2011, Journal of The Royal Society Interface.

[22]  K. Sharkey Deterministic epidemiological models at the individual level , 2008, Journal of mathematical biology.

[23]  Ursula Klingmüller,et al.  An individual-based model for collective cancer cell migration explains speed dynamics and phenotype variability in response to growth factors , 2017, npj Systems Biology and Applications.

[24]  Alexander G Fletcher,et al.  Mathematical modeling of monoclonal conversion in the colonic crypt. , 2012, Journal of theoretical biology.

[25]  Sarah Filippi,et al.  A framework for parameter estimation and model selection from experimental data in systems biology using approximate Bayesian computation , 2014, Nature Protocols.

[26]  Alex James,et al.  Collective Cell Behaviour with Neighbour-Dependent Proliferation, Death and Directional Bias , 2016, Bulletin of mathematical biology.

[27]  Andrew R. Francis,et al.  Using Approximate Bayesian Computation to Estimate Tuberculosis Transmission Parameters From Genotype Data , 2006, Genetics.

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

[29]  Matthew J Simpson,et al.  Quantifying spatial structure in experimental observations and agent-based simulations using pair-correlation functions. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

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

[32]  S. McCue,et al.  A Bayesian Computational Approach to Explore the Optimal Duration of a Cell Proliferation Assay , 2017, Bulletin of Mathematical Biology.

[33]  Vipin Kumar,et al.  Techniques for grounding agent-based simulations in the real domain: a case study in experimental autoimmune encephalomyelitis , 2012 .

[34]  Ruth E. Baker,et al.  Multiple types of data are required to identify the mechanisms influencing the spatial expansion of melanoma cell colonies , 2013, BMC Systems Biology.

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

[36]  Q. Sang,et al.  Human Mesenchymal Stem Cells Are Resistant to Paclitaxel by Adopting a Non-Proliferative Fibroblastic State , 2015, PloS one.

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

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

[39]  G. Sanguinetti,et al.  Cox process representation and inference for stochastic reaction–diffusion processes , 2016, Nature Communications.

[40]  U. Dieckmann,et al.  Population Growth in Space and Time , 2003 .

[41]  Muhammad H. Zaman,et al.  A Computational Model for Collective Cellular Motion in Three Dimensions: General Framework and Case Study for Cell Pair Dynamics , 2013, PloS one.

[42]  J. Lechner,et al.  Establishment and characterization of a human prostatic carcinoma cell line (PC-3). , 1979, Investigative urology.

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