Electrical coupling between model midbrain dopamine neurons: effects on firing pattern and synchrony

Interactions among diameter, myelination, and the Na/K pump affect axonal resilience to high-frequency spiking

The reliability of spike propagation in axons is determined by complex interactions among ionic currents, ion pumps, and morphological properties. We use compartment-based modeling to reveal that interactions of diameter, myelination, and the Na/K pump determine the reliability of high-frequency spike propagation. By acting as a “reservoir” of nodal Na⁺ influx, myelinated compartments efficiently increase propagation reliability. Although spike broadening was thought to oppose fast spiking, its effect on spike propagation is complicated, depending on the balance of Na⁺ channel inactivation gate recovery, Na⁺ influx, and axial charge. Our findings suggest that slow Na⁺ removal influences axonal resilience to high-frequency spike propagation and that different strategies may be required to overcome this constraint in different neurons.

Axons reliably conduct action potentials between neurons and/or other targets. Axons have widely variable diameters and can be myelinated or unmyelinated. Although the effect of these factors on propagation speed is well studied, how they constrain axonal resilience to high-frequency spiking is incompletely understood. Maximal firing frequencies range from ∼1 Hz to >300 Hz across neurons, but the process by which Na/K pumps counteract Na⁺ influx is slow, and the extent to which slow Na⁺ removal is compatible with high-frequency spiking is unclear. Modeling the process of Na⁺ removal shows that large-diameter axons are more resilient to high-frequency spikes than are small-diameter axons, because of their slow Na⁺ accumulation. In myelinated axons, the myelinated compartments between nodes of Ranvier act as a “reservoir” to slow Na⁺ accumulation and increase the reliability of axonal propagation. We now find that slowing the activation of K⁺ current can increase the Na⁺ influx rate, and the effect of minimizing the overlap between Na⁺ and K⁺ currents on spike propagation resilience depends on complex interactions among diameter, myelination, and the Na/K pump density. Our results suggest that, in neurons with different channel gating kinetic parameters, different strategies may be required to improve the reliability of axonal propagation.

Synchronicity: The Role of Midbrain Dopamine in Whole-Brain Coordination

Midbrain dopamine seems to play an outsized role in motivated behavior and learning. Widely associated with mediating reward-related behavior, decision making, and learning, dopamine continues to generate controversies in the field. While many studies and theories focus on what dopamine cells encode, the question of how the midbrain derives the information it encodes is poorly understood and comparatively less addressed. Recent anatomical studies suggest greater diversity and complexity of afferent inputs than previously appreciated, requiring rethinking of prior models. Here, we elaborate a hypothesis that construes midbrain dopamine as implementing a Bayesian selector in which individual dopamine cells sample afferent activity across distributed brain substrates, comprising evidence to be evaluated on the extent to which stimuli in the on-going sensorimotor stream organizes distributed, parallel processing, reflecting implicit value. To effectively generate a temporally resolved phasic signal, a population of dopamine cells must exhibit synchronous activity. We argue that synchronous activity across a population of dopamine cells signals consensus across distributed afferent substrates, invigorating responding to recognized opportunities and facilitating further learning. In framing our hypothesis, we shift from the question of how value is computed to the broader question of how the brain achieves coordination across distributed, parallel processing. We posit the midbrain is part of an “axis of agency” in which the prefrontal cortex (PFC), basal ganglia (BGS), and midbrain form an axis mediating control, coordination, and consensus, respectively.

Calcium-activated SK channels control firing regularity by modulating sodium channel availability in midbrain dopamine neurons

