Rate coding and spike-time variability in cortical neurons with two types of threshold dynamics

Hierarchies of description enable understanding of cognitive phenomena in terms of neuron activity

One objective of neuroscience is to understand a wide range of specific cognitive processes in terms of neuron activity. The huge amount of observational data about the brain makes achieving this objective challenging. Different models on different levels of detail provide some insight, but the relationship between models on different levels is not clear. Complex computing systems with trillions of components like transistors are fully understood in the sense that system features can be precisely related to transistor activity. Such understanding could not involve a designer simultaneously thinking about the ongoing activity of all the components active in the course of carrying out some system feature. Brain modeling approaches like dynamical systems are inadequate to support understanding of computing systems, because their use relies on approximations like treating all components as more or less identical. Understanding computing systems needs a much more sophisticated use of approximation, involving creation of hierarchies of description in which the higher levels are more approximate, with effective translation between different levels in the hierarchy made possible by using the same general types of information processes on every level. These types are instruction and data read/write. There are no direct resemblances between computers and brains, but natural selection pressures have resulted in brain resources being organized into modular hierarchies and in the existence of two general types of information processes called condition definition/detection and behavioral recommendation. As a result, it is possible to create hierarchies of description linking cognitive phenomena to neuron activity, analogous with but qualitatively different from the hierarchies of description used to understand computing systems. An intuitively satisfying understanding of cognitive processes in terms of more detailed brain activity is then possible.

Bridging the Molecular-Cellular Gap in Understanding Ion Channel Clustering

The clustering of many voltage-dependent ion channel molecules at unique neuronal membrane sites such as axon initial segments, nodes of Ranvier, or the post-synaptic density, is an active process mediated by the interaction of ion channels with scaffold proteins and is of immense importance for electrical signaling. Growing evidence indicates that the density of ion channels at such membrane sites may affect action potential conduction properties and synaptic transmission. However, despite the emerging importance of ion channel density for electrical signaling, how ion channel-scaffold protein molecular interactions lead to cellular ion channel clustering, and how this process is regulated are largely unknown. In this review, we emphasize that voltage-dependent ion channel density at native clustering sites not only affects the density of ionic current fluxes but may also affect the conduction properties of the channel and/or the physical properties of the membrane at such locations, all changes that are expected to affect action potential conduction properties. Using the concrete example of the prototypical Shaker voltage-activated potassium channel (Kv) protein, we demonstrate how insight into the regulation of cellular ion channel clustering can be obtained when the molecular mechanism of ion channel-scaffold protein interaction is known. Our review emphasizes that such mechanistic knowledge is essential, and when combined with super-resolution imaging microscopy, can serve to bridge the molecular-cellular gap in understanding the regulation of ion channel clustering. Pressing questions, challenges and future directions in addressing ion channel clustering and its regulation are discussed.

Past and Future of Analog-Digital Modulation of Synaptic Transmission

Action potentials (APs) are generally produced in response to complex summation of excitatory and inhibitory synaptic inputs. While it is usually considered as a digital event, both the amplitude and width of the AP are significantly impacted by the context of its emission. In particular, the analog variations in subthreshold membrane potential determine the spike waveform and subsequently affect synaptic strength, leading to the so-called analog-digital modulation of synaptic transmission. We review here the numerous evidence suggesting context-dependent modulation of spike waveform, the discovery analog-digital modulation of synaptic transmission in invertebrates and its recent validation in mammals. We discuss the potential roles of analog-digital transmission in the physiology of neural networks.

An Adenosine A2A Receptor Antagonist Improves Multiple Symptoms of Repeated Quinpirole-Induced Psychosis

Obsessive-compulsive disorder (OCD) is a neuropsychiatric disorder characterized by the repeated rise of concerns (obsessions) and repetitive unwanted behavior (compulsions). Although selective serotonin reuptake inhibitors (SSRIs) is the first-choice drug, response rates to SSRI treatment vary between symptom dimensions. In this study, to find a therapeutic target for SSRI-resilient OCD symptoms, we evaluated treatment responses of quinpirole (QNP) sensitization-induced OCD-related behaviors in mice. SSRI administration rescued the cognitive inflexibility, as well as hyperactivity in the lateral orbitofrontal cortex (lOFC), while no improvement was observed for the repetitive behavior. D₂ receptor signaling in the central striatum (CS) was involved in SSRI-resistant repetitive behavior. An adenosine A_2A antagonist, istradefylline, which rescued abnormal excitatory synaptic function in the CS indirect pathway medium spiny neurons (MSNs) of sensitized mice, alleviated both of the QNP-induced abnormal behaviors with only short-term administration. These results provide a new insight into therapeutic strategies for SSRI-resistant OCD symptoms and indicate the potential of A_2A antagonists as a rapid-acting anti-OCD drug.

Dopamine Neurons Change the Type of Excitability in Response to Stimuli

The dynamics of neuronal excitability determine the neuron's response to stimuli, its synchronization and resonance properties and, ultimately, the computations it performs in the brain. We investigated the dynamical mechanisms underlying the excitability type of dopamine (DA) neurons, using a conductance-based biophysical model, and its regulation by intrinsic and synaptic currents. Calibrating the model to reproduce low frequency tonic firing results in N-methyl-D-aspartate (NMDA) excitation balanced by γ-Aminobutyric acid (GABA)-mediated inhibition and leads to type I excitable behavior characterized by a continuous decrease in firing frequency in response to hyperpolarizing currents. Furthermore, we analyzed how excitability type of the DA neuron model is influenced by changes in the intrinsic current composition. A subthreshold sodium current is necessary for a continuous frequency decrease during application of a negative current, and the low-frequency "balanced" state during simultaneous activation of NMDA and GABA receptors. Blocking this current switches the neuron to type II characterized by the abrupt onset of repetitive firing. Enhancing the anomalous rectifier Ih current also switches the excitability to type II. Key characteristics of synaptic conductances that may be observed in vivo also change the type of excitability: a depolarized γ-Aminobutyric acid receptor (GABAR) reversal potential or co-activation of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) leads to an abrupt frequency drop to zero, which is typical for type II excitability. Coactivation of N-methyl-D-aspartate receptors (NMDARs) together with AMPARs and GABARs shifts the type I/II boundary toward more hyperpolarized GABAR reversal potentials. To better understand how altering each of the aforementioned currents leads to changes in excitability profile of DA neuron, we provide a thorough dynamical analysis. Collectively, these results imply that type I excitability in dopamine neurons might be important for low firing rates and fine-tuning basal dopamine levels, while switching excitability to type II during NMDAR and AMPAR activation may facilitate a transient increase in dopamine concentration, as type II neurons are more amenable to synchronization by mutual excitation.

