Computational modeling of spike generation in serotonergic neurons of the dorsal raphe nucleu

We consider here a single-compartment model of these neurons which is capable of describing many of the known features of spike generation, particularly the slow rhythmic pacemaking activity often observed in these cells in a variety of species. Included in the model are ten kinds of voltage dependent ion channels as well as calcium-dependent potassium current. Calcium dynamics includes buffering and pumping. In sections 3-9, each component is considered in detail and parameters estimated from voltage clamp data where possible. In the next two sections simplified versions of some components are employed to explore the effects of various parameters on spiking, using a systematic approach, ending up with the following eleven components: a fast sodium current $I_{Na}$, a delayed rectifier potassium current $I_{KDR}$, a transient potassium current $I_A$, a low-threshold calcium current $I_T$, two high threshold calcium currents $I_L$ and $I_N$, small and large conductance potassium currents $I_{SK}$ and $I_{BK}$, a hyperpolarization-activated cation current $I_H$, a leak current $I_{Leak}$ and intracellular calcium ion concentration $Ca_i$. Attention is focused on the properties usually associated with these neurons, particularly long duration of action potential, pacemaker-like spiking and the ramp-like return to threshold after a spike. In some cases the membrane potential trajectories display doublets or have kinks or notches as have been reported in some experimental studies. The computed time courses of $I_A$ and $I_T$ during the interspike interval support the generally held view of a competition between them in influencing the frequency of spiking. Spontaneous spiking could be obtained with small changes in a few parameters from their values with driven spiking.Comment: The abstract has been truncate

Dendritic spike induction of postsynaptic cerebellar LTP

The architecture of parallel fiber (PF) axons contacting cerebellar Purkinje neurons (PNs) retains spatial information over long distances. PF synapses can trigger local dendritic calcium spikes, but whether and how this calcium signal leads to plastic changes that decode the PF input organization is unknown. By combining voltage and calcium imaging, we show that PF-elicited calcium signals, mediated by voltage-gated calcium channels, increase non-linearly during high-frequency bursts of electrically constant calcium spikes because they locally and transiently saturate the endogenous buffer. We demonstrate that these non-linear calcium signals, independently of NMDA or metabotropic glutamate receptor activation, can induce PF long-term potentiation (LTP). Two-photon imaging in coronal slices revealed that calcium signals inducing LTP can be observed by stimulating either the PF or the ascending fiber pathway. We propose that local dendritic calcium spikes, evoked by synaptic potentials, provide a unique mechanism to spatially decode PF signals into cerebellar circuitry changes

A universal model for spike-frequency adaptation

Spike-frequency adaptation is a prominent feature of neural dynamics. Among other mechanisms, various ionic currents modulating spike generation cause this type of neural adaptation. Prominent examples are voltage-gated potassium currents (M-type currents), the interplay of calcium currents and intracellular calcium dynamics with calcium-gated potassium channels (AHP-type currents), and the slow recovery from inactivation of the fast sodium current. While recent modeling studies have focused on the effects of specific adaptation currents, we derive a universal model for the firing-frequency dynamics of an adapting neuron that is independent of the specific adaptation process and spike generator. The model is completely defined by the neuron's onset f-I curve, the steady-state f-I curve, and the time constant of adaptation. For a specific neuron, these parameters can be easily determined from electrophysiological measurements without any pharmacological manipulations. At the same time, the simplicity of the model allows one to analyze mathematically how adaptation influences signal processing on the single-neuron level. In particular, we elucidate the specific nature of high-pass filter properties caused by spike-frequency adaptation. The model is limited to firing frequencies higher than the reciprocal adaptation time constant and to moderate fluctuations of the adaptation and the input current. As an extension of the model, we introduce a framework for combining an arbitrary spike generator with a generalized adaptation current

Modular analysis of the control of flagellar Ca2+-spike trains produced by CatSper and CaV channels in sea urchin sperm

Intracellular calcium ([Ca2+]i) is a basic and ubiquitous cellular signal controlling a wide variety of biological processes. A remarkable example is the steering of sea urchin spermatozoa towards the conspecific egg by a spatially and temporally orchestrated series of [Ca2+]i spikes. Although this process has been an experimental paradigm for reproduction and sperm chemotaxis studies, the composition and regulation of the signalling network underlying the cytosolic calcium fluctuations are hitherto not fully understood. Here, we used a differential equations model of the signalling network to assess which set of channels can explain the characteristic envelope and temporal organisation of the [Ca2+]i-spike trains. The signalling network comprises an initial membrane hyperpolarisation produced by an Upstream module triggered by the egg-released chemoattractant peptide, via receptor activation, cGMP synthesis and decay. Followed by downstream modules leading to intraflagellar pH (pHi), voltage and [Ca2+]i fluctuations. The Upstream module outputs were fitted to kinetic data on cGMP activity and early membrane potential changes measured in bulk cell populations. Two candidate modules featuring voltage-dependent Ca2+-channels link these outputs to the downstream dynamics and can independently explain the typical decaying envelope and the progressive spacing of the spikes. In the first module, [Ca2+]i-spike trains require the concerted action of a classical CaV-like channel and a potassium channel, BK (Slo1), whereas the second module relies on pHi-dependent CatSper dynamics articulated with voltage-dependent neutral sodium-proton exchanger (NHE). We analysed the dynamics of these two modules alone and in mixed scenarios. We show that the [Ca2+]i dynamics observed experimentally after sustained alkalinisation can be reproduced by a model featuring the CatSper and NHE module but not by those including the pH-independent CaV and BK module or proportionate mixed scenarios. We conclude in favour of the module containing CatSper and NHE and highlight experimentally testable predictions that would corroborate this conclusion