Dopamine neurons in the substantia nigra pars compacta and ventral tegmental area regulate behaviours such as reward-related learning, and motor control. Dysfunction of these neurons is implicated in Schizophrenia, addiction to drugs, and Parkinson’s disease. While some dopamine neurons fire single spikes at regular intervals, others fire irregular single spikes interspersed with bursts. Pharmacological inhibition of calcium-activated potassium (SK) channels increases the variability in their firing pattern, sometimes also increasing the number of spikes fired in bursts, indicating that SK channels play an important role in maintaining dopamine neuron firing regularity and burst firing. However, the exact mechanisms underlying these effects are still unclear. Here, we develop a biophysical model of a dopamine neuron incorporating ion channel stochasticity that enabled the analysis of availability of ion channels in multiple states during spiking. We find that decreased firing regularity is primarily due to a significant decrease in the AHP that in turn resulted in a reduction in the fraction of available voltage-gated sodium channels due to insufficient recovery from inactivation. Our model further predicts that inhibition of SK channels results in a depolarisation of action potential threshold along with an increase in its variability.

Mesocorticolimbic dopamine functioning in primary psychopathy: A source of within-group heterogeneity

Despite similar emotional deficiencies, primary psychopathic individuals can be situated on a continuum that spans from controlled to disinhibited. The constructs on which primary psychopaths are found to diverge, such as self-control, cognitive flexibility, and executive functioning, are crucially regulated by dopamine (DA). As such, the goal of this review is to examine which specific alterations in the meso-cortico-limbic DA system and corresponding genes (e.g., TH, DAT, COMT, DRD2, DRD4) might bias development towards a more controlled or disinhibited expression of primary psychopathy. Based on empirical data, it is argued that primary psychopathy is generally related to a higher tonic and population activity of striatal DA neurons and lower levels of D2-type DA receptors in meso-cortico-limbic projections, which may boost motivational drive towards incentive-laden goals, dampen punishment sensitivity, and increase future reward-expectancy. However, increasingly higher levels of DA activity in the striatum (moderate versus pathological elevations), lower levels of DA functionality in the prefrontal cortex, and higher D1-to-D2-type receptor ratios in meso-cortico-limbic projections may lead to increasingly disinhibited and impetuous phenotypes of primary psychopathy. Finally, in order to provide a more coherent view on etiological mechanisms, we discuss interactions between DA and serotonin that are relevant for primary psychopathy.

Understanding dopamine and reinforcement learning: the dopamine reward prediction error hypothesis

A number of recent advances have been achieved in the study of midbrain dopaminergic neurons. Understanding these advances and how they relate to one another requires a deep understanding of the computational models that serve as an explanatory framework and guide ongoing experimental inquiry. This intertwining of theory and experiment now suggests very clearly that the phasic activity of the midbrain dopamine neurons provides a global mechanism for synaptic modification. These synaptic modifications, in turn, provide the mechanistic underpinning for a specific class of reinforcement learning mechanisms that now seem to underlie much of human and animal behavior. This review describes both the critical empirical findings that are at the root of this conclusion and the fantastic theoretical advances from which this conclusion is drawn.

Spike frequency adaptation is developmentally regulated in substantia nigra pars compacta dopaminergic neurons

Dopaminergic neurons of the substantia nigra pars compacta play a key role in the modulation of basal ganglia and provide a reward-related teaching signal essential for adaptative motor control. They are generally considered as a homogenous population despite several chemical and electrophysiological heterogeneities, which could underlie different preferential patterns of activity and/or different roles. Using whole-cell patch-clamp recordings in juvenile rat brain slices, we observed that the evoked activity of dopaminergic neurons displays variable spike frequency adaptation patterns. The intensity of spike frequency adaptation decreased during post-natal development. The adaptation was associated with an increase in the initial firing frequency due to faster kinetics of the afterhyperpolarization component of the spike. Adaptation was enhanced when small conductance calcium-activated potassium (SK) channels were blocked with bath application of apamine. Lastly, spike frequency adaptation of the evoked discharge was associated with more irregularity in the spontaneous firing pattern. Altogether these results show a developmental heterogeneity and electrophysiological maturation of substantia nigra dopaminergic neurons.

