Impact of Chloride Channel on Spiking Patterns of Morris-Lecar Model

In this paper, we study the complicated dynamics of general Morris-Lecar model with the impact of Cl − fluctuations on firing patterns of this neuron model. After adding Cl − channel in the original Morris-Lecar model, the dynamics of the original model such as its bifurcations of equilibrium points would be changed and they occurred at different values compared to the primary model. We discover these qualitative changes in the point of dynamical systems and neuroscience. We will conduct the co-dimension two bifurcations analysis with respect to different control parameters to explore the complicated behaviors for this new neuron model.


Introduction
Computational neuroscience uses knowledge of biology and combines with mathematical modeling to simulate the fluctuations of neuron cells and biological physiological characteristics of them [1] [2] [3] [4]. The progress and advances in computational neuroscience could help scientists to better understand the performance of brain and neuron cells and better fight with diseases related to neuron cells such as Parkinson's and depression [5]. Due to the complexity of nerve systems, it is impossible to fully understand the various phenomena in neuroscience only using the simple linear modeling methods, since it does not meet all neuronal properties and complicated behaviors of neurons. Thus, we need to use nonlinear methods and models which catch all of the properties of neurons using different experimental data which brings biologist and specifically neuroscientists and applied science people all together to develop appropriate models to explore all these complex dynamics of neurons [1] [5] [6] [7] [8].
The Morris-Lecar model is a two-dimensional biological neuron model which produces action potentials and neuronal fluctuations; in dynamical system point of view, oscillatory behaviors related to its ionic channels in giant barnacle muscle fibers [9] [10]. Kathleen Morris and Harold Lecar proposed this simple two-dimensional model to generate spikes in 1981 [7] [11] [12]. The Morris-Lecar model describes the electrical activities of neurons with a system of two nonlinear ordinary differential equations and includes different channels. This model reduces the four-dimensional Hodgkin-Huxley model to be a two-dimensional system of ordinary differential equations while keeping the major neuronal properties of generating action potential but through simpler mathematical model [11] [12] [13]. The general Morris-Lecar model includes three channels: a potassium channel, a leak and a calcium channel and has the following form [7] ( where V demonstrates membrane potential, and n the activation variable of the persistent K + current, so it is a two-dimensional vector ( )

Morris-Lecar Model with the Impact of Cl − Fluctuations on Firing Patterns
In this study, we will discover the dynamics of Morris-Lecar model with the impact of Cl − channel that has the following form:   Figure 1). In this paper, we focus on this novel Morris-Lecar model for a certain range of parameters value which demonstrates different types of local bifurcations such as Hopf bifurcation and homoclinic bifurcation. We use rigorous analytic dynamical systems tools to study the complicated neuronal behaviors of this model such as normal form theory. We conduct bifurcation analysis with respect to different parameters to display the effect of various biological parameters on spiking dynamics of the system. Moreover, we are interested to discover the co-dimension two bifurcations such as Bautin or generalized Hopf and Bogdanov-Takens and cusp bifurcations and we present the normal form of these bifurcations.
its stability and becomes an unstable focus as we change the parameter for bifurcation, and the stable focus which is an attractor becomes a limit cycle that this limit cycle in phase space is a closed curve [1] [14] [15]. Using the parameters values corresponding to the Hopf case in Table 1 with 0.25 φ = , we could get to the first Hopf bifurcation but unlike [12], this bifurcation happens for larger amount of app I .
We have demonstrated the time series for Hopf case in Figure 2 and Figure 3. As we can see the neuron for firing a spike needs a larger amount of app I than the model original Morris-Lecar Model in Figure 2 of [12].
And for other values of app I near the Hopf bifurcation, we can see the changing in the behaviors of solutions. The neurons that show the subcritical Hopf bifurcation are resonators that exhibit bistability and the neurons that show supercritical Hopf bifurcations are resonators that exhibit monostability [1].
Continuing further, we compare the trajectories of two types of Hopf bifurcation. For model (1), we have both types of Hopf bifurcation. For supercritical Hope bifurcation that occurs for 208 app I ≈ a small limit cycle is born at the bifurcation point. This limit cycle will grow as we increase further the injected current (the amplitude of limit cycle increases with further increase of parameter). In Figure 4 and Figure 5, we can easily distinguish the trajectories of these two types of Hopf bifurcation, before and after the bifurcation points.     In a supercritical Hopf bifurcation, we have an equilibrium point that one limit cycle bifurcates from it. It means, for the parameter values before the parameter bifurcation, the limit cycle is degenerate and its amplitude is zero, and the amplitude extends as the parameter enhances further. But, in subcritical Hopf bifurcation, we have an unstable limit cycle around the equilibrium point, and a stable limit cycle surrounds it. This unstable limit cycle is as separatrix cycle  that separates two stable states and we could see the bistability that is not available in supercritical Hopf bifurcation and as a result only mono stability occurs in this case. The unstable limit cycle in subcritical Hopf bifurcation, degenerates to the equilibrium point, and the equilibrium points becomes unstable during this bifurcation. This type of bifurcation creates catastrophic changes in the behavior as we have suddenly changed the behavior from stable focus to oscillations with large amplitude The bifurcation diagram for model (1) has been displayed in Figure 6 using the parameters values from Table 1 (1) with the parameters similar to the Hopf case in Table 1.
Thus, there should exist an unstable limit cycle, bifurcating from the equilibrium and it indicates a subcritical Hopf Bifurcation for the first point and for the second point we have supercritical Hopf bifurcation a very small limit cycle is bifurcated at the bifurcation point. This limit cycle will grow with the further increasing of the injected current.
For the first point, the eigenvalues are