Coding Properties of Three Intrinsically Distinct Retinal Ganglion Cells under Periodic Stimuli: A Computational Study

As the sole output neurons in the retina, ganglion cells play significant roles in transforming visual information into spike trains, and then transmitting them to the higher visual centers. However, coding strategies that retinal ganglion cells (RGCs) adopt to accomplish these processes are not completely clear yet. To clarify these issues, we investigate the coding properties of three types of RGCs (repetitive spiking, tonic firing, and phasic firing) by two different measures (spike-rate and spike-latency). Model results show that for periodic stimuli, repetitive spiking RGC and tonic RGC exhibit similar spike-rate patterns. Their spike- rates decrease gradually with increased stimulus frequency, moreover, variation of stimulus amplitude would change the two RGCs' spike-rate patterns. For phasic RGC, it activates strongly at medium levels of frequency when the stimulus amplitude is low. While if high stimulus amplitude is applied, phasic RGC switches to respond strongly at low frequencies. These results suggest that stimulus amplitude is a prominent factor in regulating RGCs in encoding periodic signals. Similar conclusions can be drawn when analyzes spike-latency patterns of the three RGCs. More importantly, the above phenomena can be accurately reproduced by Hodgkin's three classes of neurons, indicating that RGCs can perform the typical three classes of firing dynamics, depending on the distinctions of ion channel densities. Consequently, model results from the three RGCs may be not specific, but can also applicable to neurons in other brain regions which exhibit part(s) or all of the Hodgkin's three excitabilities.

Modeling the Influence of Ion Channels on Neuron Dynamics in Drosophila

Voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of synaptic input patterns. Drosophila and other invertebrates provide valuable model systems for investigating ion channel kinetics and their impact on firing properties. Despite the increasing importance of Drosophila as a model system, few computational models of its ion channel kinetics have been developed. In this study, experimentally observed biophysical properties of voltage gated ion channels from the fruitfly Drosophila melanogaster are used to develop a minimal, conductance based neuron model. We investigate the impact of the densities of these channels on the excitability of the model neuron. Changing the channel densities reproduces different in situ observed firing patterns and induces a switch from integrator to resonator properties. Further, we analyze the preference to input frequency and how it depends on the channel densities and the resulting bifurcation type the system undergoes. An extension to a three dimensional model demonstrates that the inactivation kinetics of the sodium channels play an important role, allowing for firing patterns with a delayed first spike and subsequent high frequency firing as often observed in invertebrates, without altering the kinetics of the delayed rectifier current.

Dynamics of on-off neural firing patterns and stochastic effects near a sub-critical Hopf bifurcation

On-off firing patterns, in which repetition of clusters of spikes are interspersed with epochs of subthreshold oscillations or quiescent states, have been observed in various nervous systems, but the dynamics of this event remain unclear. Here, we report that on-off firing patterns observed in three experimental models (rat sciatic nerve subject to chronic constrictive injury, rat CA1 pyramidal neuron, and rabbit blood pressure baroreceptor) appeared as an alternation between quiescent state and burst containing multiple period-1 spikes over time. Burst and quiescent state had various durations. The interspike interval (ISI) series of on-off firing pattern was suggested as stochastic using nonlinear prediction and autocorrelation function. The resting state was changed to a period-1 firing pattern via on-off firing pattern as the potassium concentration, static pressure, or depolarization current was changed. During the changing process, the burst duration of on-off firing pattern increased and the duration of the quiescent state decreased. Bistability of a limit cycle corresponding to period-1 firing and a focus corresponding to resting state was simulated near a sub-critical Hopf bifurcation point in the deterministic Morris-Lecar (ML) model. In the stochastic ML model, noise-induced transitions between the coexisting regimes formed an on-off firing pattern, which closely matched that observed in the experiment. In addition, noise-induced exponential change in the escape rate from the focus, and noise-induced coherence resonance were identified. The distinctions between the on-off firing pattern and stochastic firing patterns generated near three other types of bifurcations of equilibrium points, as well as other viewpoints on the dynamics of on-off firing pattern, are discussed. The results not only identify the on-off firing pattern as noise-induced stochastic firing pattern near a sub-critical Hopf bifurcation point, but also offer practical indicators to discriminate bifurcation types and neural excitability types.

Changes of Firing Rate Induced by Changes of Phase Response Curve in Bifurcation Transitions

We study dynamical mechanisms responsible for changes of the firing rate during four different bifurcation transitions in the two-dimensional Hindmarsh-Rose (2DHR) neuron model: the saddle node on an invariant circle (SNIC) bifurcation to the supercritical Andronov-Hopf (AH) one, the SNIC bifurcation to the saddle-separatrix loop (SSL) one, the AH bifurcation to the subcritical AH (SAH) one, and the SSL bifurcation to the AH one. For this purpose, we study slopes of the firing rate curve with respect to not only an external input current but also temperature that can be interpreted as a timescale in the 2DHR neuron model. These slopes are mathematically formulated with phase response curves (PRCs), expanding the firing rate with perturbations of the temperature and external input current on the one-dimensional space of the phase [Formula: see text] in the 2DHR oscillator. By analyzing the two different slopes of the firing rate curve with respect to the temperature and external input current, we find that during changes of the firing rate in all of the bifurcation transitions, the calculated slope with respect to the temperature also changes. This is largely dependent on changes in the PRC size that is also related to the slope with respect to the external input current. Furthermore, we find phase transition–like switches of the firing rate with a possible increase of the temperature during the SSL-to-AH bifurcation transition.

Density of voltage-gated potassium channels is a bifurcation parameter in pyramidal neurons

Several types of intrinsic dynamics have been identified in brain neurons. Type 1 excitability is characterized by a continuous frequency-stimulus relationship and, thus, an arbitrarily low frequency at threshold current. Conversely, Type 2 excitability is characterized by a discontinuous frequency-stimulus relationship and a nonzero threshold frequency. In previous theoretical work we showed that the density of Kv channels is a bifurcation parameter, such that increasing the Kv channel density in a neuron model transforms Type 1 excitability into Type 2 excitability. Here we test this finding experimentally, using the dynamic clamp technique on Type 1 pyramidal cells in rat cortex. We found that increasing the density of slow Kv channels leads to a shift from Type 1 to Type 2 threshold dynamics, i.e., a distinct onset frequency, subthreshold oscillations, and reduced latency to first spike. In addition, the action potential was resculptured, with a narrower spike width and more pronounced afterhyperpolarization. All changes could be captured with a two-dimensional model. It may seem paradoxical that an increase in slow K channel density can lead to a higher threshold firing frequency; however, this can be explained in terms of bifurcation theory. In contrast to previous work, we argue that an increased outward current leads to a change in dynamics in these neurons without a rectification of the current-voltage curve. These results demonstrate that the behavior of neurons is determined by the global interactions of their dynamical elements and not necessarily simply by individual types of ion channels.