Conflicting effects of excitatory synaptic and electric coupling on the dynamics of square-wave bursters

Using two-cell and 50-cell networks of square-wave bursters, we studied how excitatory coupling of individual neurons affects the bursting output of the network. Our results show that the effects of synaptic excitation vs. electrical coupling are distinct. Increasing excitatory synaptic coupling generally increases burst duration. Electrical coupling also increases burst duration for low to moderate values, but at sufficiently strong values promotes a switch to highly synchronous bursts where further increases in electrical or synaptic coupling have a minimal effect on burst duration. These effects are largely mediated by spike synchrony, which is determined by the stability of the in-phase spiking solution during the burst. Even when both coupling mechanisms are strong, one form (in-phase or anti-phase) of spike synchrony will determine the burst dynamics, resulting in a sharp boundary in the space of the coupling parameters. This boundary exists in both two cell and network simulations. We use these results to interpret the effects of gap-junction blockers on the neuronal circuitry that underlies respiration.

Calcium-activated non-selective cation currents are involved in generation of tonic and bursting activity in dopamine neurons of the substantia nigra pars compacta

Nigral dopamine neurons are transiently activated by high frequency glutamatergic inputs relaying reward-predicting sensory information. The tonic firing pattern of dopamine cells responds to these inputs with a transient burst of spikes that requires NMDA receptors. Here, we show that NMDA receptor activation further excites the cell by recruiting a calcium-activated non-selective cation current (ICAN) capable of generating a plateau potential. Burst firing in vitro is eliminated after blockade of ICAN with flufenamic acid, 9-phenanthrol, or intracellular BAPTA. ICAN is likely to be mediated by a transient receptor potential (TRP) channel, and RT-PCR was used to confirm expression of TRPM2 and TRPM4mRNA in substantia nigra pars compacta.We propose that ICAN is selectively activated during burst firing to boost NMDA currents and allow plateau potentials. This boost mechanism may render DA cells vulnerable to excitotoxicity.

Convergent processing of both positive and negative motivational signals by the VTA dopamine neuronal populations

Dopamine neurons in the ventral tegmental area (VTA) have been traditionally studied for their roles in reward-related motivation or drug addiction. Here we study how the VTA dopamine neuron population may process fearful and negative experiences as well as reward information in freely behaving mice. Using multi-tetrode recording, we find that up to 89% of the putative dopamine neurons in the VTA exhibit significant activation in response to the conditioned tone that predict food reward, while the same dopamine neuron population also respond to the fearful experiences such as free fall and shake events. The majority of these VTA putative dopamine neurons exhibit suppression and offset-rebound excitation, whereas ∼25% of the recorded putative dopamine neurons show excitation by the fearful events. Importantly, VTA putative dopamine neurons exhibit parametric encoding properties: their firing change durations are proportional to the fearful event durations. In addition, we demonstrate that the contextual information is crucial for these neurons to respectively elicit positive or negative motivational responses by the same conditioned tone. Taken together, our findings suggest that VTA dopamine neurons may employ the convergent encoding strategy for processing both positive and negative experiences, intimately integrating with cues and environmental context.

Mechanisms explaining transitions between tonic and phasic firing in neuronal populations as predicted by a low dimensional firing rate model

Several firing patterns experimentally observed in neural populations have been successfully correlated to animal behavior. Population bursting, hereby regarded as a period of high firing rate followed by a period of quiescence, is typically observed in groups of neurons during behavior. Biophysical membrane-potential models of single cell bursting involve at least three equations. Extending such models to study the collective behavior of neural populations involves thousands of equations and can be very expensive computationally. For this reason, low dimensional population models that capture biophysical aspects of networks are needed. The present paper uses a firing-rate model to study mechanisms that trigger and stop transitions between tonic and phasic population firing. These mechanisms are captured through a two-dimensional system, which can potentially be extended to include interactions between different areas of the nervous system with a small number of equations. The typical behavior of midbrain dopaminergic neurons in the rodent is used as an example to illustrate and interpret our results. The model presented here can be used as a building block to study interactions between networks of neurons. This theoretical approach may help contextualize and understand the factors involved in regulating burst firing in populations and how it may modulate distinct aspects of behavior.

