The assembly of ionic currents in a thalamic neuron. II. The stability and state diagrams

In the previous model of a thalamic neuron (R.M. Rose & J.L. Hindmarsh, Proc. R. Soc. Lond. B237, 267-288 (1989], which we referred to as the z-model, the burst response was terminated by the slow activation of a subthreshold outward current. In this paper we show that similar results can be obtained if the burst response is terminated by slow inactivation of the subthreshold inward current, Isa. We illustrate the use of this new model, which we refer to as the ha-model, by using it to explain the response of a thalamic neuron to a double ramp current. The main aim of the paper is to show how the stability and state diagrams introduced previously can be used to explain various types of firing pattern of thalamic and other neurons. We show that increasing the threshold for the fast action potentials leads to low threshold spikes of increased amplitude. Also, addition of a second subthreshold inward current adds a new stability region, which enables us to explain the origin of plateau potentials. In addition, various types of subthreshold oscillation are produced by relocating a previously stable equilibrium point in an unstable region. Finally, we predict a sequence of responses to current steps from different levels of background current that extends the burst, rest, tonic sequence of thalamic neurons. The stability and state diagrams therefore provide us with a useful way of explaining further properties of thalamic neurons and appear to have further applications to other mammalian neurons.     

State-Dependent Effects of Na Channel Noise on Neuronal Burst Generation

We explore the effects of stochastic sodium (Na) channel activation on the variability and dynamics of spiking and bursting in a model neuron. The complete model segregates Hodgin-Huxley-type currents into two compartments, and undergoes applied current-dependent bifurcations between regimes of periodic bursting, chaotic bursting, and tonic spiking. Noise is added to simulate variable, finite sizes of the population of Na channels in the fast spiking compartment.During tonic firing, Na channel noise causes variability in interspike intervals (ISIs). The variance, as well as the sensitivity to noise, depend on the model's biophysical complexity. They are smallest in an isolated spiking compartment; increase significantly upon coupling to a passive compartment; and increase again when the second compartment also includes slow-acting currents. In this full model, sufficient noise can convert tonic firing into bursting.During bursting, the actions of Na channel noise are state-dependent. The higher the noise level, the greater the jitter in spike timing within bursts. The noise makes the burst durations of periodic regimes variable, while decreasing burst length duration and variance in a chaotic regime. Na channel noise blurs the sharp transitions of spike time and burst length seen at the bifurcations of the noise-free model. Close to such a bifurcation, the burst behaviors of previously periodic and chaotic regimes become essentially indistinguishable.We discuss biophysical mechanisms, dynamical interpretations and physiological implications. We suggest that noise associated with finite populations of Na channels could evoke very different effects on the intrinsic variability of spiking and bursting discharges, depending on a biological neuron's complexity and applied current-dependent state. We find that simulated channel noise in the model neuron qualitatively replicates the observed variability in burst length and interburst interval in an isolated biological bursting neuron.     

Vasopressin induces depolarization and state-dependent firing patterns in rat thalamic paraventricular nucleus neurons in vitro

The thalamic midline paraventricular nucleus (PVT) is prominently innervated by vasopressin-immunoreactive neurons from the suprachiasmatic nucleus (SCN), site of the brain's biological clock. Using patch-clamp recordings in slice preparations taken from Wistar rats during the subjective day, we examined 90 PVT neurons for responses to bath-applied AVP (0.5-2 microM; 1-3 min). In current clamp at resting membrane potentials (-65 +/- 1 mV), PVT neurons displayed low-threshold spikes (LTSs) and burst firing patterns. In 50% of cells tested, AVP induced a slowly rising, prolonged membrane depolarization and tonic firing, returning to burst firing upon recovery. AVP modulated hyperpolarization-activated LTSs by decreasing the time to the initial sodium spike at the onset of LTS, also increasing the duration of the afterdepolarization. Responses were blockable with a V(1a) receptor antagonist (Manning compound). Under voltage clamp, AVP induced a TTX-resistant, slowly rising, and prolonged (approximately 15 min) inward current (<40 pA). Current-voltage relationship (I-V) analyses of the AVP responses revealed a decrease in membrane conductance to 73.1 +/- 6.2% of control, with net AVP current reversing at -106 +/- 4 mV, and decreased inward rectification at negative potentials. These observations are consistent with an AVP-induced closure of an inwardly rectifying potassium conductance. On the basis of these in vitro observations, we suggest that the SCN vasopressinergic innervation of PVT is excitatory in nature, possibly releasing AVP with circadian rhythmicity and contributing to state-dependent firing patterns in PVT neurons over the sleep-wake cycle.     

Ionic mechanisms underlying spontaneous CA1 neuronal firing in Ca2+-free solution

Hippocampal CA1 neurons exposed to zero-[Ca(2+)] solutions can generate periodic spontaneous synchronized activity in the absence of synaptic function. Experiments using hippocampal slices showed that, after exposure to zero-[Ca(2+)](0) solution, CA1 pyramidal cells depolarized 5-10 mV and started firing spontaneous action potentials. Spontaneous single neuron activity appeared in singlets or was grouped into bursts of two or three action potentials. A 16-compartment, 23-variable cable model of a CA1 pyramidal neuron was developed to study mechanisms of spontaneous neuronal bursting in a calcium-free extracellular solution. In the model, five active currents (a fast sodium current, a persistent sodium current, an A-type transient potassium current, a delayed rectifier potassium current, and a muscarinic potassium current) are included in the somatic compartment. The model simulates the spontaneous bursting behavior of neurons in calcium-free solutions. The mechanisms underlying several aspects of bursting are studied, including the generation of triplet bursts, spike duration, burst termination, after-depolarization behavior, and the prolonged inactive period between bursts. We show that the small persistent sodium current can play a key role in spontaneous CA1 activity in zero-calcium solutions. In particular, it is necessary for the generation of an after-depolarizing potential and prolongs both individual bursts and the interburst interval.     

A model of a thalamic neuron

We modify our recent three equilibrium-point model of neuronal bursting by a means of a small deformation of the nullclines in the x-y phase plane to give a model that can have as many as five equilibrium points. In this model the middle stable equilibrium point (e.p.) is separated from the outer stable and unstable e.ps by two saddle points. If the system is started at rest at the middle stable e.p. it has the following complex properties: A short suprathreshold current pulse switches the model from a silent state to a bursting state, or to give a single burst, depending on the choice of parameters. A subthreshold depolarizing current step gives a passive response at rest, but if the model is either constantly hyperpolarized or constantly depolarized, then the same current step gives different active responses. At a hyperpolarized level this consists of a burst response that shows refractoriness. At a depolarized level it consists of tonic firing with a linear frequency--current relationship. Hyperpolarization from rest is followed by post-inhibitory rebound. The model responds in a unique and characteristic way to an applied current ramp. These properties are very similar to those that have been recently recorded intracellularly from neurons in the mammalian thalamus. In the x-y phase plane our models of the repetitively firing neuron, the bursting neuron and the thalamic neuron form a progression of models in which the y nullcline in the subthreshold region is deformed once to give the burst neuron model, and a second time to give the thalamic neuron model. Each deformation can be interpreted as corresponding to the inclusion of a slow inward current in the model. As these currents are included so the associated firing properties increase in complexity.     

Ion channel density and threshold dynamics of repetitive firing in a cortical neuron model

Modifying the density and distribution of ion channels in a neuron (by natural up- and down-regulation, by pharmacological intervention or by spontaneous mutations) changes its activity pattern. In the present investigation, we analyze how the impulse patterns are regulated by the density of voltage-gated channels in a model neuron, based on voltage clamp measurements of hippocampal interneurons. At least three distinct oscillatory patterns, associated with three distinct regions in the Na-K channel density plane, were found. A stability analysis showed that the different regions are characterized by saddle-node, double-orbit, and Hopf bifurcation threshold dynamics, respectively. Single strongly graded action potentials occur in an area outside the oscillatory regions, but less graded action potentials occur together with repetitive firing over a considerable range of channel densities. The presently found relationship between channel densities and oscillatory behavior may be relevance for understanding principal spiking patterns of cortical neurons (regular firing and fast spiking). It may also be of relevance for understanding the action of pharmacological compounds on brain oscillatory activity.     

Spike-Threshold Adaptation Predicted by Membrane Potential Dynamics In Vivo

Neurons encode information in sequences of spikes, which are triggered when their membrane potential crosses a threshold. In vivo, the spiking threshold displays large variability suggesting that threshold dynamics have a profound influence on how the combined input of a neuron is encoded in the spiking. Threshold variability could be explained by adaptation to the membrane potential. However, it could also be the case that most threshold variability reflects noise and processes other than threshold adaptation. Here, we investigated threshold variation in auditory neurons responses recorded in vivo in barn owls. We found that spike threshold is quantitatively predicted by a model in which the threshold adapts, tracking the membrane potential at a short timescale. As a result, in these neurons, slow voltage fluctuations do not contribute to spiking because they are filtered by threshold adaptation. More importantly, these neurons can only respond to input spikes arriving together on a millisecond timescale. These results demonstrate that fast adaptation to the membrane potential captures spike threshold variability in vivo.     

Reliability of spike and burst firing in thalamocortical relay cells

The reliability and precision of the timing of spikes in a spike train is an important aspect of neuronal coding. We investigated reliability in thalamocortical relay (TCR) cells in the acute slice and also in a Morris-Lecar model with several extensions. A frozen Gaussian noise current, superimposed on a DC current, was injected into the TCR cell soma. The neuron responded with spike trains that showed trial-to-trial variability, due to amongst others slow changes in its internal state and the experimental setup. The DC current allowed to bring the neuron in different states, characterized by a well defined membrane voltage (between ?80 and ?50 mV) and by a specific firing regime that on depolarization gradually shifted from a predominantly bursting regime to a tonic spiking regime. The filtered frozen white noise generated a spike pattern output with a broad spike interval distribution. The coincidence factor and the Hunter and Milton measure were used as reliability measures of the output spike train. In the experimental TCR cell as well as the Morris-Lecar model cell the reliability depends on the shape (steepness) of the current input versus spike frequency output curve. The model also allowed to study the contribution of three relevant ionic membrane currents to reliability: a T-type calcium current, a cation selective h-current and a calcium dependent potassium current in order to allow bursting, investigate the consequences of a more complex current-frequency relation and produce realistic firing rates. The reliability of the output of the TCR cell increases with depolarization. In hyperpolarized states bursts are more reliable than single spikes. The analytically derived relations were capable to predict several of the experimentally recorded spike features.     

The potassium A-current,low firing rates and rebound excitation in Hodgkin-Huxley models

It is widely believed, following the work of Connor and Stevens (1971,J. Physiol. Lond. 214, 31–53) that the ability to fire action potentials over a wide frequency range, especially down to very low rates, is due to the transient, potassium A-current (I _A). Using a reduction of the classical Hodgkin-Huxley model, we study the effects ofI _A on steady firing rate, especially in the near-threshold regime for the onset of firing. A minimum firing rate of zero corresponds to a homoclinic bifurcation of periodic solutions at a critical level of stimulating current. It requires that the membrane's steady-state current-voltage relation be N-shaped rather than monotonic. For experimentally based genericI _A parameters, the model does not fire at arbitrarily low rates, although it can for the more atypicalI _A parameters given by Connor and Stevens for the crab axon. When theI _A inactivation rate is slow, we find that the transient potassium current can mediate more complex firing patterns, such as periodic bursting in some parameter regimes. The number of spikes per burst increases asg _A decreases and as inactivation rate decreases. We also study howI _A affects properties of transient voltage responses, such as threshold and firing latency for anodal break excitation. We provide mathematical explanations for several of these dynamic behaviors using bifurcation theory and averaging methods.     

LTS and FS inhibitory interneurons, short-term synaptic plasticity, and cortical circuit dynamics

Somatostatin-expressing, low threshold-spiking (LTS) cells and fast-spiking (FS) cells are two common subtypes of inhibitory neocortical interneuron. Excitatory synapses from regular-spiking (RS) pyramidal neurons to LTS cells strongly facilitate when activated repetitively, whereas RS-to-FS synapses depress. This suggests that LTS neurons may be especially relevant at high rate regimes and protect cortical circuits against over-excitation and seizures. However, the inhibitory synapses from LTS cells usually depress, which may reduce their effectiveness at high rates. We ask: by which mechanisms and at what firing rates do LTS neurons control the activity of cortical circuits responding to thalamic input, and how is control by LTS neurons different from that of FS neurons? We study rate models of circuits that include RS cells and LTS and FS inhibitory cells with short-term synaptic plasticity. LTS neurons shift the RS firing-rate vs. current curve to the right at high rates and reduce its slope at low rates; the LTS effect is delayed and prolonged. FS neurons always shift the curve to the right and affect RS firing transiently. In an RS-LTS-FS network, FS neurons reach a quiescent state if they receive weak input, LTS neurons are quiescent if RS neurons receive weak input, and both FS and RS populations are active if they both receive large inputs. In general, FS neurons tend to follow the spiking of RS neurons much more closely than LTS neurons. A novel type of facilitation-induced slow oscillations is observed above the LTS firing threshold with a frequency determined by the time scale of recovery from facilitation. To conclude, contrary to earlier proposals, LTS neurons affect the transient and steady state responses of cortical circuits over a range of firing rates, not only during the high rate regime; LTS neurons protect against over-activation about as well as FS neurons.     

Channel density regulation of firing patterns in a cortical neuron model

Modifying the density and distribution of ion channels in a neuron (by natural up- and downregulation or by pharmacological intervention or by spontaneous mutations) changes its activity pattern. In this investigation we analyzed how the impulse patterns are regulated by the density of voltage-gated channels in a neuron model based on voltage-clamp measurements of hippocampal interneurons. At least three distinct oscillatory patterns, associated with three distinct regions in the Na-K channel density plane, were found. A stability analysis showed that the different regions are characterized by saddle-node, double-orbit, and Hopf-bifurcation threshold dynamics, respectively. Single, strongly graded action potentials occur in an area outside the oscillatory regions, but less graded action potentials occur together with repetitive firing over a considerable range of channel densities. The relationship found here between channel densities and oscillatory behavior may partly explain the difference between the principal spiking patterns previously described for crab axons (class 1 and 2) and cortical neurons (regular firing and fast spiking).     

Ion Channel Density Regulates Switches between Regular and Fast Spiking in Soma but Not in Axons

The threshold firing frequency of a neuron is a characterizing feature of its dynamical behaviour, in turn determining its role in the oscillatory activity of the brain. Two main types of dynamics have been identified in brain neurons. Type 1 dynamics (regular spiking) shows a continuous relationship between frequency and stimulation current (f-I_stim) and, thus, an arbitrarily low frequency at threshold current; Type 2 (fast spiking) shows a discontinuous f-I_stim relationship and a minimum threshold frequency. In a previous study of a hippocampal neuron model, we demonstrated that its dynamics could be of both Type 1 and Type 2, depending on ion channel density. In the present study we analyse the effect of varying channel density on threshold firing frequency on two well-studied axon membranes, namely the frog myelinated axon and the squid giant axon. Moreover, we analyse the hippocampal neuron model in more detail. The models are all based on voltage-clamp studies, thus comprising experimentally measurable parameters. The choice of analysing effects of channel density modifications is due to their physiological and pharmacological relevance. We show, using bifurcation analysis, that both axon models display exclusively Type 2 dynamics, independently of ion channel density. Nevertheless, both models have a region in the channel-density plane characterized by an N-shaped steady-state current-voltage relationship (a prerequisite for Type 1 dynamics and associated with this type of dynamics in the hippocampal model). In summary, our results suggest that the hippocampal soma and the two axon membranes represent two distinct kinds of membranes; membranes with a channel-density dependent switching between Type 1 and 2 dynamics, and membranes with a channel-density independent dynamics. The difference between the two membrane types suggests functional differences, compatible with a more flexible role of the soma membrane than that of the axon membrane.     

成年蜜蜂脑神经细胞的培养和电生理特征

为了研究杀虫剂等对蜜蜂毒性作用的神经机制,需在体外建立成年蜜蜂脑神经细胞的分离培养和电生理记录技术并研究其正常电生理特征,而对成年蜜蜂脑神经细胞的分离培养和电生理特性的研究报道甚少。我们采用酶解和机械吹打相结合的方法获得了数量较多且活力较好的成年意大利蜜蜂Apis mellifera脑神经细胞,并用全细胞膜片钳技术研究了成年意大利蜜蜂脑神经细胞对电流和电压刺激的反应,获得了成年意蜂脑神经细胞的基本电生理特征以及钠电流和钾电流的特性。全细胞电流钳的记录结果表明,在体外培养条件下,细胞无自发放电发生,注射电流后仅引起细胞单次放电,引起细胞放电的阈电流平均为60.8±63 pA; 细胞动作电位产生的阈电位平均为−27.4±2.3 mV。用全细胞电压钳记录了神经细胞的钠电流和钾电流。钠电流的分离是在电压刺激下通过阻断钾通道和钙通道实现。细胞的内向钠电流在指令电压为−40～−30 mV左右激活,−10 mV达峰值,钠通道的稳态失活电压V_1/2为−58.4 mV; 外向钾电流成份至少包括较小的快速失活钾电流和和较大的缓慢失活钾电流（占总钾电流的80%）,其半激活膜电位V_1/2为3.86 mV,无明显的稳态失活。结果提示缓慢失活钾电流的特征可能是细胞单次放电的机制之一。     

Processing of calling songs by an L-shaped neuron in the prothoracic ganglion of the female cricket, Acheta domesticus

ABSTRACT. An L-shaped auditory intemeuron (LI) has been recorded from extracellularly and intracellularly, and identified morphologically (by Lucifer yellow or cobalt injection) in the prothoracic ganglion of mature female Acheta domesticus. The morphology of the LI is very similar to ascending, prothoracic acoustic interneurons that are most sensitive to higher carrier frequencies in both A. domesticus and other gryllid species. Its terminations in the brain are similar to ascending acoustic interneurons found in other gryllids. The LI neuron is most sensitive to 4–5 kHz model calling songs (CSs), the main carrier frequency of the natural call. Thresholds to high frequencies (8–15 kHz) are 15–20 dB higher. Increasing CS intensities of up to 15 dB above threshold at 4–5 kHz result in increased firing rates by the LI. More than 15 dB increase in intensity causes saturation with little increase in spiking rate until the intensity surpasses 80 dB. In response to 70 dB or higher stimulus intensities, the LI responds to the second and third CS syllables with one or two spikes, pauses, and then produces a burst of nerve impulses with the same or greater latency than for lower intensity stimuli. In response to CS syllables of changing duration (10–30 ms) this neuron responds with a rather constant duration burst of impulses. Syllable periods of the CS stimuli were accurately encoded by the LI. Progressively stronger injection of hyperpolarizing current reduces, and ultimately stops spiking of the LI in response to CS stimuli. More intense stimulation with reduced hyperpolarization shows an initial spike, pause and burst of spikes. Intracellular recording from axonal regions of the neuron shows large spikes, small EPSPs and a developing hyperpolarization through the response to a CS chirp. Inhibitory input to the LI is demonstrated at 4.5, 8 and 16 kHz. This probably explains the specialized response characteristics of the LI which enhanced its encoding of CS syllable period.     

Ionic currents influencing spontaneous firing and pacemaker frequency in dopamine neurons of the ventrolateral periaqueductal gray and dorsal raphe nucleus (vlPAG/DRN): A voltage-clamp and computational modelling study

Dopamine (DA) neurons of the ventrolateral periaqueductal gray (vlPAG) and dorsal raphe nucleus (DRN) fire spontaneous action potentials (APs) at slow, regular patterns in vitro but a detailed account of their intrinsic membrane properties responsible for spontaneous firing is currently lacking. To resolve this, we performed a voltage-clamp electrophysiological study in brain slices to describe their major ionic currents and then constructed a computer model and used simulations to understand the mechanisms behind autorhythmicity in silico. We found that vlPAG/DRN DA neurons exhibit a number of voltage-dependent currents activating in the subthreshold range including, a hyperpolarization-activated cation current (I_H), a transient, A-type, potassium current (I_A), a background, ‘persistent’ (I_NaP) sodium current and a transient, low voltage activated (LVA) calcium current (I_CaLVA). Brain slice pharmacology, in good agreement with computer simulations, showed that spontaneous firing occurred independently of I_H, I_A or calcium currents. In contrast, when blocking sodium currents, spontaneous firing ceased and a stable, non-oscillating membrane potential below AP threshold was attained. Using the DA neuron model we further show that calcium currents exhibit little activation (compared to sodium) during the interspike interval (ISI) repolarization while, any individual potassium current alone, whose blockade positively modulated AP firing frequency, is not required for spontaneous firing. Instead, blockade of a number of potassium currents simultaneously is necessary to eliminate autorhythmicity. Repolarization during ISI is mediated initially via the deactivation of the delayed rectifier potassium current, while a sodium background ‘persistent’ current is essentially indispensable for autorhythmicity by driving repolarization towards AP threshold.     

Bursting as an Effective Relay Mode in a Minimal Thalamic Model

In recent years, accumulating evidence indicates that thalamic bursts are present during wakefulness and participate in information transmission as an effective relay mode with distinctive properties from the tonic activity. Thalamic bursts originate from activation of the low threshold calcium cannels via a local feedback inhibition, exerted by the thalamic reticular neurons upon the relay neurons. This article, examines if this simple mechanism is sufficient to explain the distinctive properties of thalamic bursting as an effective relay mode. A minimal model of thalamic circuit composed of a retinal spike train, a relay neuron and a reticular neuron is simulated to generate the tonic and burst firing modes. The integrate-and-fire-or-burst model is used to simulate the neurons. After discriminating the burst events with criteria based on inter-spike-intervals, statistical indices show that the bursts of the minimal model are stereotypic events. The relation between the rate of bursts and the parameters of the input spike train demonstrates marked nonlinearities. Burst response is shown to be selective to spike-silence-spike sequences in the input spike train. Moreover, burst events represent the input more reliably than the tonic spike in a considerable range of the parameters of the model. In conclusion, many of the distinctive properties of thalamic bursts such as stereotypy, nonlinear dependence on the sensory stimulus, feature selectivity and reliability are reproducible in the minimal model. Furthermore, the minimal model predicts that while the bursts are more frequent in the spike train of the off-center X relay neurons (corresponding to off-center X retinal ganglion cells), they are more reliable when generated by the on-center ones (corresponding to on-center X ganglion cells).     

Variability of bursting patterns in a neuron model in the presence of noise

Spiking and bursting patterns of neurons are characterized by a high degree of variability. A single neuron can demonstrate endogenously various bursting patterns, changing in response to external disturbances due to synapses, or to intrinsic factors such as channel noise. We argue that in a model of the leech heart interneuron existing variations of bursting patterns are significantly enhanced by a small noise. In the absence of noise this model shows periodic bursting with fixed numbers of interspikes for most parameter values. As the parameter of activation kinetics of a slow potassium current is shifted to more hyperpolarized values of the membrane potential, the model undergoes a sequence of incremental spike adding transitions accumulating towards a periodic tonic spiking activity. Within a narrow parameter window around every spike adding transition, spike alteration of bursting is deterministically chaotic due to homoclinic bifurcations of a saddle periodic orbit. We have found that near these transitions the interneuron model becomes extremely sensitive to small random perturbations that cause a wide expansion and overlapping of the chaotic windows. The chaotic behavior is characterized by positive values of the largest Lyapunov exponent, and of the Shannon entropy of probability distribution of spike numbers per burst. The windows of chaotic dynamics resemble the Arnold tongues being plotted in the parameter plane, where the noise intensity serves as a second control parameter. We determine the critical noise intensities above which the interneuron model generates only irregular bursting within the overlapped windows.     

Exploring neuronal bistability at the depolarization block

Many neurons display bistability-coexistence of two firing modes such as bursting and tonic spiking or tonic spiking and silence. Bistability has been proposed to endow neurons with richer forms of information processing in general and to be involved in short-term memory in particular by allowing a brief signal to elicit long-lasting changes in firing. In this paper, we focus on bistability that allows for a choice between tonic spiking and depolarization block in a wide range of the depolarization levels. We consider the spike-producing currents in two neurons, models of which differ by the parameter values. Our dopaminergic neuron model displays bistability in a wide range of applied currents at the depolarization block. The Hodgkin-Huxley model of the squid giant axon shows no bistability. We varied parameter values for the model to analyze transitions between the two parameter sets. We show that bistability primarily characterizes the inactivation of the Na(+) current. Our study suggests a connection between the amount of the Na(+) window current and the length of the bistability range. For the dopaminergic neuron we hypothesize that bistability can be linked to a prolonged action of antipsychotic drugs.     

Firing patterns in the adaptive exponential integrate-and-fire model

For simulations of large spiking neuron networks, an accurate, simple and versatile single-neuron modeling framework is required. Here we explore the versatility of a simple two-equation model: the adaptive exponential integrate-and-fire neuron. We show that this model generates multiple firing patterns depending on the choice of parameter values, and present a phase diagram describing the transition from one firing type to another. We give an analytical criterion to distinguish between continuous adaption, initial bursting, regular bursting and two types of tonic spiking. Also, we report that the deterministic model is capable of producing irregular spiking when stimulated with constant current, indicating low-dimensional chaos. Lastly, the simple model is fitted to real experiments of cortical neurons under step current stimulation. The results provide support for the suitability of simple models such as the adaptive exponential integrate-and-fire neuron for large network simulations.     

Interactive responses of a thalamic neuron to formalin induced lasting pain in behaving mice

Thalamocortical (TC) neurons are known to relay incoming sensory information to the cortex via firing in tonic or burst mode. However, it is still unclear how respective firing modes of a single thalamic relay neuron contribute to pain perception under consciousness. Some studies report that bursting could increase pain in hyperalgesic conditions while others suggest the contrary. However, since previous studies were done under either neuropathic pain conditions or often under anesthesia, the mechanism of thalamic pain modulation under awake conditions is not well understood. We therefore characterized the thalamic firing patterns of behaving mice in response to nociceptive pain induced by inflammation. Our results demonstrated that nociceptive pain responses were positively correlated with tonic firing and negatively correlated with burst firing of individual TC neurons. Furthermore, burst properties such as intra-burst-interval (IntraBI) also turned out to be reliably correlated with the changes of nociceptive pain responses. In addition, brain stimulation experiments revealed that only bursts with specific bursting patterns could significantly abolish behavioral nociceptive responses. The results indicate that specific patterns of bursting activity in thalamocortical relay neurons play a critical role in controlling long-lasting inflammatory pain in awake and behaving mice.