Synaptic current at the rat ganglionic synapse and its interactions with the neuronal voltage-dependent currents

Activity-induced spontaneous spikes in GABAergic neurons suppress seizure discharges: an implication of computational modeling

Background

Epilepsy, a prevalent neurological disorder, appears self-termination. The endogenous mechanism for seizure self-termination remains to be addressed in order to develop new strategies for epilepsy treatment. We aim to examine the role of activity-induced spontaneous spikes at GABAergic neurons as an endogenous mechanism in the seizure self-termination.

Methods and Results

Neuronal spikes were induced by depolarization pulses at cortical GABAergic neurons from temporal lobe epilepsy patients and mice, in which some of these neurons fired activity-induced spontaneous spikes. Neural networks including excitatory and inhibitory neurons were computationally constructed, and their functional properties were based on our studies from whole-cell recordings. With the changes in the portion and excitability of inhibitory neurons that generated activity-induced spontaneous spike, the efficacies to suppress synchronous seizure activity were analyzed, such as its onset time, decay slope and spike frequency. The increases in the proportion and excitability of inhibitory neurons that generated activity-induced spontaneous spikes effectively suppressed seizure activity in neural networks. These factors synergistically strengthened the efficacy of seizure activity suppression.

Conclusion

Our study supports a notion that activity-induced spontaneous spikes in GABAergic neurons may be an endogenous mechanism for seizure self-termination. A potential therapeutic strategy for epilepsy is to upregulate the cortical inhibitory neurons that generate activity-induced spontaneous spikes.

Changes in cationic selectivity of the nicotinic channel at the rat ganglionic synapse: a role for chloride ions?

The permeability of the nicotinic channel (nAChR) at the ganglionic synapse has been examined, in the intact rat superior cervical ganglion in vitro, by fitting the Goldman current equation to the synaptic current (EPSC) I–V relationship. Subsynaptic nAChRs, activated by neurally-released acetylcholine (ACh), were thus analyzed in an intact environment as natively expressed by the mature sympathetic neuron. Postsynaptic neuron hyperpolarization (from −40 to −90 mV) resulted in a change of the synaptic potassium/sodium permeability ratio (P_K/P_Na) from 1.40 to 0.92, corresponding to a reversible shift of the apparent acetylcholine equilibrium potential, E_ACh, by about +10 mV. The effect was accompanied by a decrease of the peak synaptic conductance (g_syn) and of the EPSC decay time constant. Reduction of [Cl⁻]_o to 18 mM resulted in a change of P_K/P_Na from 1.57 (control) to 2.26, associated with a reversible shift of E_ACh by about −10 mV. Application of 200 nM αBgTx evoked P_K/P_Na and g_syn modifications similar to those observed in reduced [Cl⁻]_o. The two treatments were overlapping and complementary, as if the same site/mechanism were involved. The difference current before and after chloride reduction or toxin application exhibited a strongly positive equilibrium potential, which could not be explained by the block of a calcium component of the EPSC. Observations under current-clamp conditions suggest that the driving force modification of the EPSC due to P_K/P_Na changes represent an additional powerful integrative mechanism of neuron behavior. A possible role for chloride ions is suggested: the nAChR selectivity was actually reduced by increased chloride gradient (membrane hyperpolarization), while it was increased, moving towards a channel preferentially permeable for potassium, when the chloride gradient was reduced.

The nicotinic activation of the denervated sympathetic neuron of the rat

Nicotinic responses to endogenous acetylcholine and to exogenously applied agonists have been studied in the intact or denervated rat sympathetic neuron in vitro, by using the two-microelectrode voltage-clamp technique. Preganglionic denervation resulted in progressive decrease of the synaptic current (excitatory postsynaptic current, EPSC) amplitude, which disappeared within 24 h. These effects were accompanied by changes in ion selectivity of the nicotinic channel (nAChR). The extrapolated EPSC null potential (equilibrium potential for acetylcholine action, E(Syn)) shifted from a mean value of -15.9+/-0.7 mV, in control, to -7.4+/-1.6 mV, in denervated neurons, indicating a decrease of the permeability ratio for the main components of the synaptic current (P(K)/P(Na)) from 1.56 to 1.07. The overall properties of AChRs were investigated by applying dimethylphenylpiperazinium or cytisine and by examining the effects of endogenous ACh, diffusing within the ganglion after preganglionic tetanization in the presence of neostigmine. The null potentials of these macrocurrents (equilibrium potential for dimethylphenylpiperazinium action, E(DMPP); and equilibrium potential for diffusing acetylcholine, E(ACh), respectively) were evaluated by applying voltage ramps and from current-voltage plots. In normal neurons, E(Syn) (-15.9+/-0.7 mV) was significantly different from E(DMPP) (-26.1+/-1.0) and E(ACh) (-31.1+/-3.3); following denervation, nerve-evoked currents displayed marked shifts in their null potentials (E(Syn)=-7.4+/-1.6 mV), whereas the amplitude and null potential of the agonist-evoked macrocurrents were unaffected by denervation and its duration (E(DMPP)=-26.6+/-1.2 mV). It is suggested that two populations of nicotinic receptors, synaptic and extrasynaptic, are present on the neuron surface, and that only the synaptic type displays sensitivity to denervation.

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

