A modelling approach towards Epidermal homoeostasis control

Gernot Schaller,Michael Meyer-Hermann
DOI: https://doi.org/10.1016/j.jtbi.2007.03.023
2007-04-11
Abstract:In order to grasp the features arising from cellular discreteness and individuality, in large parts of cell tissue modelling agent-based models are favoured. The subclass of off-lattice models allows for a physical motivation of the intercellular interaction rules. We apply an improved version of a previously introduced off-lattice agent-based model to the steady-state flow equilibrium of skin. The dynamics of cells is determined by conservative and drag forces,supplemented with delta-correlated random forces. Cellular adjacency is detected by a weighted Delaunay triangulation. The cell cycle time of keratinocytes is controlled by a diffusible substance provided by the dermis. Its concentration is calculated from a diffusion equation with time-dependent boundary conditions and varying diffusion coefficients. The dynamics of a nutrient is also taken into account by a reaction-diffusion equation. It turns out that the analysed control mechanism suffices to explain several characteristics of epidermal homoeostasis formation. In addition, we examine the question of how {\em in silico} melanoma with decreased basal adhesion manage to persist within the steady-state flow-equilibrium of the <a class="link-external link-http" href="http://skin.Interestingly" rel="external noopener nofollow">this http URL</a>, even for melanocyte cell cycle times being substantially shorter than for keratinocytes, tiny stochastic effects can lead to completely different outcomes. The results demonstrate that the understanding of initial states of tumour growth can profit significantly from the application of off-lattice agent-based models in computer simulations.
Biological Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to understand and simulate the formation mechanism of epidermal homoeostasis. Specifically, the author studies the flow equilibrium of the skin in a stable state by improving the off - lattice agent - based model. The following are the main issues involved in the paper: 1. **Cell discreteness and individuality**: In order to capture the characteristics caused by cell discreteness and individuality, the author uses an agent - based model. This type of model allows for physically - motivated inter - cellular interaction rules. 2. **Control mechanism of epidermal homoeostasis**: The author attempts to explain multiple characteristics of the formation of epidermal homoeostasis, including how the cell - cycle time of keratinocytes is regulated by diffusible substances provided by the dermis. The concentration of this substance is calculated by a diffusion equation with time - varying boundary conditions and different diffusion coefficients. 3. **Dynamics of nutrients**: The dynamic changes of nutrients are considered through reaction - diffusion equations. 4. **Simulation of melanoma**: Study how melanoma cells with reduced basal adhesion can persist under steady - state flow equilibrium. In particular, it is explored that even if the cell - cycle time of melanocytes is significantly shorter than that of keratinocytes, small random effects may lead to completely different results. 5. **Understanding the initial stage of tumor growth**: Improve the understanding of the initial state of tumor growth by applying the off - lattice agent - based model through computer simulation. ### Formula summary 1. **Inter - cellular interaction**: - Use the Johnson - Kendall - Roberts (JKR) model to describe the elastic, adhesive and viscous non - isotropic damping forces between cells: \[ A(t)\dot{\mathbf{x}}=\mathbf{b} \] where \(\mathbf{x}\) is the position vector of \(N\) cells, \(A(t)\) is the friction matrix, and \(\mathbf{b}(t)\) contains all non - viscous forces. 2. **Total cell binding strength**: - The total cell binding strength \(\Sigma_i(t)\) can be expressed as: \[ \Sigma_i(t)=\sum_{j\in NN(i)}\sigma_{ij}(t)A_{ij}(t) \] where \(\sigma_{ij}\) is the adhesive energy density between cells \(i\) and \(j\), and \(A_{ij}\) is the contact area. 3. **Change in cell elastic modulus**: - When a cell enters the keratinization or necrosis state, the elastic modulus \(E_i\) changes with time: \[ \dot{E}_i = - 2\alpha E_i \] 4. **Cell volume growth**: - In the G1 phase, the cell volume grows at a constant rate \(r_V\): \[ \dot{R}=\frac{r_V}{4\pi R^2} \] 5. **Effect of local mobile water concentration**: - If the local mobile water concentration is lower than the critical value \(U_{\text{crit}}\), the cell will directly continue the cell cycle; otherwise, the cell will enter the G0 phase. Through these models and formulas, the author aims to better understand the formation mechanism of epidermal homoeostasis and explore the behavior of melanoma cells in the early stage.