No Arabic abstract
In this paper we introduce a new mathematical model for the active contraction of cardiac muscle, featuring different thermo-electric and nonlinear conductivity properties. The passive hyperelastic response of the tissue is described by an orthotropic exponential model, whereas the ionic activity dictates active contraction incorporated through the concept of orthotropic active strain. We use a fully incompressible formulation, and the generated strain modifies directly the conductivity mechanisms in the medium through the pull-back transformation. We also investigate the influence of thermo-electric effects in the onset of multiphysics emergent spatiotemporal dynamics, using nonlinear diffusion. It turns out that these ingredients have a key role in reproducing pathological chaotic dynamics such as ventricular fibrillation during inflammatory events, for instance. The specific structure of the governing equations suggests to cast the problem in mixed-primal form and we write it in terms of Kirchhoff stress, displacements, solid pressure, electric potential, activation generation, and ionic variables. We also propose a new mixed-primal finite element method for its numerical approximation, and we use it to explore the properties of the model and to assess the importance of coupling terms, by means of a few computational experiments in 3D.
A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport, and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold, and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models, and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.
Self-regulation of living tissue as an example of self-organization phenomena in active fractal systems of biological, ecological, and social nature is under consideration. The characteristic feature of these systems is the absence of any governing center and, thereby, their self-regulation is based on a cooperative interaction of all the elements. The paper develops a mathematical theory of a vascular network response to local effects on scales of individual units of peripheral circulation. First, it formulates a model for the self-processing of information about the cellular tissue state and cooperative interaction of blood vessels governing redistribution of blood flow over the vascular network. Mass conservation (conservation of blood flow as well as transported biochemical compounds) plays the key role in implementing these processes. The vascular network is considered to be of the tree form and the blood vessels are assumed to respond individually to an activator in blood flowing though them. Second, the constructed governing equations are analyzed numerically. It is shown that at the first approximation the blood perfusion rate depends locally on the activator concentration in the cellular tissue, which is due to the hierarchical structure of the vascular network. Then the distinction between the reaction threshold of individual vessels and that of the vascular network as a whole is demonstrated. In addition, the nonlocal component of the dependence of the blood perfusion rate on the activator concentration is found to change its form as the activator concentration increases.
The growth of several biological tissues is known to be controlled in part by local geometrical features, such as the curvature of the tissue interface. This control leads to changes in tissue shape that in turn can affect the tissues evolution. Understanding the cellular basis of this control is highly significant for bioscaffold tissue engineering, the evolution of bone microarchitecture, wound healing, and tumour growth. While previous models have proposed geometrical relationships between tissue growth and curvature, the role of cell density and cell vigor remains poorly understood. We propose a cell-based mathematical model of tissue growth to investigate the systematic influence of curvature on the collective crowding or spreading of tissue-synthesising cells induced by changes in local tissue surface area during the motion of the interface. Depending on the strength of diffusive damping, the model exhibits complex growth patterns such as undulating motion, efficient smoothing of irregularities, and the generation of cusps. We compare this model with in-vitro experiments of tissue deposition in bioscaffolds of different geometries. By accounting for the depletion of active cells, the model is able to capture both smoothing of initial substrate geometry and tissue deposition slowdown as observed experimentally.
Identifying the mechanism of intercellular feedback regulation is critical for the basic understanding of tissue growth control in organisms. In this paper, we analyze a tissue growth model consisting of a single lineage of two cell types regulated by negative feedback signalling molecules that undergo spatial diffusion. By deriving the fixed points for the uniform steady states and carrying out linear stability analysis, phase diagrams are obtained analytically for arbitrary parameters of the model. Two different generic growth modes are found: blow-up growth and final-state controlled growth which are governed by the non-trivial fixed point and the trivial fixed point respectively, and can be sensitively switched by varying the negative feedback regulation on the proliferation of the stem cells. Analytic expressions for the characteristic time scales for these two growth modes are also derived. Remarkably, the trivial and non-trivial uniform steady states can coexist and a sharp transition occurs in the bistable regime as the relevant parameters are varied. Furthermore, the bi-stable growth properties allows for the external control to switch between these two growth modes. In addition, the condition for an early accelerated growth followed by a retarded growth can be derived. These analytical results are further verified by numerical simulations and provide insights on the growth behavior of the tissue. Our results are also discussed in the light of possible realistic biological experiments and tissue growth control strategy. Furthermore, by external feedback control of the concentration of regulatory molecules, it is possible to achieve a desired growth mode, as demonstrated with an analysis of boosted growth, catch-up growth and the design for the target of a linear growth dynamic.
Several biological tissues undergo changes in their geometry and in their bulk material properties by modelling and remodelling processes. Modelling synthesises tissue in some regions and removes tissue in others. Remodelling overwrites old tissue material properties with newly formed, immature tissue properties. As a result, tissues are made up of different patches, i.e., adjacent tissue regions of different ages and different material properties, within evolving boundaries. In this paper, generalised equations governing the spatio-temporal evolution of such tissues are developed within the continuum model. These equations take into account nonconservative, discontinuous surface mass balance due to creation and destruction of material at moving interfaces, and bulk balance due to tissue maturation. These equations make it possible to model patchy tissue states and their evolution without explicitly maintaining a record of when/where resorption and formation processes occurred. The time evolution of spatially averaged tissue properties is derived systematically by integration. These spatially-averaged equations cannot be written in closed form as they retain traces that tissue destruction is localised at tissue boundaries. The formalism developed in this paper is applied to bone tissues, which exhibit strong material heterogeneities due to their slow mineralisation and remodelling processes. Evolution equations are proposed in particular for osteocyte density and bone mineral density. Effective average equations for bone mineral density (BMD) and tissue mineral density (TMD) are derived using a mean-field approximation. The error made by this approximation when remodelling patchy tissue is investigated. The specific time signatures of BMD or TMD during remodelling events may provide a way to detect these events occurring at lower, unseen spatial resolutions from microCT scans.