A CMOS IC-based multisite measuring system for stimulation and recording in neural preparations in vitro

In this report, we describe the system integration of a complementary metal oxide semiconductor (CMOS) integrated circuit (IC) chip, capable of both stimulation and recording of neurons or neural tissues, to investigate electrical signal propagation within cellular networks in vitro. The overall system consisted of three major subunits: a 5.0 × 5.0 mm CMOS IC chip, a reconfigurable logic device (field-programmable gate array, FPGA), and a PC. To test the system, microelectrode arrays (MEAs) were used to extracellularly measure the activity of cultured rat cortical neurons and mouse cortical slices. The MEA had 64 bidirectional (stimulation and recording) electrodes. In addition, the CMOS IC chip was equipped with dedicated analog filters, amplification stages, and a stimulation buffer. Signals from the electrodes were sampled at 15.6 kHz with 16-bit resolution. The measured input-referred circuitry noise was 10.1 μ V root mean square (10 Hz to 100 kHz), which allowed reliable detection of neural signals ranging from several millivolts down to approximately 33 μ Vpp. Experiments were performed involving the stimulation of neurons with several spatiotemporal patterns and the recording of the triggered activity. An advantage over current MEAs, as demonstrated by our experiments, includes the ability to stimulate (voltage stimulation, 5-bit resolution) spatiotemporal patterns in arbitrary subsets of electrodes. Furthermore, the fast stimulation reset mechanism allowed us to record neuronal signals from a stimulating electrode around 3 ms after stimulation. We demonstrate that the system can be directly applied to, for example, auditory neural prostheses in conjunction with an acoustic sensor and a sound processing system.

Afferent inputs to cortical fast-spiking interneurons organize pyramidal cell network oscillations at high-gamma frequencies (60-200 Hz)

High-gamma activity, ranging in frequency between ∼60 Hz and 200 Hz, has been observed in local field potential, electrocorticography, EEG and magnetoencephalography signals during cortical activation, in a variety of functional brain systems. The origin of these signals is yet unknown. Using computational modeling, we show that a cortical network model receiving thalamic input generates high-gamma responses comparable to those observed in local field potential recorded in monkey somatosensory cortex during vibrotactile stimulation. These high-gamma oscillations appear to be mediated mostly by an excited population of inhibitory fast-spiking interneurons firing at high-gamma frequencies and pacing excitatory regular-spiking pyramidal cells, which fire at lower rates but in phase with the population rhythm. The physiological correlates of high-gamma activity, in this model of local cortical circuits, appear to be similar to those proposed for hippocampal ripples generated by subsets of interneurons that regulate the discharge of principal cells.

The past, present, and future of real-time control in cellular electrophysiology

For over 60 years, real-time control has been an important technique in the study of excitable cells. Two such control-based technologies are reviewed here. First, voltage-clamp methods revolutionized the study of excitable cells. In this family of techniques, membrane potential is controlled, allowing one to parameterize a powerful class of models that describe the voltage-current relationship of cell membranes simply, flexibly, and accurately. Second, dynamic-clamp methods allow the addition of new, "virtual" membrane mechanisms to living cells. Dynamic clamp allows researchers unprecedented ways of testing computationally based hypotheses in biological preparations. The review ends with predictions of how control-based technologies will be improved and adapted for new uses in the near future.

Seizure-induced alterations in fast-spiking basket cell GABA currents modulate frequency and coherence of gamma oscillation in network simulations

Gamma frequency oscillations have been proposed to contribute to memory formation and retrieval. Fast-spiking basket cells (FS-BCs) are known to underlie development of gamma oscillations. Fast, high amplitude GABA synapses and gap junctions have been suggested to contribute to gamma oscillations in FS-BC networks. Recently, we identified that, apart from GABAergic synapses, FS-BCs in the hippocampal dentate gyrus have GABAergic currents mediated by extrasynaptic receptors. Our experimental studies demonstrated two specific changes in FS-BC GABA currents following experimental seizures [Yu et al., J. Neurophysiol. 109, 1746 (2013)]: increase in the magnitude of extrasynaptic (tonic) GABA currents and a depolarizing shift in GABA reversal potential (E(GABA)). Here, we use homogeneous networks of a biophysically based model of FS-BCs to examine how the presence of extrasynaptic GABA conductance (g(GABA-extra)) and experimentally identified, seizure-induced changes in g(GABA-extra) and E(GABA) influence network activity. Networks of FS-BCs interconnected by fast GABAergic synapses developed synchronous firing in the dentate gamma frequency range (40-100 Hz). Systematic investigation revealed that the biologically realistic range of 30 to 40 connections between FS-BCs resulted in greater coherence in the gamma frequency range when networks were activated by Poisson-distributed dendritic synaptic inputs rather than by homogeneous somatic current injections, which were balanced for FS-BC firing frequency in unconnected networks. Distance-dependent conduction delay enhanced coherence in networks with 30-40 FS-BC interconnections while inclusion of gap junctional conductance had a modest effect on coherence. In networks activated by somatic current injections resulting in heterogeneous FS-BC firing, increasing g(GABA-extra) reduced the frequency and coherence of FS-BC firing when E(GABA) was shunting (-74 mV), but failed to alter average FS-BC frequency when E(GABA) was depolarizing (-54 mV). When FS-BCs were activated by biologically based dendritic synaptic inputs, enhancing g(GABA-extra) reduced the frequency and coherence of FS-BC firing when E(GABA) was shunting and increased average FS-BC firing when E(GABA) was depolarizing. Shifting E(GABA) from shunting to depolarizing potentials consistently increased network frequency to and above high gamma frequencies (>80 Hz). Since gamma oscillations may contribute to learning and memory processing [Fell et al., Nat. Neurosci. 4, 1259 (2001); Jutras et al., J. Neurosci. 29, 12521 (2009); Wang, Physiol. Rev. 90, 1195 (2010)], our demonstration that network oscillations are modulated by extrasynaptic inhibition in FS-BCs suggests that neuroactive compounds that act on extrasynaptic GABA receptors could impact memory formation by modulating hippocampal gamma oscillations. The simulation results indicate that the depolarized FS-BC GABA reversal, observed after experimental seizures, together with enhanced spillover extrasynaptic GABA currents are likely to promote generation of focal high frequency activity associated with epileptic networks.

