Mixed mode synchronization and network bursting of neurons with post-inhibitory rebound
|
|
- Mark Wilkins
- 5 years ago
- Views:
Transcription
1 Recent experimental and theoretical studies have disclosed a distinct role of CPGs in generation of adaptive and coordinated motor activity of animals An important feature of CPGs is their ability to produce various types of rhythmic bursting activity, what causes flexible and adaptive locomotion of an organism. To roarxiv: v1 [nlin.cd] 5 Dec 14 Mixed mode synchronization and network bursting of neurons with post-inhibitory rebound Roman Nagornov, 1, a) Grigory Osipov, 1 Maxim Komarov, 1, 2 Arkady Pikovsky, 1, 2 and Andrey Shilnikov 1, 3 1) Department of Control Theory, Nizhni Novgorod State University, Gagarin Av. 23, 695, Nizhni Novgorod, Russia 2) Department of Physics and Astronomy, University of Potsdam, Karl-Liebknecht-Str 24/25, Potsdam, Germany 3) Neuroscience Institute and Department of Mathematics and Statistics, Georgia State University, Piedmont Str., Atlanta, GA 333, USA (Dated: 22 September 18) We study activity rhythms produced by a half-center oscillator a pair of reciprocally coupled neurons that are capable of producing post-inhibitory rebounds. Network (coupling-induced) burstings possessing two time scales, one of fast spiking and another of slow bursting and quiescent periods, are shown to exhibit nontrivial synchronization properties. We consider several network configurations composed of both endogenously bursting, tonic spiking and quiescent neurons, as well as various mixed combinations. We show that synchronization at low frequency of bursting can be accompanied by complex, sometimes chaotic, patters of fast spikes. This study is focused on the mechanisms of rhythmogenesis and robustness of anti-phase bursting in half-center-oscillators (HCOs) consisting of two reciprocally inhibitory coupled neurons. There is a growing body of experimental evidence that a HCO is a universal building block for larger neural networks, including central pattern generators (CPGs) controlling a variety of locomotion behaviors in spineless animals and mammals. It remains unclear how CPGs achieve the level of robustness and stability observed in nature. There has been a vastly growing consensus in the community of neurophysiologists and computational researchers that some basic structural and functional elements are likely shared by CPGs of both invertebrate and vertebrate animals. In this study we consider several configurations of HCOs including coupled endogenous bursters, tonic spiking, and quiescent neurons, that become network bursters only when coupled by fast inhibitory synapses through the mechanism of post-inhibitory-rebound (PIR). The goal is better understanding the PIR mechanism as a key component for robust anti-phase bursting in generic HCOs. I. INTRODUCTION Synchronization of coupled oscillators is a fundamental phenomenon in nonlinear systems that has been observed in a wide range of diverse applications 1. The mathematical concept of synchronization, first introduced and developed for periodic oscillators 2, has further been generalized for other aperiodic systems, including ones with a) Electronic mail: nagornov.r@gmail.com chaotic dynamics. In life sciences, of a keen interest is synchronization or phase locking among oscillators with multiple time scales. They may include mixed-mode and slow-fast relaxation-type oscillators 3, whose synergetic interaction can give rise to the onset of a variety of synchronization patterns 4 7. In neuroscience, a plethora of rhythmic motor behaviors with diverse time scales, such as heartbeat, respiration, chewing, and locomotion on land and in water are produced and governed by neural networks called Central Pattern Generators (CPGs) 8,9. A CPG is a microcircuit of neurons that is able to autonomously generate an array of polyrhythmic bursting patterns, underlying various motor behaviors. Endogenous (self-sustained) bursting and network (coupling-induced) bursting are composite oscillatory behaviors, featuring active phases during which a neuron or a group of neurons generates trains of fast action potentials, which are alternated with long interburst intervals during which it remains inactive or quiescent, until a new cycle of bursting occurs. In this paper we examine synchronization of bursting patterns emerging through interactions of two interneurons coupled reciprocally by fast inhibitory synapses. This study has been driven by two major motivations: first, a general one concerning questions on synchronization of mixed-mode oscillators. The second is a neuroscience related one, aimed at understanding of intrinsic mechanisms of rhythmogenesis in CPGs composed, often symmetrically, of such small networks of interneurons, as outlined below. It is still unclear how CPGs can achieve the level of synergy and robustness to produce a plethora of rhythmic patterns observed in nature.
2 2 bustly govern motor patterns, CPGs are able to flexibly adjust their oscillatory properties (such as bursts duration, frequency of spiking, phase relations of bursts) due to feedback from sensory inputs, for example, in response to changes of environment 9,1. Up to a certain extent, the flexibility of CPG behaviors may be attributed to its multistability (of several coexistent attractors representing different bursting rhythms in a phase space of the dynamical system) allowing for fast switching between operating modes 7,13. From the theoretical point of view, a CPG is modeled as a small network of coupled oscillatory, or quiescent, interneurons, each described by a system of nonlinear ordinary differential or difference equations (dynamical system) 11,14. The study of CPGs allows one to develop a general understanding of synchronization patterns in mixed-mode oscillators, applicable to systems of different physical and biological origin. There is a growing body of experimental evidence that a universal building block of most identified CPGs is a half-center oscillator (HCO) 15. A HCO is a pair of two reciprocally inhibitory interneurons bursting in alternation. Such a pair can be composed of endogenously bursting interneurons, as well as of endogenously tonic spiking or quiescent ones that start bursting only when they are coupled in a network. Theoretical studies 16,17 have indicated that formation of an anti-phase bursting rhythm is always based on some slow subsystem dynamics. The slow subsystem is strongly associated with the slow membrane current, such as persistent sodium current or slow calcium-dependent current; following 18 we term currents associated with slow-varying concentrations and gating variables as slow ones. There are three basic mechanisms to generate alternating bursting in the HCO: release, escape, and post-inhibitory rebound (PIR). The first mechanism is typical for endogenously bursting neurons 4,6. The other two mechanisms underlie network bursting in HCOs comprised of neurons which are depolarized or hyperpolarized quiescent interneurons, respectively, in isolation 9,19 22 The PIR mechanism uses reciprocal inhibition to keep network bursting. It triggers an onset of a single or a series of action potentials in the post-synaptic neuron after it has been prolongedly hyperpolarized and abruptly released from inhibition generated by the pre-synaptic neuron during an active, tonic spiking phase of bursting. After that, the neurons of the HCO swap their opposite roles to repeat the PIR mechanism. The PIR promotes the action potential generation after a period of sufficiently strong hyper-polarizing (inhibiting) input, as illustrated in Fig. 1. PIR is often caused by a lowthreshold activated calcium current in neurons and their biophysically plausible models. This paper examine the PIR contribution to formation, synchronization and robustness of bursting rhythms in inhibitory coupled neurons. Here, we employ a modification of the Hodgkin-Huxley-type model introduced in 18, to plausibly describe the PIR mechanism. Depending on its parameters, the model is known to produce an array of typical neuronal activities such as excitable dynamics, periodic spiking and bursting. Below we first examine the conditions that confirm and reproduce the PIR mechanism in the neurons. We will argue that the PIR is a pivotal component that promotes an alternating bursting rhythm in the HCO made of intrinsically spiking and excitable neurons. Finally, we will show that the PIR mechanism makes anti-phase network bursting even more robust and easier to initiate in the case of endogenously bursting neurons. The paper is organized as follows: in Section II we introduce the neuronal model, and discuss its dynamical properties in Section III. Section IV is focused on synchronization properties of a pair of endogenous bursters, while Section IV examines anti-phase bursting onset in coupled tonic-spiking neurons. Sections VI and VII discusses the PIR-mechanisms of the HCO, and synergetic dynamics of mixed neurons: individually tonic spiking and hyperpolarized quiescent. I ext (µa/cm 2 ) 2. g Ca =1.75 ms/cm 2 g Ca =.1 ms/cm FIG. 1. Abrupt release of a hyperpolarized pulse of the external current I ext (bottom panel) triggering a post-inhibitory rebound of bursting in the quiescent neuron (top panel). The parameter g Ca stated in the legend denotes maximal conductance of the current I T. When the parameter g Ca is relatively small, the effect of the I T current is negligible what causes absence of PIR in the dynamics of the neuron (dashed green curve in the top panel). II. BASIC MODEL In this study we use a neuronal model of the Hodgkin- Huxley type that has been introduced in 18. The basic equations describing the dynamics of the membrane potentials of the neurons read as follows (prime denotes the
3 3 time derivative, i = 1, 2): C m V i = I (i) ext f i (V i ) N j=1,j i I syn (V i, V j ), f i (V i ) = I (i) leak + I(i) Na + I(i) K + I(i) K[Ca] + I(i) T. (1) The variable V i (t) describes evolution of the membrane potential of the i-th neuron. The first two terms on the right-hand side of Eq. (1) govern intrinsic dynamics of the neuron: I (i) ext stands for a constant external current applied to the neuron, while term f i (V i ) represents the sum of intrinsic ionic currents. The sum of synaptic currents I syn (V i, V j ) describes the coupling N interac- j=1,j i tions between the neurons. The full details of current representation are given in the Appendix below. There are two principal control parameters in this system: the external current, I (i) ext, governs the activity type in an isolated, uncoupled neuron. The second parameter, the maximal conductance G of the synaptic current given by Eq. (13), controls the coupling strength, and thus the network dynamics of the coupled neurons. The low-threshold activated calcium current I T, modeled by Eq. (7) and regulated by the maximal conductance g Ca, causes the PIR mechanism in the neurons (1), this effect is illustrated in Fig. 1. Rebound bursts can be triggered by injection of a hyper-polarizing pulse of external current I (i) ext (bottom panel in Fig. 1). If the impact of I T current is strong enough, the neuron generates bursts of spikes, after recovering from the hyperpolarized quiescent state (top panel in Fig. 1). III. DYNAMICS OF AN ISOLATED NEURON Let us first consider the dynamics of an isolated neuron at I syn = in Eq. (1). We treat the external current I ext as the primary bifurcation parameter in both computational and biological experiments, which allows one to examine and plausibly model various neuron activities and transitions between them. For the same purpose, we employ interspike intervals (ISIs), which are distances between adjacent spikes generated by the neuron, as measurements characterizing the neuron s dynamical regimes. Figure 2 serves as an illustration of how the ISIs changes with variations of I ext. First we describe the neuron dynamics at large I ext values. For I ext > 4 µa/cm 2 the neuron produces tonic spiking activity (corresponding to a stable periodic orbit in the phase space of the model (1)). Here, the value of ISI (which is the period of the stable orbit) is inversely proportional to the I ext -value. The period decreases with increase of I ext, and vice-versa. As one can see from Fig. 2, near a critical value I ext 4 µa/cm 2 the ISI shows an unbounded growth, which is an indication of a bifurcation of the stable periodic orbit with an I ext (µa/cm 2 ) FIG. 2. Bifurcation diagram for an isolated neuron: plotting inter-spike intervals (ISI) against the external current I ext (this horizontal axis is broken in the middle to highlight the ranges of nontrivial behaviors) reveals windows of bursting (.14 µa/cm 2 < I ext <.5 µa/cm 2 ), quiescence (.5 µa/cm 2 < I ext < 4 µa/cm 2 ), and tonic spiking (I ext > 4 µa/cm 2 ) activity in the neuron. The low ISI branch corresponds to short time intervals between fast spikes, and the top branch represent long interburst intervals between consecutive spike trains. arbitrarily long period. Detailed examinations indicate the occurrence of a homoclinic saddle-node bifurcation 23 underlying the transition from tonic spiking activity to hyperpolarized quiescence, represented by a stable equilibrium state at low values of the membrane potential V. This stable equilibrium state (a node with real and negative characteristic exponents) persists within the parameter window.8 µa/cm 2 < I ext < 4 µa/cm 2, where the neuron remains ready for PIRs. Small perturbations of the quiescent state of the neuron have no pronounced effects. Relatively strong perturbations can trigger a spike, after which the neuron comes back to the over-damped quiescence state. With decreases of I ext from the threshold I ext 4 µa/cm 2, the steady state becomes a focus. The stable focus loses the stability through a supercritical Andronov-Hopf bifurcation at I ext.8 µa/cm 2. These oscillations are not seen in Fig. 2 because they are under the threshold of spike generation. The amplitude of the oscillations increases rapidly as I ext is further decreased. At I ext.5 µa/cm 2 the shape of the periodic orbit in the phase space changes via the addition of turns, the so-called spike-adding mechanism 24,25. Such a periodic orbit is associated with bursting activity, with turns corresponding to spikes within a burst. Bursting that includes two time scales can be recognized in Fig. 2 through two characteristic branches: the bottom branch corresponds to small ISI values due to fast spiking, while the top branch is due to long interburst intervals. More details on geometry of bursting can be found in 14,26. Decreasing the external current I ext makes 7
4 4 the bursting neuron more depolarized and increases the number of spikes per burst. Having reached some maximal value of spikes per burst at I ext.136 µa/cm 2, decreasing I ext gives rise to bursts with fewer spikes. The value I ext.14 µa/cm 2 corresponds to the occurrence of another saddle-node bifurcation of equilibria since the time interval between two consecutive bursts becomes arbitrarily large. For I ext <.14 µa/cm 2 the neuron remains quiescent and excitable. The examination of bifurcations have been repeated for several values of the maximal conductance, g Ca, regulating the low-threshold Ca 2+ current responsible for PIR in the neuron. We have found that this current bears an insignificant affect on the intrinsic dynamics of the individual neurons, mainly because the membrane potential does not decrease below the threshold value to activate the current. As we will see below, this is not the case for coupled neurons, where the parameter g Ca has a pronounced effect on the collective dynamics. IV. SYNCHRONIZATION OF TWO BURSTING NEURONS In this section and sections V-VII below, we explore a repertoire of rhythmic bursting outcomes generated by a HCO of two neurons, coupled reciprocally by fast, nondelayed inhibitory synapses. It is known that such coupled, endogenous bursters can produce a range of synchronous rhythmic outcomes with various fixed phaselags due to overlapping spike interactions, see 6,27 and references therein. We will first examine the cooperative dynamics in the network of two coupled neurons (N = 2 in Eqs. (1)) with different I (1,2) ext within the interval [.14 µa/cm 2,.8 µa/cm 2 ]. This range of I ext corresponds to endogenous bursting in both neurons. Due to the difference, I ext = I (2) ext I ext, (1) the temporal and quantitative characteristics of endogenous bursters such as their period, duty cycles, the spike numbers per burst, are different. The strength of coupling is quantified by the maximal conductance, G, of the inhibitory synaptic current I syn. In what follows, we show that while I ext remains relatively small, increasing the coupling strength shall give rise to an onset of synchrony between the neurons with the same bursting period. However, this will no longer be true for larger I ext values. A. Nearly identical bursters While I ext <.5 µa/cm2, synchronous bursting can already occur at a relatively weak coupling. Here, the burst periods of both neurons become equal, however the spike number per burst in the neurons may not be the same. This is shown in Fig. 3, which depicts the dependence of the mean phase-lag between oscillations of the low-threshold Ca 2+ -currents in the neurons, on ϕ g Ca =1.3 ms/cm 2 g Ca =1.75 ms/cm 2 g Ca =2.3 ms/cm FIG. 3. Mean phase lag, ϕ plotted against the coupling strength G for different values of the maximal conductance g Ca of the low-threshold Ca 2+ -current (larger g Ca values promote stronger PIRs in the neurons); here I (1) ext =.2 µa/cm 2 and I (2) ext =.15 µa/cm 2. the coupling strength G for several values of the conductance g Ca. Note that the post-inhibitory rebounds in the neurons under consideration are due to the slowest, low-threshold Ca 2+ -current I T, whose magnitude is controlled by the maximal conductance g Ca. It is worth noticing that PIRs occur more reliably with increasing g Ca. Frequency(kHz) I =.2 µa/cm 2 I =.5 µa/cm 2 I =.1 µa/cm 2 I =.15 µa/cm 2 I =.2 µa/cm FIG. 4. Bifurcation diagram for antiphase synchronization regimes: averaged frequency of bursting is plotted against the coupling strength for several increasing I ext values. Double overlapping branches are the indication of bursting dichotomy with two slow different frequencies in the two coupled neurons. Note a large plateau of the pronounced 4:3 frequency locking at I ext =.2 µa/cm 2, collapsing into anti-phase synchrony locked at an 1:1 ratio at higher G values. Nevertheless, for a broad range of initial condition, the
5 5 neurons will predominantly burst in anti- phase rather than in-phase. This implies that the HCO can be bistable. This observation may shed light on properties of multifunctional CPGs composed of coupled neurons including coupled HCOs, which can switch between multiple functional states associated with different locomotion types. While in-phase bursting appears to be atypical for neurons with fast inhibitory coupling like in our case, the coexistence of anti-phase and in-phase regimes have been reported elsewhere 4,6. However, in-phase synchrony can prevail whenever both HCO neurons are driven externally by another inhibiting burster or HCO 28,29. Next, we discuss the effects that the PIR activity of the neurons may have on antiphase bursting. Figure 3 discloses how the coupling threshold for the onset of synchronous bursting depends on mismatches between the neurons of the HCO, as well as on the phase-locked frequency of slow I T -oscillations. The increasing phases, defined on modulo 2π, of the neurons are reset at the very beginning of each burst. Figure 3 shows the evolution of the phase lag, ϕ, between the oscillations by the slow gating variables, h (1,2) T, for inactivation of the lowthreshold Ca 2+ -current, which is is plotted against the coupling constant G for various g Ca -values. We found in most cases that the phase lag tends eventually to π, i.e. the anti-phase bursting occurs in the HCO, with the coupling strength above a critical value. Depending on I ext, the actual ϕ-value may deviate from π. One can see from Fig. 3 that increasing g Ca (promoting stronger post-inhibitory rebound activity) leads to quite symmetric antiphase bursting even at small values of the coupling strength. This is an explicit manifestation of the contribution of the slow low-threshold Ca 2+ -current in fostering the antiphase synchronization between the bursting neurons in the HCO. B. Bursters far from identical Increasing I ext between coupled oscillators leads to an array of pronounced synchronization effects. Figure 4 represents the slow (bursting) frequencies plotted against coupling strength, G, at several I ext values. One observes the presence of overlapping branches at low G- values. This indicates that the coupled neurons generate bursting activities at different frequencies, until the coupling is increased over a threshold value at which both branches merge. This threshold value becomes higher with increases of I ext, as the individual neurons become more and more different. We note from this figure that for the largest case I ext =.2 µa/cm 2, the HCO neurons become locked at the 4:3 frequency ratio, prior to the occurrence of the ultimate 1:1 frequency locking at higher coupling values. Figure 5 illustrates transitions to synchrony for the largest case I ext =.2 µa/cm 2. The left panels demonstrate the established antiphase bursting in voltage traces produced by the HCO at various values of the coupling h T 2 h T 2 h T 2 h T 2.25 spikes spikes...12 h T 1.25 spikes spikes...14 h T 1.2 spikes spikes...12 h T 1.2 spikes spikes...12 h T 1 FIG. 5. Left column shows progressions of phase locked voltage traces of two coupled neurons at I ext =.2 µa/cm 2, I (1) ext =.2 µa/cm 2 for G =.2 ms/cm 2 ;.445 ms/cm 2 ;.464 ms/cm 2 ;.48 ms/cm 2. Right panels demonstrate the Lissajous curves drawn by slow variables h (1) T and h (2) T, which correspond to (top) quasi-periodic dynamics at G =.2; a 4:3-frequency locking regime at G =.445 ms/cm 2 ; an 1:1 chaotic locking at G =.464 ms/cm 2 ; (bottom) an 1:1 periodic locking at G =.48 ms/cm 2. Circles and squares mark spike events in bursting. strength, G, below and above the threshold. The right panels represent the so-called Lissajous curves, which are parametrically traced down by the slow variables h (1) T and h (2) T of both neurons. These curves help to interpret the corresponding types of frequency locking including quasiperiodic dynamics. On these curves we mark timing of individual spikes. It is worth noticing that such a large I ext makes the phases of slow oscillations shift closer to π/2, rather than to π. With a π/2 phase shift, active phases of bursting in the neurons partially overlap, showing the emergence of complex dynamics caused by spike interactions. To understand these dynamics, we evaluated the distribution of numbers pf spikes per burst for several coupling strengths. Our findings are presented in
6 6 Fig. 6. It is clear from this figure that the transition to the 1:1 locking occurs at G.45mS/cm 2. Moreover, our simulation indicated that there is a small window, G ms/cm 2, of hysteresis occurring at the transition. Figure 6 reveals a pecular feature of the synchronous state, existing at strong coupling G [.45 ms/cm 2,.485 ms/cm 2 ]: while the spike number in bursts generated by neuron 2 remains nearly constant, the number of spikes in bursts by neuron 2 deviates quite strongly before it becomes fixed with stronger coupling. From the two bottom insets in Fig. 5 one can observe that depending on coupling strengths the 1:1 locking features chaotic and periodic modulations via the slow gating variables of the neurons. In the chaotic case, the number of spikes per burst in either neuron or both neurons will vary while the alternating bursting in the HCO persist locked at an 1:1-ratio within a windows [.46 ms/cm 2,.47 ms/cm 2 ]. This is an explicit manifestation of the phenomena of chaotic phase synchronization 5. spikes/burst FIG. 6. Bifurcation diagram showing spike number per burst generated by the neurons plotted against the coupling constant G; other parameters are same as in Fig. 5. Value G.45 ms/cm 2 is a threshold towards the 1:1-frequency locking; G.464 ms/cm 2 corresponds to phase slipping in the chaotic 1:1 locking state of alternating bursting. V. ANTIPHASE BURSTING OUT OF SPIKING NEURONS In this section we study emergence of alternating bursting in a HCO made of two coupled neurons that are tonically spiking in isolation at I (1) ext = 5 µa/cm 2 with relatively small I ext =.2 µa/cm 2. Our goal is to reveal how increasing the inhibitory coupling strength transforms such tonic spikers into network bursters, as illustrated in Fig. 7. Weak coupling gives rise to periodic antiphase tonic spiking in the network at G =.2 ms/cm 2 (Fig. 7(a)). With further increases of the coupling strength, the HCO (a) (c) (b) FIG. 7. Voltage oscillatory activity generated by the HCO at g Ca = 1 ms/cm 2 and I ext 1 = 5, and I ext =.2 µa/cm 2 : (a) antiphase spiking at G =.2 ms/cm 2 ; (b) chaotic spiking activity at G =.3 ms/cm 2 ; (c) forced sub-threshold oscillations in neuron 2 due to fast spiking in neuron 1 at G = 2 ms/cm 2, (d) network bursting at G = 4 ms/cm 2. neurons evolve through a window of irregular, unreliable spiking at G =.3 ms/cm 2, towards an asymmetric network dynamics. In this robust regime, either neuron generates tonic spiking activity, while the other one is forced to produce sub-threshold oscillations of quite a large amplitude, such as ones shown in Figs. 7(b,c) at at G = 2 ms/cm 2. A dramatic increase in the inhibitory coupling strength beyond G 3 ms/cm 2 finally forces the neurons to begin bursting in alternation. Such bursting activity is often referred to as a network bursting, which is the result of strongly reciprocal interactions of two individually tonic spiking neurons. Figure 2 provides an explanation why inhibitory coupling must be so strong to induce network bursting in the given model. There is a wide gap between the I ext parameter values corresponding to bursting (at the low end) and to tonic spiking (on the high end). Within this gap, the neurons remain hyperpolarized quiescent, and hence the network bursting can only occur through PIR mechanisms. A strong flux of inhibitory current is required to originate from the presynaptic, tonic spiking neuron in order to make the postsynaptic neuron drive over the gap to become a network burster. Then, the neurons swap the roles of the driving and driven neurons and so forth. This assertion is further supported by Fig. 8 plotting inter-spike intervals (ISIs) in the driven postsynaptic neurons (neurons 1 and 2 in Insets (a) and (b), resp.) against the coupling strength for three values of g Ca. We can see the corresponding gap in the right panel, within which neuron 2, through strongly inhibitory force, generates large sub-threshold oscillations, partially because of the network asymmetry due to I ext =.2 µa/cm 2. (d)
7 7 Figure 8 also discloses the qualitative role of the slow low-threshold Ca 2+ -current I T ; namely, increasing the corresponding rebound parameter g Ca lowers the threshold of coupling-induced bursts, and shrinks the parameter interval of hyperpolarized quiescence in the neurons g Ca =1 ms/cm 2 g Ca =1.3 ms/cm 2 g Ca =2.3 ms/cm (a) g Ca =1 ms/cm 2 g Ca =1.3 ms/cm 2 g Ca =2.3 ms/cm FIG. 8. Bifurcation diagram of inter-spike intervals for coupled neuron 1 (a) and neuron 2 (b) plotted against the coupling constant G for three increasing values of g Ca; here I ext 1 = 5 µa/cm 2 and I ext =.2 µa/cm 2. (b) which the HCO robustly produces anti-phase bursting. Figure 9 also demonstrates that the PIR mechanism of network bursting becomes more reliable with increasing g Ca, which lowers the threshold value of inhibitory coupling. Regularization of PIR network bursting is quite sensitive to variations of the coupling strength. By inspecting the top ISI-branch in Fig. 9 one can identify the stability windows of G-values corresponding to plateau-like ISI intervals alternating with the windows of fluctuations of ISI values. Fluctuations are not marginal about 1% off the mean value. The PIR network bursting is illustrated in Fig. 1(a). Fig. 1(b) reveals a quasi-periodic modulation of the slow gating variables; here the periodic Lissajous curve has a shape of figure-eight. h T 2.6 spikes spikes VI. PIR MECHANISM FOR ANTIPHASE BURSTING h T 1 (a) (b) g Ca =1.3 ms/cm 2 g Ca =1.75 ms/cm 2 g Ca =2.3 ms/cm (a) 18 1 g Ca =1.3 ms/cm 2 g Ca =1.75 ms/cm 2 14 g Ca =2.3 ms/cm (b) FIG. 1. (a) Antiphase network bursting through the PIR mechanism is filled out densely by alternating voltage traces without quiescent gaps unlike the case of coupled endogenous bursters through the release mechanism. (b) Lissajous curve of a figure-eight shape traced down by the slow gating variables of the models at G = 1.24 ms/cm 2 and g Ca = 1.75 ms/cm 2 ; other parameters the same as in Fig. 9). Circles and squares in panel (b) mark occurrences of spikes. FIG. 9. Bifurcation diagram of ISIs evolutions for coupled neuron 1 (a) and neuron 2 (b) plotted against the coupling constant G; here I (1) ext = 2 µa/cm 2 and I ext =.2 µa/cm 2 correspond to hyperpolarized quiescent neurons in isolation. In this section we examine the emergence of network antiphase bursting through the PIR mechanism. Here, it is imperative that both neurons remain hyperpolarized within a parameter window quiescence.5 µa/cm 2 < I ext < 4 µa/cm 2, as seen from Fig. 2. It is also necessary for PIR network bursting that the initial states of the neurons must be different: one tonic spiking and one hyperpolarized quiescent. An alternative is an application of negative pulse of current to trigger a PIR in a targeted neuron. In addition to the above constrains, the coupling strength must exceed a certain threshold, see Fig. 9. It is seen from this plot that relatively weak coupling cannot initiate PIR network bursting. Rather some subthreshold oscillations are generated in the post-synaptic neuron, which cease as soon as the pre-synaptic neuron ends its active spiking phase and becomes hyperpolarized quiescent as well. For each set of the parameters there is a threshold for the coupling strength beyond VII. COUPLING OF NEURONS IN DIFFERENT MODES In this section we consider the network dynamics of coupled neurons operating in different modes. Specifically, the neurons are set close to quiescence-spiking and quiescence-bursting transitions, respectively. As Fig. 2 suggests, we choose the following values of the external drive: I (1) ext = 4.1 µa/cm 2 and I (2) ext = 3.7 µa/cm 2. This means that in isolation neuron 1 produces tonic spiking activity, while neuron 2 is hyperpolarized and ready for the PIR mechanism. Our network simulations of two such neurons are summarized in Figs. 11 and 12. We can conclude that while coupling remains weak, neuron 1 tonically spikes (Fig. 12(a)). Above a threshold, G 1.75 ms/cm 2, neuron 2 begins spiking as well. Fig. 12(b)) depicts alternating spiking at G = 1.8 ms/cm 2 with the locking ratio 2:1, which means two spikes in neuron 1 are followed by a single spike in neuron 2. The network begins to burst in antiphase, first with irregularly varying spike numbers, as G exceeds 2 ms/cm 2, as illustrated in Fig. 12(c). Further increase in the coupling strength above G = 2.5 ms/cm 2, regu-
8 8 larizes network bursting that becomes stable and similar to the pattern shown in Fig FIG. 11. Bifurcation diagram of ISI plotted against the coupling parameter G indicates the threshold beyond which the network made of the tonically spiking neuron 1 and the hyperpolarized quiescent neuron 2 begins generating anti-phase bursting; g Ca = 1.75 ms/cm (a) (b) FIG. 12. Established voltages traces generated by the HCO neurons: (a) tonic spiking in neuron 1 and inhibition induced subthreshold oscillations in neuron 2 at G = 1 ms/cm 2 ; (b) periodic antiphase spiking with a 2:1 locked ratio at G = 1.8 ms/cm 2 ; (c) chaotic antiphase bursting near the threshold at G = 2 ms/cm 2. Figure 13 represents the bifurcation diagram for the network of an endogenous burster neuron 2 at I (2) ext =.4 µa/cm 2, coupled with neuron 1 hyperpolarized at I (1) ext = 1 µa/cm 2. Because the quiescent neuron 1 initially remains below the synaptic threshold, anti-phase bursting in the network may only start when the PIR mechanism is induced by the endogenous burster. This explains a qualitative resemblance of the diagram in Fig. 13 with that in Fig. 11 for the HCOP made of tonicspiking and quiescent neurons (c) VIII. CONCLUSIONS This study focused on the mechanisms of rhythmogenesis of anti-phase bursting in HCO networks consisting of two reciprocally inhibitory coupled neurons. Such HCOs are primary building blocks for larger neural networks, including CPGs controlling a plethora of locomotion behaviors in spineless animals and mammals. There is a growing consensus in the neuroscience community that CPGs share some universal principles of their functioning. In this study, for CPG neurons we used the biophysically plausible Hodgkin-Huxley-type model. Its main feature is post-inhibitory rebound dynamics: a quiescent neuron is able to produce a single or a series of spikes after it has been quickly released from inhibition by another pre-synaptic neuron, or by a hyperpolarized pulse of external current. The PIR mechanism allows a pair of such neurons to stably generate anti-phase bursting in a network initiated externally by inhibitory perturbation(s). We considered several configurations of HCOs including coupled endogenous bursters. We also discussed HCOs involving tonic spiking and quiescent neurons that become network bursters only when coupled by fast inhibitory synapses. In our examination of synchronization properties of bursting we found that, in all considered cases, the network can reliable achieve synchrony in antiphase bursting. We also described some partial configurations leading to incomplete synchronization where the neurons become partially synchronized, through slowvarying currents. Meanwhile, cross-correlations in their fast voltage dynamics are not always obvious. These fast voltage dynamics may give rise to the emergence of synergistically complex states, including chaos in the neural network. We found that while enhancing the PIR mechanism does not lead to drastic changes in the dynamics of the individual neurons, it causes significant modulation of network dynamics. Specifically, we found that the windows of anti-phase bursting rhythms can be extended in the parameter space of the network, when increasing the PIR mechanisms in individual neurons become more prominent. This suggests that the PIR is a key component for robust and stable anti-phase bursting in HCOs. In the future, we plan to examine specific networks consisting of several HCOs, which have been identified in swim CPGs of several sea mollusks ACKNOWLEDGMENTS FIG. 13. Bifurcation diagram of ISI plotted against the coupling parameter G indicates the threshold beyond which the network made of the bursting neuron 2 and the hyperpolarized quiescent neuron 1 begins generating anti-phase bursting; g Ca = 1.75 ms/cm 2. The research is supported by the grant (the agreement of August 27, 13 N 2.B between The Ministry of education and science of the Russian Federation and Lobachevsky State University of Nizhni Novgorod, Sections 2-4) and by the Russian Science Foundation (Project No , Sections 5-7). A.P. thanks
9 9 A. Politi for interesting discussions and acknowledges the G.Galilei Institute for Theoretical Physics (Italy) for the hospitality and the INFN for partial support during the completion of this work. M.K. thanks Alexander von Humboldt foundation for support. A.S. also acknowledges the support from NSF grants DMS-9591, RFFI , and RSF grant We thank A. Kelley for helpful suggestions. IX. APPENDIX: CONDUCTANCE BASED MODEL The model in this study is adopted from Ref. 18. The dynamics of the membrane potential, V is governed by the following equation: C m V = I ext I T I leak I Na I K I K[Ca] I syn ; (2) here, C m is the specific membrane capacity, I ext is the external current in na, I T is the slow low-threshold Ca 2+ - current, I leak is the leakage current, I Na is the Na + - current, I K is the K + -current, I K[Ca] is the slow Ca 2+ - activated K + -current which moderates bursting activity in the model, and I syn is the synaptic current from other neurons. Leak current I leak is given by I leak = g Na (V E leak ) (3) with E leak = 78 mv being reversal potential for leak current, and maximal conductance g L =.5 ms/cm 2. Dynamics of fast Na + -current I Na = g Na m 3 h(v E Na ) is described by the following equations: m =.32(13 V ).28(V 4) (1 m) e.25(13 V ) 1 e.2(v 4) 1 m h =.128 e 17 V 18 (1 h) 4 e.2(v 4) + 1 h, (4) where g Na = ms/cm 2 is the maximal conductance of Na + -current, E Na = 5 mv is the Nernst equation potential for Na + -current, m and h are the gating variables describing activation and inactivation of the current. Dynamics of fast K + -current I K is described by I K = g K n 4 (V E K ), n =.32(15 V ) e.2(15 V ) 1 (1 n).5e (1 V ) 4 n (5) with n being the gating activation variable; here E K = 95 mv and g K = 1 ms/cm 2. Dynamics of intraneuronular concentration of calcium ions [Ca] is described [Ca] = ki T 2F d K T [Ca], (6) [Ca] + K d where the first term is a inflow through thin membrane due low-threshold Ca 2+ -current, and the second is a contribution of Ca 2+ ion-pump, F = 96, 469 C/mol is the Faraday constant, d = 1 µm is the membrane thickness, k =.1 is the dimension constant for [I T ] = µa/cm 2 and [Ca], that has dimension of millimole, K T = 1 4 mm/ms 1 and K d = 1 4 mm are the speed constants from the Michaelis Menten approximation of ion-pump kinetics, they are taken the way that one can observe fast output of calcium from a neuron in millisecond time scale. Dynamics of slow low-threshold Ca 2+ -current I T is modeled as follows: I T = g Ca m 2 T h T i (V E Ca ), (7) with g Ca = 1.75mS/cm 2, and the gating activation m and inactivation h are given by with m = m m (V ) τ m (V ), h = h h (V ), (8) τ h (V ) 1 m (V ) =, 1 + e V τ m (V ) =.44 + (e V e V ), 1 h (V ) =, 1 + e V τ m (V ) = (e V e V ). (9) These constants were taken at temperature 36 o C and extraneuronular concentration calcium concentration [Ca] = 2 mm. Equilibrium potential E Ca depending on intraneuronular concentration of Ca 2+ are found from the Nernst equation: E Ca = k R T ( ) 2F ln [Ca] ; (1) [Ca] here R = J K/mol o, T = o K, dimensionless constant k = for E Ca is measured in millivolts, extraneuronular concentration of calcium ions is [Ca] = 2mM. Dynamics of Ca 2+ -activated K + -current is determined through the following equations: and I K[Ca] = g K[Ca] m 2 K[Ca] (V i E K ), where m = m m ([Ca]), τ m ([Ca]) τ m ([Ca]) = m ([Ca]) = 1 48[Ca] 2 +.3, 48[Ca] 2 48[Ca] 2 +.3, (11) (12) with g K[Ca] = 1mS/cm 2, E K = 95 mv m being the gating activation variable. Synaptic currents are models using the fast-threshold modulation paradigm 3 : I syn (V i, V j ) = GS(V j θ syn )(V i E syn ), (13)
10 1 where G is the maximal conductance of synaptic current flowing from pre-synaptic j-th neuron into the postsynaptic i-th neuron. For inhibitory coupling we set E syn = 8 mv ; the synaptic activity function S(V (j)) is given by 1 S(V ) =, (14) 1 + e (V θ) with the synaptic threshold θ = mv set in a middle of fast spikes. 1 A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in Nonlinear Sciences, Cambridge Nonlinear Science Series (Cambridge University Press, 3). 2 A. Shilnikov, L. Shilnikov, and D. Turaev, On some mathematical topics in classical synchronization: a tutorial, J Bifurcations and Chaos 14, (4). 3 M. Desroches, J. Guckenheimer, B. Krauskopf, C. Kuehn, H. M. Osinga, and M. Wechselberger, Mixed-mode oscillations with multiple time scales, SIAM Rev. 54, (12). 4 S. Jalil, I. Belykh, and A. Shilnikov, Fast reciprocal inhibition can synchronize bursting neurons, Physical Review E 81 (1), 1.113/PhysRevE I. Omelchenko, M. Rosenblum, and A. Pikovsky, Synchronization of slow-fast systems, European Physical Journal-Special Topics 191, 3 14 (1). 6 S. Jalil, I. Belykh, and A. Shilnikov, Spikes matter for phaselocked bursting in inhibitory neurons, Physical Review E 85 (12), 1.113/PhysRevE J. Wojcik, J. Schwabedal, R. Clewley, and A. Shilnikov, Key bifurcation of bursting polyrhythms in 3-cell central pattern generators, PLoS ONE 9, e92918 (14). 8 A. Selverston, ed., Model Neural Networks and Behavior (Springer, Berlin, 1985). 9 E. Marder and R. Calabrese, Principles of rhythmic motor pattern generation, Physiological Reviews 76, (1996). 1 S. Gillner and P. Wallen, Central pattern generators for locomotion, with special references to vertebrates, Ann. Rev. Neurosci. 8, (1985). 11 M. I. Rabinovich, P. Varona, A. I. Selverston, and H. D. I. Abarbanel, Dynamical principles in neuroscience, Rev. Mod. Phys. 78, (6). 12 K. M. Shaw, D. N. Lyttle, J. P. Gill, M. J. Cullins, J. M. Mc- Manus, H. Lu, P. J. Thomas, and H. J. Chiel, The significance of dynamical architecture for adaptive responses to mechanical loads during rhythmic behavior, Journal of Computational Neuroscience, 1 27 (14). 13 J. T. C. Schwabedal, A. B. Neiman, and A. L. Shilnikov, Robust design of polyrhythmic neural circuits, Phys. Rev. E 9, (14). 14 E. M. Izhikevich, Dynamical Systems in Neuroscience (MIT Press, Cambridge, mass., 7). 15 A. A. V. Hill, S. Van Hooser, and R. L. Calabrese, Half-center oscillators underlying rhythmic movements, in The Handbook of Brain Theory and Neural Networks, edited by M. A. Arbib (The MIT Press, 3). 16 X.-J. Wang and J. Rinzel, Alternating and synchronous rhythms in reciprocally inhibitory model neurons, Ann. Rev. Neurosci. 8, (1985). 17 S. Daun, J. E. Rubin, and I. A. Rybak, Control of oscillation periods and phase durations in half-center central pattern generators: a comparative mechanistic analysis, Journal of Computational Neuroscience 27, 3 36 (9). 18 A. Destexhe, D. Contreras, T. Sejnowski, and M. Steriade, A Model of Spindle Rhythmicity in the Isolated Thalamic Reticular Nucleus, Journal of Neurophysiology 72, (1994). 19 D. Perkel and B. Mulloney, Motor Pattern Production in Reciprocally Inhibitory Neurons Exhibiting Postinhibitory Rebound, Science 185, (1974). F. Skinner, N. Kopell, and E. Marder, Mechanisms for oscillation and frequency control in reciprocally inhibitory model neural networks, Journal of Computational Neuroscience 1, (1994). 21 J. Angstadt, J. Grassmann, K. Theriault, and S. Levasseur, Mechanisms of postinhibitory rebound and its modulation by serotonin in excitatory swim motor neurons of the medicinal leech, Journal of Comparative Physiology A-Neuroethology, Sensory, Neural and Behavioral Physiology 191, (5). 22 N. Kopell and G. Ermentrout, Mechanisms of phase-locking and frequency control in pairs of coupled neural oscillators, Handbook on Dynamical Systems 2, 3 54 (2). 23 L. Shilnikov, A. Shilnikov, D. Turaev, and L. Chua, Methods of Qualitative Theory in Nonlinear Dynamics, Parts I and II (World Scientific Publ., 1998,1). 24 P. Channell, G. Cymbalyuk, and A. Shilnikov, Origin of bursting through homoclinic spike adding in a neuron model, Phys.Rev.Lett. 98, (7). 25 P. Channell, I. Fuwape, A. Neiman, and A. Shilnikov, Variability of bursting patterns in a neuron model in the presence of noise. Journal Computational Neuroscience 27, (9). 26 A. Shilnikov, Complete dynamical analysis of a neuron model, Nonlinear Dynamics 68, (12). 27 S. Jalil, D. Allen, J. Youker, and A. Shilnikov, Toward robust phase-locking in melibe swim central pattern generator models, Chaos 23, 4615 (13). 28 I. Belykh and A. Shilnikov, When weak inhibition synchronizes strongly desynchronizing networks of bursting neurons. Phys Rev Lett 11, 7812 (8). 29 A. Shilnikov, R. Gordon, and I. Belykh, Polyrhythmic synchronization in bursting networking motifs, Chaos 18 (8), 1.163/ N. Kopell and D. Somers, Rapid synchronization through fast threshold modulation, Biol. Cybern. 68, 5 (1993).
Order parameter for bursting polyrhythms in multifunctional central pattern generators
Order parameter for bursting polyrhythms in multifunctional central pattern generators Jeremy Wojcik, Robert Clewley, and Andrey Shilnikov Neuroscience Institute and Department of Mathematics and Statistics,
More informationWhen Transitions Between Bursting Modes Induce Neural Synchrony
International Journal of Bifurcation and Chaos, Vol. 24, No. 8 (2014) 1440013 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414400136 When Transitions Between Bursting Modes Induce
More informationSpikes matter for phase-locked bursting in inhibitory neurons
Spikes matter for phase-locked bursting in inhibitory neurons Sajiya Jalil, Igor Belykh, and Andrey Shilnikov Department of Mathematics and Statistics, Georgia State University, 3 Pryor Street, Atlanta,
More informationMultistability in Bursting Patterns in a Model of a Multifunctional Central Pattern Generator.
Georgia State University ScholarWorks @ Georgia State University Mathematics Theses Department of Mathematics and Statistics 7-15-2009 Multistability in Bursting Patterns in a Model of a Multifunctional
More informationKey Bifurcations of Bursting Polyrhythms in 3-Cell Central Pattern Generators
Georgia State University ScholarWorks @ Georgia State University Neuroscience Institute Faculty Publications Neuroscience Institute 4-04 Key Bifurcations of Bursting Polyrhythms in 3-Cell Central Pattern
More informationNumerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a Mathematical Model of Snail Neurons
International Journal of Theoretical and Mathematical Physics 2015, 5(5): 145-150 DOI: 10.5923/j.ijtmp.20150505.08 Numerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a
More informationA Mechanism of Co-Existence of Bursting and Silent Regimes of Activities of a Neuron
Georgia State University ScholarWorks @ Georgia State University Physics and Astronomy Theses Department of Physics and Astronomy 8-3-2007 A Mechanism of Co-Existence of Bursting and Silent Regimes of
More informationAn Introductory Course in Computational Neuroscience
An Introductory Course in Computational Neuroscience Contents Series Foreword Acknowledgments Preface 1 Preliminary Material 1.1. Introduction 1.1.1 The Cell, the Circuit, and the Brain 1.1.2 Physics of
More informationDynamical Systems in Neuroscience: The Geometry of Excitability and Bursting
Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting Eugene M. Izhikevich The MIT Press Cambridge, Massachusetts London, England Contents Preface xv 1 Introduction 1 1.1 Neurons
More informationThe Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms
The Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms Jiawei Zhang Timothy J. Lewis Department of Mathematics, University of California, Davis Davis, CA 9566,
More information2 1. Introduction. Neuronal networks often exhibit a rich variety of oscillatory behavior. The dynamics of even a single cell may be quite complicated
GEOMETRIC ANALYSIS OF POPULATION RHYTHMS IN SYNAPTICALLY COUPLED NEURONAL NETWORKS J. Rubin and D. Terman Dept. of Mathematics; Ohio State University; Columbus, Ohio 43210 Abstract We develop geometric
More informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More informationNetwork Oscillations Generated by Balancing Graded Asymmetric Reciprocal Inhibition in Passive Neurons
The Journal of Neuroscience, April 1, 1999, 19(7):2765 2779 Network Oscillations Generated by Balancing Graded Asymmetric Reciprocal Inhibition in Passive Neurons Yair Manor, 1 Farzan Nadim, 1 Steven Epstein,
More informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 12, W.W. Lytton, Hodgkin-Huxley Model) 9.1 Passive and active membrane models In the previous lecture we have considered a passive
More informationEvaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons
International Journal of Theoretical and Mathematical Physics 2016, 6(3): 99-103 DOI: 10.5923/j.ijtmp.20160603.02 Evaluating Bistability in a Mathematical Model of Circadian Pacemaker Neurons Takaaki Shirahata
More informationPhase Response Properties of Half-Center. Oscillators
Phase Response Properties of Half-Center Oscillators Jiawei Calvin Zhang Timothy J. Lewis Department of Mathematics, University of California, Davis Davis, CA 95616, USA June 17, 212 Abstract We examine
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More informationTopics in Neurophysics
Topics in Neurophysics Alex Loebel, Martin Stemmler and Anderas Herz Exercise 2 Solution (1) The Hodgkin Huxley Model The goal of this exercise is to simulate the action potential according to the model
More informationConsider the following spike trains from two different neurons N1 and N2:
About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in
More informationA Mathematical Study of Electrical Bursting in Pituitary Cells
A Mathematical Study of Electrical Bursting in Pituitary Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Collaborators on
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More informationCHAPTER 1. Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting
CHAPTER 1 Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting John Guckenheimer Mathematics Department, Cornell University Ithaca, NY 14853 Joseph H. Tien Center for Applied Mathematics,
More informationA model of slow plateau-like oscillations based upon the fast Na current in a window mode
Neurocomputing 38}40 (2001) 159}166 A model of slow plateau-like oscillations based upon the fast Na current in a window mode G.S. Cymbalyuk*, R.L. Calabrese Biology Department, Emory University, 1510
More informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationLinearization of F-I Curves by Adaptation
LETTER Communicated by Laurence Abbott Linearization of F-I Curves by Adaptation Bard Ermentrout Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A. We show that negative
More informationBasic elements of neuroelectronics -- membranes -- ion channels -- wiring
Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wires -- signal propagation -- processing
More informationΝευροφυσιολογία και Αισθήσεις
Biomedical Imaging & Applied Optics University of Cyprus Νευροφυσιολογία και Αισθήσεις Διάλεξη 5 Μοντέλο Hodgkin-Huxley (Hodgkin-Huxley Model) Response to Current Injection 2 Hodgin & Huxley Sir Alan Lloyd
More informationPhase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable
Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a
More informationBursting Oscillations of Neurons and Synchronization
Bursting Oscillations of Neurons and Synchronization Milan Stork Applied Electronics and Telecommunications, Faculty of Electrical Engineering/RICE University of West Bohemia, CZ Univerzitni 8, 3064 Plzen
More informationEffects of synaptic conductance on the voltage distribution and firing rate of spiking neurons
PHYSICAL REVIEW E 69, 051918 (2004) Effects of synaptic conductance on the voltage distribution and firing rate of spiking neurons Magnus J. E. Richardson* Laboratory of Computational Neuroscience, Brain
More informationStochastic Oscillator Death in Globally Coupled Neural Systems
Journal of the Korean Physical Society, Vol. 52, No. 6, June 2008, pp. 19131917 Stochastic Oscillator Death in Globally Coupled Neural Systems Woochang Lim and Sang-Yoon Kim y Department of Physics, Kangwon
More informationInternal and external synchronization of self-oscillating oscillators with non-identical control parameters
Internal and external synchronization of self-oscillating oscillators with non-identical control parameters Emelianova Yu.P., Kuznetsov A.P. Department of Nonlinear Processes, Saratov State University,
More informationMANY scientists believe that pulse-coupled neural networks
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 3, MAY 1999 499 Class 1 Neural Excitability, Conventional Synapses, Weakly Connected Networks, and Mathematical Foundations of Pulse-Coupled Models Eugene
More informationSynchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014
Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,
More informationarxiv:nlin/ v1 [nlin.cd] 13 May 2006
On transition to bursting via deterministic chaos arxiv:nlin/0605030v1 [nlin.cd] 13 May 2006 Georgi S. Medvedev September 23, 2018 Abstract We study statistical properties of the irregular bursting arising
More informationIN THIS turorial paper we exploit the relationship between
508 IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 3, MAY 1999 Weakly Pulse-Coupled Oscillators, FM Interactions, Synchronization, Oscillatory Associative Memory Eugene M. Izhikevich Abstract We study
More information9 Generation of Action Potential Hodgkin-Huxley Model
9 Generation of Action Potential Hodgkin-Huxley Model (based on chapter 2, W.W. Lytton, Hodgkin-Huxley Model) 9. Passive and active membrane models In the previous lecture we have considered a passive
More informationSynaptic dynamics. John D. Murray. Synaptic currents. Simple model of the synaptic gating variable. First-order kinetics
Synaptic dynamics John D. Murray A dynamical model for synaptic gating variables is presented. We use this to study the saturation of synaptic gating at high firing rate. Shunting inhibition and the voltage
More informationCORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE. by Pamela Reitsma. B.S., University of Maine, 2007
CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE by Pamela Reitsma B.S., University of Maine, 27 Submitted to the Graduate Faculty of the Department of Mathematics in partial
More informationChapter 48 Neurons, Synapses, and Signaling
Chapter 48 Neurons, Synapses, and Signaling Concept 48.1 Neuron organization and structure reflect function in information transfer Neurons are nerve cells that transfer information within the body Neurons
More informationConductance-Based Integrate-and-Fire Models
NOTE Communicated by Michael Hines Conductance-Based Integrate-and-Fire Models Alain Destexhe Department of Physiology, Laval University School of Medicine, Québec, G1K 7P4, Canada A conductance-based
More informationDynamical phase transitions in periodically driven model neurons
Dynamical phase transitions in periodically driven model neurons Jan R. Engelbrecht 1 and Renato Mirollo 2 1 Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA 2 Department
More informationNeurophysiology. Danil Hammoudi.MD
Neurophysiology Danil Hammoudi.MD ACTION POTENTIAL An action potential is a wave of electrical discharge that travels along the membrane of a cell. Action potentials are an essential feature of animal
More information6.3.4 Action potential
I ion C m C m dφ dt Figure 6.8: Electrical circuit model of the cell membrane. Normally, cells are net negative inside the cell which results in a non-zero resting membrane potential. The membrane potential
More informationThe Spike Response Model: A Framework to Predict Neuronal Spike Trains
The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology
More information1 Introduction and neurophysiology
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 (2009) 535-549 Copyright c 2009 Watam Press http://www.watam.org ASYMPTOTIC DYNAMICS OF THE SLOW-FAST HINDMARSH-ROSE
More informationChapter 9. Nerve Signals and Homeostasis
Chapter 9 Nerve Signals and Homeostasis A neuron is a specialized nerve cell that is the functional unit of the nervous system. Neural signaling communication by neurons is the process by which an animal
More informationMathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued
Mathematical Foundations of Neuroscience - Lecture 3. Electrophysiology of neurons - continued Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland
More informationControl and Integration. Nervous System Organization: Bilateral Symmetric Animals. Nervous System Organization: Radial Symmetric Animals
Control and Integration Neurophysiology Chapters 10-12 Nervous system composed of nervous tissue cells designed to conduct electrical impulses rapid communication to specific cells or groups of cells Endocrine
More informationSelf-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire Neuron Model Based on Small World Networks
Commun. Theor. Phys. (Beijing, China) 43 (2005) pp. 466 470 c International Academic Publishers Vol. 43, No. 3, March 15, 2005 Self-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire
More informationOverview Organization: Central Nervous System (CNS) Peripheral Nervous System (PNS) innervate Divisions: a. Afferent
Overview Organization: Central Nervous System (CNS) Brain and spinal cord receives and processes information. Peripheral Nervous System (PNS) Nerve cells that link CNS with organs throughout the body.
More informationStrange Nonchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley Neuron
Journal of the Korean Physical Society, Vol. 57, o. 1, July 2010, pp. 23 29 Strange onchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley euron Woochang Lim and Sang-Yoon Kim Department of Physics,
More informationLecture 11 : Simple Neuron Models. Dr Eileen Nugent
Lecture 11 : Simple Neuron Models Dr Eileen Nugent Reading List Nelson, Biological Physics, Chapter 12 Phillips, PBoC, Chapter 17 Gerstner, Neuronal Dynamics: from single neurons to networks and models
More informationThe Mixed States of Associative Memories Realize Unimodal Distribution of Dominance Durations in Multistable Perception
The Mixed States of Associative Memories Realize Unimodal Distribution of Dominance Durations in Multistable Perception Takashi Kanamaru Department of Mechanical Science and ngineering, School of Advanced
More informationarxiv: v3 [q-bio.nc] 18 Dec 2018
A biophysical model explains the spontaneous bursting behavior in the developing retina Dora Matzakou-Karvouniari 1, Lionel Gil 2, Elaine Orendorff 3, Olivier Marre 3, Serge Picaud 3, and Bruno Cessac
More informationLocalized activity patterns in excitatory neuronal networks
Localized activity patterns in excitatory neuronal networks Jonathan Rubin Amitabha Bose February 3, 2004 Abstract. The existence of localized activity patterns, or bumps, has been investigated in a variety
More informationCanonical Neural Models 1
Canonical Neural Models 1 Frank Hoppensteadt 1 and Eugene zhikevich 2 ntroduction Mathematical modeling is a powerful tool in studying fundamental principles of information processing in the brain. Unfortunately,
More informationAction Potentials and Synaptic Transmission Physics 171/271
Action Potentials and Synaptic Transmission Physics 171/271 Flavio Fröhlich (flavio@salk.edu) September 27, 2006 In this section, we consider two important aspects concerning the communication between
More informationEmergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity
Emergence of resonances in neural systems: the interplay between adaptive threshold and short-term synaptic plasticity Jorge F. Mejias 1,2 and Joaquín J. Torres 2 1 Department of Physics and Center for
More informationSingle neuron models. L. Pezard Aix-Marseille University
Single neuron models L. Pezard Aix-Marseille University Biophysics Biological neuron Biophysics Ionic currents Passive properties Active properties Typology of models Compartmental models Differential
More informationNeurons and Nervous Systems
34 Neurons and Nervous Systems Concept 34.1 Nervous Systems Consist of Neurons and Glia Nervous systems have two categories of cells: Neurons, or nerve cells, are excitable they generate and transmit electrical
More informationAll-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model
All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model S. A. Sadegh Zadeh, C. Kambhampati International Science Index, Mathematical and Computational Sciences waset.org/publication/10008281
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/319/5869/1543/dc1 Supporting Online Material for Synaptic Theory of Working Memory Gianluigi Mongillo, Omri Barak, Misha Tsodyks* *To whom correspondence should be addressed.
More informationReduction of Conductance Based Models with Slow Synapses to Neural Nets
Reduction of Conductance Based Models with Slow Synapses to Neural Nets Bard Ermentrout October 1, 28 Abstract The method of averaging and a detailed bifurcation calculation are used to reduce a system
More informationSynchrony in Neural Systems: a very brief, biased, basic view
Synchrony in Neural Systems: a very brief, biased, basic view Tim Lewis UC Davis NIMBIOS Workshop on Synchrony April 11, 2011 components of neuronal networks neurons synapses connectivity cell type - intrinsic
More informationBiological Modeling of Neural Networks
Week 4 part 2: More Detail compartmental models Biological Modeling of Neural Networks Week 4 Reducing detail - Adding detail 4.2. Adding detail - apse -cable equat Wulfram Gerstner EPFL, Lausanne, Switzerland
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best Koch, Chapters 6, 8, 9 See also Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913
More informationMembrane Potentials, Action Potentials, and Synaptic Transmission. Membrane Potential
Cl Cl - - + K + K+ K + K Cl - 2/2/15 Membrane Potentials, Action Potentials, and Synaptic Transmission Core Curriculum II Spring 2015 Membrane Potential Example 1: K +, Cl - equally permeant no charge
More informationNervous Systems: Neuron Structure and Function
Nervous Systems: Neuron Structure and Function Integration An animal needs to function like a coherent organism, not like a loose collection of cells. Integration = refers to processes such as summation
More informationNeurons, Synapses, and Signaling
Chapter 48 Neurons, Synapses, and Signaling PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions
More informationACTION POTENTIAL. Dr. Ayisha Qureshi Professor MBBS, MPhil
ACTION POTENTIAL Dr. Ayisha Qureshi Professor MBBS, MPhil DEFINITIONS: Stimulus: A stimulus is an external force or event which when applied to an excitable tissue produces a characteristic response. Subthreshold
More informationSpontaneous voltage oscillations and response dynamics of a Hodgkin-Huxley type model of sensory hair cells
Journal of Mathematical Neuroscience (2011) 1:11 DOI 10.1186/2190-8567-1-11 RESEARCH OpenAccess Spontaneous voltage oscillations and response dynamics of a Hodgkin-Huxley type model of sensory hair cells
More informationIntroduction and the Hodgkin-Huxley Model
1 Introduction and the Hodgkin-Huxley Model Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference:
More informationPhase Response Properties and Phase-Locking in Neural Systems with Delayed Negative-Feedback. Carter L. Johnson
Phase Response Properties and Phase-Locking in Neural Systems with Delayed Negative-Feedback Carter L. Johnson Faculty Mentor: Professor Timothy J. Lewis University of California, Davis Abstract Oscillatory
More informationMATH 3104: THE HODGKIN-HUXLEY EQUATIONS
MATH 3104: THE HODGKIN-HUXLEY EQUATIONS Parallel conductance model A/Prof Geoffrey Goodhill, Semester 1, 2009 So far we have modelled neuronal membranes by just one resistance (conductance) variable. We
More informationPatterns of Synchrony in Neural Networks with Spike Adaptation
Patterns of Synchrony in Neural Networks with Spike Adaptation C. van Vreeswijky and D. Hanselz y yracah Institute of Physics and Center for Neural Computation, Hebrew University, Jerusalem, 9194 Israel
More informationTitle. Author(s)Fujii, Hiroshi; Tsuda, Ichiro. CitationNeurocomputing, 58-60: Issue Date Doc URL. Type.
Title Neocortical gap junction-coupled interneuron systems exhibiting transient synchrony Author(s)Fujii, Hiroshi; Tsuda, Ichiro CitationNeurocomputing, 58-60: 151-157 Issue Date 2004-06 Doc URL http://hdl.handle.net/2115/8488
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More informationSingle-Compartment Neural Models
Single-Compartment Neural Models BENG/BGGN 260 Neurodynamics University of California, San Diego Week 2 BENG/BGGN 260 Neurodynamics (UCSD) Single-Compartment Neural Models Week 2 1 / 18 Reading Materials
More informationDESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS
Letters International Journal of Bifurcation and Chaos, Vol. 8, No. 8 (1998) 1733 1738 c World Scientific Publishing Company DESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS I. P.
More informationA Neuromorphic VLSI System for Modeling the Neural Control of Axial Locomotion
A Neuromorphic VLSI System for Modeling the Neural Control of Axial Locomotion Girish N. Patel girish@ece.gatech.edu Edgar A. Brown ebrown@ece.gatech.edu Stephen P. De Weerth steved@ece.gatech.edu School
More informationNEURONS, SENSE ORGANS, AND NERVOUS SYSTEMS CHAPTER 34
NEURONS, SENSE ORGANS, AND NERVOUS SYSTEMS CHAPTER 34 KEY CONCEPTS 34.1 Nervous Systems Are Composed of Neurons and Glial Cells 34.2 Neurons Generate Electric Signals by Controlling Ion Distributions 34.3
More informationMathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells
Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University
More informationMEMBRANE POTENTIALS AND ACTION POTENTIALS:
University of Jordan Faculty of Medicine Department of Physiology & Biochemistry Medical students, 2017/2018 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Review: Membrane physiology
More informationSynchronization in two neurons: Results for a two-component dynamical model with time-delayed inhibition
Physica D 4 (998 47 7 Synchronization in two neurons: Results for a two-component dynamical model with time-delayed inhibition Gilles Renversez Centre de Physique Théorique, UPR 76 du Centre National de
More informationEffects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell
Effects of Betaxolol on Hodgkin-Huxley Model of Tiger Salamander Retinal Ganglion Cell 1. Abstract Matthew Dunlevie Clement Lee Indrani Mikkilineni mdunlevi@ucsd.edu cll008@ucsd.edu imikkili@ucsd.edu Isolated
More information/639 Final Solutions, Part a) Equating the electrochemical potentials of H + and X on outside and inside: = RT ln H in
580.439/639 Final Solutions, 2014 Question 1 Part a) Equating the electrochemical potentials of H + and X on outside and inside: RT ln H out + zf 0 + RT ln X out = RT ln H in F 60 + RT ln X in 60 mv =
More informationElectrodiffusion Model of Electrical Conduction in Neuronal Processes
i. ISMS OF CONDlTlONlff PLASTICITY ditecl by Charks B. Woody, Daniel L. Alk~n, and James b. McGauyh (Plenum Publishing Corporation, 19 23 Electrodiffusion Model of Electrical Conduction in Neuronal Processes
More informationCoordination of Cellular Pattern-Generating Circuits that Control Limb Movements: The Sources of Stable Differences in Intersegmental Phases
The Journal of Neuroscience, April 15, 2003 23(8):3457 3468 3457 Coordination of Cellular Pattern-Generating Circuits that Control Limb Movements: The Sources of Stable Differences in Intersegmental Phases
More informationThis script will produce a series of pulses of amplitude 40 na, duration 1ms, recurring every 50 ms.
9.16 Problem Set #4 In the final problem set you will combine the pieces of knowledge gained in the previous assignments to build a full-blown model of a plastic synapse. You will investigate the effects
More informationSynchronization of Elliptic Bursters
SIAM REVIEW Vol. 43,No. 2,pp. 315 344 c 2001 Society for Industrial and Applied Mathematics Synchronization of Elliptic Bursters Eugene M. Izhikevich Abstract. Periodic bursting behavior in neurons is
More informationModeling neural oscillations
Physiology & Behavior 77 (2002) 629 633 Modeling neural oscillations G. Bard Ermentrout*, Carson C. Chow Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, USA Received 30 July
More informationPlasticity and learning in a network of coupled phase oscillators
PHYSICAL REVIEW E, VOLUME 65, 041906 Plasticity and learning in a network of coupled phase oscillators Philip Seliger, Stephen C. Young, and Lev S. Tsimring Institute for onlinear Science, University of
More informationSynchronization and Phase Oscillators
1 Synchronization and Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Synchronization
More informationStochastic resonance in the absence and presence of external signals for a chemical reaction
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 Stochastic resonance in the absence and presence of external signals for a chemical reaction Lingfa Yang, Zhonghuai Hou, and Houwen Xin
More informationSingle-Cell and Mean Field Neural Models
1 Single-Cell and Mean Field Neural Models Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 The neuron
More informationNonlinear systems, chaos and control in Engineering
Nonlinear systems, chaos and control in Engineering Module 1 block 3 One-dimensional nonlinear systems Cristina Masoller Cristina.masoller@upc.edu http://www.fisica.edu.uy/~cris/ Schedule Flows on the
More informationarxiv: v1 [q-bio.nc] 13 Feb 2018
Gain control with A-type potassium current: I A as a switch between divisive and subtractive inhibition Joshua H Goldwyn 1*, Bradley R Slabe 2, Joseph B Travers 3, David Terman 2 arxiv:182.4794v1 [q-bio.nc]
More information