Mechanochemical Models for Calcium Waves in Embryonic Epithelia

• DOI: 10.1007/s10013-022-00579-y
• Publication Date: 2022-09-14T16:02:53Z
• Authors: Katerina Kaouri; Paul E. Méndez; Ricardo Ruiz-Baier
• Keywords: 0301 basic medicine; Finite element method; Cell Mechanics and Extracellular Matrix Interactions; Biomedical Engineering; Organic chemistry; FOS: Medical engineering; Mechanics; Mathematical analysis; Quantum mechanics; 03 medical and health sciences; Engineering; Context (archaeology); Biochemistry, Genetics and Molecular Biology; FOS: Mathematics; Hyperelastic Modeling; Mathematics - Numerical Analysis; Tissues and Organs (q-bio.TO); Biomechanical Modeling of Arterial Tissues; Physics; Life Sciences; Paleontology; Quantitative Biology - Tissues and Organs; Viscoelasticity; Geology; Numerical Analysis (math.NA); Cell Biology; FOS: Earth and related environmental sciences; 3. Good health; Chemistry; Biology and Pathology of Keratins and Related Disorders; FOS: Biological sciences; Physical Sciences; Nonlinear system; 92C15, 65M60, 35K57, 74L15; Thermodynamics; Calcium; Bifurcation; Mathematics

AbstractIn embryogenesis, epithelial cells acting as individual entities or as coordinated aggregates in a tissue, exhibit strong coupling between mechanical responses to internally or externally applied stresses and chemical signalling. One of the most important chemical signals in this process is calcium. This mechanochemical coupling and intercellular communication drive the coordination of morphogenetic movements which are characterised by drastic changes in the concentration of calcium in the tissue. In this paper we extend the recent mechanochemical model in Kaouri et al. (J. Math. Biol.78, 2059–2092, 2019), for an epithelial continuum in one dimension, to a more realistic multi-dimensional case. The resulting parametrised governing equations consist of an advection-diffusion-reaction system for calcium signalling coupled with active-stress linear viscoelasticity and equipped with pure Neumann boundary conditions. We implement a finite element method in perturbed saddle-point form for the simulation of this complex multiphysics problem. Special care is taken in the treatment of the stress-free boundary conditions for the viscoelasticity in order to eliminate rigid motions from the space of admissible displacements. The stability and solvability of the continuous weak formulation is shown using fixed-point theory. Guided by the bifurcation analysis of the one-dimensional model, we analyse the behaviour of the system as two bifurcation parameters vary: the level of IP3concentration and the strength of the mechanochemical coupling. We identify the parameter regions giving rise to solitary waves and periodic wavetrains of calcium. Furthermore, we demonstrate the nucleation of calcium sparks into synchronous calcium waves coupled with deformation. This model can be employed to gain insights into recent experimental observations in the context of embryogenesis, but also in other biological systems such as cancer cells, wound healing, keratinocytes, or white blood cells.