  We present a computational model to study the spatiotemporal dynamics of the
epidermis homeostasis 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 homeostatic 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 the psoriasis as a parallel homeostasis of
normal and psoriatic keratinocytes originated from a shared stem-cell niche
environment and predicts two homeostatic modes of the psoriasis: a disease mode
and a quiescent mode. Interconversion between the two modes can be controlled
by interactions between psoriatic stem cells and the immune system and by the
normal and psoriatic stem cells competing for growth niches. The prediction of
a quiescent state potentially explains the efficacy of the 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.
