Ionic mechanisms underlying tonic and phasic firing behaviors in retinal ganglion cells

In the retina, the firing behaviors that ganglion cells exhibit when exposed to light stimuli are very important due to the significant roles they play in encoding the visual information. However, the detailed mechanisms, especially the intrinsic properties that generate and modulate these firing behaviors is not completely clear yet. In this study, 2 typical firing behaviors—i.e., tonic and phasic activities, which are widely observed in retinal ganglion cells (RGCs)—are investigated. A modified computational model was developed to explore the possible ionic mechanisms that underlie the generation of these 2 firing patterns. Computational results indicate that the generation of tonic and phasic activities may be attributed to the collective actions of 2 kinds of adaptation currents, i.e., an inactivating sodium current and a delayed-rectifier potassium current. The concentration of magnesium ions has crucial but differential effects in the modulation of tonic and phasic firings, when the model neuron is driven by N-methyl-D-aspartate (NMDA) -type synaptic input instead of constant current injections. The proposed model has robust features that account for the ionic mechanisms underlying the tonic and phasic firing behaviors, and it may also be used as a good candidate for modeling some other firing patterns in RGCs.     

Modeling temporal behavior of postnatal cat retinal ganglion cells

During development, mammalian retinal ganglion cells (RGCs) go through marked ontogenetic changes with respect to their excitable membrane properties. Voltage-clamp studies conducted in our laboratory have shown that the amplitude, voltage-dependence and kinetics of activation and inactivation (where present) of Na(+), K(+) and Ca(2+) conductances all exhibit developmental changes during a time when the firing patterns of mammalian ganglion cells shift from being transient to being predominantly sustained in nature. In order to better understand the contribution of each conductance to the generation of spikes and spiking patterns, we have developed a model based on our experimental data. For simplicity, we have initially used experimental data obtained from postnatal ganglion cells. At this age the ontogenetic changes observed in the characteristics of the various ionic currents are complete. Utilizing the methods adopted by Hodgkin and Huxley for the giant squid axon, we have determined rate equations for the activation and inactivation properties of the I(A), I(K dr), I(Na), I(Ca L), I(Ca N), and I(leak) currents in postnatal cat RGCs. Combining these with a simplified model of the calcium-activated potassium current (I(KCa)), we have solved and analysed the resulting differential equations. While spikes and spiking patterns resembling experimental data could be obtained from a model in which [Ca(2+)i] was averaged across the whole cell, more accurate simulations were obtained when the diffusion of intracellular Ca(2+) was modeled spatially. The resulting spatial calcium gradients were more effective in gating I(KCa), and our simulations more accurately matched the recorded amplitude and shape of individual spikes as well as the frequency of maintained discharges observed in mammalian postnatal RGCs.     

A model for repetitive firing in neurons

This paper describes a model for the generation of repetitive firing patterns in single neurons to be used as a module in large-scale network simulation studies. The model is based on the combination of extended versions of Hill's model for accomodation and of Kernell's model for adaptation. Both digital computer and electronic circuit realizations of the model are presented. The model is shown to produce strength-duration curves for accomodation which are compatible with available data from real neurons. Both “high ceiling” and “low ceiling” cell types can be matched by adjusting parameters in the model. An equation relating steady-state firing rate to amplitude of applied steady current is presented which includes the accumulation of potassium conductance changes with repetitive firing. The occurence of phasic and tonic responses to step stimulation is mapped in the parameter space of the model. Several representative response patterns to irregular inputs are presented.     

Visual pattern discrimination by population retinal ganglion cells’ activities during natural movie stimulation

In the visual system, neurons often fire in synchrony, and it is believed that synchronous activities of group neurons are more efficient than single cell response in transmitting neural signals to down-stream neurons. However, whether dynamic natural stimuli are encoded by dynamic spatiotemporal firing patterns of synchronous group neurons still needs to be investigated. In this paper we recorded the activities of population ganglion cells in bullfrog retina in response to time-varying natural images (natural scene movie) using multi-electrode arrays. In response to some different brief section pairs of the movie, synchronous groups of retinal ganglion cells (RGCs) fired with similar but different spike events. We attempted to discriminate the movie sections based on temporal firing patterns of single cells and spatiotemporal firing patterns of the synchronous groups of RGCs characterized by a measurement of subsequence distribution discrepancy. The discrimination performance was assessed by a classification method based on Support Vector Machines. Our results show that different movie sections of the natural movie elicited reliable dynamic spatiotemporal activity patterns of the synchronous RGCs, which are more efficient in discriminating different movie sections than the temporal patterns of the single cells’ spike events. These results suggest that, during natural vision, the down-stream neurons may decode the visual information from the dynamic spatiotemporal patterns of the synchronous group of RGCs’ activities.     