A computational study of spike time reliability in two types of threshold dynamics

Spike time reliability (STR) refers to the phenomenon in which repetitive applications of a frozen copy of one stochastic signal to a neuron trigger spikes with reliable timing while a constant signal fails to do so. Observed and explored in numerous experimental and theoretical studies, STR is a complex dynamic phenomenon depending on the nature of external inputs as well as intrinsic properties of a neuron. The neuron under consideration could be either quiescent or spontaneously spiking in the absence of the external stimulus. Focusing on the situation in which the unstimulated neuron is quiescent but close to a switching point to oscillations, we numerically analyze STR treating each spike occurrence as a time localized event in a model neuron. We study both the averaged properties as well as individual features of spike-evoking epochs (SEEs). The effects of interactions between spikes is minimized by selecting signals that generate spikes with relatively long interspike intervals (ISIs). Under these conditions, the frequency content of the input signal has little impact on STR. We study two distinct cases, Type I in which the f-I relation (f for frequency, I for applied current) is continuous and Type II where the f-I relation exhibits a jump. STR in the two types shows a number of similar features and differ in some others. SEEs that are capable of triggering spikes show great variety in amplitude and time profile. On average, reliable spike timing is associated with an accelerated increase in the "action" of the signal as a threshold for spike generation is approached. Here, "action" is defined as the average amount of current delivered during a fixed time interval. When individual SEEs are studied, however, their time profiles are found important for triggering more precisely timed spikes. The SEEs that have a more favorable time profile are capable of triggering spikes with higher precision even at lower action levels.

Toggling between gamma-frequency activity and suppression of cell assemblies

Gamma (30–80 Hz) rhythms in hippocampus and neocortex resulting from the interaction of excitatory and inhibitory cells (E- and I-cells), called Pyramidal-Interneuronal Network Gamma (PING), require that the I-cells respond to the E-cells, but don't fire on their own. In idealized models, there is a sharp boundary between a parameter regime where the I-cells have weak-enough drive for PING, and one where they have so much drive that they fire without being prompted by the E-cells. In the latter regime, they often de-synchronize and suppress the E-cells; the boundary was therefore called the “suppression boundary” by Börgers and Kopell (2005). The model I-cells used in the earlier work by Börgers and Kopell have a “type 1” phase response, i.e., excitatory input always advances them. However, fast-spiking inhibitory basket cells often have a “type 2” phase response: Excitatory input arriving soon after they fire delays them. We study the effect of the phase response type on the suppression transition, under the additional assumption that the I-cells are kept synchronous by gap junctions. When many E-cells participate on a given cycle, the resulting excitation advances the I-cells on the next cycle if their phase response is of type 1, and this can result in suppression of more E-cells on the next cycle. Therefore, strong E-cell spike volleys tend to be followed by weaker ones, and vice versa. This often results in erratic fluctuations in the strengths of the E-cell spike volleys. When the phase response of the I-cells is of type 2, the opposite happens: strong E-cell spike volleys delay the inhibition on the next cycle, therefore tend to be followed by yet stronger ones. The strengths of the E-cell spike volleys don't oscillate, and there is a nearly abrupt transition from PING to ING (a rhythm involving I-cells only).

Moving beyond Type I and Type II neuron types

In 1948, Hodgkin delineated different classes of axonal firing.  This has been mathematically translated allowing insight and understanding to emerge.  As such, the terminology of ‘Type I’ and ‘Type II’ neurons is commonplace in the Neuroscience literature today.  Theoretical insights have helped us realize that, for example, network synchronization depends on whether neurons are Type I or Type II.  Mathematical models are precise with analyses (considering Type I/II aspects), but experimentally, the distinction can be less clear.  On the other hand, experiments are becoming more sophisticated in terms of distinguishing and manipulating particular cell types but are limited in terms of being able to consider network aspects simultaneously.   Although there is much work going on mathematically and experimentally, in my opinion it is becoming common that models are either superficially linked with experiment or not described in enough detail to appreciate the biological context.  Overall, we all suffer in terms of impeding our understanding of brain networks and applying our understanding to neurological disease.  I suggest that more modelers become familiar with experimental details and that more experimentalists appreciate modeling assumptions. In other words, we need to move beyond our comfort zones.

Identifying type I excitability using dynamics of stochastic neural firing patterns

The stochastic firing patterns are simulated near saddle-node bifurcation on an invariant cycle corresponding to type I excitability in stochastic Morris-Lecar model. In absence of external periodic signal, the stochastic firing manifests continuous distribution in ISI histogram (ISIH), whose amplitude at first increases sharply and then decreases exponentially. In presence of the external periodic signal, stochastic firing patterns appear as two cases of integer multiple firing with multiple discrete peaks in ISIH. One manifests perfect exponential decay in all peaks and the other imperfect exponential decay except a lower first peak. These stochastic firing patterns simulated with or without external periodic signal can be demonstrated in the experiments on rat hippocampal CA1 pyramidal neurons. The exponential decay laws in the multiple peaks are also acquired using probability analysis method. The perfect decay law is determined by the independent characteristic within the firing while the imperfect decay law is from the inhibitory effect. In addition, the stochastic firing patterns corresponding to type I excitability are compared to those of type II excitability. The results not only reveal the dynamics of stochastic firing patterns with or without external signal corresponding to type I excitability, but also provide practical indicators to availably identify type I excitability.

Inferring single neuron properties in conductance based balanced networks

Balanced states in large networks are a usual hypothesis for explaining the variability of neural activity in cortical systems. In this regime the statistics of the inputs is characterized by static and dynamic fluctuations. The dynamic fluctuations have a Gaussian distribution. Such statistics allows to use reverse correlation methods, by recording synaptic inputs and the spike trains of ongoing spontaneous activity without any additional input. By using this method, properties of the single neuron dynamics that are masked by the balanced state can be quantified. To show the feasibility of this approach we apply it to large networks of conductance based neurons. The networks are classified as Type I or Type II according to the bifurcations which neurons of the different populations undergo near the firing onset. We also analyze mixed networks, in which each population has a mixture of different neuronal types. We determine under which conditions the intrinsic noise generated by the network can be used to apply reverse correlation methods. We find that under realistic conditions we can ascertain with low error the types of neurons present in the network. We also find that data from neurons with similar firing rates can be combined to perform covariance analysis. We compare the results of these methods (that do not requite any external input) to the standard procedure (that requires the injection of Gaussian noise into a single neuron). We find a good agreement between the two procedures.

