Whereas development of mathematical models describing the endocrine system as a whole remains a challenging problem, visible progress has been demonstrated in modeling its subsystems, or axes. Models of hormonal axes portray only the most essential interactions between the hormones and can be described by low-order systems of differential equations. This paper analyzes the properties of a novel model of a hypothalamic-pituitary axis, portraying the interactions in a chain of a release hormone (secreted by the hypothalamus), a tropic hormone (produced by the pituitary gland) and an effector hormone (secreted by a target gland). This model, unlike previously published ones, captures two prominent features of neurohormonal regulation systems, namely, the pulsatile (episodic) production of the release hormone and a complex non-cyclic feedback mechanism that maintains the involved hormone concentrations within certain biological limits. At the same time, the discussed model is analytically tractable; in particular, the existence of a so-called 1-cycle featured by a single pulse over one period is proven mathematically.

The 2017 Nobel Prize in Physiology or Medicine awarded for discoveries of molecular mechanisms controlling the circadian rhythm has called attention to the challenging area of nonlinear dynamics that deals with synchronization and entrainment of oscillations. Biological circadian clocks keep time in living organisms, orchestrating hormonal cycles and other periodic rhythms. The periodic oscillations of circadian pacemakers are self-sustained; at the same time, they are entrainable by external periodic signals that adjust characteristics of autonomous oscillations. Whereas modeling of biological oscillators is a well-established research topic, mathematical analysis of entrainment, i.e. the nonlinear phenomena imposed by periodic exogenous signals, remains an open problem. Along with sustained periodic rhythms, periodically forced oscillators can exhibit various “irregular” behaviors, such as quasiperiodic or chaotic trajectories. This paper presents an overview of the mathematical models of circadian rhythm with respect to endocrine regulation, as well as biological background. Dynamics of the human endocrine system, comprising numerous glands and hormones operating under neural control, are highly complex. Therefore, only endocrine subsystems (or axes) supporting certain biological functions are usually studied. Low-order dynamical models that capture the essential characteristics and interactions between a few hormones can than be derived. Goodwin's oscillator often serves as such a model and is widely regarded as a prototypical biological oscillator. A comparative analysis of forced dynamics arising in two versions of Goodwin's oscillator is provided in the present paper: the classical continuous oscillator and a more recent impulsive one, capturing e.g. pulsatile secretion of hormones due to neural regulation. The main finding of this study is that, while the continuous oscillator is always forced to a periodic solution by a sufficiently large exogenous signal amplitude, the impulsive one commonly exhibits a quasiperiodic or chaotic behavior due to non-smooth dynamics in entrainment.