Non-Trivial Dynamics in the FizHugh-Rinzel Model and Non-Homogeneous Oscillatory-Excitable Reaction-Diffusions Systems

This article focuses on the qualitative analysis of complex dynamics arising in a few mathematical models in neuroscience context. We first discuss the dynamics arising in the three-dimensional FitzHugh-Rinzel (FHR) model and then illustrate those arising in a class of non-homogeneous FitzHugh-Nagumo (Nh-FHN) reaction-diffusion systems. FHR and Nh-FHN models can be used to generate relevant complex dynamics and wave-propagation phenomena in neuroscience context. Such complex dynamics include canards, mixed-mode oscillations (MMOs), Hopf-bifurcations and their spatially extended counterpart. Our article highlights original methods to characterize these complex dynamics and how they emerge in ordinary differential equations and spatially extended models.

Mathematical models in GnRH research

Mathematical modelling is an indispensable tool in modern biosciences, enabling quantitative analysis and integration of biological data, transparent formulation of our understanding of complex biological systems, and efficient experimental design based on model predictions. This review article provides an overview of the impact that mathematical models had on GnRH research. Indeed, over the last 20 years mathematical modelling has been used to describe and explore the physiology of the GnRH neuron, the mechanisms underlying GnRH pulsatile secretion, and GnRH signalling to the pituitary. Importantly, these models have contributed to GnRH research via novel hypotheses and predictions regarding the bursting behaviour of the GnRH neuron, the role of kisspeptin neurons in the emergence of pulsatile GnRH dynamics, and the decoding of GnRH signals by biochemical signalling networks. We envisage that with the advent of novel experimental technologies, mathematical modelling will have an even greater role to play in our endeavour to understand the complex spatiotemporal dynamics underlying the reproductive neuroendocrine system.

Mathematical modeling approaches of cellular endocrinology within the hypothalamo-pituitary-gonadal axis

The reproductive neuroendocrine axis, or hypothalamo-pituitary-gonadal (HPG) axis, is a paragon of complex biological system involving numerous cell types, spread over several anatomical levels communicating through entangled endocrine feedback loops. The HPG axis exhibits remarkable dynamic behaviors on multiple time and space scales, which are an inexhaustible source of studies for mathematical and computational biology. In this review, we will describe a variety of modeling approaches of the HPG axis from a cellular endocrinology viewpoint. We will in particular investigate the questions raised by some of the most striking features of the HPG axis: (i) the pulsatile secretion of hypothalamic and pituitary hormones, and its counterpart, the cell signaling induced by frequency-encoded hormonal signals, and (ii) the dual, gametogenic and glandular function of the gonads, which relies on the tight control of the somatic cell populations ensuring the proper maturation and timely release of the germ cells.

Mathematical Modelling of Endocrine Systems

Hormone rhythms are ubiquitous and essential to sustain normal physiological functions. Combined mathematical modelling and experimental approaches have shown that these rhythms result from regulatory processes occurring at multiple levels of organisation and require continuous dynamic equilibration, particularly in response to stimuli. We review how such an interdisciplinary approach has been successfully applied to unravel complex regulatory mechanisms in the metabolic, stress, and reproductive axes. We discuss how this strategy is likely to be instrumental for making progress in emerging areas such as chronobiology and network physiology. Ultimately, we envisage that the insight provided by mathematical models could lead to novel experimental tools able to continuously adapt parameters to gradual physiological changes and the design of clinical interventions to restore normal endocrine function.

Perturbing the Hypothalamic-Pituitary-Adrenal Axis: A Mathematical Model for Interpreting PTSD Assessment Tests

We use a dynamical systems model of the hypothalamic-pituitary-adrenal (HPA) axis to understand the mechanisms underlying clinical protocols used to probe patient stress response. Specifically, we address dexamethasone (DEX) and ACTH challenge tests, which probe pituitary and adrenal gland responses, respectively. We show that some previously observed features and experimental responses can arise from a bistable mathematical model containing two steady-states, rather than relying on specific and permanent parameter changes due to physiological disruption. Moreover, we show that the timing of a perturbation relative to the intrinsic oscillation of the HPA axis can affect challenge test responses. Conventional mechanistic hypotheses supported and refuted by the challenge tests are reexamined by varying parameters in our mathematical model associated with these hypotheses. We show that (a) adrenal hyposensitivity can give rise to the responses seen in ACTH challenge tests and (b) enhanced cortisol-mediated suppression of the pituitary in subjects with PTSD is not necessary to explain the responses observed in DEX stress tests. We propose a new two-stage DEX/external stressor protocol to more clearly distinguish between the conventional hypothesis of enhanced suppression of the pituitary and bistable dynamics hypothesized in our model.

The Effects of GABAergic Polarity Changes on Episodic Neural Network Activity in Developing Neural Systems

Early in development, neural systems have primarily excitatory coupling, where even GABAergic synapses are excitatory. Many of these systems exhibit spontaneous episodes of activity that have been characterized through both experimental and computational studies. As development progress the neural system goes through many changes, including synaptic remodeling, intrinsic plasticity in the ion channel expression, and a transformation of GABAergic synapses from excitatory to inhibitory. What effect each of these, and other, changes have on the network behavior is hard to know from experimental studies since they all happen in parallel. One advantage of a computational approach is that one has the ability to study developmental changes in isolation. Here, we examine the effects of GABAergic synapse polarity change on the spontaneous activity of both a mean field and a neural network model that has both glutamatergic and GABAergic coupling, representative of a developing neural network. We find some intuitive behavioral changes as the GABAergic neurons go from excitatory to inhibitory, shared by both models, such as a decrease in the duration of episodes. We also find some paradoxical changes in the activity that are only present in the neural network model. In particular, we find that during early development the inter-episode durations become longer on average, while later in development they become shorter. In addressing this unexpected finding, we uncover a priming effect that is particularly important for a small subset of neurons, called the “intermediate neurons.” We characterize these neurons and demonstrate why they are crucial to episode initiation, and why the paradoxical behavioral change result from priming of these neurons. The study illustrates how even arguably the simplest of developmental changes that occurs in neural systems can present non-intuitive behaviors. It also makes predictions about neural network behavioral changes that occur during development that may be observable even in actual neural systems where these changes are convoluted with changes in synaptic connectivity and intrinsic neural plasticity.