The low KM-phosphodiesterase inhibitor denbufylline enhances neuronal excitability in guinea pig hippocampus in vitro

The actions of the phosphodiesterase inhibitor denbufylline on the excitability of hippocampal neurons were investigated by means of extracellular and intracellular recordings. Denbufylline, which has been shown to selectively inhibit a low KM, Ca2+/calmodulin-independent phosphodiesterase isozyme, concentration-dependently increased the amplitude of the extracellularly recorded CAI population spike evoked by electrical stimulation of the Schaffer collateral/commissural pathway. Concentration-response-curves yielded an EC50 for denbufylline of 0.76 M. In comparison, the nonselective phosphodiesterase inhibitor 3-isobutyl-lmethylxanthine (IBMX) also produced an increase in the amplitude of the population spike. From the concentration-response-curve, which was steeper than that of denbufylline, an EC50 for IBMX of 1.04 M was obtained. However, despite their similar EC50 values, denbufylline was found to be significantly more potent at lower concentrations (<- 300 nM) than IBMX. Intracellular recordings from CAI pyramidal cells revealed postsynaptic actions of denbufylline (300 nM) as indicated by a small drug-induced depolarization (2 – 5 mV) associated with an increase in membrane input resistance by 10–20%. In addition, denbufylline blocked the accommodation of trains of action potentials evoked by the injection of depolarizing current pulses. The results suggest i) that accumulation of adenosine-3,5-monophosphate (CAMP) in the postsynaptic cell and/or in the presynaptic terminal produced by blockade of phosphodiesterases leads to enhanced synaptic transmission in the CAI area of the hippocampus and ii) that a low KM, Ca 2+/calmodulin-independent cAMP-phosphodiesterase is an important component involved in the regulation of the intracellular cAMP level at synapses of central nervous system neurons

Fast non-negative deconvolution for spike train inference from population calcium imaging

Calcium imaging for observing spiking activity from large populations of neurons are quickly gaining popularity. While the raw data are fluorescence movies, the underlying spike trains are of interest. This work presents a fast non-negative deconvolution filter to infer the approximately most likely spike train for each neuron, given the fluorescence observations. This algorithm outperforms optimal linear deconvolution (Wiener filtering) on both simulated and biological data. The performance gains come from restricting the inferred spike trains to be positive (using an interior-point method), unlike the Wiener filter. The algorithm is fast enough that even when imaging over 100 neurons, inference can be performed on the set of all observed traces faster than real-time. Performing optimal spatial filtering on the images further refines the estimates. Importantly, all the parameters required to perform the inference can be estimated using only the fluorescence data, obviating the need to perform joint electrophysiological and imaging calibration experiments.Comment: 22 pages, 10 figure

Slow synaptic transmission in frog sympathetic ganglia

Bullfrog ganglia contain two classes of neurone, B and C cells, which receive different inputs and exhibit different slow synaptic potentials. B cells, to which most effort has been directed, possess slow and late slow EPSPs. The sEPSP reflects a muscarinic action of acetylcholine released from boutons on B cells, whereas the late sEPSP is caused by a peptide (similar to teleost LHRH) released from boutons on C cells. During either sEPSP there is a selective reduction in two slow potassium conductances, designated 'M' and 'AHP'. The M conductance is voltage dependent and the AHP conductance is calcium dependent. Normally they act synergistically to prevent repetitive firing of action potentials during maintained stimuli. Computer stimulation of the interactions of these conductances with the other five voltage-dependent conductances present in the membrane allows a complete reconstruction of the effects of slow synaptic transmission on electrical behaviour

Spike frequency adaptation affects the synchronization properties of networks of cortical oscillators

Oscillations in many regions of the cortex have common temporal characteristics with dominant frequencies centered around the 40 Hz (gamma) frequency range and the 5–10 Hz (theta) frequency range. Experimental results also reveal spatially synchronous oscillations, which are stimulus dependent (Gray&Singer, 1987;Gray, König, Engel, & Singer, 1989; Engel, König, Kreiter, Schillen, & Singer, 1992). This rhythmic activity suggests that the coherence of neural populations is a crucial feature of cortical dynamics (Gray, 1994). Using both simulations and a theoretical coupled oscillator approach, we demonstrate that the spike frequency adaptation seen in many pyramidal cells plays a subtle but important role in the dynamics of cortical networks. Without adaptation, excitatory connections among model pyramidal cells are desynchronizing. However, the slow processes associated with adaptation encourage stable synchronous behavior