g inh as Bifurcation Parameter
In this section, we consider inh g and app I as free parameters. If we continue to increase app I , the period of limit cycle becomes 2-times and then 2 2 and … Here we say that the system undergoes the period-doubling bifurcation. As we do continuation of Hopf point, the trajectory repeats itself while we passed through PD point and we can see the trajectory is going around the loop twice before coming back to PD point. For period-doubling bifurcation a new limit cycle creates from the primary limit cycle, and the period of this new limit cycle is two times the previous limit cycle (see Figure 7 and Figure 8).  If we increase app I further, the limit cycle that is born from supercritical Hopf bifurcation undergoes another period-doubling bifurcation. The Period-doubling can continue to infinity. But, in this procedure, the values of app I at which this bifurcation occurs become closer than before.
For localization of Branch Points of cycles, a special algorithm which is of Newton type Continuation Method has been used to show these singularities. Also, the equilibrium curve can have two generic codimension one bifurcation limit point and Hopf point and plus continuation of equilibrium curve gives us Branch Point that in Figure 9, Figure 10, it has been shown by BPC.    The complete bifurcation diagram with respect to inh g has been demonstrated in Figure 11.
In Figure 11, codimension two bifurcations Bogdanov-Takens and Cusp bifurcations have been demonstrated. Bogdanov-Takens bifurcation is a codimension two bifurcation that has the following normal form x a x γ = +  . Figure 11. Complete bifurcation continuation.  So that we can easily write the normal form for this equilibrium point and we can find the bifurcation set and saddle-node curves.

g exc as Bifurcation Parameter
In this section, we consider exc g and app I as free parameters. At first, we assume exc g as bifurcation parameter (see Figure 12). In Figure 12     Also, at 0.071067341 inh g = and considering exc g as bifurcation, parameter with changing the period we have the result in Figure 13 and Figure 14.

The Homoclinic Case
In this section, we study the homoclinic bifurcation using the parameter values from Table 1.
When a saddle homoclinic bifurcation occurs a limit cycle disappears or born. A saddle point has two stable and unstable submanifolds. When we change the parameters for bifurcation, these two submanifolds collide and create one homoclinic orbit as we see clearly in Figure 15. After and before the moment of bifurcation we do not have this homoclinic orbit.

