Abstract
AbstractCentre-based or cell-centre models are a framework for the computational study of multicellular systems with widespread use in cancer modelling and computational developmental biology. At the core of these models are the numerical method used to update cell positions and the force functions that encode the pairwise mechanical interactions of cells. For the latter, there are multiple choices that could potentially affect both the biological behaviour captured, and the robustness and efficiency of simulation. For example, available open-source software implementations of centre-based models rely on different force functions for their default behaviour and it is not straightforward for a modeller to know if these are interchangeable. Our study addresses this problem and contributes to the understanding of the potential and limitations of three popular force functions from a numerical perspective. We show empirically that choosing the force parameters such that the relaxation time for two cells after cell division is consistent between different force functions results in good agreement of the population radius of a two-dimensional monolayer relaxing mechanically after intense cell proliferation. Furthermore, we report that numerical stability is not sufficient to prevent unphysical cell trajectories following cell division, and consequently, that too large time steps can cause geometrical differences at the population level.
Publisher
Springer Science and Business Media LLC
Subject
Computational Theory and Mathematics,General Agricultural and Biological Sciences,Pharmacology,General Environmental Science,General Biochemistry, Genetics and Molecular Biology,General Mathematics,Immunology,General Neuroscience
Reference75 articles.
1. Ascher UM, Petzold LR (1998) Computer methods for ordinary differential equations and differential-algebraic equations. SIAM
2. Atkins P, de Paula J, Keeler J (2018) Atkins’ physical chemistry 11e. Oxford University Press, Oxford
3. Atwell K (2016) Investigating the interplay between cellular mechanics and decision-making in the c. elegans germ line. PhD thesis, University of Oxford
4. Atwell K, Qin Z, Gavaghan D, Kugler H, Hubbard EJA, Osborne JM (2015) Mechano-logical model of c. elegans germ line suggests feedback on the cell cycle. Development 142(22):3902–3911
5. Byrd RH, Lu P, Nocedal J, Zhu C (1995) A limited memory algorithm for bound constrained optimization. SIAM J Sci Comput 16(5):1190–1208
Cited by
11 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献