Synchronous regimes in ensembles of coupled Bonhoeffer-van der Pol oscillators

We study synchronous behavior in ensembles of locally coupled nonidentical Bonhoeffer-van der Pol oscillators. We show that, in a chain of N elements not less than 2;{N-1}, different coexisting regimes of global synchronization are possible, and we investigate wave-induced synchronous regimes in a chain and in a lattice of coupled nonidentical Bonhoeffer-van der Pol oscillators.

Chemical transmission between dopaminergic neuron pairs

Midbrain dopaminergic (DAergic) neurons play a major regulatory role in in goal-directed behavior and reinforcement learning. DAergic neuron activity, and therefore spatiotemporal properties of dopamine release, precisely encodes reward signals. Neuronal activity is shaped both by external afferences and local interactions (chemical and electrical transmissions). Numerous hints suggest the existence of chemical interactions between DAergic neurons, but direct evidence and characterization are still lacking. Here, we show, using dual patch-clamp recordings in rat brain slices, a widespread bidirectional chemical transmission between DAergic neuron pairs. Hyperpolarizing postsynaptic potentials were partially mediated by D2-like receptors, and entirely resulted from the inhibition of the hyperpolarization-activated depolarizing current (Ih). These results constitute the first evidence in paired recordings of a chemical transmission relying on conductance decrease in mammals. In addition, we show that chemical transmission and electrical synapses frequently coexist within the same neuron pair and dynamically interact to shape DAergic neuron activity.

From the Golgi–Cajal mapping to the transmitter-based characterization of the neuronal networks leading to two modes of brain communication: Wiring and volume transmission

After Golgi-Cajal mapped neural circuits, the discovery and mapping of the central monoamine neurons opened up for a new understanding of interneuronal communication by indicating that another form of communication exists. For instance, it was found that dopamine may be released as a prolactin inhibitory factor from the median eminence, indicating an alternative mode of dopamine communication in the brain. Subsequently, the analysis of the locus coeruleus noradrenaline neurons demonstrated a novel type of lower brainstem neuron that monosynaptically and globally innervated the entire CNS. Furthermore, the ascending raphe serotonin neuron systems were found to globally innervate the forebrain with few synapses, and where deficits in serotonergic function appeared to play a major role in depression. We propose that serotonin reuptake inhibitors may produce antidepressant effects through increasing serotonergic neurotrophism in serotonin nerve cells and their targets by transactivation of receptor tyrosine kinases (RTK), involving direct or indirect receptor/RTK interactions. Early chemical neuroanatomical work on the monoamine neurons, involving primitive nervous systems and analysis of peptide neurons, indicated the existence of alternative modes of communication apart from synaptic transmission. In 1986, Agnati and Fuxe introduced the theory of two main types of intercellular communication in the brain: wiring and volume transmission (WT and VT). Synchronization of phasic activity in the monoamine cell clusters through electrotonic coupling and synaptic transmission (WT) enables optimal VT of monoamines in the target regions. Experimental work suggests an integration of WT and VT signals via receptor-receptor interactions, and a new theory of receptor-connexin interactions in electrical and mixed synapses is introduced. Consequently, a new model of brain function must be built, in which communication includes both WT and VT and receptor-receptor interactions in the integration of signals. This will lead to the unified execution of information handling and trophism for optimal brain function and survival.

Electrical synapses in basal ganglia