The mechanism of ethanol action on midbrain dopaminergic neuron firing: a dynamic-clamp study of the role of I(h) and GABAergic synaptic integration

Hyperpolarization-activated and cyclic nucleotide-gated (HCN) channels are expressed in dopaminergic (DA) neurons of the ventral tegmental area (VTA) as well as in DA and GABAergic neurons of the substantia nigra (SN). The excitation of DA neurons induced by ethanol has been proposed to result from its enhancing HCN channel current, I(h). Using perforated patch-clamp recordings in rat midbrain slices, we isolated I(h) in these neurons by voltage clamp. We showed that ethanol reversibly increased the amplitude and accelerated the activation kinetics of I(h) and caused a depolarizing shift in its voltage dependence. Using dynamic-clamp conductance injection, we injected artificial I(h) and fluctuating GABAergic synaptic conductance inputs into neurons following block of intrinsic I(h). This demonstrated directly a major role of I(h) in promoting rebound spiking following phasic inhibition, which was enhanced as the kinetics and amplitude of I(h) were changed in the manner induced by ethanol. Similar effects of ethanol were observed on I(h) and firing rate in non-DA, putatively GABAergic interneurons, indicating that in addition to its direct effects on firing, ethanol will produce large changes in the inhibition and disinhibition (via GABAergic interneurons) converging on DA neurons. Thus the overall effects of ethanol on firing of DA cells of the VTA and SN in vivo, and hence on phasic dopamine release in the striatum, appear to be determined substantially by its action on I(h) in both DA cells and GABAergic interneurons.

Cellularly-driven differences in network synchronization propensity are differentially modulated by firing frequency

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.

Synchronization of the firing of neurons in the brain is related to many cognitive functions, such as recognizing faces, discriminating odors, and coordinating movement. It is therefore important to understand what properties of neuronal networks promote synchrony of neural firing. One measure that is often used to determine the contribution of individual neurons to network synchrony is called the phase response curve (PRC). PRCs describe how the timing of neuronal firing changes depending on when input, such as a synaptic signal, is received by the neuron. A characteristic of PRCs that has previously not been well understood is that they change dramatically as the neuron's firing frequency is modulated. This effect carries potential significance, since cognitive functions are often associated with specific frequencies of network activity in the brain. We showed computationally that the frequency dependence of PRCs can be explained by the relative timing of ionic membrane currents with respect to the time between spike firings. Our simulations also showed that the frequency dependence of neuronal PRCs leads to frequency-dependent changes in network synchronization that can be different for different neuron types. These results further our understanding of how synchronization is generated in the brain to support various cognitive functions.

Spike-timing-dependent plasticity and reliability optimization: the role of neuron dynamics

Plastic changes in synaptic efficacy can depend on the time ordering of presynaptic and postsynaptic spikes. This phenomenon is called spike-timing-dependent plasticity (STDP). One of the most striking aspects of this plasticity mechanism is that the STDP windows display a great variety of forms in different parts of the nervous system. We explore this issue from a theoretical point of view. We choose as the optimization principle the minimization of conditional entropy or maximization of reliability in the transmission of information. We apply this principle to two types of postsynaptic dynamics, designated type I and type II. The first is characterized as being an integrator, while the second is a resonator. We find that, depending on the parameters of the models, the optimization principle can give rise to a wide variety of STDP windows, such as antisymmetric Hebbian, predominantly depressing or symmetric with one positive region and two lateral negative regions. We can relate each of these forms to the dynamical behavior of the different models. We also propose experimental tests to assess the validity of the optimization principle.

Rate maintenance and resonance in the entorhinal cortex

Throughout the brain, neurons encode information in fundamental units of spikes. Each spike represents the combined thresholding of synaptic inputs and intrinsic neuronal dynamics. Here, we address a basic question of spike train formation: how do perithreshold synaptic inputs perturb the output of a spiking neuron? We recorded from single entorhinal principal cells in vitro and drove them to spike steadily at ∼5 Hz (theta range) with direct current injection, then used a dynamic-clamp to superimpose strong excitatory conductance inputs at varying rates. Neurons spiked most reliably when the input rate matched the intrinsic neuronal firing rate. We also found a striking tendency of neurons to preserve their rates and coefficients of variation, independently of input rates. As mechanisms for this rate maintenance, we show that the efficacy of the conductance inputs varied with the relationship of input rate to neuronal firing rate, and with the arrival time of the input within the natural period. Using a novel method of spike classification, we developed a minimal Markov model that reproduced the measured statistics of the output spike trains and thus allowed us to identify and compare contributions to the rate maintenance and resonance. We suggest that the strength of rate maintenance may be used as a new categorization scheme for neuronal response and note that individual intrinsic spiking mechanisms may play a significant role in forming the rhythmic spike trains of activated neurons; in the entorhinal cortex, individual pacemakers may dominate production of the regional theta rhythm.

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.

All activity of the brain is manifested in electrical oscillatory patterns, shaped by the firing dynamics of the many neurons forming the brain networks. The underlying mechanisms of the firing pattern in the single neurons are still not fully understood. The distribution and identity of different channel types have been suggested as critical factors. We have suggested that the density of channels in the membrane is a fundamental complementary mechanism. In a hippocampal soma membrane model study we have shown that altering the ion channel densities can cause the membrane to switch between two qualitatively different firing patterns. Here we extend the analysis to two axon membranes. Unexpectedly, both show that channel density alterations do not cause switches between different firing behaviours. We believe that this is an important property of axon membranes, explaining their limited flexibility.

Diversity of intrinsic frequency encoding patterns in rat cortical neurons--mechanisms and possible functions

Extracellular recordings of single neurons in primary and secondary somatosensory cortices of monkeys in vivo have shown that their firing rate can increase, decrease, or remain constant in different cells, as the external stimulus frequency increases. We observed similar intrinsic firing patterns (increasing, decreasing or constant) in rat somatosensory cortex in vitro, when stimulated with oscillatory input using conductance injection (dynamic clamp). The underlying mechanism of this observation is not obvious, and presents a challenge for mathematical modelling. We propose a simple principle for describing this phenomenon using a leaky integrate-and-fire model with sinusoidal input, an intrinsic oscillation and Poisson noise. Additional enhancement of the gain of encoding could be achieved by local network connections amongst diverse intrinsic response patterns. Our work sheds light on the possible cellular and network mechanisms underlying these opposing neuronal responses, which serve to enhance signal detection.