Processing of behaviorally relevant temporal parameters of acoustic stimuli by single neurons in the superior olivary nucleus of the leopard frog

Summary Response characteristics of 130 single neurons in the superior olivary nucleus of the northern leopard frog (Rana pipiens pipiens) were examined to determine their selectivity to various behaviorally relevant temporal parameters [rise-fall time, duration, and amplitude modulation (AM) rate of acoustic signals. Response functions were constructed with respect to each of these variables. Neurons with different temporal firing patterns such as tonic, phasic or phasic-burst firing patterns, participated in time domain analysis in specific manners. Phasic neurons manifested preferences for signals with short rise-fall times, thus possessing low-pass response functions with respect to this stimulus parameter; conversely, tonic and phasic-burst units were non-selective and possessed all-pass response functions. A distinction between temporal firing patterns was also observed for duration coding. Whereas phasic units showed no change in the mean spike count with a change in stimulus duration (i.e., all-pass duration response functions), tonic and phasic-burst units gave higher mean spike counts with an increase in stimulus duration (i.e., primary-like high-pass response functions). Phasic units manifested greater response selectivity for AM rate than did tonic or phasic-burst units, and many phasic units were tuned to a narrow range of modulation rates (i.e., band-pass). The results suggest that SON neurons play an important role in the processing of complex acoustic patterns; they perform extensive computations on AM rate as well as other temporal parameters of complex sounds. Moreover, the response selectivities for rise-fall time, duration, and AM rate could often be shown to contribute to the differential responses to complex synthetic and natural sounds.Abbreviations SON superior olivary nucleus - DMN dorsal medullary nucleus - TS torus semicircularis - FTC frequency threshold curve - BF best excitatory frequency - PAM pulsatile amplitude modulation - SAM sinusoidal amplitude modulation - SQAM square-wave amplitude modulation - MTF modulation transfer function - PSTH peri-stimulus time histogram     

Activity-mediated accumulation of potassium induces a switch in firing pattern and neuronal excitability type

During normal neuronal activity, ionic concentration gradients across a neuron’s membrane are often assumed to be stable. Prolonged spiking activity, however, can reduce transmembrane gradients and affect voltage dynamics. Based on mathematical modeling, we investigated the impact of neuronal activity on ionic concentrations and, consequently, the dynamics of action potential generation. We find that intense spiking activity on the order of a second suffices to induce changes in ionic reversal potentials and to consistently induce a switch from a regular to an intermittent firing mode. This transition is caused by a qualitative alteration in the system’s voltage dynamics, mathematically corresponding to a co-dimension-two bifurcation from a saddle-node on invariant cycle (SNIC) to a homoclinic orbit bifurcation (HOM). Our electrophysiological recordings in mouse cortical pyramidal neurons confirm the changes in action potential dynamics predicted by the models: (i) activity-dependent increases in intracellular sodium concentration directly reduce action potential amplitudes, an effect typically attributed solely to sodium channel inactivation; (ii) extracellular potassium accumulation switches action potential generation from tonic firing to intermittently interrupted output. Thus, individual neurons may respond very differently to the same input stimuli, depending on their recent patterns of activity and/or the current brain-state.     

Rhythmic activities of hypothalamic magnocellular neurons: Autocontrol mechanisms

Electrophysiological recordings in lactating rats show that oxytocin (OT) and vasopressin (AVP) neurons exhibit specific patterns of activities in relation to peripheral stimuli: periodic bursting firing for OT neurons during suckling, phasic firing for AVP neurons during hyperosmolarity (systemic injection of hypertonic saline). These activities are autocontrolled by OT and AVP released somato-dentritically within the hypothalamic magnocellular nuclei. In vivo, OT enhances the amplitude and frequency of bursts, an effect accompanied with an increase in basal firing rate. However, the characteristics of firing change as facilitation proceeds: the spike patterns become very irregular with clusters of spikes spaced by long silences; the firing rate is highly variable and clearly oscillates before facilitated bursts. This unstable behaviour dramatically decreases during intense tonic activation which temporarily interrupts bursting, and could therefore be a prerequisite for bursting. In vivo, the effects of AVP depend on the initial firing pattern of AVP neurons: AVP excites weakly active neurons (increasing duration of active periods and decreasing silences), inhibits highly active neurons, and does not affect neurons with intermediate phasic activity. AVP brings the entire population of AVP neurons to discharge with a medium phasic activity characterised by periods of firing and silence lasting 20–40 s, a pattern shown to optimise the release of AVP from the neurohypophysis. Each of the peptides (OT or AVP) induces an increase in intracellular Ca²⁺ concentration, specifically in the neurons containing either OT or AVP respectively. OT evokes the release of Ca²⁺ from IP3-sensitive intracellular stores. AVP induces an influx of Ca²⁺ through voltage-dependent Ca²⁺ channels of T-, L- and N-types. We postulate that the facilitatory autocontrol of OT and AVP neurons could be mediated by Ca²⁺ known to play a key role in the control of the patterns of phasic neurons.     

Tonic and phasic receptor neurons in the vertebrate olfactory epithelium

Olfactory receptor neurons (ORNs) respond to odorants with characteristic patterns of action potentials that are relevant for odor coding. Prolonged odorant exposures revealed three populations of dissociated toad ORNs, which were mimicked by depolarizing currents: tonic (TN, displaying sustained firing, 49% of 102 cells), phasic (PN, exhibiting brief action potential trains, 36%) and intermediate neurons (IN, generating trains longer than PN, 15%). We studied the biophysical properties underlying the differences between TNs and PNs, the most extreme cases among ORNs. TNs and PNs possessed similar membrane capacitances (approximately 4 pF), but they differed in resting potential (-82 versus -64 mV), input resistance (4.2 versus 2.9 G(Omega)) and unspecific current, I(u) (TNs: 0 < I(u) 1 pA/pF). Firing behavior did not correlate with differences in voltage-gated conductances. We developed a mathematical model that accurately simulates tonic and phasic patterns. Whole cell recordings from rat ORNs in fragments (approximately 4 mm(2)) of olfactory epithelium showed that such a tissue normally contains tonic and phasic receptor neurons, suggesting that this feature is common across a wide range of vertebrates. Our findings show that the individual passive electrical properties can govern the firing patterns of ORNs.     

Gap junction permeability modulated by dopamine exerts effects on spatial and temporal correlation of retinal ganglion cells’ firing activities

Synchronized activities among retinal ganglion cells (RGCs) via gap junctions can be increased by exogenous dopamine (DA). During DA application, single neurons’ firing activities become more synchronized with its adjacent neighbors. One intriguing question is how the enhanced spatial synchronization alters the temporal firing structure of single neurons. In the present study, firing activities of bullfrog’s dimming detectors in response to binary pseudo-random checker-board flickering were recorded via a multi-channel recording system. DA was applied in the retina to modulate synchronized activities between RGCs, and the effect of DA on firing activities of single neurons was examined. It was found that, during application of DA, synchronized activities between single neuron and its neighboring neurons was enhanced. At the meantime, the temporal structures of single neuron spike train changed significantly, and the temporal correlation in single neuron’s response decreased. The pharmacological study results indicated that the activation of D1 receptor might have effects on gap junction permeability between RGCs. Our results suggested that the dopaminergic pathway participated in the modulation of spatial and temporal correlation of RGCs’ firing activities, and may exert critical effects on visual information processing in the retina.     

Modelling intrinsic electrophysiological properties of ON and OFF retinal ganglion cells

ON and OFF retinal ganglion cells (RGCs) display differences in their intrinsic electrophysiology: OFF cells maintain spontaneous activity in the absence of any input, exhibit subthreshold membrane potential oscillations, rebound excitation and burst firing; ON cells require excitatory input to drive their activity and display none of the aforementioned phenomena. The goal of this study was to identify and characterize ionic currents that explain these intrinsic electrophysiological differences between ON and OFF RGCs. A mathematical model of the electrophysiological properties of ON and OFF RGCs was constructed and validated using published patch-clamp data from isolated intact mouse retina. The model incorporates three ionic currents hypothesized to play a role in generating behaviors that are different between ON and OFF RGCs. These currents are persistent Na⁺, I _NaP, hyperpolarization-activated, I _h, and low voltage activated Ca^2 +, I _T, currents. Using computer simulations of Hodgkin-Huxley type neuron with a single compartment model we found two distinct sets of I _NaP, I _h, I _T conductances that correspond to ON and OFF RGCs populations. Simulations indicated that special properties of I _T explain the differences in intrinsic electrophysiology between ON and OFF RGCs examined here. The modelling shows that the maximum conductance of I _T is higher in OFF than in ON cells, in agreement with recent experimental data.     

The development of intrinsic excitability in mouse retinal ganglion cells

Activity-dependent refinement of synaptic connections occurs throughout the developing nervous system, including the visual system. Retinal ganglion cells (RGCs) overproduce synapses then refine them in an activity-dependent manner that segregates RGC connections into multicellular patterns, such as eye-specific regions and retinotopic maps. Ferrets additionally segregate ON and OFF retinogeniculate pathways in an activity-dependent manner. It was unknown whether differences in ON versus OFF intrinsic and spontaneous activity occur in postnatal mouse. The work reported here measured the intrinsic properties and spontaneous activity of morphologically identified postnatal mouse RGCs, and tested the hypothesis that mouse ON and OFF RGCs develop differences in spontaneous activity. We found developmental changes in resting potential, action potential threshold, depolarization to threshold, action potential width, action potential patterns, and maximal firing rates. These results are consistent with the maturation of the intrinsic properties of RGCs extending through the first three postnatal weeks. However, there were no differences among mouse ON, OFF, and multistratified RGCs in intrinsic excitability, spontaneous synaptic drive or spontaneous action potential patterns. The absence of differences between ON and OFF activity patterns is unlike the differences that arise in ferrets. In contrast to the ferret, the ON and OFF target neurons in the mouse are organized in a random pattern, not layers. This supports the hypothesis that the absence of systematic differences in activity results in the nonlayered distribution of retinogeniculate connections.     

Period doubling of calcium spike firing in a model of a Purkinje cell dendrite

Recordings from cerebellar Purkinje cell dendrites have revealed that in response to sustained current injection, the cell firing pattern can move from tonic firing of Ca²⁺ spikes to doublet firing and even to quadruplet firing or more complex firing. These firing patterns are not modified substantially if Na⁺ currents are blocked. We show that the experimental results can be viewed as a slow transition of the neuronal dynamics through a period-doubling bifurcation. To further support this conclusion and to understand the underlying mechanism that leads to doublet firing, we develop and study a simple, one-compartment model of Purkinje cell dendrite. The neuron can also exhibit quadruplet and chaotic firing patterns that are similar to the firing patterns that some of the Purkinje cells exhibit experimentally. The effects of parameters such as temperature, applied current, and potassium reversal potential in the model resemble their effects in experiments. The model dynamics involve three time scales. Ca²⁺- dependent K⁺ currents, with intermediate time scales, are responsible for the appearance of doublet firing, whereas a very slow hyperpolarizing current transfers the neuron from tonic to doublet firing. We use the fast-slow analysis to separate the effects of the three time scales. Fast-slow analysis of the neuronal dynamics, with the activation variable of the very slow, hyperpolarizing current considered as a parameter, reveals that the transitions occurs via a cascade of period-doubling bifurcations of the fast and intermediate subsystem as this slow variable increases. We carry out another analysis, with the Ca²⁺ concentration considered as a parameter, to investigate the conditions for the generation of doublet firing in systems with one effective variable with intermediate time scale, in which the rest state of the fast subsystem is terminated by a saddle-node bifurcation. We find that the scenario of period doubling in these systems can occur only if (1) the time scale of the intermediate variable (here, the decay rate of the calcium concentration) is slow enough in comparison with the interspike interval of the tonic firing at the transition but is not too slow and (2) there is a bistability of the fast subsystem of the spike-generating variables.     

Phasic Firing in Vasopressin Cells: Understanding Its Functional Significance through Computational Models

Vasopressin neurons, responding to input generated by osmotic pressure, use an intrinsic mechanism to shift from slow irregular firing to a distinct phasic pattern, consisting of long bursts and silences lasting tens of seconds. With increased input, bursts lengthen, eventually shifting to continuous firing. The phasic activity remains asynchronous across the cells and is not reflected in the population output signal. Here we have used a computational vasopressin neuron model to investigate the functional significance of the phasic firing pattern. We generated a concise model of the synaptic input driven spike firing mechanism that gives a close quantitative match to vasopressin neuron spike activity recorded in vivo, tested against endogenous activity and experimental interventions. The integrate-and-fire based model provides a simple physiological explanation of the phasic firing mechanism involving an activity-dependent slow depolarising afterpotential (DAP) generated by a calcium-inactivated potassium leak current. This is modulated by the slower, opposing, action of activity-dependent dendritic dynorphin release, which inactivates the DAP, the opposing effects generating successive periods of bursting and silence. Model cells are not spontaneously active, but fire when perturbed by random perturbations mimicking synaptic input. We constructed one population of such phasic neurons, and another population of similar cells but which lacked the ability to fire phasically. We then studied how these two populations differed in the way that they encoded changes in afferent inputs. By comparison with the non-phasic population, the phasic population responds linearly to increases in tonic synaptic input. Non-phasic cells respond to transient elevations in synaptic input in a way that strongly depends on background activity levels, phasic cells in a way that is independent of background levels, and show a similar strong linearization of the response. These findings show large differences in information coding between the populations, and apparent functional advantages of asynchronous phasic firing.     

Respiratory rhythmicity in the cat.

Brain stem respiratory neuron activity in the cat was studied in relation to efferent outflow (phrenic discharge) under the influence of several forcing inputs: 1) CO2 tension: hypocapnia produces disappearance of firing in some neurons, and conversion of respiratory-modulated to continuous (tonic) firing in others. 2) Lung inflation: during the Bruer-Hering reflex, some neurons have "classical" responses and others have "paradoxical" responses (i.e., opposite in direction to peripheral discharge). 3) Electrical stimulation: stimulus trains to the pneumotaxic center region (rostral lateral pons) produce phase-switching, whose threshold is: a) sharp (indicating action of positive-feedback mechanisms), and b) dependent on timing of stimulus delivery (indicating continuous excitability changes during each respiratory phase). Auto- and crosscorrelation analysis revealed the existence of short-term interactions between: a) medullary inspiratory (I) neurons and phrenic motoneurons; b) pairs of medullary I neurons; c) medullary I neurons and expiratory (E) neurons. A model of the respiratory oscillator is presented, in which the processes of conversion of tonic to phasic activity and switching of the respiratory phases are explained by recurrent excitatory and inhibitory loops.     