The basal ganglia (BG) provide a major integrative system of the forebrain involved in the organization of goal-directed behaviour. Pathological alteration of BG function leads to major motor and cognitive impairments such as observed in Parkinson's disease. Recent advances in BG research stress the role of neural oscillations and synchronization in the normal and pathological function of BG. As demonstrated in several brain structures, these patterns of neural activity can emerge from electrically coupled neuronal networks. This review aims at addressing the presence, functionality and putative role of electrical synapses in BG, with a particular emphasis on the striatum and the substantia nigra pars compacta (SNc), two main BG nuclei in which the existence and functional properties of neuronal coupling are best documented.

Somato-dendritic mechanisms underlying the electrophysiological properties of hypothalamic magnocellular neuroendocrine cells: a multicompartmental model study

Magnocellular neuroendocrine cells (MNCs) of the hypothalamus synthesize the neurohormones vasopressin and oxytocin, which are released into the blood and exert a wide spectrum of actions, including the regulation of cardiovascular and reproductive functions. Vasopressin- and oxytocin-secreting neurons have similar morphological structure and electrophysiological characteristics. A realistic multicompartmental model of a MNC with a bipolar branching structure was developed and calibrated based on morphological and in vitro electrophysiological data in order to explore the roles of ion currents and intracellular calcium dynamics in the intrinsic electrical MNC properties. The model was used to determine the likely distributions of ion conductances in morphologically distinct parts of the MNCs: soma, primary dendrites and secondary dendrites. While reproducing the general electrophysiological features of MNCs, the model demonstrates that the differential spatial distributions of ion channels influence the functional expression of MNC properties, and reveals the potential importance of dendritic conductances in these properties.

On reaction dynamics at dopamine synapses

Dopamine neurons play a key role in normal and pathological cognitive processes as well as in the effects of certain drugs of addiction. Models of the synapses of such neurons include transporter mechanisms and reaction dynamics. This article focuses attention on the fundamental reaction that converts tyrosine to DOPA, which involves a cofactor. The Michaelis-Menten formula for the rate of an enzymatic reaction is modified by the presence of cofactors, which may be either essential or nonessential. In the essential case, the reaction rate is found to depend on the relative magnitudes of the concentrations of the primary enzyme and the cofactor. The case of a nonessential cofactor is more complex and it is shown for the first time how this leads to reaction rate formulas that depend explicitly on the concentrations of the enzyme and cofactor. The extended Michaelis-Menten formulas are applied to the aforementioned reaction with tyrosine hydroxylase as enzyme and biopterin as cofactor. The results are useful in constructing accurate models of dopamine synapses.

Network mechanism for burst generation

We report on the mechanism of burst generation by populations of intrinsically spiking neurons, when a certain threshold in coupling strength is exceeded. These ensembles synchronize at relatively low coupling strength and lose synchronization at stronger coupling via spatiotemporal intermittency. The latter transition triggers fast repetitive spiking, which results in synchronized bursting. We present evidence that this mechanism is generic for various network topologies from regular to small-world and scale-free ones, different types of coupling and neuronal model.

Just a phase they're going through: the complex interaction of intrinsic high-threshold bursting and gap junctions in the generation of thalamic alpha and theta rhythms

Rhythms in the alpha frequency band (8-13 Hz) are a defining feature of the human EEG during relaxed wakefulness and are known to be influenced by the thalamus. In the early stages of sleep and in several neurological and psychiatric conditions alpha rhythms are replaced by slower activity in the theta (3-7 Hz) band. Of particular interest is how these alpha and theta rhythms are generated at the cellular level. Recently we identified a subset of thalamocortical (TC) neurons in the lateral geniculate nucleus (LGN) which exhibit rhythmic high-threshold (>-55 mV) bursting at approximately 2-13 Hz and which are interconnected by gap junctions (GJs). These cells combine to generate a locally synchronized continuum of alpha and theta oscillations, thus providing direct evidence that the thalamus can act as an independent pacemaker of alpha and theta rhythms. Interestingly, GJ coupled pairs of TC neurons can exhibit both in-phase and anti-phase synchrony and will often spontaneously alternate between these two states. This dictates that the local field oscillation amplitude is not simply linked to the extent of cell recruitment into a single synchronized neuronal assembly but also to the degree of destructive interference between dynamic, spatially overlapping, competing anti-phase groups of continuously bursting neurons. Thus, the waxing and waning of thalamic alpha/theta rhythms should not be assumed to reflect a wholesale increase and reduction, respectively, in underlying neuronal synchrony. We argue that these network dynamics might have important consequences for relating changes in the amplitude of EEG alpha and theta rhythms to the activity of thalamic networks.