Inhibition of adult rat retinal ganglion cells by D1-type dopamine receptor activation

The spike output of neural pathways can be regulated by modulating output neuron excitability and/or their synaptic inputs. Dopaminergic interneurons synapse onto cells that route signals to mammalian retinal ganglion cells, but it is unknown whether dopamine can activate receptors in these ganglion cells and, if it does, how this affects their excitability. Here, we show D(1a) receptor-like immunoreactivity in ganglion cells identified in adult rats by retrogradely transported dextran, and that dopamine, D(1)-type receptor agonists, and cAMP analogs inhibit spiking in ganglion cells dissociated from adult rats. These ligands curtailed repetitive spiking during constant current injections and reduced the number and rate of rise of spikes elicited by fluctuating current injections without significantly altering the timing of the remaining spikes. Consistent with mediation by D(1)-type receptors, SCH-23390 [R-(+)-7-chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine] reversed the effects of dopamine on spikes. Contrary to a recent report, spike inhibition by dopamine was not precluded by blocking I(h). Consistent with the reduced rate of spike rise, dopamine reduced voltage-gated Na(+) current (I(Na)) amplitude, and tetrodotoxin, at doses that reduced I(Na) as moderately as dopamine, also inhibited spiking. These results provide the first direct evidence that D(1)-type dopamine receptor activation can alter mammalian retinal ganglion cell excitability and demonstrate that dopamine can modulate spikes in these cells by a mechanism different from the presynaptic and postsynaptic means proposed by previous studies. To our knowledge, our results also provide the first evidence that dopamine receptor activation can reduce excitability without altering the temporal precision of spike firing.

Generalization of the event-based Carnevale-Hines integration scheme for integrate-and-fire models

An event-based integration scheme for an integrate-and-fire neuron model with exponentially decaying excitatory synaptic currents and double exponential inhibitory synaptic currents has been introduced by Carnevale and Hines. However, the integration scheme imposes nonphysiological constraints on the time constants of the synaptic currents, which hamper its general applicability. This letter addresses this problem in two ways. First, we provide physical arguments demonstrating why these constraints on the time constants can be relaxed. Second, we give a formal proof showing which constraints can be abolished. As part of our formal proof, we introduce the generalized Carnevale-Hines lemma, a new tool for comparing double exponentials as they naturally occur in many cascaded decay systems, including receptor-neurotransmitter dissociation followed by channel closing. Through repeated application of the generalized lemma, we lift most of the original constraints on the time constants. Thus, we show that the Carnevale-Hines integration scheme for the integrate-and-fire model can be employed for simulating a much wider range of neuron and synapse types than was previously thought.

Causes of transient instabilities in the dynamic clamp

The dynamic clamp is a widely used method for integrating mathematical models with electrophysiological experiments. This method involves measuring the membrane voltage of a cell, using it to solve computational models of ion channel dynamics in real-time, and injecting the calculated current(s) back into the cell. Limitations of this technique include those associated with single electrode current clamping and the sampling effects caused by the dynamic clamp. In this study, we show that the combination of these limitations causes transient instabilities under certain conditions. Through physical experiments and simulations, we show that dynamic clamp instability is directly related to the sampling delay and the maximum simulated conductance being injected. It is exaggerated by insufficient electrode series resistance and capacitance compensation. Increasing the sampling rate of the dynamic clamp system increases dynamic clamp stability; however, this improvement, is constrained by how well the electrode series resistance and capacitance are compensated. At present, dynamic clamp sampling rates are justified solely on the temporal dynamics of the models being simulated; here we show that faster rates increase the stable range of operation for the dynamic clamp system. In addition, we show that commonly accepted levels of resistance compensation nevertheless significantly compromise the stability of a dynamic clamp system.

Integration of broadband conductance input in rat somatosensory cortical inhibitory interneurons: an inhibition-controlled switch between intrinsic and input-driven spiking in fast-spiking cells

Quantitative understanding of the dynamics of particular cell types when responding to complex, natural inputs is an important prerequisite for understanding the operation of the cortical network. Different types of inhibitory neurons are connected by electrical synapses to nearby neurons of the same type, enabling the formation of synchronized assemblies of neurons with distinct dynamical behaviors. Under what conditions is spike timing in such cells determined by their intrinsic dynamics and when is it driven by the timing of external input? In this study, we have addressed this question using a systematic approach to characterizing the input-output relationships of three types of cortical interneurons (fast spiking [FS], low-threshold spiking [LTS], and nonpyramidal regular-spiking [NPRS] cells) in the rat somatosensory cortex, during fluctuating conductance input designed to mimic natural complex activity. We measured the shape of average conductance input trajectories preceding spikes and fitted a two-component linear model of neuronal responses, which included an autoregressive term from its own output, to gain insight into the input-output relationships of neurons. This clearly separated the contributions of stimulus and discharge history, in a cell-type dependent manner. Unlike LTS and NPRS cells, FS cells showed a remarkable switch in dynamics, from intrinsically driven spike timing to input-fluctuation-controlled spike timing, with the addition of even a small amount of inhibitory conductance. Such a switch could play a pivotal role in the function of FS cells in organizing coherent gamma oscillations in the local cortical network. Using both pharmacological perturbations and modeling, we show how this property is a consequence of the particular complement of voltage-dependent conductances in these cells.

Biophysical basis for three distinct dynamical mechanisms of action potential initiation