Neurons and dynamical models of spike generation display two different classes of threshold behavior: type 1 [firing frequency vs. current (f-I) relationship is continuous at threshold] and type 2 (discontinuous f-I). With steady current or conductance stimulation, regular-spiking (RS) pyramidal neurons and fast-spiking (FS) inhibitory interneurons in layer 2/3 of somatosensory cortex exhibit type 1 and type 2 threshold behaviors, respectively. We compared the postsynaptic firing variability of type 1 RS and type 2 FS cells, during naturalistic, fluctuating conductance input. In RS neurons, increasing the level of independently random, shunting inhibition caused a monotonic increase in spike reliability, whereas in FS interneurons, there was an optimum level of shunting inhibition for achieving the most reliable spike generation and the most precise spike-time encoding. This was observed over a range of different degrees of synchrony, or correlation, in the input. RS cells displayed a progressive rise in spike jitter during natural-like transient burst inputs, whereas for FS cells, jitter was mostly kept low. Furthermore, RS cells showed encoding of the input level in the spike shape, whereas FS cells did not. These differences between the two cell types are consistent with a role of RS neurons as rate-coding integrators, and a role of FS neurons as resonators controlling the coherence of synchronous firing.

Automated classification of evoked quantal events

We provide both theoretical and computational improvements to the analysis of synaptic transmission data. Theoretically, we demonstrate the correlation structure of observations within evoked postsynaptic potentials (EPSP) are consistent with multiple random draws from a common autoregressive moving-average (ARMA) process of order (2, 2). We use this observation and standard time series results to construct a statistical hypothesis testing procedure for determining whether a given trace is an EPSP. Computationally, we implement this method in R, a freeware statistical language, which reduces the amount of time required for the investigator to classify traces into EPSPs or non-EPSPs and eliminates investigator subjectivity from this classification. In addition, we provide a computational method for calculating common functionals of EPSPs (peak amplitude, decay rate, etc.). The methodology is freely available over the internet. The automated procedure to index the quantal characteristics greatly facilitates determining if any one or multiple parameters are changing due to experimental conditions. In our experience, the software reduces the time required to perform these analyses from hours to minutes.

Synaptic and somatic effects of axotomy in the intact, innervated rat sympathetic neuron

A biophysical description of the axotomized rat sympathetic neuron is reported, obtained by the two-electrode voltage-clamp technique in mature, intact superior cervical ganglia in vitro. Multiple aspects of neuron functioning were tested. Synaptic conductance activated by the whole presynaptic input decreased to 29% of the control value (0.92 muS per neuron) 1 day after axotomy and to 18% after 3 days. Despite the decrease in amplitude of the macroscopic current, miniature excitatory postsynaptic current (mEPSC) mean conductance, acetylcholine (ACh) equilibrium potential, and EPSC decay time constant were unaffected. Synaptic efficacy was tested during paired-pulse or maintained stimulation (5, 10, and 15 Hz, 10-s duration). Quantal release in axotomized neurons was preserved during the tetanus despite the reduction of the initial EPSC amplitude, suggesting that ACh secretion depended on the number of surviving synapses; each of them exhibited dynamic behavior during trains similar to that of normal synapses. Facilitation of EPSC amplitude was noted in 2-day axotomized neurons during the first few impulses in the train. Voltage-dependent potassium currents (the delayed I(KD) and the transient I(A)) exhibited an early drastic decrease in peak amplitude; these effects persisted 7 days after axotomy. Marked changes in I(A) kinetics occurred after injury: the steady-state inactivation curve shifted by up to +17 mV toward positive potentials and the voltage sensitivity of inactivation removal became steeper. I(A) impairment was reflected in a reduced inward threshold charge for discharge and reduced spike repolarization rate. Synaptic and somatic data were applied in a mathematical model to describe the progressive decrease in the safety factor, and the eventual failure of ganglionic transmission after axotomy.

