- Poster Presentation
- Open Access
Mixed mode oscillations in a gonadotropin-releasing hormone (GnRH) neuron model
BMC Neuroscience volume 11, Article number: P161 (2010)
Mixed mode oscillations (MMOs) are complex patterns resulting from an inter-mixing of large- and small-amplitude oscillations. Such patterns have been widely observed in experiments as well as mathematical models  and are currently an active area of mathematical research. In neurons MMOs often consist of a mixture of spikes and subthreshold oscillations. Our GnRH work leads to a seven-dimensional soma model aimed at investigating the firing properties of GnRH neurons across puberty. The results of our model are supported by published electrophysiological data of Liu and Herbison . The model exhibits robust subthreshold oscillations in addition to action potentials, an observation validated by experimental evidence. Previously identified mechanisms giving rise to MMOs in mathematical models include canard-based mechanisms such as folded node singularities and singular Hopf bifurcations [3, 4]. MMOs arising out of these mechanisms have characteristic patterns of variation in oscillation amplitudes; in the case of folded nodes there are also bounds for the maximum number of subthreshold oscillations per cycle. The MMOs in our system do not follow these patterns but are instead “exotic” MMOs. In this work we use geometric singular perturbation theory to analyze the source of these MMOs in the model and compare to known mechanisms; we also explore the possible physiological function of MMOs in GnRH cells.
Bifurcation analysis of our system was performed with the help of the simulation and continuation tool XPPAUT. Our analysis suggests that MMOs in our model persist near a saddle-node of periodic orbits and the geometry further points to a tangency between an unstable manifold and a family of periodic orbits. In order to facilitate a more detailed analysis of this phenomenon we reduced to a three dimensional fast-slow system with two slow variables that exhibits qualitatively similar MMOs. We note here that MMOs in three or higher dimensional systems have been observed to occur robustly when two of the variables evolve on a much slower time-scale than the rest [3–5] unlike in two dimensional systems where they occur over an exponentially small parameter interval. The bifurcation analysis of our reduced model provides insights into the geometry of stable and unstable invariant manifolds of our system and their role in the generation of the exotic MMOs.
Brons M, Kaper TJ, Rotstein HG: Introduction to focus issue: mixed mode oscillations: experiment, computation, and analysis. Chaos. 2008, 18: 015101-0151014. 10.1063/1.2903177.
Liu X, Herbison AE: Small-conductance calcium-activated potassium channels control excitability and firing dynamics in gonadotropin-releasing hormone (GnRH) neurons. Endocrinology. 2008, 149 (7): 3598-3604. 10.1210/en.2007-1631.
Wechselberger M: Existence and bifurcation of canards in R3 in the case of a folded node. SIAM Journal of Applied Dynamical Systems. 2005, 4 (1): 101-139. 10.1137/030601995.
Guckenheimer J: Singular hopf bifurcation in systems with two slow variables. SIAM Journal of Applied Dynamical Systems. 2008, 7 (4): 1355-1377. 10.1137/080718528.
Rubin J, Wechselberger M: Giant squid-hidden canard: the 3d geometry of the hodgkin-huxley model. Biological Cybernetics. 2007, 97: 5-32. 10.1007/s00422-007-0153-5.
This work was supported by HD-060818.
About this article
Cite this article
Banerjee, S., Best, J. & Suter, K. Mixed mode oscillations in a gonadotropin-releasing hormone (GnRH) neuron model. BMC Neurosci 11, P161 (2010) doi:10.1186/1471-2202-11-S1-P161
- Periodic Orbit
- Hopf Bifurcation
- Invariant Manifold
- Unstable Manifold
- Bifurcation Analysis