No Arabic abstract
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.
In this paper, we derive an effective macroscale description suitable to describe the growth of biological tissue within a porous tissue-engineering scaffold. As in our recent work (Holden textit{et al.} A multiphase multiscale model for nutrient limited tissue growth, The ANZIAM Journal, 2018, doi:10.1017/S1446181118000044) the underlying tissue dynamics is described as a multiphase mixture, thereby naturally accommodating features such as interstitial growth and active cell motion. Via a linearisation of the underlying multiphase model (whose nonlinearity poses significant challenge for such analyses), we obtain, by means of multiple-scales homogenisation, a simplified macroscale model that nevertheless retains explicit dependence on both the microscale scaffold structure and the tissue dynamics. The model we obtain comprises Darcy flow, and differential equations for the volume fraction of cells within the scaffold and the concentration of nutrient, required for growth. These are coupled to underlying Stokes-type cell problems that provide permeability tensors to parameterise the macroscale description. In Holden textit{et al.}, the cell problems retain macroscale dependence, posing significant computational challenges; here, we obtain a decoupled system whereby the quasi-steady cell-problems may be solved separately from the macroscale description, thereby greatly reducing the complexity associated with fully-coupled multiscale descriptions. Moreover, we indicate how the formulation is influenced by a set of alternative microscale boundary conditions.S
To maintain bone mass during bone remodelling, coupling is required between bone resorption and bone formation. This coordination is achieved by a network of autocrine and paracrine signalling molecules between cells of the osteoclast lineage and cells of the osteoblastic lineage. Mathematical modelling of signalling between cells of both lineages can assist in the interpretation of experimental data, clarify signalling interactions and help develop a deeper understanding of complex bone diseases. In this paper, we further develop a mathematical model of bone cell interactions by Pivonka et al. (2008) to include the proliferation of precursor osteoblasts into the model. This inclusion is important to be able to account for Wnt signalling, believed to play an important role in anabolic responses of bone. We show that an increased rate of differentiation to precursor cells or an increased rate of proliferation of precursor osteoblasts themselves both result in increased bone mass. However, modelling these different processes separately enables the new model to represent recent experimental discoveries such as the role of Wnt signalling in bone biology and the recruitment of osteoblast progenitor cells by transforming growth factor beta. Finally, we illustrate the power of the new models capabilities by applying the model to prostate cancer metastasis to bone. In the bone microenvironment, prostate cancer cells are believed to release some of the same signalling molecules used to coordinate bone remodelling (i.e. Wnt and PTHrP), enabling the cancer cells to disrupt normal signalling and coordination between bone cells. This disruption can lead to either bone gain or bone loss. We demonstrate that the new computational model developed here is capable of capturing some key observations made on the evolution of the bone mass due to metastasis of prostate cancer to the bone microenvironment
Bone remodelling maintains the functionality of skeletal tissue by locally coordinating bone-resorbing cells (osteoclasts) and bone-forming cells (osteoblasts) in the form of Bone Multicellular Units (BMUs). Understanding the emergence of such structured units out of the complex network of biochemical interactions between bone cells is essential to extend our fundamental knowledge of normal bone physiology and its disorders. To this end, we propose a spatio-temporal continuum model that integrates some of the most important interaction pathways currently known to exist between cells of the osteoblastic and osteoclastic lineage. This mathematical model allows us to test the significance and completeness of these pathways based on their ability to reproduce the spatio-temporal dynamics of individual BMUs. We show that under suitable conditions, the experimentally-observed structured cell distribution of cortical BMUs is retrieved. The proposed model admits travelling-wave-like solutions for the cell densities with tightly organised profiles, corresponding to the progression of a single remodelling BMU. The shapes of these spatial profiles within the travelling structure can be linked to the intrinsic parameters of the model such as differentiation and apoptosis rates for bone cells. In addition to the cell distribution, the spatial distribution of regulatory factors can also be calculated. This provides new insights on how different regulatory factors exert their action on bone cells leading to cellular spatial and temporal segregation, and functional coordination.
The primary exchange units in the human placenta are terminal villi, in which fetal capillary networks are surrounded by a thin layer of villous tissue, separating fetal from maternal blood. To understand how the complex spatial structure of villi influences their function, we use an image-based theoretical model to study the effect of tissue metabolism on the transport of solutes from maternal blood into the fetal circulation. For solute that is taken up under first-order kinetics, we show that the transition between flow-limited and diffusion-limited transport depends on two new dimensionless parameters defined in terms of key geometric quantities, with strong solute uptake promoting flow-limited transport conditions. We present a simple algebraic approximation for solute uptake rate as a function of flow conditions, metabolic rate and villous geometry. For oxygen, accounting for nonlinear kinetics using physiological parameter values, our model predicts that villous metabolism does not significantly impact oxygen transfer to fetal blood, although the partitioning of fluxes between the villous tissue and the capillary network depends strongly on the flow regime.
This paper proposes a model for the growth two interacting populations of cells that do not mix. The dynamics is driven by pressure and cohesion forces on the one hand and proliferation on the other hand. Following earlier works on the single population case, we show that the model approximates a free boundary Hele Shaw type model that we characterise using both analytical and numerical arguments.