Comparison of the inhibition of Renshaw cells during subthreshold and suprathreshold conditions using anatomically and physiologically realistic models

Inhibitory synaptic inputs to Renshaw cells are concentrated on the soma and the juxtasomatic dendrites. In the present study, we investigated whether this proximal bias leads to more effective inhibition under different neuronal operating conditions. Using compartmental models based on detailed anatomical measurements of intracellularly stained Renshaw cells, we compared the inhibition produced by glycine/gamma-aminobutyric acid-A (GABA(A)) synapses when distributed with a proximal bias to the inhibition produced when the same synapses were distributed uniformly (i.e., with no regional bias). The comparison was conducted in subthreshold and suprathreshold conditions. The latter were mimicked by voltage clamping the soma to -55 mV. The voltage clamp reduces nonlinear interactions between excitatory and inhibitory synapses. We hypothesized that for electrotonically compact cells such as Renshaw cells, the strength of the inhibition would become much less dependent on synaptic location in suprathreshold conditions. This hypothesis was not confirmed. The inhibition produced when inhibitory inputs were proximally distributed was always stronger than when the same inputs were uniformly distributed. In fact, the relative effectiveness of proximally distributed inhibitory inputs over uniformly distributed synapses was greater in suprathreshold conditions than that in subthreshold conditions. The somatic voltage clamp minimized saturation of inhibitory driving potentials. Because this effect was greatest near the soma, the current produced by more distal synapses suffered a greater loss because of saturation. Conversely, in subthreshold conditions, the effectiveness of proximal synapses was substantially reduced at high levels of background synaptic activity because of saturation. Our results suggest glycine/GABA(A) synapses on Renshaw cells are strategically distributed to block the powerful excitatory drive produced by recurrent collaterals from motoneurons.

Voltage- and activity-dependent chloride conductance controls the resting status of the intact rat sympathetic neuron

Remarkable activity dependence was uncovered in the chloride conductance that operates in the subthreshold region of membrane potential, by using the two-microelectrode voltage-clamp technique in the mature and intact rat sympathetic neuron. Both direct and synaptic neuron tetanization (15 Hz, 10-s duration to saturate the response) resulted in a long-lasting (not less than 15 min) increase of cell input conductance (+70-150% 10 min after tetanus), accompanied by the onset of an inward current with the same time course. Both processes developed with similar properties in the postganglionic neuron when presynaptic stimulation was performed under current- or voltage-clamp conditions and were unaffected by external calcium on direct stimulation. The posttetanic effects were sustained by gCl increase because both conductance and current modifications were blocked by 0.5 mM Anthracene-9-carboxylic acid (a chloride channel blocker) but were unaffected by TEACl or cesium chloride treatments. The chloride channel properties were modified by stimulation: their voltage sensitivity and rate of closure in response to hyperpolarization strongly increased. The voltage dependence of the three major conductances governing the cell subthreshold status (gCl, gK, and gL) was evaluated over the -40/-110 mV membrane potential range in unstimulated neurons and compared with previous results in stimulated neurons. A drastic difference between the voltage-conductance profiles was observed, exclusively sustained by gCl increase. The chloride channel thus hosts an intrinsic mechanism, a memory of previous neuron activity, which makes the chloride current a likely candidate for natural controller of the balance between opposite resting currents and thus of membrane potential level.

Postsynaptic variability of firing in rat cortical neurons: the roles of input synchronization and synaptic NMDA receptor conductance