The origin of spontaneous synchronized burst in cultured neuronal networks based on multi-electrode arrays

Many neural networks in mammalian central nervous system (CNS) fire single spike and complex spike burst. In fact, the conditions for triggering burst are not well understood. In the paper multi-electrode arrays (MEA) are used to record the spontaneous electrophysiological activities of cultured rat hippocampal neuronal network for a long time. After about 3 weeks culture, a transition from single spike to burst is observed in several networks. All of these spikes fire quickly before burst begins. The firing rate during the burst is lower than that just before the burst, but differences of inter-spike intervals (ISIs) between two firing patterns are not clear. Moreover, the electrical activities on neighboring electrodes show strong synchrony during the burst activities. In a word, the generation of the burst requires that network should have a sufficient level of excitation as well as a balance of synaptic inhibition.

Transient high-frequency firing in a coupled-oscillator model of the mesencephalic dopaminergic neuron

Dopaminergic neurons of the midbrain fire spontaneously at rates <10/s and ordinarily will not exceed this range even when driven with somatic current injection. When driven at higher rates, these cells undergo spike failure through depolarization block. During spontaneous bursting of dopaminergic neurons in vivo, bursts related to reward expectation in behaving animals, and bursts generated by dendritic application of N-methyl-d-aspartate (NMDA) agonists, transient firing attains rates well above this range. We suggest a way such high-frequency firing may occur in response to dendritic NMDA receptor activation. We have extended the coupled oscillator model of the dopaminergic neuron, which represents the soma and dendrites as electrically coupled compartments with different natural spiking frequencies, by addition of dendritic AMPA (voltage-independent) or NMDA (voltage-dependent) synaptic conductance. Both soma and dendrites contain a simplified version of the calcium-potassium mechanism known to be the mechanism for slow spontaneous oscillation and background firing in dopaminergic cells. The compartments differ only in diameter, and this difference is responsible for the difference in natural frequencies. We show that because of its voltage dependence, NMDA receptor activation acts to amplify the effect on the soma of the high-frequency oscillation of the dendrites, which is normally too weak to exert a large influence on the overall oscillation frequency of the neuron. During the high-frequency oscillations that result, sodium inactivation in the soma is removed rapidly after each action potential by the hyperpolarizing influence of the dendritic calcium-dependent potassium current, preventing depolarization block of the spike mechanism, and allowing high-frequency spiking.

Ultrastructure of jaw muscle spindle afferents within the rat trigeminal mesencephalic nucleus

This study examined the ultrastructures of neuronal elements within trigeminal mesencephalic nucleus by labeling masseteric mesencephalic neurons and masseter motoneurons with injection of horseradish peroxidase into masseteric muscle. Of eight horseradish peroxidase-labeled muscle spindle afferents examined, four terminals showed synaptic contact with labeled dendrites of masseteric motoneurons, two with labeled somata, and the remaining two with unlabeled dendrites. A few of the labeled dendrites showed intimate contact with the somata of the trigeminal mesencephalic nucleus neurons. These results provide morphological evidence of synaptic contact of recurring masseteric muscle spindle afferents with the trigeminal mesencephalic nucleus somata and also suggest the presence of electrical synapses between the somata of the trigeminal mesencephalic nucleus neurons and dendrites of jaw-closing motoneurons.

Electrical synapses between dopaminergic neurons of the substantia nigra pars compacta

