Understanding the relationship between cell death and tissue shrinkage via a stochastic agent-based model.

Cell death, a process which can occur both naturally and in response to insult, is both a complex and diverse phenomenon. Under some circumstances, dying cells actively contract and cause their neighbors to rearrange and maintain tissue integrity. Under other circumstances, dying cells leave behind gaps, which results in tissue separation. A better understanding of how the cellular scale features of cell death manifest on the population scale has implications ranging from morphogenesis to tumor response to treatment. However, the mechanistic relationship between cell death and population scale shrinkage is not well understood, and computational methods for studying these relationships are not well established. Here we propose a mechanically robust agent-based cell model designed to capture the implications of cell death on the population scale. In our agent-based model, algorithmic rules applied on the cellular level emerge on the population scale where their effects are quantified. To better quantify model uncertainty and parameter interactions, we implement a recently developed technique for conducting a variance-based sensitivity analysis on the stochastic model. From this analysis and subsequent investigation, we find that cellular scale shrinkage has the largest influence of all model parameters tested, and that by adjusting cellular scale shrinkage population shrinkage varies widely even across simulations which contain the same fraction of dying cells. We anticipate that the methods and results presented here are a starting point for significant future investigation toward quantifying the implications of different mechanisms of cell death on population and tissue scale behavior.

[1]  Stefano Tarantola,et al.  Global Sensitivity Analysis: An Introduction , 2005 .

[2]  J F Kerr,et al.  Shrinkage necrosis: A distinct mode of cellular death , 1971, The Journal of pathology.

[3]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[4]  T. Stylianopoulos,et al.  Biphasic modeling of brain tumor biomechanics and response to radiation treatment. , 2016, Journal of biomechanics.

[5]  David J. Williams,et al.  Changes in the stiffness of human mesenchymal stem cells with the progress of cell death as measured by atomic force microscopy. , 2014, Journal of biomechanics.

[6]  C. Linder,et al.  Modeling tumor growth with peridynamics , 2017, Biomechanics and modeling in mechanobiology.

[7]  S. Elmore Apoptosis: A Review of Programmed Cell Death , 2007, Toxicologic pathology.

[8]  S. Silling Reformulation of Elasticity Theory for Discontinuities and Long-Range Forces , 2000 .

[9]  C. Zeiss,et al.  The Apoptosis-Necrosis Continuum: Insights from Genetically Altered Mice , 2003, Veterinary pathology.

[10]  B. Zhivotovsky,et al.  Death through a tragedy: mitotic catastrophe , 2008, Cell Death and Differentiation.

[11]  J. Tinsley Oden,et al.  Analysis and numerical solution of stochastic phase‐field models of tumor growth , 2015 .

[12]  L. Covarrubias,et al.  Interdigital cell death function and regulation: New insights on an old programmed cell death model , 2011, Development, growth & differentiation.

[13]  Will Usher,et al.  SALib: An open-source Python library for Sensitivity Analysis , 2017, J. Open Source Softw..

[14]  Alen Alexanderian,et al.  Efficient Computation of Sobol' Indices for Stochastic Models , 2016, SIAM J. Sci. Comput..

[15]  F J Vermolen,et al.  A multi-agent cell-based model for wound contraction. , 2016, Journal of biomechanics.

[16]  R. Lehoucq,et al.  Peridynamic Theory of Solid Mechanics , 2010 .

[17]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[18]  C. Bortner,et al.  Apoptotic volume decrease and the incredible shrinking cell , 2002, Cell Death and Differentiation.

[19]  Erdogan Madenci,et al.  An adaptive dynamic relaxation method for quasi-static simulations using the peridynamic theory , 2010 .

[20]  J D Humphrey,et al.  Perspectives on biological growth and remodeling. , 2011, Journal of the mechanics and physics of solids.

[21]  Alexander G. Fletcher,et al.  Chaste: A test-driven approach to software development for biological modelling , 2009, Comput. Phys. Commun..

[22]  Paola Annoni,et al.  Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index , 2010, Comput. Phys. Commun..

[23]  R. Baskar,et al.  Cancer and Radiation Therapy: Current Advances and Future Directions , 2012, International journal of medical sciences.

[24]  David John Littlewood,et al.  Roadmap for Peridynamic Software Implementation , 2015 .

[25]  Arthur Konnerth,et al.  Apoptotic Force and Tissue Dynamics During Drosophila Embryogenesis , 2008 .