Neurons in the functioning cortex fire erratically, with highly variable intervals between spikes. How much irregularity comes from the process of postsynaptic integration and how much from fluctuations in synaptic input? We have addressed these questions by recording the firing of neurons in slices of rat visual cortex in which synaptic receptors are blocked pharmacologically, while injecting controlled trains of unitary conductance transients, to electrically mimic natural synaptic input. Stimulation with a Poisson train of fast excitatory (AMPA-type) conductance transients, to simulate independent inputs, produced much less variability than encountered in vivo. Addition of NMDA-type conductance to each unitary event regularized the firing but lowered the precision and reliability of spikes in repeated responses. Independent Poisson trains of GABA-type conductance transients (reversing at the resting potential), which simulated independent activity in a population of presynaptic inhibitory neurons, failed to increase timing variability substantially but increased the precision of responses. However, introduction of synchrony, or correlations, in the excitatory input, according to a nonstationary Poisson model, dramatically raised timing variability to in vivo levels. The NMDA phase of compound AMPA-NMDA events conferred a time-dependent postsynaptic variability, whereby the reliability and precision of spikes degraded rapidly over the 100 msec after the start of a synchronous input burst. We conclude that postsynaptic mechanisms add significant variability to cortical responses but that substantial synchrony of inputs is necessary to explain in vivo variability. We suggest that NMDA receptors help to implement a switch from precise firing to random firing during responses to concerted inputs.

Synaptic stimulation of nicotinic receptors in rat sympathetic ganglia is followed by slow activation of postsynaptic potassium or chloride conductances

Two slow currents have been described in rat sympathetic neurons during and after tetanization of the whole preganglionic input. Both effects are mediated by nicotinic receptors activated by native acetylcholine (ACh). A first current, indicated as IAHPsyn, is calcium dependent and voltage independent, and is consistent with an IAHP-type potassium current sustained by calcium ions accompanying the nicotinic synaptic current. The conductance activated by a standard synaptic train was approximately 3.6 nS per neuron; it was detected in isolation in 14 out of a 52-neuron sample. A novel current, IADPsyn, was described in 42/52 of the sample as a post-tetanic inward current, which increased in amplitude with increasing membrane potential negativity and exhibited a null-point close to the holding potential and the cell momentary chloride equilibrium potential. IADPsyn developed during synaptic stimulation and decayed thereafter according to a single exponential (mean tau = 148.5 ms) in 18 neurons or according to a two-exponential time course (tau = 51.8 and 364.9 ms, respectively) in 19 different neurons. The mean peak conductance activated was approximately 20 nS per neuron. IADPsyn was calcium independent, it was affected by internal and external chloride concentration, but was insensitive to specific blockers (anthracene-9-carboxylic acid, 9AC) of the chloride channels open in the resting neuron. It is suggested that gADPsyn represents a specific chloride conductance activatable by intense nicotinic stimulation; in some neurons it is even associated with single excitatory postsynaptic potentials (EPSCs). Both IAHP and IADPsyn are apparently devoted to reduce neuronal excitability during and after intense synaptic stimulation.

Nicotinic EPSCs in intact rat ganglia feature depression except if evoked during intermittent postsynaptic depolarization

The involvement of the postsynaptic membrane potential level in controlling synaptic strength at the ganglionic synapse was studied by recording nicotinic fast synaptic currents (EPSCs) from neurons in the intact, mature rat superior cervical ganglion, using the two-electrode voltage-clamp technique. EPSCs were evoked by 0.05-Hz supramaximal stimulation of the preganglionic sympathetic trunk over long periods; their peak amplitude (or synaptic charge transfer) over time appeared to depend on the potential level of the neuronal membrane where the nicotinic receptors are embedded. EPSC amplitude remained constant (n = 6) only if ACh was released within repeated depolarizing steps of the postganglionic neuron, which constantly varied between -50 and -20 mV in consecutive 10-mV steps, whereas it decreased progressively by 45% (n = 9) within 14 min when the sympathetic neuron was held at constant membrane potential. Synaptic channel activation, channel ionic permeation and depolarization of the membrane in which the nicotinic receptor is localized must occur simultaneously to maintain constant synaptic strength at the ganglionic synapse during low-rate stimulation (0.03-1 Hz). Different posttetanic (20 Hz for 10 s) behaviors were observed depending on the mode of previous stimulation. In the neuron maintained at constant holding potential during low-rate stimulation, the depressed EPSC showed posttetanic potentiation, recovering approximately 23% of the mean pretetanic values (n = 10). The maximum effect was immediate in 40% of the neurons tested and developed over a 3- to 6-min period in the others; thereafter potentiation vanished within 40 min of 0.05-Hz stimulation. In contrast, no statistically significant synaptic potentiation was observed when EPSC amplitudes were kept constant by repeated -50/-20-mV command cycles (n = 12). It is suggested that, under these conditions, posttetanic potentiation could represent an attempt at recovering the synaptic strength lost during inappropriate functioning of the ganglionic synapse.

