Complex oscillations in the delayed Fitzhugh-Nagumo equation

Motivated by the dynamics of neuronal responses, we analyze the dynamics of the Fitzhugh-Nagumo slow-fast system with delayed self-coupling. This system provides a canonical example of a canard explosion for sufficiently small delays. Beyond this regime, delays significantly enrich the dynamics, leading to mixed-mode oscillations, bursting and chaos. These behaviors emerge from a delay-induced subcritical Bogdanov-Takens instability arising at the fold points of the S-shaped critical manifold. Underlying the transition from canard-induced to delay-induced dynamics is an abrupt switch in the nature of the Hopf bifurcation

A Mechanism of Co-Existence of Bursting and Silent Regimes of Activities of a Neuron

The co-existence of bursting activity and silence is a common property of various neuronal models. We describe a novel mechanism explaining the co-existence of and the transition between these two regimes. It is based on the specific homoclinic and Andronov-Hopf bifurcations of the hyper- and depolarized steady states that determine the co-existence domain in the parameter space of the leech heart interneuron models: canonical and simplified. We found that a sub-critical Andronov-Hopf bifurcation of the hyperpolarized steady state gives rise to small amplitude sub-threshold oscillations terminating through the secondary homoclinic bifurcation. Near the corresponding boundary the system can exhibit long transition from bursting oscillations into silence, as well as the bi-stability where the observed regime is determined by the initial state of the neuron. The mechanism found is shown to be generic for the simplified 4D and the original 14D leech heart interneuron models

Complex Dynamics in Dedicated / Multifunctional Neural Networks and Chaotic Nonlinear Systems

We study complex behaviors arising in neuroscience and other nonlinear systems by combining dynamical systems analysis with modern computational approaches including GPU parallelization and unsupervised machine learning. To gain insights into the behaviors of brain networks and complex central pattern generators (CPGs), it is important to understand the dynamical principles regulating individual neurons as well as the basic structural and functional building blocks of neural networks. In the first section, we discuss how symbolic methods can help us analyze neural dynamics such as bursting, tonic spiking and chaotic mixed-mode oscillations in various models of individual neurons, the bifurcations that underlie transitions between activity types, as well as emergent network phenomena through synergistic interactions seen in realistic neural circuits, such as network bursting from non-intrinsic bursters. The second section is focused on the origin and coexistence of multistable rhythms in oscillatory neural networks of inhibitory coupled cells. We discuss how network connectivity and intrinsic properties of the cells affect the dynamics, and how even simple circuits can exhibit a variety of mono/multi-stable rhythms including pacemakers, half-center oscillators, multiple traveling-waves, fully synchronous states, as well as various chimeras. Our analyses can help generate verifiable hypotheses for neurophysiological experiments on central pattern generators. In the last section, we demonstrate the inter-disciplinary nature of this research through the applications of these techniques to identify the universal principles governing both simple and complex dynamics, and chaotic structure in diverse nonlinear systems. Using a classical example from nonlinear laser optics, we elaborate on the multiplicity and self-similarity of key organizing structures in 2D parameter space such as homoclinic and heteroclinic bifurcation curves, Bykov T-point spirals, and inclination flips. This is followed by detailed computational reconstructions of the spatial organization and 3D embedding of bifurcation surfaces, parametric saddles, and isolated closed curves (isolas). The generality of our modeling approaches could lead to novel methodologies and nonlinear science applications in biological, medical and engineering systems

Codimension-two homoclinic bifurcations underlying spike adding in the Hindmarsh-Rose burster

The well-studied Hindmarsh-Rose model of neural action potential is revisited from the point of view of global bifurcation analysis. This slow-fast system of three paremeterised differential equations is arguably the simplest reduction of Hodgkin-Huxley models capable of exhibiting all qualitatively important distinct kinds of spiking and bursting behaviour. First, keeping the singular perturbation parameter fixed, a comprehensive two-parameter bifurcation diagram is computed by brute force. Of particular concern is the parameter regime where lobe-shaped regions of irregular bursting undergo a transition to stripe-shaped regions of periodic bursting. The boundary of each stripe represents a fold bifurcation that causes a smooth spike-adding transition where the number of spikes in each burst is increased by one. Next, numerical continuation studies reveal that the global structure is organised by various curves of homoclinic bifurcations. In particular the lobe to stripe transition is organised by a sequence of codimension-two orbit- and inclination-flip points that occur along {\em each} homoclinic branch. Each branch undergoes a sharp turning point and hence approximately has a double-cover of the same curve in parameter space. The sharp turn is explained in terms of the interaction between a two-dimensional unstable manifold and a one-dimensional slow manifold in the singular limit. Finally, a new local analysis is undertaken using approximate Poincar\'{e} maps to show that the turning point on each homoclinic branch in turn induces an inclination flip that gives birth to the fold curve that organises the spike-adding transition. Implications of this mechanism for explaining spike-adding behaviour in other excitable systems are discussed.Comment: 32 pages, 18 figures, submitted to SIAM Journal on Applied Dynamical System

Parameter-sweeping techniques for temporal dynamics of neuronal systems: case study of Hindmarsh-Rose model

Background: Development of effective and plausible numerical tools is an imperative task for thorough studies of nonlinear dynamics in life science applications. Results: We have developed a complementary suite of computational tools for twoparameter screening of dynamics in neuronal models. We test a ‘brute-force’ effectiveness of neuroscience plausible techniques specifically tailored for the examination of temporal characteristics, such duty cycle of bursting, interspike interval, spike number deviation in the phenomenological Hindmarsh-Rose model of a bursting neuron and compare the results obtained by calculus-based tools for evaluations of an entire spectrum of Lyapunov exponents broadly employed in studies of nonlinear systems. Conclusions: We have found that the results obtained either way agree exceptionally well, and can identify and differentiate between various fine structures of complex dynamics and underlying global bifurcations in this exemplary model. Our future planes are to enhance the applicability of this computational suite for understanding of polyrhythmic bursting patterns and their functional transformations in small networks

Flipping Biological Switches: Solving for Optimal Control: A Dissertation

Switches play an important regulatory role at all levels of biology, from molecular switches triggering signaling cascades to cellular switches regulating cell maturation and apoptosis. Medical therapies are often designed to toggle a system from one state to another, achieving a specified health outcome. For instance, small doses of subpathologic viruses activate the immune system’s production of antibodies. Electrical stimulation revert cardiac arrhythmias back to normal sinus rhythm. In all of these examples, a major challenge is finding the optimal stimulus waveform necessary to cause the switch to flip. This thesis develops, validates, and applies a novel model-independent stochastic algorithm, the Extrema Distortion Algorithm (EDA), towards finding the optimal stimulus. We validate the EDA’s performance for the Hodgkin-Huxley model (an empirically validated ionic model of neuronal excitability), the FitzHugh-Nagumo model (an abstract model applied to a wide range of biological systems that that exhibit an oscillatory state and a quiescent state), and the genetic toggle switch (a model of bistable gene expression). We show that the EDA is able to not only find the optimal solution, but also in some cases excel beyond the traditional analytic approaches. Finally, we have computed novel optimal stimulus waveforms for aborting epileptic seizures using the EDA in cellular and network models of epilepsy. This work represents a first step in developing a new class of adaptive algorithms and devices that flip biological switches, revealing basic mechanistic insights and therapeutic applications for a broad range of disorders