An age-structured model of epidermis growth

The force of cell-cell adhesion in determining the outcome in a nonlocal advection diffusion model of wound healing

A model of wound healing is presented to investigate the connection of the force of cell-cell adhesion to the sensing radius of cells in their spatial environment. The model consists of a partial differential equation with nonlocal advection and diffusion terms, describing the movement of cells in a spatial environment. The model is applied to biological wound healing experiments to understand incomplete wound closure. The analysis demonstrates that for each value of the force of adhesion parameter, there is a critical value of the sensing radius above which complete wound healing does not occur.

Stability analysis of an age-structured model of cervical cancer cells and HPV dynamics

Stability analysis of an autonomous epidemic model of an age-structured sub-populations of susceptible, infected, precancerous and cancer cells and unstructured sub-population of human papilloma virus (HPV) (SIPCV epidemic model) aims to gain an insight into the features of cervical cancer disease. The model considers the immune functional response of organism to the virus population growing by the HPV-density dependent death rate, while the death rates of infected, precancerous and cancerous cells do not depend on the HPV quantity because the immune system of organism does not respond to its own cells. Interaction between susceptible cells and HPV is described by the Lotka-Voltera incidence rate and leads to the growth of infected cells. Some of infected cells become precancerous cells, and the other apoptosis when viruses leave infected cells and are ready to infect new susceptible cells. Precancerous cells partially become cancer cells with the density-dependent saturated rate. Conditions of existence of the endemic equilibrium of system were obtained. It was proved that this equilibrium is always locally asymptotically stable whenever it exists. We obtained: (i) the conditions of cancer tumor localization (asymptotically stable dynamical regimes), (ii) outbreak of cancer cell population (that may correspond to metastasis), (iii) outbreak of dysplasia (precancerous cells) which induces the outbreak of cancer cells (that may correspond to metastasis). In cases (ii), (iii) the conditions of existence of endemic equilibrium do not hold. All cases are illustrated by numerical experiments.

Stability Analysis of an Age-Structured SEIRS Model with Time Delay

This paper is concerned with the stability of an age-structured susceptible–exposed– infective–recovered–susceptible (SEIRS) model with time delay. Firstly, the traveling wave solution of system can be obtained by using the method of characteristic. The existence and uniqueness of the continuous traveling wave solution is investigated under some hypotheses. Moreover, the age-structured SEIRS system is reduced to the nonlinear autonomous system of delay ODE using some insignificant simplifications. It is studied that the dimensionless indexes for the existence of one disease-free equilibrium point and one endemic equilibrium point of the model. Furthermore, the local stability for the disease-free equilibrium point and the endemic equilibrium point of the infection-induced disease model is established. Finally, some numerical simulations were carried out to illustrate our theoretical results.

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.

A Dynamic Model for Prediction of Psoriasis Management by Blue Light Irradiation

Clinical investigations prove that blue light irradiation reduces the severity of psoriasis vulgaris. Nevertheless, the mechanisms involved in the management of this condition remain poorly defined. Despite the encouraging results of the clinical studies, no clear guidelines are specified in the literature for the irradiation scheme regime of blue light-based therapy for psoriasis. We investigated the underlying mechanism of blue light irradiation of psoriatic skin, and tested the hypothesis that regulation of proliferation is a key process. We implemented a mechanistic model of cellular epidermal dynamics to analyze whether a temporary decrease of keratinocytes hyper-proliferation can explain the outcome of phototherapy with blue light. Our results suggest that the main effect of blue light on keratinocytes impacts the proliferative cells. They show that the decrease in the keratinocytes proliferative capacity is sufficient to induce a transient decrease in the severity of psoriasis. To study the impact of the therapeutic regime on the efficacy of psoriasis treatment, we performed simulations for different combinations of the treatment parameters, i.e., length of treatment, fluence (also referred to as dose), and intensity. These simulations indicate that high efficacy is achieved by regimes with long duration and high fluence levels, regardless of the chosen intensity. Our modeling approach constitutes a framework for testing diverse hypotheses on the underlying mechanism of blue light-based phototherapy, and for designing effective strategies for the treatment of psoriasis.

The steady state of epidermis: mathematical modeling and numerical simulations

We consider a model with age and space structure for the epidermis evolution. The model, previously presented and analyzed with respect to the suprabasal epidermis, includes different types of cells (proliferating cells, differentiated cells, corneous cells, and apoptotic cells) moving with the same velocity, under the constraint that the local volume fraction occupied by the cells is constant in space and time. Here, we complete the model proposing a mechanism regulating the cell production in the basal layer and we focus on the stationary case of the problem, i.e. on the case corresponding to the normal status of the skin. A numerical scheme to compute the solution of the model is proposed and its convergence is studied. Simulations are provided for realistic values of the parameters, showing the possibility of reproducing the structure of both "thin" and "thick" epidermis.

