- Poster presentation
- Open Access
Bifurcation structure of phase locked modes in Type I excitable cells based on phase and spike time resetting curves
BMC Neuroscience volume 13, Article number: P78 (2012)
Some neurons are endogenous bursters, i.e., they intrinsically fire with a steady period (T i ) in the absence of external stimuli. In neural networks, the neurons constantly receive stimuli, e.g., presynaptic inputs that perturb their intrinsic activity. The transient changes in the firing frequency of a neuron due to a presynaptic input received at a stimulus time (t s ), or phase φ = t s /T i , is tabulated by the phase resetting curve (PRC). The first order PRC is extracted in open loop, i.e., from isolated neurons, by tabulating the relative change in the firing period F(φ) = ∆T 1 /T i = (T 1 -T i )/T i due to a stimulus received during the current cycle at phase φ (see Fig. 1A). Equivalently, we use the spike time resetting curve (STRC) that tabulates the recovery time t r of a neuron that received an input during the current cycle at a stimulus time t s (see Fig. 1A). The assumption is that the open loop PRC/STRC can predict the phase locked modes of a network provided, among other requirements, that the open loop stimulus closely resembles the input the neuron receives in a fully connected network. For example, the PRC/STRC method can predict the 1:1 phase locked modes of a two-neuron network (see Fig. 1B). Assuming that the effect of the presynaptic input dies out by the time the neuron receives another input, we can derive (Figure 1B) the recursions that determine the phase locked mode:
which can be translated into phase recursions by using the definitions φ = ts/Ti and T1(φ) = Ti(1+F(φ)). Very close to the bifurcation point of a Type I excitable cell, the infinitesimal PRC is F(φ) = cTi(1-cos(φ)), where c is a model-dependent constant [1, 2]. We explored the bifurcation structure of phase locked modes of two coupled Type I intrinsic bursters (see a particular example of bifurcation diagram in Figure 1C). We also derived stability criteria for predicted phase locked modes. The advantage of the proposed approach is that it reveals the route to phase locked mode transition. The challenge of the method is representing efficiently a high dimensional solution in two dimensions. For this purpose, we used the method of multidimensional stacks.
Brown E, Moehlis J, Holmes P: On the Phase Reduction and Response Dynamics of Neural Oscillator Populations,. Neural Comp. 2004, 16: 673-715. 10.1162/089976604322860668.
Ermentrout G: Type I membranes, phase resetting curves, and synchrony. Neural Comp. 1996, 8: 979-1001. 10.1162/neco.19188.8.131.529.
Research supported by NSF IOS – 1054914 grant to S.A.O.