The subunit dominates the relaxation kinetics of heteromeric neuronal nicotinic receptors

The ACh-induced voltage-jump relaxation currents of the nicotinic receptors formed by pair-wise expression of the rat alpha2, alpha3, or alpha4 subunits with the beta2 or beta4 subunit in Xenopus oocytes were fitted best by the sum of two exponentials and a constant between -60 and -150 mV. As the ACh concentration approached zero, the relaxation time constants approached limiting values that should equal the single-channel burst duration at low ACh concentrations and the synaptic current decay time constants. beta4 co-expression prolonged the zero ACh concentration limits for the relaxation time constants. The fast beta4 zero ACh concentration limits ranged from 40 to 121 ms between -60 and -150 mV, and the slow beta4 zero ACh concentration limits ranged from 274 to 1039 ms. In contrast, the fast beta2 limits were 4-6 ms over the same voltage range and the slow beta2 limits were 30-53 ms. Expression with the beta4 subunit increased the voltage sensitivity of the alpha2, alpha3 and slow alpha4 relaxation time constants but not that of the fast alpha4 relaxation time constant. Reducing the temperature from 22 C to 8-9 C increased the alpha4beta2 and alpha3beta4 relaxation time constants 2.3- to 6.6-fold and reduced the fractional amplitude of the fast relaxation component. It also increased the voltage dependence of the fast alpha3beta4 relaxation time constant and decreased that of the slow time constant. The Q10 for alpha4beta2 and alpha3beta4 relaxation time constants ranged from 1.9 to 3.9 between 10 and 20 C. The beta subunit appears to have a dominant influence on the voltage-jump relaxation kinetics of heteromeric neuronal nicotinic receptors.

Selective reduction in the nicotinic acetylcholine receptor and dystroglycan at the postsynaptic apparatus of mdx mouse superior cervical ganglion

Our previous data suggested that in mouse sympathetic superior cervical ganglion (SCG) the dystrophin-dystroglycan complex may be involved in the stabilization of the nicotinic acetylcholine receptor (nAChR) clusters. Here we used SCG of dystrophic mdx mice, which express only the shorter isoforms of dystrophin (Dys), to investigate whether the lack of the full-length dystrophin (Dp427) could affect the localization of the dystroglycan and the alpha3 nAChR subunit (alpha3AChR) at the postsynaptic apparatus. We found a selective reduction in intraganglionic postsynaptic specializations immunopositive for alpha3AChR and for alpha- and beta-dystroglycan compared with the wild-type. Moreover, in mdx mice, unlike the wild-type, the disassembly of intraganglionic synapses induced by postganglionic nerve crush occurred at the slower rate and was not preceded by the loss of immunoreactivity for Dys isoforms, beta-dystroglycan, and alpha3AChR. These data indicate that the absence of Dp427 at the intraganglionic postsynaptic apparatus of mdx mouse SCG interferes with the presence of both dystroglycan and nAChR clusters at these sites and affects the rate of synapse disassembly induced by postganglionic nerve crush. Moreover, they suggest that the decrease in ganglionic nAChR may be one of the factors responsible for autonomic imbalance described in Duchenne muscular dystrophy patients.

Participation of a chloride conductance in the subthreshold behavior of the rat sympathetic neuron