A glycemia-structured population model

Structured models are population models in which the individuals are characterized with respect to the value of some variable of interest, called the structure variable. In the present paper, we propose a glycemia-structured population model, based on a linear partial differential equation with variable coefficients. The model is characterized by three rate functions: a new-adult population glycemic profile, a glycemia-dependent mortality rate and a glycemia-dependent average worsening rate. First, we formally analyze some properties of the solution, the transient behavior and the equilibrium distribution. Then, we identify the key parameters and functions of the model from real-life data and we hypothesize some plausible modifications of the rate functions to obtain a more beneficial steady-state behavior. The interest of the model is that, while it summarizes the evolution of diabetes in the population in a completely different way with respect to previously published Monte Carlo aggregations of individual-based models, it does appear to offer a good approximation of observed reality and of the features expected in the clinical setting. The model can offer insights in pharmaceutical research and be used to assess possible public health intervention strategies.

Modelling epidermis homoeostasis and psoriasis pathogenesis

We present a computational model to study the spatio-temporal dynamics of epidermis homoeostasis under normal and pathological conditions. The model consists of a population kinetics model of the central transition pathway of keratinocyte proliferation, differentiation and loss and an agent-based model that propagates cell movements and generates the stratified epidermis. The model recapitulates observed homoeostatic cell density distribution, the epidermal turnover time and the multilayered tissue structure. We extend the model to study the onset, recurrence and phototherapy-induced remission of psoriasis. The model considers psoriasis as a parallel homoeostasis of normal and psoriatic keratinocytes originated from a shared stem cell (SC) niche environment and predicts two homoeostatic modes of psoriasis: a disease mode and a quiescent mode. Interconversion between the two modes can be controlled by interactions between psoriatic SCs and the immune system and by normal and psoriatic SCs competing for growth niches. The prediction of a quiescent state potentially explains the efficacy of multi-episode UVB irradiation therapy and recurrence of psoriasis plaques, which can further guide designs of therapeutics that specifically target the immune system and/or the keratinocytes.

Towards a quantitative theory of epidermal calcium profile formation in unwounded skin

We propose and mathematically examine a theory of calcium profile formation in unwounded mammalian epidermis based on: changes in keratinocyte proliferation, fluid and calcium exchange with the extracellular fluid during these cells’ passage through the epidermal sublayers, and the barrier functions of both the stratum corneum and tight junctions localised in the stratum granulosum. Using this theory, we develop a mathematical model that predicts epidermal sublayer transit times, partitioning of the epidermal calcium gradient between intracellular and extracellular domains, and the permeability of the tight junction barrier to calcium ions. Comparison of our model’s predictions of epidermal transit times with experimental data indicates that keratinocytes lose at least 87% of their volume during their disintegration to become corneocytes. Intracellular calcium is suggested as the main contributor to the epidermal calcium gradient, with its distribution actively regulated by a phenotypic switch in calcium exchange between keratinocytes and extracellular fluid present at the boundary between the stratum spinosum and the stratum granulosum. Formation of the extracellular calcium distribution, which rises in concentration through the stratum granulosum towards the skin surface, is attributed to a tight junction barrier in this sublayer possessing permeability to calcium ions that is less than 15 nm s⁻¹ in human epidermis and less than 37 nm s⁻¹ in murine epidermis. Future experimental work may refine the presented theory and reduce the mathematical uncertainty present in the model predictions.

Active regulation of the epidermal calcium profile

A distinct calcium profile is strongly implicated in regulating the multi-layered structure of the epidermis. However, the mechanisms that govern the regulation of this calcium profile are currently unclear. It clearly depends on the relatively impermeable barrier of the stratum corneum (passive regulation) but may also depend on calcium exchanges between keratinocytes and extracellular fluid (active regulation). Using a mathematical model that treats the viable sublayers of unwounded human and murine epidermis as porous media and assumes that their calcium profiles are passively regulated, we demonstrate that these profiles are also actively regulated. To obtain this result, we found that diffusion governs extracellular calcium motion in the viable epidermis and hence intracellular calcium is the main source of the epidermal calcium profile. Then, by comparison with experimental calcium profiles and combination with a hypothesised cell velocity distribution in the viable epidermis, we found that the net influx of calcium ions into keratinocytes from extracellular fluid may be constant and positive throughout the stratum basale and stratum spinosum, and that there is a net outflux of these ions in the stratum granulosum. Hence, the calcium exchange between keratinocytes and extracellular fluid differs distinctly between the stratum granulosum and the underlying sublayers, and these differences actively regulate the epidermal calcium profile. Our results also indicate that plasma membrane dysfunction may be an early event during keratinocyte disintegration in the stratum granulosum.