Transduction of graded synaptic input into trains of all-or-none action potentials (spikes) is a crucial step in neural coding. Hodgkin identified three classes of neurons with qualitatively different analog-to-digital transduction properties. Despite widespread use of this classification scheme, a generalizable explanation of its biophysical basis has not been described. We recorded from spinal sensory neurons representing each class and reproduced their transduction properties in a minimal model. With phase plane and bifurcation analysis, each class of excitability was shown to derive from distinct spike initiating dynamics. Excitability could be converted between all three classes by varying single parameters; moreover, several parameters, when varied one at a time, had functionally equivalent effects on excitability. From this, we conclude that the spike-initiating dynamics associated with each of Hodgkin's classes represent different outcomes in a nonlinear competition between oppositely directed, kinetically mismatched currents. Class 1 excitability occurs through a saddle node on invariant circle bifurcation when net current at perithreshold potentials is inward (depolarizing) at steady state. Class 2 excitability occurs through a Hopf bifurcation when, despite net current being outward (hyperpolarizing) at steady state, spike initiation occurs because inward current activates faster than outward current. Class 3 excitability occurs through a quasi-separatrix crossing when fast-activating inward current overpowers slow-activating outward current during a stimulus transient, although slow-activating outward current dominates during constant stimulation. Experiments confirmed that different classes of spinal lamina I neurons express the subthreshold currents predicted by our simulations and, further, that those currents are necessary for the excitability in each cell class. Thus, our results demonstrate that all three classes of excitability arise from a continuum in the direction and magnitude of subthreshold currents. Through detailed analysis of the spike-initiating process, we have explained a fundamental link between biophysical properties and qualitative differences in how neurons encode sensory input.

Information is transmitted through the nervous system in the form of action potentials or spikes. Contrary to popular belief, a spike is not generated instantaneously when membrane potential crosses some preordained threshold. In fact, different neurons employ different rules to determine when and why they spike. These different rules translate into diverse spiking patterns that have been observed experimentally and replicated time and again in computational models. In this study, our aim was not simply to replicate different spiking patterns; instead, we sought to provide deeper insight into the connection between biophysics and neural coding by relating each to the process of spike initiation. We show that Hodgkin's three classes of excitability result from a nonlinear competition between oppositely directed, kinetically mismatched currents; the outcome of that competition is manifested as dynamically distinct spike-initiating mechanisms. Our results highlight the benefits of forward engineering minimal models capable of reproducing phenomena of interest and then dissecting those models in order to identify general explanations of how those phenomena arise. Furthermore, understanding nonlinear dynamical processes such as spike initiation is crucial for definitively explaining how biophysical properties impact neural coding.

Cellular mechanisms underlying stimulus-dependent gain modulation in primary visual cortex neurons in vivo

Gain modulation is a widespread neuronal phenomenon that modifies response amplitude without changing selectivity. Computational and in vitro studies have proposed cellular mechanisms of gain modulation based on the postsynaptic effects of background synaptic activation, but these mechanisms have not been studied in vivo. Here, we used intracellular recordings from cat primary visual cortex to measure neuronal gain while changing background synaptic activity with visual stimulation. We found that increases in the membrane fluctuations associated with increases in synaptic input do not obligatorily result in gain modulation in vivo. However, visual stimuli that evoked sustained changes in resting membrane potential, input resistance, and membrane fluctuations robustly modulated neuronal gain. The magnitude of gain modulation depended critically on the spatiotemporal properties of the visual stimulus. Gain modulation in vivo may thus be determined on a moment-to-moment basis by sensory context and the consequent dynamics of synaptic activation.

Thalamocortical relay fidelity varies across subthalamic nucleus deep brain stimulation protocols in a data-driven computational model

The therapeutic effectiveness of deep brain stimulation (DBS) of the subthalamic nucleus (STN) may arise through its effects on inhibitory basal ganglia outputs, including those from the internal segment of the globus pallidus (GPi). Changes in GPi activity will impact its thalamic targets, representing a possible pathway for STN-DBS to modulate basal ganglia-thalamocortical processing. To study the effect of STN-DBS on thalamic activity, we examined thalamocortical (TC) relay cell responses to an excitatory input train under a variety of inhibitory signals, using a computational model. The inhibitory signals were obtained from single-unit GPi recordings from normal monkeys and from monkeys rendered parkinsonian through arterial 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine injection. The parkinsonian GPi data were collected in the absence of STN-DBS, under sub-therapeutic STN-DBS, and under therapeutic STN-DBS. Our simulations show that inhibition from parkinsonian GPi activity recorded without DBS-compromised TC relay of excitatory inputs compared with the normal case, whereas TC relay fidelity improved significantly under inhibition from therapeutic, but not sub-therapeutic, STN-DBS GPi activity. In a heterogeneous model TC cell population, response failures to the same input occurred across multiple TC cells significantly more often without DBS than in the therapeutic DBS case and in the normal case. Inhibitory signals preceding successful TC relay were relatively constant, whereas those before failures changed more rapidly. Computationally generated inhibitory inputs yielded similar effects on TC relay. These results support the hypothesis that STN-DBS alters parkinsonian GPi activity in a way that may improve TC relay fidelity.

Characterizing synaptic conductance fluctuations in cortical neurons and their influence on spike generation

Cortical neurons are subject to sustained and irregular synaptic activity which causes important fluctuations of the membrane potential (V(m)). We review here different methods to characterize this activity and its impact on spike generation. The simplified, fluctuating point-conductance model of synaptic activity provides the starting point of a variety of methods for the analysis of intracellular V(m) recordings. In this model, the synaptic excitatory and inhibitory conductances are described by Gaussian-distributed stochastic variables, or "colored conductance noise". The matching of experimentally recorded V(m) distributions to an invertible theoretical expression derived from the model allows the extraction of parameters characterizing the synaptic conductance distributions. This analysis can be complemented by the matching of experimental V(m) power spectral densities (PSDs) to a theoretical template, even though the unexpected scaling properties of experimental PSDs limit the precision of this latter approach. Building on this stochastic characterization of synaptic activity, we also propose methods to qualitatively and quantitatively evaluate spike-triggered averages of synaptic time-courses preceding spikes. This analysis points to an essential role for synaptic conductance variance in determining spike times. The presented methods are evaluated using controlled conductance injection in cortical neurons in vitro with the dynamic-clamp technique. We review their applications to the analysis of in vivo intracellular recordings in cat association cortex, which suggest a predominant role for inhibition in determining both sub- and supra-threshold dynamics of cortical neurons embedded in active networks.

A scriptable DSP-based system for dynamic conductance injection

A variety of software and hardware systems have been developed to inject controlled electrical conductances into excitable cells, to investigate the physiological mechanisms of action potential generation. These systems face several challenges: the need to model complex conductances, including voltage-gated ion channels, synaptic conductances controlled by electrical models of entire cells or even networks of cells, to do so rapidly and stably, with precisely controlled update intervals of 20micros or less, and to present an easy and flexible interface to the user, allowing new experiments to be designed and executed easily. In this paper I describe a new software system (SM-2) which is designed to meet these requirements, and which runs on the current generation of digital-signal-processing (DSP) analog input-output I/O boards, hosted in Windows PCs. Its key innovation is its configurability by simple user-written text scripts, or "scriptability", which gives it a high flexibility of purpose, and allows non-programmers the capacity to rapidly design and use new Hodgkin-Huxley-type active conductances, conductances with arbitrary current-voltage relationships, Markov process conductance mechanisms with user-specified rate matrices, and hybrid networks of virtual cells. At the same time, the hardware platform allows this to be achieved with a fast and accurately timed input-computation-output cycle.

