Search

Cells mechanically interact with their environment to sense, for example, topography, elasticity and mechanical cues from other cells. Mechano-sensing has profound effects on cellular behaviour, including motility. The current study aims to develop a mathematical model of cellular mechano-sensing on planar elastic substrates and demonstrate the model’s predictive capabilities for the motility of individual cells in a colony.In the model, a cell is assumed to transmit an adhesion force, derived from a dynamic focal adhesion integrin density, that locally deforms a substrate, and to sense substrate deformation originating from neighbouring cells. The substrate deformation from multiple cells is expressed as total strain energy density with a spatially varying gradient. The magnitude and direction of the gradient at the cell location define the cell motion. Cell-substrate friction, partial motion randomness, and cell death and division are included.The substrate deformation by a single cell and the motility of two cells are presented for several substrate elasticities and thicknesses. The collective motility of 25 cells on a uniform substrate mimicking the closure of a circular wound of 200 μm is predicted for deterministic and random motion. Cell motility on substrates with varying elasticity and thickness is explored for four cells and 15 cells, the latter again mimicking wound closure. Wound closure by 45 cells is used to demonstrate the simulation of cell death and division during migration.The mathematical model can adequately simulate the mechanically-induced collective cell motility on planar elastic substrates. The model is suitable for extension to other cell and substrates shapes and the inclusion of chemotactic cues, offering the potential to complementin vitroandin vivostudies.

We consider the stability analysis of a two-dimensional model for post-burn contraction. The model is based on morphoelasticity for permanent deformations and combined with a chemical-biological model that incorporates cellular densities, collagen density, and the concentration of chemoattractants. We formulate stability conditions depending on the decay rate of signaling molecules for both the continuous partial differential equations-based problem and the (semi-)discrete representation. We analyze the difference and convergence between the resulting spatial eigenvalues from the continuous and semi-discrete problems.

Deep dermal wounds induce skin contraction as a result of the traction forcing exerted by (myo)fibroblasts on their immediate environment. These (myo)fibroblasts are skin cells that are responsible for the regeneration of collagen that is necessary for the integrity of skin We consider several mathematical issues regarding models that simulate traction forces exerted by (myo)fibroblasts. Since the size of cells (e.g. (myo)fibroblasts) is much smaller than the size of the domain of computation, one often considers point forces, modelled by Dirac Delta distributions on boundary segments of cells to simulate the traction forces exerted by the skin cells. In the current paper, we treat the forces that are directed normal to the cell boundary and toward the cell centre. Since it can be shown that there exists no smooth solution, at least not in H1 for solutions to the governing momentum balance equation, we analyse the convergence and quality of approximation. Furthermore, the expected finite element problems that we get necessitate to scrutinize alternative model formulations, such as the use of smoothed Dirac Delta distributions, or the so-called smoothed particle approach as well as the so-called 'hole' approach where cellular forces are modelled through the use of (natural) boundary conditions. In this paper, we investigate and attempt to quantify the conditions for consistency between the various approaches. This has resulted into error analyses in the L2-norm of the numerical solution based on Galerkin principles that entail Lagrangian basis functions. The paper also addresses well-posedness in terms of existence and uniqueness. The current analysis has been performed for the linear steady-state (hence neglecting inertia and damping) momentum equations under the assumption of Hooke's law. (C) 2022 The Author(s). Published by Elsevier B.V. on behalf of International Association for Mathematics and Computers in Simulation (IMACS).

We consider a two-dimensional biomorphoelastic model describing post-burn scar contraction. This model describes skin displacement and the development of the effective Eulerian strain in the tissue. Besides these mechanical components, signaling molecules, fibroblasts, myofibroblasts, and collagen also play a significant role in the model. We perform a sensitivity analysis for the independent parameters of the model and focus on the effects on features of the relative surface area and the total strain energy density. We conclude that the most sensitive parameters are the Poisson's ratio, the equilibrium collagen concentration, the contraction inhibitor constant, and the myofibroblast apoptosis rate. Next to these insights, we perform a sensitivity analysis where the proliferation rates of fibroblasts and myofibroblasts are not the same. The impact of this model adaptation is significant. This project is funded by the Dutch Burns Foundation under Project 17.105 The authors are grateful for the fnancial support from the Dutch Burns Foundation under Project 17.105.

To deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order $${{\mathcal {O}}}(h^2)$$ O ( h 2 ) . Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.

We consider a one-dimensional morphoelastic model describing post-burn scar contraction. Contraction can lead to a limited range of motion (contracture). Reported prevalence of burn scar contractures are 58.6% at 3-6 weeks and 20.9% at 12 months post-reconstructive surgery after burns. This model describes the displacement of the dermal layer of the skin and the development of the effective Eulerian strain in the tissue. Besides these components, the model also contains components that play a major role in the skin repair after trauma. These components are signaling molecules, fibroblasts, myofibroblasts, and collagen. We perform a sensitivity analysis for many parameters of the model and use the results for a feasibility study. In this study, we test whether the model is suitable for predicting the extent of contraction in different age groups. To this end, we conduct an extensive literature review to find parameter values. From the sensitivity analysis, we conclude that the most sensitive parameters are the equilibrium collagen concentration in the dermal layer, the apoptosis rate of fibroblasts and myofibroblasts, and the secretion rate of signaling molecules. Further, although we can use the model to simulate significant distinct contraction densities in different age groups, our results differ from what is seen in the clinic. This particularly concerns children and elderly patients. In children we see more intense contractures if the burn injury occurs near a joint, because the growth induces extra forces on the tissue. Elderly patients seem to suffer less from contractures, possibly because of excess skin. The authors are grateful for the fnancial support by the Dutch Burns Foundation under Project 17.105. This project is funded by the Dutch Burns Foundation under Project 17.105.

The existence of waiting times, before boundary motion sets in, for a diffusion–diffusion reaction equation with a discontinuous switch mechanism, is demonstrated. Limit cases of the waiting times are discussed in mathematical rigor. Further, analytic solutions for planar and circular wounds are derived. The waiting times, as predicted using these analytic solutions, are perfectly between the derived bounds. Furthermore, it is demonstrated by both physical reasoning and mathematical rigor that the movement of the boundary can be delayed once it starts moving. The proof of this assertion resides on continuity and monotonicity arguments. The theory sustains the construction of analytic solutions. The model is applied to simulation of biological processes with a threshold behavior, such as wound healing or tumor growth.

A simplified finite-element model for wound healing is proposed. The model takes into account the sequential steps of dermal regeneration, wound contraction, angiogenesis and wound closure. An innovation in the present study is the combination of the aforementioned partially overlapping processes, which can be used to deliver novel insights into the process of wound healing, such as geometry related influences, as well as the influence of coupling between the various existing subprocesses on the actual healing behavior. The model confirms the clinical observation that epidermal closure proceeds by a crawling and climbing mechanism at the early stages, and by a stratification process in layers parallel to the skin surface at the later stages. The local epidermal oxygen content may play an important role here. The model can also be used to investigate the influence of local injection of hormones that stimulate partial processes occurring during wound healing. These insights can be used to improve wound healing treatments.