Methamphetamine neurotoxicity decreases phasic, but not tonic, dopaminergic signaling in the rat striatum

Neurotoxic doses of methamphetamine (METH) are known to cause depletions in striatal dopamine (DA) tissue content. However, the effects of METH-induced insults on dopaminergic neurotransmission are not fully understood. Here, we employed fast-scan cyclic voltammetry at a carbon-fiber microelectrode in the anesthetized rat striatum to assess the effects of a neurotoxic regimen of METH on phasic and tonic modes of dopaminergic signaling and underlying mechanisms of DA release and uptake. Extracellular DA was electrically evoked by stimulation of the medial forebrain bundle mimicking tonic and phasic firing patterns for dopaminergic cells and was monitored simultaneously in both the dorsomedial and dorsolateral striatum. Kinetic analysis of evoked recordings determined parameters describing DA release and uptake. Striatal DA tissue content was quantified by high performance liquid chromatography with electrochemical detection. METH-pretreatment (four doses of 7.5 or 10.0 mg/kg s.c.) induced DA depletions of ～ 40% on average, which are reported in both striatal subregions. METH pre-treatment significantly decreased the amplitude of signals evoked by phasic, but not tonic, stimulation. Parameters for DA release and uptake were also similarly reduced by ～ 40%, consistent with effects on evoked phasic-like responses and DA tissue content. Taken together, these results suggest that METH-pretreatment selectively diminishes phasic, but not tonic, dopaminergic signaling in the dorsal striatum.     

