HPA axis differential flatness and Liouvillian study for higher resiliency investigations

In this paper, we study several existing quantitative models of the hypothalamic–pituitary–adrenal (HPA) axis from a control systems theory viewpoint, that is, we suppose that we can act on the dynamics of the HPA axis throughout some parameters, which are the system inputs. In particular, we will focus on flatness and Liouvillian properties of the considered control systems of the HPA axis. We first study the minimal three-dimensional model of Bangsgaard and Ottesen (2017, Math. Biosci., 287:24–35) and the semi-mechanistic four-dimensional model of Gupta et al. (2007, Theor. Biol. Medical Model., 4(1):8) which are shown to be flat, and then, we consider the more involved and important model proposed in Rao & Androulakis (2019, Sci. Rep., 9(1):11212; 2020, IFAC-PapersOnLine, 53(2):15858–15863), with seven states, for which we prove that for the nominal values of the parameters involved in the model, flatness no longer holds. The more involved model satisfies however a similar but weaker property than flatness: it is a Liouvillian system.