[26]  C. Linder,et al.  Modeling mechanical inhomogeneities in small populations of proliferating monolayers and spheroids , 2017, Biomechanics and Modeling in Mechanobiology.

[27]  T. Newman,et al.  Modeling cell rheology with the Subcellular Element Model , 2008, Physical biology.

[28]  E. Bernhard,et al.  How does radiation kill cells? , 1999, Current opinion in chemical biology.

[29]  S. Silling,et al.  A meshfree method based on the peridynamic model of solid mechanics , 2005 .

[30]  G. An,et al.  Agent‐based models in translational systems biology , 2009, Wiley interdisciplinary reviews. Systems biology and medicine.

[31]  A. Saltelli,et al.  Making best use of model evaluations to compute sensitivity indices , 2002 .

[32]  Dirk Drasdo,et al.  Off-Lattice Agent-Based Models for Cell and Tumor Growth: Numerical Methods, Implementation, and Applications , 2018 .

[33]  Timon Rabczuk,et al.  Dual‐horizon peridynamics , 2015, 1506.05146.

[34]  H. Steller,et al.  Shaping organisms with apoptosis , 2013, Cell Death and Differentiation.

[35]  E Kuhl,et al.  Computational modeling of chemo-bio-mechanical coupling: a systems-biology approach toward wound healing , 2016, Computer methods in biomechanics and biomedical engineering.

[36]  P. Vandenabeele,et al.  Regulated necrosis: the expanding network of non-apoptotic cell death pathways , 2014, Nature Reviews Molecular Cell Biology.

[37]  D. Drasdo,et al.  A single-cell-based model of tumor growth in vitro: monolayers and spheroids , 2005, Physical biology.

[38]  I. Sobola,et al.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[39]  Selda Oterkus,et al.  Peridynamics for the solution of multiphysics problems , 2015 .

[40]  S. Silling,et al.  Peridynamic States and Constitutive Modeling , 2007 .

[41]  Yusheng Feng,et al.  A fully coupled space-time multiscale modeling framework for predicting tumor growth. , 2017, Computer methods in applied mechanics and engineering.

[42]  Winfried Kurth,et al.  Facilitating Parameter Estimation and Sensitivity Analysis of Agent-Based Models: A Cookbook Using NetLogo and 'R' , 2014, J. Artif. Soc. Soc. Simul..

[43]  J. Freidman,et al.  Multivariate adaptive regression splines , 1991 .

[44]  Erkan Oterkus,et al.  Peridynamic Theory and Its Applications , 2013 .

[45]  M. Suzanne,et al.  The Morphogenetic Role of Apoptosis. , 2015, Current topics in developmental biology.

[46]  D. Kirschner,et al.  A methodology for performing global uncertainty and sensitivity analysis in systems biology. , 2008, Journal of theoretical biology.

[47]  C. Linder,et al.  Quantifying the relationship between cell division angle and morphogenesis through computational modeling. , 2017, Journal of theoretical biology.

[48]  Yang Cao,et al.  Sensitivity analysis of discrete stochastic systems. , 2005, Biophysical journal.

[49]  Nick Jagiella,et al.  Inferring Growth Control Mechanisms in Growing Multi-cellular Spheroids of NSCLC Cells from Spatial-Temporal Image Data , 2016, PLoS Comput. Biol..

[50]  Wei Sun,et al.  Three dimensional multi-scale modelling and analysis of cell damage in cell-encapsulated alginate constructs. , 2010, Journal of biomechanics.

[51]  C. Potten,et al.  Education and debate , 2022 .

[52]  D. Kirschner,et al.  Computational Modeling Predicts Simultaneous Targeting of Fibroblasts and Epithelial Cells Is Necessary for Treatment of Pulmonary Fibrosis , 2016, Front. Pharmacol..

[53]  Robert Jeraj,et al.  An imaging-based tumour growth and treatment response model: investigating the effect of tumour oxygenation on radiation therapy response , 2008, Physics in medicine and biology.

[54]  Johan Pallud,et al.  A Tumor Growth Inhibition Model for Low-Grade Glioma Treated with Chemotherapy or Radiotherapy , 2012, Clinical Cancer Research.

[55]  Thomas Mangeat,et al.  Apico-basal forces exerted by apoptotic cells drive epithelium folding , 2015, Nature.

[56]  Konstantina S. Nikita,et al.  Simulating growth dynamics and radiation response of avascular tumour spheroids - model validation in the case of an EMT6/Ro multicellular spheroid , 2004, Comput. Methods Programs Biomed..