Abstract:
Gonadotropin-releasing hormone (GnRH) neurons are cells in the hypothalamus that produce GnRH, one of the major hormones that controls fertility and reproduction. However, despite their importance, little is known about the mechanisms by which GnRH is produced. GnRH neurons exhibit complicated membrane potential dynamics, in the form of electrical bursting, and this bursting is closely coupled to the dynamics of intracellular calcium (Ca²⁺ ) in ways that are not yet well understood. A mathematical model has been constructed to help understand the mechanisms underlying the observed behaviours of GnRH neurons, and how electrical bursting synchronizes with transients in the cytosolic Ca²⁺ concentration ([Ca²⁺]cyt). Simulations show that the model is consistent with all of the most important experimental data in- volving TTX (tetrodotoxin), zero extracellular-[ Ca²⁺] ([Ca²⁺]e) solution, apamin, CPA (cyclopiazonic acid) and 2-APB (2-aminoethoxydiphenyl borate). The mathematical model predicted the existence of particular [Ca²⁺]cyt-activated potassium (K⁺) channels (sIAHP-UCL), which were then confirmed experimentally. In contrast to the apamin-sensitive [Ca²⁺]cyt-activated K⁺ channels (sIAHP-SK), which control both the structure of firing within bursts and the interburst intervals, sIAHP-UCL solely determines the interburst dynamics. Furthermore, we show how a fast-slow bifurcation analysis of the model can be used to understand the behaviour of the cells under a wide range of experimental conditions. The high complexity of the initial GnRH neuron model prompted us to develop a much simpler bursting model based on the FitzHugh-Nagumo (FHN) model. This simplified bursting model exhibits the similar qualitative behaviour as our initial GnRH neuron model. The use of bifurcation analysis with the FHN model simplifies the simulation of bursting, and the results of this analysis are also briefly explained.