T. Azizi, B. Alali
In neuroscience point of view, before the beginning of bifurcation, all trajectories go back to the resting state after a counter-clockwise excursion (firing spike) and we have limit cycle attractors corresponding to periodic spiking. At the moment of bifurcation, we have periodic spiking activity with the period that goes to infinity (the period of limit cycle goes to infinity). After the bifurcation the state of the system again goes back to the rest and we do not have tonic spiking.
There is another bifurcation similar to saddle Homoclinic bifurcation and in both of them we have a Homoclinic orbit that begins and finishes in the same equilibrium point. But saddle-Homoclinic is of codimension 1 and saddle-node Homoclinic is of codimension two. For saddle-node Homoclinic we must have the condition for occurrence of saddle-node bifurcation 1 0 λ = and also the saddle-node quantity 1 2 0 λ λ + < it means that 2 0 λ < . As a result, in saddle-node Homoclinic, we observe the appearance or disappearance of a stable limit cycle.

Discussion and Conclusion
In this study we discovered the biological and spiking behaviors of Morris-Lecal neuron model including a new more channel, chloride channel, and we com- In supercritical Hopf bifurcation, when we changed the control parameters (bifurcation parameters), we could see a stable equilibrium point lost its stability and a stable limit cycle bifurcated from it. In this case, the state of the system stayed near the equilibrium and we could see only small amplitude oscillations around the equilibrium point. When we changed the parameter in the opposite direction the limit cycle degenerated to a point and we reached the previous equilibrium point. We noticed that as it is expected this neuron that has undertaken supercritical Hopf bifurcations did not show an immediate spike and it remained quiescent (it displayed small amplitude sustained oscillations).
In other cases, for subcritical Hopf bifurcation, when we changed the bifurcation parameter, we noticed the immediate spike or a large amplitude jump in time series of system equations. In this case, when we changed the parameter in opposite direction, the equilibrium became stable again but the state of the system went back to its similar saddle-node bifurcation or saddle-Homoclinic bifurcations.
As we have seen for model (1), there are different types of codimension-two bifurcations which occurred when we had two bifurcation parameters. One of them was the Bogdanov-Takens bifurcation which is a codimension-two bifurcation of an equilibrium point in a two-dimensional system for which the equilibrium point has two zero eigenvalues. Near the bifurcation parameter, we have two equilibrium points, a saddle point and others that are not a saddle point which collides and disappears via a saddle-node bifurcation. For the equilibrium that was not saddle, we have seen Hopf bifurcation that produced a limit cycle. This single limit cycle degenerated to a Homoclinic orbit to the saddle and disappeared via a saddle homoclinic bifurcation. Bogdanov-Takens bifurcation can occur when an equilibrium undergoes Hopf bifurcation and saddle-node bifurcation simultaneously and also it occurs when we have two-dimensional systems. Therefore, for this model we showed that we have a transition between integrator and resonator as the neuron that undergoes saddle-node bifurcation is an integrator (without damped subthreshold oscillations) and the neuron that undergoes Hopf bifurcation is a resonator (with damped subthreshold oscillations).
Continuation of limit cycles from the Hopf point, gave us the limit point (Fold/Saddle Node) bifurcation. Limit point cycles are a set of limit cycles that bifurcate from the Hopf point, and moreover, limit point cycle (LPC) is a fold bifurcation for which two limit cycles with different periods are near the LPC point. Likewise, as we saw, the continuation of cycles from a Hopf point also gave period doubling (flip) bifurcations. Period doubling bifurcation is a bifurcation for which a new limit cycle creates from the limit cycle that exists and the period of the new limit cycle is twice the first limit cycle. Moreover, we have shown changing the type of Hopf bifurcation from supercritical to subcritical which means that the type of stability changes from monostability in supercritical to bistability in subcritical Hopf bifurcation.