Targeted Expression of Channelrhodopsin-2 to the Axon Initial Segment Alters the Temporal Firing Properties of Retinal Ganglion Cells

The axon initial segment (AIS) is essential for initiating action potentials and maintaining neuronal excitability in axon-bearing neurons in the CNS. There is increasing interest in the targeting of optogenetic tools to subcellular compartments, including the AIS, to gain precise control of neuronal activity for basic research and clinical applications. In particular, targeted expression of optogenetic tools in retinal ganglion cells (RGCs) has been explored as an approach for restoring vision after photoreceptor degeneration. Thus, understanding the effects of such targeting on spiking abilities and/or patterns is important. Here, we examined the effects of recombinant adeno-associated virus (rAAV)-mediated targeted expression of channelrhodopsin-2 (ChR2)-GFP with a Na_V channel motif in mouse RGCs. We found that this targeted expression disrupted Na_V channel clustering at the AIS and converted the spike firing patterns of RGCs from sustained to transient. Our results suggest that the clustering of membrane channels, including Na_V channels, at the AIS is important for the ability of RGCs to generate sustained spike firing. Additionally, the targeting of optogenetic tools to the AIS with the Na_V channel motif may offer a way to create transient light responses in RGCs for vision restoration.     

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.     

Electrical activity of rat retinal ganglion cells