Spatiotemporal properties of dopamine release play a major role both in striatal and nigral physiology because dopamine is released from nerve terminals and dendrites of nigrostriatal dopaminergic (DA) neurons. Pioneering work revealed gap junctional communication (assessed by dye-coupling experiments) between DA cells in the substantia nigra pars compacta (SNc). However, direct evidence of functional electrical synapses between DA neurons is still lacking. In this study, gap junctional communication between DA neurons was investigated in rat brain slices. Tracer coupling was observed in postnatal day 5 (P5) to P10 and P15-P25 rats. Dual whole-cell patch-clamp recordings revealed that 96% of DA neurons were coupled by electrical synapses in P7-P10 rats, and 20% were coupled in P17-P21 rats. These electrical synapses were mainly symmetrical and displayed strong low-pass filtering properties. When spontaneous firing activity was monitored, no significant synchronization was observed. Nevertheless, an efficient modulation of the spontaneous firing frequency of the postsynaptic cell occurred during injection of hyperpolarizing and depolarizing currents in the coupled presynaptic cell. Together, these observations demonstrate the existence of a fast communication between SNc DA neurons through electrical synapses.

A Modeling Study Suggests Complementary Roles for GABAA and NMDA Receptors and the SK Channel in Regulating the Firing Pattern in Midbrain Dopamine Neurons

Midbrain dopaminergic (DA) neurons in vivo exhibit two major firing patterns: single-spike firing and burst firing. The firing pattern expressed is dependent on both the intrinsic properties of the neurons and their excitatory and inhibitory synaptic inputs. Experimental data suggest that the activation of N-methyl-d-aspartate (NMDA) and GABAA receptors is a crucial contributor to the initiation and suppression of burst firing, respectively, and that blocking Ca2+-activated potassium SK channels can facilitate burst firing. A multi-compartmental model of a DA neuron with a branching structure was developed and calibrated based on in vitro experimental data to explore the effects of different levels of activation of NMDA and GABAA receptors as well as the modulation of the SK current on the firing activity. The simulated tonic activation of GABAA receptors was calibrated by taking into account the difference in the electrotonic properties in vivo versus in vitro. Although NMDA-evoked currents are required for burst generation in the model, currents evoked by GABAA-receptor activation can also regulate the firing pattern. For example, the model predicts that increasing the level of NMDA receptor activation can produce excessive depolarization that prevents burst firing, but a concurrent increase in the activation of GABAA receptors can restore burst firing. Another prediction of the model is that blocking the SK channel current in vivo will facilitate bursting, but not as robustly as blocking the GABAA receptors.

Alternative to hand-tuning conductance-based models: construction and analysis of databases of model neurons

Conventionally, the parameters of neuronal models are hand-tuned using trial-and-error searches to produce a desired behavior. Here, we present an alternative approach. We have generated a database of about 1.7 million single-compartment model neurons by independently varying 8 maximal membrane conductances based on measurements from lobster stomatogastric neurons. We classified the spontaneous electrical activity of each model neuron and its responsiveness to inputs during runtime with an adaptive algorithm and saved a reduced version of each neuron's activity pattern. Our analysis of the distribution of different activity types (silent, spiking, bursting, irregular) in the 8-dimensional conductance space indicates that the coarse grid of conductance values we chose is sufficient to capture the salient features of the distribution. The database can be searched for different combinations of neuron properties such as activity type, spike or burst frequency, resting potential, frequency-current relation, and phase-response curve. We demonstrate how the database can be screened for models that reproduce the behavior of a specific biological neuron and show that the contents of the database can give insight into the way a neuron's membrane conductances determine its activity pattern and response properties. Similar databases can be constructed to explore parameter spaces in multicompartmental models or small networks, or to examine the effects of changes in the voltage dependence of currents. In all cases, database searches can provide insight into how neuronal and network properties depend on the values of the parameters in the models.