The presence of a novel voltage-dependent chloride current, active in the subthreshold range of membrane potential, was detected in the mature and intact rat sympathetic neuron in vitro by using the two-microelectrode voltage-clamp technique. Hyperpolarizing voltage steps applied to a neuron held at -40/-50 mV elicited inward currents, whose initial magnitude displayed a linear instantaneous current-voltage (I-V) relationship; afterward, the currents decayed exponentially with a single voltage-dependent time constant (63.5 s at -40 mV; 10.8 s at -130 mV). The cell input conductance decreased during the command step with the same time course as the current. On returning to the holding potential, the ensuing outward currents were accompanied by a slow increase in input conductance toward the initial values; the inward charge movement during the transient ON response (a mean of 76 nC in 8 neurons stepped from -50 to -90 mV) was completely balanced by outward charge displacement during the OFF response. The chloride movements accompanying voltage modifications were studied by estimating the chloride equilibrium potential (E(Cl)) at different holding potentials from the reversal of GABA evoked currents. [Cl(-)](i) was strongly affected by membrane potential, and at steady state it was systematically higher than expected from passive ion distribution. The transient current was blocked by substitution of isethionate for chloride and by Cl(-) channel blockers (9AC and DIDS). It proved insensitive to K(+) channel blockers, external Cd(2+), intracellular Ca(2+) chelators [bis-(o-aminophenoxy)-N,N,N',N'-tetraacetic acid (BAPTA)] and reduction of [Na(+)](e). It is concluded that membrane potential shifts elicit a chloride current that reflects readjustment of [Cl(-)](i). The cell input conductance was measured over the -40/-120-mV voltage range, in control medium, and under conditions in which either the chloride or the potassium current was blocked. A mix of chloride, potassium, and leakage conductances was detected at all potentials. The leakage component was voltage independent and constant at approximately 14 nS. Conversely, gCl decreased with hyperpolarization (80 nS at -40 mV, undetectable below -110 mV), whereas gK displayed a maximum at -80 mV (55.3 nS). Thus the ratio gCl/gK continuously varied with membrane polarization (2.72 at -50 mV; 0.33 at -110 mV). These data were forced in a model of the three current components here described, which accurately simulates the behavior observed in the "resting" neuron during membrane migrations in the subthreshold potential range, thereby confirming that active K and Cl conductances contribute to the genesis of membrane potential and possibly to the control of neuronal excitability.

The kinetics of nerve-evoked quantal secretion

Current views on quantal release of neurotransmitters hold that after the vesicle migrates towards release sites (active zones), multiple protein interactions mediate the docking of the vesicle to the presynaptic membrane and the formation of a multimolecular protein complex (the 'fusion machine') which ultimately makes the vesicle competent to release a quantum in response to the action potential. Classical biophysical studies of quantal release have modelled the process by a binomial system where n vesicles (sites) competent for exocytosis release a quantum, with probability p, in response to the action potential. This is likely to be an oversimplified model. Furthermore, statistical and kinetic studies have given results which are difficult to reconcile within this framework. Here, data are presented and discussed which suggest a revision of the biophysical model. Transient silencing of release is shown to occur following the pulse of synchronous transmitter release, which is evoked by the presynaptic action potential. This points to a schema where the vesicle fusion complex assembly is a reversible, stochastic process. Asynchronous exocytosis may occur at several intermediate stages in the process, along paths which may be differentially regulated by divalent cations or other factors. The fusion complex becomes competent for synchronous release (armed vesicles) only at appropriately organized sites. The action potential then triggers (deterministically rather than stochastically) the synchronous discharge of all armed vesicles. The existence of a specific conformation for the fusion complex to be competent for synchronous evoked fusion reconciles statistical and kinetic results during repetitive stimulation and helps explain the specific effects of toxins and genetic manipulation on the synchronization of release in response to an action potential.

A model of signal processing at a mammalian sympathetic neurone

A computational model has been developed for the action potential and, more generally, the electrical behaviour of the rat sympathetic neurone. The neurone is simulated as a complex system in which five voltage-dependent conductances (gNa, gCa, gKV, gA, gKCa), one Ca2+-dependent voltage-independent conductance (gAHP) and the activating synaptic conductance coexist. The individual currents are mathematically described, based on a systematic analysis obtained for the first time in a mature and intact mammalian neurone using two-electrode voltage-clamp experiments. The simulation initiates by setting the starting values of each variable and by evaluating the holding current required to maintain the imposed membrane potential level. It is then possible to simulate current injection to reproduce either the experimental direct stimulation of the neurone or the physiological activation by the synaptic current flow. The subthreshold behaviour and the spiking activity, even during long-lasting current application, can be analysed. At every time step, the program calculates the amplitude of the individual currents and the ensuing changes; it also takes into account the accompanying K+ accumulation process in the perineuronal space and changes in Ca2+ load. It is shown that the computed time course of membrane potential must be filtered, in order to reproduce the limited bandwidth of the recording instruments, if it is to be compared with experimental measurements under current-clamp conditions. The membrane potential trajectory and single current data are written in files readable by graphic software. Finally, a screen image is obtained which displays in separate graphs the membrane potential time course, the synaptic current and the six ionic current flows. The simulated action potentials are comparable to the experimental ones as concerns overshoot amplitude and rising and falling rates. Therefore, this program is potentially helpful in investigating many aspects of neurone behaviour.