The electrical activity of rat retinal ganglion cells is described. It was found that most such cells generate tonic discharges, while cells that demonstrate a phasic type of activity are less numerous. Neirofiziologiya/Neurophysiology, Vol. 39, Nos. 4/5, pp. 382–384, July–October, 2007.     

Properties of mouse retinal ganglion cell dendritic growth during postnatal development

The property of dendritic growth dynamics during development is a subject of intense interest. Here, we investigated the dendritic motility of retinal ganglion cells (RGCs) during different developmental stages, using ex vivo mouse retina explant culture, Semliki Forest Virus transfection and time-lapse observations. The results illustrated that during development, the dendritic motility underwent a change from rapid growth to a relatively stable state, i.e., at P0 (day of birth), RGC dendrites were in a highly active state, whereas at postnatal 13 (P13) they were more stable, and at P3 and P8, the RGCs were in an intermediate state. At any given developmental stage, RGCs of different types displayed the same dendritic growth rate and extent. Since the mouse is the most popular mammalian model for genetic manipulation, this study provided a methodological foundation for further exploring the regulatory mechanisms of dendritic development.     

Intrinsic membrane properties and dynamics of medial vestibular neurons: a simulation

The vestibulo-ocular and vestibulo-spinal network provides the ability to hold gaze fixed on an object during passive head movement. Within that network, most of the second-order neurons of the medial vestibular nucleus (MVNn) compute internal representations of head movement velocity in the horizontal plane. Our previous in vitro studies of the MVNn membrane properties indicated that they may play a major role in determining the dynamic properties of these neurons independently of their connectivity. The present study investigated that hypothesis at a theoretical level. Biophysical models of type A and B MVNn were developed. Two factors were found to be important in modeling tonic and phasic firing activity: the activation of the delayed potassium current and the rate of calcium flux. In addition, the model showed that the strength of the delayed potassium current may determine the different forms of action potentials observed experimentally. These two models (type A and B cells) were examined using depolarizing stimulation, random noise, step, ramp and sinusoidal inputs. For random noise, type A cells showed stable (regular) firing frequencies, while type B cells exhibited irregular activity. With step stimulation, the models exhibited tonic and phasic firing responses, respectively. Using ramp stimulations, frequency versus current curves showed a linear response for the type B neuron model. Finally, with sinusoidal stimulation of increasing frequencies, the type A model demonstrated a decrease in sensitivity, while the type B model exhibited an increase in sensitivity. These theoretical results support the hypothesis that MVNn intrinsic membrane properties specify various types of dynamic properties amongst these cells and therefore contribute to the wide range of dynamic responses which characterize the vestibulo-ocular and vestibulo-spinal network. Received: 1 August 1996 / Accepted in revised form: 16 December 1998