Stochastic dynamics of uncoupled neural oscillators: Fokker-Planck studies with the finite element method

We have investigated the effect of the phase response curve on the dynamics of oscillators driven by noise in two limit cases that are especially relevant for neuroscience. Using the finite element method to solve the Fokker-Planck equation we have studied (i) the impact of noise on the regularity of the oscillations quantified as the coefficient of variation, (ii) stochastic synchronization of two uncoupled phase oscillators driven by correlated noise, and (iii) their cross-correlation function. We show that, in general, the limit of type II oscillators is more robust to noise and more efficient at synchronizing by correlated noise than type I.

Effects of imperfect dynamic clamp: computational and experimental results

In the dynamic clamp technique, a typically nonlinear feedback system delivers electrical current to an excitable cell that represents the actions of "virtual" ion channels (e.g., channels that are gated by local membrane potential or by electrical activity in neighboring biological or virtual neurons). Since the conception of this technique, there have been a number of different implementations of dynamic clamp systems, each with differing levels of flexibility and performance. Embedded hardware-based systems typically offer feedback that is very fast and precisely timed, but these systems are often expensive and sometimes inflexible. PC-based systems, on the other hand, allow the user to write software that defines an arbitrarily complex feedback system, but real-time performance in PC-based systems can be deteriorated by imperfect real-time performance. Here, we systematically evaluate the performance requirements for artificial dynamic clamp knock-in of transient sodium and delayed rectifier potassium conductances. Specifically, we examine the effects of controller time step duration, differential equation integration method, jitter (variability in time step), and latency (the time lag from reading inputs to updating outputs). Each of these control system flaws is artificially introduced in both simulated and real dynamic clamp experiments. We demonstrate that each of these errors affect dynamic clamp accuracy in a way that depends on the time constants and stiffness of the differential equations being solved. In simulations, time steps above 0.2ms lead to catastrophic alteration of spike shape, but the frequency-current relationship is much more robust. Latency (the part of the time step that occurs between measuring membrane potential and injecting re-calculated membrane current) is a crucial factor as well. Experimental data are substantially more sensitive to inaccuracies than simulated data.

Optimal time scale for spike-time reliability: theory, simulations, and experiments

Use of spike timing to encode information requires that neurons respond with high temporal precision and with high reliability. Fast fluctuating stimuli are known to result in highly reproducible spike times across trials, whereas constant stimuli result in variable spike times. Here, we first studied mathematically how spike-time reliability depends on the rapidness of aperiodic stimuli. Then, we tested our theoretical predictions in computer simulations of neuron models (Hodgkin-Huxley and modified quadratic integrate-and-fire), as well as in patch-clamp experiments with real neurons (mitral cells in the olfactory bulb and pyramidal cells in the neocortex). As predicted by our theory, we found that for firing frequencies in the beta/gamma range, spike-time reliability is maximal when the time scale of the input fluctuations (autocorrelation time) is in the range of a few milliseconds (2-5 ms), coinciding with the time scale of fast synapses, and decreases substantially for faster and slower inputs. Finally, we comment how these findings relate to mechanisms causing neuronal synchronization.

Temporal patterning of saccadic eye movement signals

Electrical microstimulation is used widely in experimental neurophysiology to examine causal links between specific brain areas and their behavioral functions and is used clinically to treat neurological and psychiatric disorders in patients. Typically, microstimulation is applied to local brain regions as a train of equally spaced current pulses. We were interested in the sensitivity of a neural circuit to a train of variably spaced pulses, as is observed in physiological spike trains. We compared the effect of fixed, decelerating, accelerating, and randomly varying microstimulation patterns on the likelihood and metrics of eye movements evoked from the frontal eye field of monkeys, while holding the mean interpulse interval constant. Our results demonstrate that the pattern of microstimulation pulses strongly influences the probability of evoking a saccade, as well as the metrics of the saccades themselves. Specifically, the pattern most closely resembling physiological spike trains (accelerating pattern) was most effective at evoking a saccade, three times more so than the least effective decelerating pattern. A saccade-triggered average of effective random trains confirmed the positive relationship between accelerating rate and efficacy. These results have important implications for the use of electrical microstimulation in both experimental and clinical settings and suggest a means to study the role of temporal pattern in the encoding of behavioral and cognitive functions.

Coding with spike shapes and graded potentials in cortical networks

In cortical neurones, analogue dendritic potentials are thought to be encoded into patterns of digital spikes. According to this view, neuronal codes and computations are based on the temporal patterns of spikes: spike times, bursts or spike rates. Recently, we proposed an 'action potential waveform code' for cortical pyramidal neurones in which the spike shape carries information. Broader somatic action potentials are reliably produced in response to higher conductance input, allowing for four times more information transfer than spike times alone. This information is preserved during synaptic integration in a single neurone, as back-propagating action potentials of diverse shapes differentially shunt incoming postsynaptic potentials and so participate in the next round of spike generation. An open question has been whether the information in action potential waveforms can also survive axonal conduction and directly influence synaptic transmission to neighbouring neurones. Several new findings have now brought new light to this subject, showing cortical information processing that transcends the classical models.

Density dependent neurodynamics

The dynamics of a neural network depends on density parameters at (at least) two different levels: the subcellular density of ion channels in single neurons, and the density of cells and synapses at a network level. For the Frankenhaeuser-Huxley (FH) neural model, the density of sodium (Na) and potassium (K) channels determines the behaviour of a single neuron when exposed to an external stimulus. The features of the onset of single neuron oscillations vary qualitatively among different regions in the channel density plane. At a network level, the density of neurons is reflected in the global connectivity. We study the relation between the two density levels in a network of oscillatory FH neurons, by qualitatively distinguishing between three regions, where the mean network activity is (1) spiking, (2) oscillating with enveloped frequencies, and (3) bursting, respectively. We demonstrate that the global activity can be shifted between regions by changing either the density of ion channels at the subcellular level, or the connectivity at the network level, suggesting that different underlying mechanisms can explain similar global phenomena. Finally, we model a possible effect of anaesthesia by blocking specific inhibitory ion channels.