JN AJP citation statistics
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 96: 2564-2577, 2006. First published July 26, 2006; doi:10.1152/jn.00446.2006
0022-3077/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
96/5/2564    most recent
00446.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van Drongelen, W.
Right arrow Articles by Ramirez, J.-M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van Drongelen, W.
Right arrow Articles by Ramirez, J.-M.

Role of Persistent Sodium Current in Bursting Activity of Mouse Neocortical Networks In Vitro

Wim van Drongelen1,2, Henner Koch3, Frank P. Elsen1, Hyong C. Lee1, Ana Mrejeru3, Erin Doren1, Charles J. Marcuccilli1, Mark Hereld2,4, Rick L. Stevens2,4 and Jan-Marino Ramirez3

1Department of Pediatrics, 2Computation Institute, and 3Department of Anatomy and Organismal Biology, The University of Chicago, Chicago; and 4Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, Illinois

Submitted 28 April 2006; accepted in final form 20 July 2006


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Most types of electrographic epileptiform activity can be characterized by isolated or repetitive bursts in brain electrical activity. This observation is our motivation to determine mechanisms that underlie bursting behavior of neuronal networks. Here we show that the persistent sodium (NaP) current in mouse neocortical slices is associated with cellular bursting and our data suggest that these cells are capable of driving networks into a bursting state. This conclusion is supported by the following observations. 1) Both low concentrations of tetrodotoxin (TTX) and riluzole reduce and eventually stop network bursting while they simultaneously abolish intrinsic bursting properties and sensitivity levels to electrical stimulation in individual intrinsically bursting cells. 2) The sensitivity levels of regular spiking neurons are not significantly affected by riluzole or TTX at the termination of network bursting. 3) Propagation of cellular bursting in a neuronal network depended on excitatory connectivity and disappeared on bath application of CNQX (20 µM) + CPP (10 µM). 4) Voltage-clamp measurements show that riluzole (20 µM) and very low concentrations of TTX (50 nM) attenuate NaP currents in the neural membrane within a 1-min interval after bath application of the drug. 5) Recordings of synaptic activity demonstrate that riluzole at this concentration does not affect synaptic properties. 6) Simulations with a neocortical network model including different types of pyramidal cells, inhibitory interneurons, neurons with and without NaP currents, and recurrent excitation confirm the essence of our experimental observations that NaP conductance can be a critical factor sustaining slow population bursting.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Epilepsy is a common neurological disease characterized by chronic seizures, affecting nearly 1% of the world population. Seizure activity in the neocortex plays a major role in children with intractable epilepsy and a better understanding of neocortical physiology is critical for the development of a rational approach to anticonvulsant therapy. One of the major characteristics of most types of epileptiform behavior is a rhythmic pattern of electrical discharges of the neuronal population. In principle these discharges may be caused by intrinsic neuronal mechanisms, network function, or—most likely—both network and cellular properties. Increasing evidence suggests that epileptiform activity originates within the neocortex (Timofeev and Steriade 2004Go). When considering network bursting and different cell types known to occur in neocortex, the intrinsically bursting cell type (Connors and Gutnick 1990Go; Gray and McCormick 1996Go; Steriade et al. 1998Go) seems especially interesting. Intrinsically bursting neurons, in particular the so-called fast rapid bursting neurons, have been hypothesized to play an important role in epileptogenesis and various forms of neocortical activity (Blatow et al. 2003Go; Cunningham et al. 2004Go; Timofeev and Steriade 2004Go).

Several studies indicated that intrinsic properties such as the persistent sodium currents may play a role in bursting and oscillatory behavior (e.g., Mantegazza et al. 1998Go; Timofeev et al. 2000aGo, 2004Go; Traub et al. 2003Go). In animal models it was established that intrinsically bursting neurons are increased in seizing tissue (Jacobs et al. 1999Go; Sanabria et al. 2002Go; Steriade et al. 2001Go; Timofeev et al. 2000bGo). Indeed, recent evidence indicates that this cell type occurs in neocortex of pediatric patients with epilepsy (Foehring and Wyler 1990Go; Van Drongelen et al. 2003bGo). The simplest scenario is a unidirectional path from cell to network, where a bursting neuron paces the activity of a neuronal population. This "pacemaker hypothesis" was first raised in 1949, when Bremer postulated that neurons with intrinsic pacemaker or bursting properties contribute to the generation of neocortical EEG activity (Bremer 1949Go). However, a better understanding of how neuronal networks generate rhythmic activity led to a major modification of a pure "pacemaker hypothesis" because several studies indicated that bursting and synaptic properties interact in concert (e.g., Steriade et al. 1998Go, 2001Go; Traub et al. 2001Go, 2003Go). Unfortunately, the close interaction between these properties makes it extremely difficult or even impossible to determine the relative contribution of these properties in functional neuronal networks. The finding that the strength of synaptic and intrinsic bursting properties is not fixed but dynamically modulated by endogenously released neuromodulators or by manipulations such as deafferentation further complicates the question to what extent intrinsic and synaptic properties contribute to rhythm generation (Foehring et al. 2002Go; Peña and Ramirez 2004Go; Topolnik et al. 2003Go).

The difficulty of determining the relative contribution of synaptic and bursting properties to the generation of a given activity pattern was overcome in some small invertebrate neuronal networks that were particularly amenable to a rigorous cellular and systems level analysis (Marder and Calabrese 1996Go). The increased difficulty of dissecting synaptic and bursting properties in mammalian networks created considerable uncertainty as to how intrinsic and synaptic properties interact to generate rhythmic activity. This is an issue of general interest because bursting neurons exist almost everywhere in the mammalian nervous system (Ramirez et al. 2004Go) including the thalamus, basal ganglia, spinal cord (Darbon et al. 2004Go), medulla (Peña and Ramirez 2004Go), hypothalamus, olfactory bulb, ventral tegmentum, hippocampus (Sanabria et al. 2001Go), and neocortex (Connors and Gutnick 1990Go; Dégenètais et al. 2003Go; Istvan and Zarzecki 1994Go; Nowak et al. 2003Go). Although some researchers hypothesize that these neurons are primarily responsible for the generation of rhythmic activity in areas where bursting neurons are found (for a review see Arshavsky 2003Go), other studies claim that bursting neurons or pacemaker activity plays no obligatory role in the generation of a given rhythmic activity. Prominent examples are the generation of the respiratory rhythm (Del Negro et al. 2005Go), the slow neocortical oscillations that might underlie sleep (Sanchez-Vives and McCormick 2000Go), and the slow oscillations generated in the hippocampus (Staley et al. 2001Go), all of which are believed to be primarily the result of emergent network (i.e., synaptic) properties.

In the present study we examined the role of bursting neurons in neocortical networks by selectively blocking persistent sodium (NaP) currents of pyramidal neurons in brain slices obtained from the neocortex of neonatal mice (P8–P13), while observing whether this blockade affects the generation of spontaneously generated slow network oscillations (0.1–1 Hz) described previously by Van Drongelen et al. (2003a)Go. The spontaneously active isolated neocortical network enabled us to pharmacologically manipulate cellular and network properties, a process that is more difficult under in vivo conditions. Our data suggest that blockade of NaP currents with low concentrations of tetrodotoxin (TTX) or riluzole simultaneously reduces bursting activity of intrinsically bursting pyramidal neurons and spontaneous bursting of neocortical networks. Riluzole is known to have a variety of effects in a variety of systems at higher concentrations (O’Neill et al. 1997Go; Peña and Tapia 2000Go; Siniscalchi et al. 1997Go; Stefani et al. 1997Go; Zona et al. 1998Go). For example, riluzole causes a rapid and sustained plateau increase in Ca2+ with an EC50 of 150 µM in human osteosarcoma cells (Jan et al. 2002Go). At 40–50 µM riluzole accelerates the activation kinetics of calcium channels without affecting voltage dependency of activation (Huang et al. 1997Go). However, riluzole seems to selectively block the voltage-gated NaP current at concentrations of 0.5–20 µM (Del Negro et al. 2002Go; Spadoni et al. 2002Go; Urbani and Belluzi 2000Go). In voltage-clamp studies we showed that low concentrations of TTX or riluzole reduce the NaP current. Further, in a combined experimental and modeling approach, we confirmed that reduction of the NaP current is associated with a block of cellular and population bursting activity. Overall our findings support the hypothesis that the NaP current and intrinsically bursting neurons play a critical role in the generation of rhythmic bursting behavior of neocortical networks.


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Neocortical slices

Neonatal (P8–P13) male or female CD-1 mice (n = 105) were deeply anesthetized and decapitated at the C3/C4 spinal level and the forebrain was isolated in ice-cold artificial cerebrospinal fluid (ACSF). One hemisphere was then glued onto an agar block with its rostral end up and mounted onto a vibrating tissue slicer. Coronal slices (500 µm for current-clamp and 350 µm for voltage-clamp experiments) of the somatosensory cortex were transferred into a recording chamber and submerged under a stream of ACSF (temperature, 30°C; flow rate, 10 ml/min) containing (in mM) 118 NaCl, 3 KCl, 1.5 CaCl2, 1 MgSO4, 25 NaHCO3, 1 NaH2PO4, and 30 D-glucose equilibrated with carbogen (95% O2-5% CO2). In all current-clamp and extracellular experiments, the potassium concentration was routinely raised from 3 to 5 mM over 30 min to obtain spontaneous rhythmic activity close to 0.2 Hz.

Extracellular and intracellular current-clamp recording

To obtain a signal containing multiunit action potential (AP) activity, extracellular signals were amplified 10,000-fold and filtered between 0.25 and 1.5 kHz (Fig. 1B, bottom trace). To facilitate detection of bursts, this signal was rectified and integrated by using an electronic integrator with a time constant of 50 ms. By applying this procedure we obtained an index for multiunit AP activity in which bursts were easily distinguished (Fig. 1B, top trace). Both the extracellular recording and the integration were performed with home-built circuitry.


Figure 1
View larger version (13K):
[in this window]
[in a new window]
 
FIG. 1. Recording techniques. A: extracellular and intracellular electrodes are placed in close proximity in the neocortical slice. B: population activity is recorded extracellularly, and a rectifier plus integrator is used to emphasize burst activity of the network. C: current-clamp measurement shows the activity pattern of individual pyramidal neurons associated with the population bursts. D: voltage-clamp measurement was applied to evaluate sodium currents in the cell membrane and to measure postsynaptic current.

 
Intracellular current-clamp recordings were obtained from cortical neurons with the blind-patch technique. The patch electrodes were manufactured from filamented borosilicate glass tubes (Clarke GC 150TF), filled with a solution containing (in mM/l): 140 K-gluconic acid, 1 CaCl2 · 6H2O, 10 EGTA, 2 MgCl2 · 6H2O, 4 Na2ATP, and 10 HEPES. The intracellular pipettes contained biocytin (4.5 mg/ml) to allow for identification of neuron location and morphology. Electrodes with a positive pressure of 35–50 mmHg were penetrated deep into the slice in 10-µm steps by using a piezo-driven micromanipulator (Böhm, Neustadt, Germany). As the electrode approaches the cell, the measured electrode resistance increased. To obtain a Giga seal, we removed the positive pressure and applied negative suction. Recordings were low-pass filtered (0–2 kHz, Bessel four-pole filter, –3 dB) (Fig. 1C). Because the extracellular and intracellular electrodes were in close proximity to each other (within 150 µm), the extracellular signal was used as an index for the activity of the surrounding population (network) of neurons. After recording, the slices were placed in a paraformaldehyde solution for subsequent staining procedures.

Voltage-clamp recordings

Pyramidal neurons in layer 5 of the cortex were visually identified (Zeiss Axioskop 2 FS microscope with IR-DIC). Whole cell patch-clamp recordings (Fig. 1D) were obtained with a sample frequency of 10 kHz and a low-pass filter setting of 2 kHz. Recordings were made with unpolished patch electrodes, manufactured from borosilicate glass pipettes with filament (Warner Instruments G150F-4). The electrodes had a resistance of 3–5 M{Omega} when filled with the whole cell patch-clamp pipette solution containing (in mM) 110 CsCl, 30 TEA-Cl, 1 CaCl2, 10 EGTA, 2 MgCl2, 4 Na2ATP, and 10 HEPES (pH 7.2). The patch-clamp experiments were performed with a patch-clamp amplifier (AxoPatch 200B), a digitizing interface (Digidata 1322A), and the software program pClamp 9.2 (Axon Instruments).

Neurons located at least three to four cell layers (about 80–150 µm) caudal from the rostral surface of the slice were recorded under visual control. Neurons located directly at the slice surface were not examined because they were more likely to be severely damaged during the preparation than were neurons located deeper within the slice. Current-response traces were recorded with either off- or on-line leak subtraction (P/4 protocol), eliminating the linear leak current and residual capacity currents. The 2-mV liquid junction potential was manually subtracted with the amplifier’s pipette offset regulator immediately before establishing the patch-clamp configuration. The series resistance was always 80% compensated and regularly corrected throughout the experiments. We emphasize that whole cell voltage-clamp recordings from neurons embedded in a functional network are accompanied by difficult space-clamp control. This could lead to incorrect values for current amplitudes. Thus recordings with obvious space-clamp problems (Armstrong and Gilly 1992Go) were discarded. Poor space clamping was indicated by rebound spikes (rapid, fast inactivating inward currents, which were induced by steps from depolarizing test potentials to the former holding potential) or an increase in the delay to onset of an inward current with increasing magnitude of test pulse. Steps to higher test potentials were typically associated with a reduction in delay to current onset. We also discarded neurons with insufficiently blocked K+ currents, evident in outward currents typically commencing at voltage steps to 10 mV.

Using the conventional patch-clamp technique, we pharmacologically isolated both the voltage-activated calcium and the voltage-activated sodium currents by intracellular blockade of voltage-activated potassium currents with 110 mM CsCl and 30 mM tetraethylammonium (TEA) chloride. As shown in a previous study (Elsen and Ramirez 1998Go), this blockade was sufficient for the investigation of the maximal inward current amplitude. In addition, voltage-activated calcium currents were blocked by extracellular bath application of 200 µM cadmium chloride (CdCl2) to isolate the NaP current.

Synaptic input measurements

To examine potential synaptic effects of riluzole at a concentration of 20 µM, we eliminated AP-evoked synaptic events and added 1 µM TTX to the ACSF. Miniature excitatory and inhibitory postsynaptic currents (mEPSCs and mIPSCs, respectively) were recorded and digitized for a period of 4 min per data file. Signals were sampled with a frequency of 10 kHz (filtered with 2 kHz) at a holding potential (Vh) of –60 mV and stored for subsequent analysis. The recording electrodes had a resistance of 3–5 M{Omega} when filled with the whole cell patch-clamp pipette solution containing (in mM) 138 Cs-Me-sulfonate, 2 NaCl, 1 CaCl2, 10 EGTA, 4 Na2ATP, 0.3 Na-GTP, and 10 HEPES (pH 7.2). Consulting equilibrium potential to holding potential relationships, all upward deflections had to be chloride-conducting, inhibitory currents, whereas all downward deflection had to be mEPSCs.

To ensure steady drug (riluzole) concentration equilibrium in the recording chamber, a 4-min preapplication of riluzole-containing solution was carried out before the start of the actual 4-min data recording that was used for the later analysis.

Experiments

Rhythmically active slices were obtained 30 min after increasing the level of potassium to 5 mM. This preparation allowed us to study both network activity and cellular activity (Fig. 1, B and C). To investigate the role of the NaP current, in the generation of network bursting and cellular activity, we applied either riluzole (20 µM) or a very low concentration TTX (50 nM) to the ACSF. Typically the network bursting decreased and stopped after attenuating the NaP current (Fig. 2). In some of the experiments, the neurons were decoupled from the network by blocking non–N-methyl-D-aspartate (NMDA) and NMDA glutamatergic receptors using 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX, 20 µM) and (±)-3-(2-carboxypiperazin-4-yl)-propyl-1-phosphonic acid (CPP, 10 µM); adding this combination of drugs in the ACSF stopped network activity within a few minutes.


Figure 2
View larger version (9K):
[in this window]
[in a new window]
 
FIG. 2. Integrated network activity in a slice where riluzole was added to the artificial cerebrospinal fluid (ACSF).

 
To explore the effectiveness of riluzole and a low concentration of TTX and the role of persistent sodium in cellular and network activity, we carried out the following experiments.
  1. )Riluzole (n = 13) or TTX (n = 8) was added to the ACSF and network bursting activity was observed and quantified for 10 min after drug application.
  2. )The neuronal response to current injection (50–150 pA) was determined before and after bath application of riluzole (n = 12) or low-concentration TTX (n = 8) and every 30 s thereafter for 10 min. Responses of regular spiking (RS)–type neurons were evaluated only during absence of overt network activity. The intrinsically bursting (IB)–type neurons were isolated from the network to confirm that the bursting was intrinsic and not network driven. The response evoked by the current injection was quantified by the number of action potentials in 1 s after onset of injection.
  3. )By means of the voltage-clamp technique, the persistent sodium conductivity before and after addition of riluzole (n = 5) or a low concentration of TTX (n = 7) was determined in pyramidal neurons after they were isolated from the network.
  4. )In five experiments we also determined the effect of riluzole on calcium currents in the presence of TTX (1 µM).
  5. )The effect of riluzole on postsynaptic activity was recorded in voltage-clamp mode (n = 9).

The extracellular and intracellular current-clamp recordings were performed in a different setup from that for the voltage-clamp measurements. Because of differences in the tubing of each setup, there was a 60-s difference in latency between drug delivery to the ACSF system and arrival in the bath. The drug delivery difference is compensated for in the graphs of GoGoFig. 5 by shifting the graphs relative to each other.


Figure 3
View larger version (20K):
[in this window]
[in a new window]
 
FIG. 3. Diagram of the computational model of neocortex. Superficial pyramidal cells (S) and deep pyramidal cells (D) are located in different layers. Interneurons (I) are inhibitory cells located in a 3rd layer. A: excitatory connections between pyramidal cells and between pyramidal cells and interneurons. Interneurons of each type are interconnected with gap junctions (R). Parameters associated with the excitatory connections indicated by 1–14 can be found in Table A1 in the APPENDIX. B: inhibitory connectivity between basket cells (B) and the soma of the pyramidal neurons and between chandelier cells (C) and the initial segment of the pyramidal neurons. Inhibitory connections indicated by a–h and disinhibition I–IX correspond to Table A1 in the APPENDIX.

 

Figure 4
View larger version (41K):
[in this window]
[in a new window]
 
FIG. 4. Simulation results for single-cell activity showing the patterns evoked by current injections in regular spiking (RS) and intrinsically bursting (IB) types of the cell models included in the network simulation. Top trace in each panel shows the action potential activity; the bottom trace indicates onset and level of current injected in the cell soma. To show similar responses, the level of current injection was adapted for each cell type. AD: responses of superficial and deep pyramidal neurons. EH: activity from 3 types of basket cell and the axo-axonic interneuron (chandelier cell). Except for the bursting cells (A and C) that included a persistent sodium (NaP) current, all neuron types had regular spiking characteristics. Compared with the pyramidal neurons, the smaller basket cells (EG) and the chandelier cell (H) spiked at a lower threshold of current. Trace in I represents 15 s of simulated network activity with NaP current of 22.5 S/m2. Network activity was presented as the integrated power of the simulated extracellular signal. In this signal, bursts were detected above a threshold of 0.01 arbitrary units (AUs).

 

Figure 5
View larger version (59K):
[in this window]
[in a new window]
 
FIG. 5. A: effect of riluzole (20 µM) on network bursting frequency. Extracellular network activity was recorded for ≤10 min after drug application and the number of bursts was recorded in 50-s intervals. B: effect of riluzole (20 µM) on the action potential firing rate of RS (n = 8, open triangle) and IB neurons (n = 4, open square). Rates were determined for 30-s intervals. C: effect of riluzole (20 µM) on the amplitude of NaP current. Maximal NaP amplitude was determined every 10 s by a single voltage step from –80 to –20 mV; 100% was determined by the average from the last 3 measurements before onset of riluzole application. D: effect of low-concentration tetrodotoxin (TTX, 50 nM) on network bursting frequency. Extracellular network activity was recorded for ≤10 min after drug application and the number of bursts was recorded in 50-s intervals. E: effect of low-concentration TTX (50 nM) on the action potential firing rate of RS (n = 4, open triangle) and IB neurons (n = 4, open square). Rates were determined for 30-s intervals. F: effect of low-concentration TTX (50 nM) on the amplitude of NaP current. Maximal NaP amplitude was determined every 10 s by a single voltage step from –80 to –20 mV; 100% was determined by the average from the last 3 measurements before onset of TTX application. Note: Because cellular + network activity was recorded with a different system from the voltage-clamp data, there is a 60-s difference in latency between drug administration and arrival in the bath. For this reason there is a 60-s shift of the graphs in C and F relative to the others.

 
Data analysis

Data were analyzed by using MiniAnalysis 5.41 (Synaptosoft) and statistical analysis was performed with Prism 4.03 (Graphpad). Significance between data values was assessed with the Student’s t-test and assumed when P ≤ 0.05. Voltage-clamp data were analyzed off-line with the software program ClampFit 9.2 (Axon Instruments). Quantitative data are given as means ± SE.

A computerized preanalysis of each 4-min recording was performed to tag all events that qualified as miniature excitatory postsynaptic currents (mEPSCs) according to the user-defined, program-specific detection parameters. However, an additional manual user-performed screen through each data trace was necessary to eliminate computer errors. During this eyeballed-screening process, wrongly tagged artificial events were unmarked and missed real mEPSC events were tagged. To determine the representative amplitude and decay time of the recorded mEPSCs in each data trace, an average mEPSC from 50 single mEPSCs was created as follows: For each 4-min data trace, the computer randomly picked 75 single events. Those events were superimposed and aligned at 50% rise time. Single events that did not line up within a rise-time window of 4 ms were manually eliminated. In addition, the decay phase of each event had to steadily decay without any interruption by additional events or noise for 50 ms. This procedure was followed until 50 events were included into the average. Computerized single-exponential fit through the decay phase (10–90%) of the average mEPSC (GoGoGoFig. 9B) delivered a representative amplitude and decay time constant ({tau}). This analysis was not performed on the inhibitory component because of a low frequency of occurrence of mIPSCs before and after adding riluzole.


Figure 6
View larger version (20K):
[in this window]
[in a new window]
 
FIG. 6. Neuronal activity of the population (top traces in each panel) usually showed a relationship with the individual neuron (bottom traces). In RS neurons (A panels) and IB neurons (B panels), we found 3 types of activity associated with the population burst: excitation (A1, B1), a combination of excitation followed by a phase of inhibition (A2, B2), or inhibition (A3, B3). A subclass of bursting cells remained active after the network effect was removed by bath application of CPP and CNQX.

 

Figure 7
View larger version (23K):
[in this window]
[in a new window]
 
FIG. 7. Responses of representative RS (A1) and IB (A2, A3) neurons to current injection (100 pA, bottom traces) before and after riluzole (A1, A2) or a low-concentration TTX (A3) was added to the ACSF. Application of riluzole or TTX inhibited the intrinsic bursting properties of the IB-type cells; firing of the RS was not significantly altered. B: instantaneous spike frequency plotted against time elapsed since stimulus onset for the neurons is shown in A. Lines represent the frequency before adding riluzole, the dotted lines after adding riluzole measured immediately after the network stopped bursting. At the cessation of network bursting, the IB neurons exhibited regular spiking activity like the RS-type cell. Losing the bursting property was associated with a reduction of action potentials per second. Histograms in C show the overall change in response expressed in the number of action potentials in the first second after onset of current injection of (C1) RS (n = 8) and (C2) IB neurons (n = 4) before adding riluzole (CTRL) and the effect of riluzole after the network stopped bursting (RIL). Effects for low concentrations of TTX on RS (n = 4) and IB (n = 4) type cells were identical. For the RS neurons, no significant difference was observed; the difference for the IB type neuron was significant (P < 0.05).

 

Figure 8
View larger version (40K):
[in this window]
[in a new window]
 
FIG. 8. A: voltage-activated sodium current response traces to different test potentials (inset: holding potential: –80 mV, 10-mV steps from –90 to 20 mV, step duration: 150 ms). We measured the persistent sodium (NaP) current amplitudes at the steady-state level as indicated by the filled circle. Figure shows cortical pyramidal neuron with attached patch pipette. B: current–voltage relationship diagram (IV curve) of persistent sodium (NaP) current (n = 43) and NaP after bath application of 20 µM riluzole. Average current amplitude values are plotted against the respective test potentials. C: normalized current amplitudes of NaP, NaP + riluzole (20 µM, n = 5), and NaP + TTX (50 nM, n = 7). All absolute values were normalized to the respective maximal NaP current amplitude (I/Imax). D: inactivation test protocol: example of NaP current response traces to a test potential of –15 mV (peak voltage range) after different holding potentials (10 s) from –20 to –120 mV (10-mV steps). Still activatable NaP current amplitudes were measured at the steady-state level as indicated by the open circle. E: average activation and inactivation curves for NaP. Activation parameters (filled circles) are plotted as normalized conductance values (right y-axis) against the respective test potentials (Vt). Inactivation parameters (open circles) are plotted as normalized current amplitudes from 28 NaP inactivation measurements (left y-axis) against the respective Vt. Sigmoidal Boltzmann fits through both sets of data points deliver the activation parameters V50 = –46.2 ± 1.4 mV and slope = 9.9 ± 1.3) and inactivation parameters V50 = –65.4 ± 2.3 mV and slope = –19.0 ± 2.5 for NaP. F: activation parameter comparison between NaP, NaP + riluzole (20 µM), and NaP + TTX (50 nM). Sigmoidal Boltzmann fits through the data sets deliver the following activation parameters: NaP + riluzole: V50 = –49.6 ± 2.8 mV and slope = 4.8 ± 2.8; and NaP + TTX: V50 = –48.6 ± 3.0 mV and slope = 5.6 ± 4.1.

 

Figure 9
View larger version (25K):
[in this window]
[in a new window]
 
FIG. 9. Synaptic activity before and after adding riluzole (20 µM). First the action potential generation mechanism was disabled with TTX (1 µM), and the control synaptic activity was measured; then the effect of riluzole was determined. B: example miniature excitatory postsynaptic current traces (mEPSCs) during control conditions (TTX) and after additional application of riluzole (RIL). Traces presented here are examples of average current traces from 50 single events of a 4-min recording during the respective conditions. C: comparison of the decay time constants ({tau}) and amplitudes of mEPSCs during control (TTX) and riluzole (RIL) conditions. Time constants and amplitudes were determined by a single exponential fit through the decay phase of the average mEPSCs (see METHODS). Changes for {tau} and amplitude were not significantly different.

 
Model

The purpose of the computational approach in this study was to explain the relationship between persistent sodium conductivity and emergent network bursting behavior. The neural network in the computational model representing neocortex includes excitatory and inhibitory cell populations with a multicompartmental representation of each cell type. The network connectivity is based on histological and physiological work (DeFelipe et al. 2002Go; Feldmeyer and Sakman 2000Go; Hellwig 2000Go; Krimer and Goldman-Rakic 2001Go; Mountcastle 1997Go; Nieuwenhuys 1994Go) and includes the essence of neocortical microcircuitry.

  1. )A large fraction of excitatory cell types relative to inhibitory interneurons (a ratio of approximately 4:1).
  2. )Excitatory connectivity ending on the dendritic portion, whereas important components of the inhibition are located on the initial segment, soma, and dendrite (DeFelipe et al. 2002Go).
  3. )Recurrent excitation and reciprocal connection between inhibitory cells. A unique neocortical feature is the presence of axo-axonic interneurons (the chandelier cells) that do not make synaptic interconnections between themselves (DeFelipe et al. 2002Go; Somogyi et al. 1998Go).
  4. )Direct electrical contact between inhibitory interneurons (Amitai et al. 2002Go).

The spatiotemporal relationship between the neuronal elements is determined by the geometry of the network, that is, realistic cell dimensions and spacing, associated with realistic conduction speeds. Details of the model were published previously (Van Drongelen et al. 2004Go, 2005Go); further parameters are summarized in the APPENDIX.

In short, the excitatory network (Fig. 3A) consists of superficial pyramidal cells from layers 2/3 (five compartments, Table 1) and deep pyramidal cells from layers 5/6 (seven compartments, Table 1). The inhibitory cells (I, Fig. 3A) receive input from both types of pyramidal neurons. Gap junctions (R, Fig. 3A) between inhibitory cells show nearest-neighbor connectivity (e.g., Amitai et al. 2002Go). The implemented inhibitory interneurons are three types of basket cells and the chandelier cell (B, C, respectively, in Fig. 3B). The basket cell types inhibit the pyramidal cell soma, whereas the chandelier cell directly inhibits the initial segment (Fig. 3B). The neurons are placed in a three-dimensional (3D) grid and connected to each other with connectivity rules derived from connections in the mammalian cortex. In a short range, the synaptic connectivity decreases with distance between source and target elements (e.g., Feldmeyer and Sakman 2000Go; Hellwig 2000Go; Krimer and Goldman-Rakic 2001Go; Nieuwenhuys 1994Go). Excitatory synaptic connections were simulated by an alpha function with a time constant of 1–3 ms for excitatory connectivity; inhibitory synaptic activity was represented by a dual-exponential function with time constants of 1 and 7 ms.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Overview of cell compartments, compartment size, associated voltage-sensitive ion channels, synaptic channels, and gap junctions

 
For fast sodium and potassium currents, implemented on the soma and initial segment of the model cells, we used conductivity parameters from Hodgkin and Huxley (1952)Go; these values are justified because the parameters reported in several modeling studies of cortical neurons differ only slightly with respect to activation and inactivation coefficients and associated time constants (Destexhe and Paré 1999Go; Fleidervish and Gutnick 1996Go; Golomb and Amitai 1997Go; Traub et al. 1994Go). Dynamics of the persistent sodium channels located on the soma of a subset of the superficial pyramidal cells reflect the measurements reported in this study. The computational model included bursting and RS-type pyramidal neurons (Fig. 4, AD) and interneurons that fired at lower thresholds (Fig. 4, EH). The bursting properties of individual neurons were quantified by the number of action potentials (nAP) in the burst and the burst frequency (BF). A compound cellular bursting parameter (CB) reflecting burst intensity was calculated as CB = nAP x BF.

The computational model is implemented in the parallel GENESIS simulator (Bower and Beeman 1998Go) to facilitate a high degree of detail, large-scale simulations, and multiple parameter searches. A scalable version of the model runs on the Jazz computing cluster (Argonne National Laboratory), which enabled us to perform large sets of parameter searches within a reasonable time period. [For those interested, the current version of the pGENESIS scripts are available; please contact the first author for details.] Neuronal parameters were determined with a brute-force parameter search (see the APPENDIX) in which intrinsic and evoked cellular behaviors were used as a target. A subsequent brute-force parameter search was used to determine network behavior. The compartmental and cellular activity consists of complete knowledge of the local membrane potential and current. The extracellular activity was obtained as a weighted sum of currents generated by the model neurons (Nunez 1981Go). Where needed, network bursting (NB) was quantified by using the procedure described below. To obtain an estimate of the power, a 150-s epoch of the extracellular signal was first squared; this signal was then integrated with the same filter used for the recorded network activity (a filter with a 50-ms time constant). In these integrated traces, bursts were clearly visible and a threshold separating ongoing activity and bursts was visually determined (Fig. 4I). Finally the total area above the threshold served as an index for the bursting power in the trace. Although selection of the threshold was not critical for the qualitative result (reduced bursting on reduction of NaP current), a threshold at 0.01 arbitrary units (AUs; Fig. 4I) for all traces was used to generate the NB burst index shown in Fig. 10, A and B.


Figure 10
View larger version (27K):
[in this window]
[in a new window]
 
FIG. 10. A: network bursting (NB, line) and cellular bursting (CB, dashed line) as a function of the maximum conductivity of persistent sodium (GNaP). B: linear regression between NB and CB was characterized by the line NB = 154.2CB + 53.9; correlation coefficient = 0.89 (P < 0.01). C: network activity in the computational network. Traces represent the raw network activity (top trace) and integrated activity (bottom trace), the latter of which is the simulated equivalent of the recording in Fig. 2. Computational results show that bursting activity of the network stops after reducing the NaP conductivity in the time interval indicated by the horizontal bar.

 

 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Network activity

To obtain stable baseline bursting activity of 0.23 ± 0.08 Hz, 5 mM of [K+] was used (Van Drongelen et al. 2003aGo). Network activity expressed in bursts per second decreased drastically after adding 20 µM riluzole (n = 13) or 50 nM TTX (n = 5) to the bath (Figs. 2 and 5, A and D). In both cases there was a 100-s latency between adding the drug and a reduction in the burst activity level. All network bursting activity in the slice disappeared 280.0 ± 49.0 s after administration of 50 nM TTX and 403.8 ± 40.2 s after 20 µM riluzole was added to the ACSF. In both TTX and riluzole application, the overall reduction in network bursting was significantly reduced for latencies ≥150 s (P < 0.05).

Cellular activity

Simultaneous extracellular recordings from population bursts and intracellular recordings from individual neurons revealed that all neurons were either excited or inhibited in phase with population activity; that is, neurons received synaptic inputs from the network. Typical patterns (Van Drongelen et al. 2003aGo) were excitation, excitation followed by inhibition, and inhibition as shown for RS-type cells in Fig. 6A and for IB-type neurons in Fig. 6B. The membrane potentials were 61.0 ± 2.7 and 66.4 ± 5.1 mV for the RS- and IB-type neurons, respectively.

Attenuation of the NaP current affects the intrinsic excitability of neocortical pyramidal neurons; however, the effect depended on cell type (Fig. 7). Examples of the effect of riluzole on RS and IB cells are shown in Fig. 7, A1 and A2. When network bursting stopped, the RS-type neurons were not significantly affected by riluzole (n = 8, P = 0.26) (Fig. 7C1); the average number of action potentials per second (AP/s) was 10.7 ± 0.8 before and 9.4 ± 1.3 after riluzole application. However, 10 min after network cessation, the neural response evoked by current injection into RS neurons became significantly reduced (P < 0.05, 3.0 ± 1.2 AP/s). By contrast, bursting activity in IB neurons completely disappeared and became regular spiking within 5–6 min after riluzole application (Fig. 7, A2 and B2). As a consequence of the loss of bursts the number of APs generated by the current injection was significantly decreased (n = 4: P < 0.05, AP/s 9.2 ± 1.0 before adding riluzole and 4.4 ± 0.8 after the application of riluzole, Fig. 7C2). The difference between activity levels of the RS- and IB-type neurons under control conditions was not significant. The time course of the activity levels of both RS- and IB-type neurons in the presence of riluzole is depicted in Fig. 5B.

To address the possibility that these results were caused by the known unspecific effects of riluzole, we repeated these experiments with TTX at low concentrations. The response in RS-type neurons (n = 4) was not significantly affected (P = 0.55) within 5–6 min of bath-applied TTX (Fig. 7C1). By contrast, the bursting in IB-type neurons was blocked within 5–6 min (Fig. 7, A3 and B3), significantly decreasing the frequency of APs (n = 4: P < 0.01, AP/s 11.1 ± 2.8 before adding TTX and 4.7 ± 2.2 after the network stopped bursting, Fig. 7C2). The time course of the activity levels of both RS- and IB-type neurons is depicted in Fig. 5E.

Voltage-clamp data

We recorded in the voltage-clamp mode from 51 neurons in layer 5 of the neocortex. The average cell capacity was 107 ± 14 pF and cell resistance was 365 ± 92 M{Omega}. After blockade of voltage-activated calcium (bath: 200 µM CdCl2) and potassium (intracellular: 110 mM CsCl, 30 mM TEA-Cl) currents, we evoked voltage-activated NaP currents with specific voltage-step protocols (Fig. 8, A and D, insets). The NaP amplitudes were measured at the steady-state level toward the end of the voltage steps (bullet in Fig. 8A) and were plotted against the respective test potentials (average IV curve, Fig. 8B; normalized curves, Fig. 8C). The peak voltage of the NaP currents was found at –10 mV, with an amplitude of 333 ± 34 pA.

To determine the activation properties, we calculated the NaP conductance (g) to compensate for the changing driving force at respective test potentials [g = I/(VtEq), where Vt is test potential and Eq is equilibrium potential]. The NaP conductance values were normalized and the average was plotted against the respective test potentials (Fig. 8F). A sigmoidal Boltzmann-curve fit 1/[(1 + e–(V–V50)/slope] for the normalized NaP amplitudes as a function of the test potential V delivered the activation parameters: V50 = –46.2 ± 1.4 mV and slope = 9.9 ± 1.3. Both low TTX and riluzole significantly reduced the NaP current (Fig. 8B) and resulted in a lower level of activation (Fig. 8F).

To determine the inactivation properties of the NaP currents, we measured the NaP current amplitudes at the peak voltage (–15 mV) immediately after different prolonged (10-s) test potentials (Fig. 8D). The remaining, activatable NaP amplitudes were normalized (I/Imax) and the average was plotted against the respective test potentials (Fig. 8E). A sigmoidal Boltzmann-fit delivered the inactivation parameters: V50 = –65.4 ± 2.3 mV and slope = –19.0 ± 2.5. No significant effects of riluzole on Ca2+ currents were measured (n = 5).

Synaptic properties

To investigate possible synaptic effects of riluzole, we compared miniature excitatory postsynaptic currents (mEPSCs) under control conditions [with TTX (1 µM) to disable spike activity] with mEPSCs after additional application of 20 µM riluzole (Fig. 9). The difference of the frequency of mEPSCs, 0.9 ± 0.5 Hz under control conditions and 1.3 ± 0.6 Hz after adding riluzole, was not significant. For the mIPSCs the frequency was very low both under control conditions (0.05 Hz) and after adding riluzole (0.04 Hz). In five cells using 50 single events of a 4-min recording for each neuron, we found that neither the mEPSC amplitudes (TTX: 19.8 ± 1.1 pA, TTX + riluzole: 20.4 ± 1.2 pA), nor the mEPSC decay time constants (TTX: 4.6 ± 0.4 ms, TTX + riluzole: 4.2 ± 0.2 ms) were significantly different from each other (Fig. 9, B and C). Because of these low numbers for the inhibitory events, a meaningful averaging procedure and subsequent quantitative analysis was not feasible.

Simulated activity

The purpose of the simulations is to relate the persistent sodium (NaP) current to bursting both at the cellular and network levels; for details of the modeling procedure see METHODS and the APPENDIX. Although we found that the NaP current is not strictly a persistent current because it does inactivate, its voltage-clamp characteristic shows a fairly large window in which these channels can be activated (Fig. 8E). This property indicates that a small, "persistent" sodium current is present during the resting state of the neurons, causing the membrane to slowly depolarize up to the threshold for generating action potentials. This small depolarizing current causes neurons to burst at regular intervals. At small offset depolarizations (or current injections) the neurons burst more frequently, whereas at higher levels of depolarization the neurons show beating (continuous spiking) behavior identical to the activity pattern described by Butera et al. (1999)Go.

The model network that included the bursting pyramidal cell type immediately settled into a pattern of bursting network activity. This activity was used to relate neuronal and network bursting. In Fig. 10A, the burst index CB of the neuronal response to a 220-pA current injection and the burst index of the network NB (see METHODS for the definition of these indices) are both plotted against the maximum conductance of persistent sodium (GNaP). In the absence of sufficient NaP conductance, the cellular response disappeared and spontaneous network bursts were absent as well. The relationship between the cellular and network bursting indices is depicted in Fig. 10B, where the line represents linear regression. The correlation coefficient r = 0.89 between both the cellular and network bursting indices was significant (P < 0.01). The trace shown in Fig. 10C was obtained by reducing the maximum conductance of NaP from 65 to 0.065 S/m2 over a time interval of 12 min. This decrease was done in a logarithmic fashion to mimic the penetration of riluzole or TTX into the slice. By following this procedure we gradually attenuated the effect of NaP current similar to the riluzole or TTX experiments shown in Fig. 5. A typical example of a recorded trace of network activity in a slice is shown in Fig. 2; in this experiment riluzole was used to attenuate bursting activity. Both in the recorded (Fig. 2) and simulated traces (Fig. 10C) network activity initially decreases and stops within 10 min.


 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Our data suggest that the NaP current is important both for the generation bursting activity in a subset of neocortical neurons and for the generation of bursting at the network level. In the presence of low concentrations of TTX or the neuroprotective drug riluzole, the spontaneous network bursting activity in neocortical slices of neonatal mice was suppressed (Fig. 5, A and D). To our knowledge this is the first study in which the role of riluzole is quantified and modeled at the level of

  1. )network activity (Fig. 5)
  2. )single-neuron activity (Figs. 5 and 7)
  3. )intrinsic (persistent sodium) membrane currents (Fig. 8)

Although we emphasize that our modeling was not an exact-fitting procedure to the recorded data, the experimental findings were essentially confirmed in the computational model (Fig. 10).

Our study also confirmed previous experimental and computational studies that indicate that the Nap conductance plays a major role in providing intrinsic excitation to cortical networks (Agrawal et al. 2001Go; Bazhenov et al. 2002Go, 2004Go). Combined in vivo recordings and computational models indicate that the interplay between the hyperpolarization-activated depolarizing current (IH), Ca2+-sensitive potassium current, and a NaP current could organize paroxysmal oscillations (Timofeev et al. 2004Go). The NaP current plays an important role in activating intrinsic bursting in neocortical neurons (Guatteo et al. 1996Go) and a number of studies demonstrated that the number of intrinsically bursting neurons is significantly increased after neocortical trauma or conditions that enhance the propensity of epilepsy (Sanabria et al. 2002Go; Topolnik et al. 2003Go). Steriade and coworkers demonstrated that intrinsically bursting neurons are frequently found in vivo (20–25% of the total cortical neuronal population) and these authors proposed that these bursting neurons play a critical role in the generation of seizure activity (Timofeev and Steriade 2004Go; Timofeev et al. 2000bGo). Indeed, intrinsically bursting neurons are rhythmically active during Lennox–Gastaut-type seizures, specifically during spike-wave (SW) discharges, polyspike-wave (PSW) complexes, and fast runs (Timofeev and Steriade 2004Go). Mechanistically, bursting in these neurons could have a substantial impact on local networks that initiate seizures focally and in related structures as proposed by Steriade and coworkers (Steriade et al. 1998Go; Timofeev and Steriade 2004Go). Such neurons could also play a role in triggering ripples that in turn would facilitate seizure generation, not only in the neocortex, but also in the hippocampus (Bragin et al. 1999aGo,bGo,c, 2002Go; Fisher et al. 1992Go; Grenier et al. 2001Go; Timofeev and Steriade 2004Go; Traub et al. 2001Go). The cellular bursting property together with a modest level of synchrony and recurrent excitation can lead to a seizurelike pattern (Van Drongelen et al. 2003aGo, 2005Go). Our current result, suggesting that intrinsic cellular bursting properties associated with the NaP current are a critical component of synchronous activity in cortical neuronal populations, agrees well with the reported anticonvulsant action of persistent sodium blockers in some patients with epilepsy (Köhling 2002Go).

In addition to the potential role of NaP in epilepsy, the bursting properties associated with this persistent current may very well drive physiological oscillations observed in sleep and memory formation (e.g., Agrawal et al. 2001Go; Bazhenov et al. 2002Go; Timofeev et al. 2000bGo). The rhythmic network bursting as described in the present study could represent an in vitro form of slow sleep oscillation. Such oscillations are known to develop into SW/PSW seizures in vivo (Steriade and Amzica 2003Go; Timofeev and Steriade 2004Go). We emphasize, however, that in a functional network, the persistent sodium channels not only function in form of intrinsic bursting but also could amplify currents flowing through glutamate-activated dendritic channels (Crill 1999Go; Schwindt and Crill 1999Go; Timofeev et al. 2000bGo). The involvement of the NaP current in the generation of intrinsic bursting could also explain the finding that the propensity to burst is a "flexible" property (Steriade 2004Go). The NaP current and bursting neurons are a target of various neuromodulators (Peña and Ramirez 2002Go; Staiger et al. 2002Go).

From a theoretical perspective, network bursting activity can be caused by cellular oscillations arising from intrinsic membrane conductance, network properties, or a combination of both. Current modeling results and earlier simulation studies demonstrate that networks of RS-type pyramidal cells can oscillate at an appropriate strength of connectivity but that intrinsic bursting neuronal activity is critical for sustaining slow bursting network activity (Traub et al. 2003Go; Van Drongelen et al. 2004Go, 2005Go). In our study, low concentrations of TTX or riluzole inhibited the NaP current (Fig. 8), while leaving synaptic function intact (Fig. 9). At the cellular level, we found that reducing the NaP current through low TTX or riluzole can dramatically change the firing properties of IB neurons (Fig. 7). The bursting behavior of these neurons was altered and they converted to a pattern similar to that of regular firing, with a reduced output of action potentials per injected current pulse (Fig. 7C2). Similar experiments with regular firing (RS-type) neurons showed that their firing pattern was unchanged at the cessation of network bursting (Fig. 7, A1, B1, and C1). However, current injections ≥10 min after cessation of network bursting showed that all cell types had a significantly decreased activity level. Our findings suggest that the NaP current plays a role in enhancing excitability of the majority of neocortical neurons but that the intrinsic bursting property is more sensitive to the blockade of NaP than the generation of regular spiking activity. Such a differential sensitivity could explain why intrinsic bursting was blocked first. However, we have not attempted to investigate the reason for the observed differential effect; the focus of the present study was to use this differential effect as a pharmacological tool to determine factors involved in reducing network activity that finally leads to the cessation of network bursting.

The computational approach confirmed the essence of our experimental observation and showed that NaP caused periodic bursting at the cellular level (Fig. 4A), identical to the pattern reported in Butera et al. (1999)Go. At the network level this feature was associated with a bursting pattern, although without the simple periodicity (Fig. 10C). This finding together with the reduced network bursting at attenuated NaP currents suggests that the bursting neurons in a network recruit surrounding neurons at irregular intervals. In this state of the network the neurons are not acting as simple pacemakers that synchronize other units into their intrinsic periodic rhythm; rather, they function merely as units that occasionally can drive or recruit a critical mass of neurons resulting in a burst of activity, whereas the network in its turn codetermines the state of these bursters—clearly a bidirectional process, as similarly discussed by Ramirez et al. (2004)Go. Therefore we propose that the intrinsic (NaP current) and network property (excitatory connectivity) together sustain the rhythmic network behavior that we observed in the neocortical slices.


 APPENDIX
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The overall approach to our modeling effort in this study was to determine acceptable cellular behavior first and then to obtain emergent network behavior. To determine the part of parameter space that would create the desired cellular behavior, we first determined the voltage sensitivity of the sodium and potassium channels by using the literature values originally reported by Hodgkin and Huxley (1952)Go and our recorded data for NaP (Fig. 8). The compartment sizes (Table 1) were based on cellular dimensions and a procedure to collapse cellular parts into single cylindrical compartments (Bush and Sjenowski 1993Go; Van Drongelen et al. 2004Go, 2005Go). Because the model’s surface areas represent an estimate, the maximum levels of conductivity were determined in a brute-force parameter search of 58,644 simulations. For regular and fast spiking cells, the target was to obtain spiking at different levels of current injection (Fig. 4) and no spontaneous activity at a normal resting potential. The cells, including NaP current, were adjusted to show bursting behavior during rest and moderate current injection levels. At higher levels of current injection the behavior transitioned into beating behavior (Butera et al. 1999Go). We evaluated the following parameter ranges; the values we selected for the simulations reported here are indicated in parentheses.

  1. )Membrane resistance: 0.1–5.0 {Omega}m2 (5 and 0.1 for the axon/initial segment)
  2. )Axial intracellular resistance: 0.3–10 {Omega}m (2 and 1 for the axon/initial segment)
  3. )Maximum fast sodium conductance in the soma: 300–3,000 S/m2 (2000)
  4. )Maximum persistent sodium conductance in the soma: 0.3–350 S/m2 (65)
  5. )Maximum potassium conductance in the soma: 6–1,000 S/m2 (650)
  6. )Maximum fast sodium conductance in the initial segment/axon: 30–100,000 S/m2 (6000)
  7. )Maximum potassium conductance in the initial segment/axon: 6–100,000 S/m2 (2000)

Membrane capacitance was 0.01 F/m2. Gap junctions had a resistance of 200 M{Omega}, a value representing five gap junctions with a 1-nS conductivity for each (Traub et al. 2005Go). In spite of the common belief that one may easily simulate any type of activity pattern with the rich parameter set of biophysically realistic models, the large number of parameters is more an obstruction than a facilitator for finding the part of parameter space associated with acceptable behavior; in our example above we found about 0.05% of the simulations generated the desired, realistic behavior. This set was not distributed over the parameter space but encompassed a small area showing desired spontaneous bursting; this area was located around the parameter choice indicated in parentheses in the list above.

As the next step in defining the network, we connected the cells according to the diagrams in Fig. 3. By evaluating a wide range of synaptic strength parameters for the individual synapse model in a brute-force parameter search, we adjusted the synaptic strength to obtain desynchronized network bursting behavior (Van Drongelen et al. 2005Go). The synaptic connectivity parameters based on values reported for mammalian cortex (Van Drongelen et al. 2004Go) are summarized in Table A1. The reference numbers of the different connections in column 1 of Table A1 correspond with those in Fig. 3. The area of connectivity for each connection type was determined by the destination area around the cell and the hole in the destination area that was excluded from making synaptic contacts (columns 24, Table A1). In this destination area, the probability for a synapse to occur declined exponentially with horizontal distance between source (presynaptic) and target (postsynaptic) cells (Probability and Exponential in Table A1). The conduction velocity of the connection and neuronal interdistance determined the delay between spike generation and synaptic effect (column 7, Table A1). Synaptic conductance was determined by a dual exponential function and maximum conductance gmax. In the network gmax was weighed with a factor W, which was determined by the distance between source and target neurons according to the exponential function: W = (max W – min W) exp(d/3D), where max W, min W, and the spatial three-dimensional constant 3D are summarized in the last three columns of Table A1.


View this table:
[in this window]
[in a new window]
 
TABLE A1. Overview of the connectivity parameters in the simulation

 

 GRANTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Part of this work was supported by the Falk Foundation, the Linn family, and by the U.S. Department of Energy under Contract W-31-109-ENG-38.


 ACKNOWLEDGMENTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Dr. K. E. Hecox for support.

Present address of Drs. F. P. Elsen and C. J. Marcuccilli: Department of Neurology, Medical College of Wisconsin, Milwaukee, WI 53226.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: W. van Drongelen, Department of Pediatrics, C372, MC 3055, 5841 S. Maryland Ave., Chicago, IL 60637-1470 (E-mail: wvandron{at}peds.bsd.uchicago.edu)


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Agrawal N, Hamam BN, Magistretti J, Alonso A, and Ragsdale DS. Persistent sodium channel activity mediates subthreshold membrane potential oscillations and low-threshold spikes in rat entorhinal cortex layer V neurons. Neuroscience 102: 53–64, 2001.[CrossRef][Web of Science][Medline]

Amitai Y, Gibson JR, Beierlein M, Patrick SL, Ho AM, Connors BW, and Golomb D. The spatial dimensions of electrically coupled networks of interneurons in the neocortex. J Neurosci 22: 4142–4152, 2002.[Abstract/Free Full Text]

Armstrong CM and Gilly WF. Access resistance and space clamp problems associated with whole cell patch clamping. Methods Enzymol 207: 100–122, 1992.[Web of Science][Medline]

Arshavsky YI. Cellular and network properties in the functioning of the nervous system: from central pattern generators to cognition. Brain Res Brain Res Rev 41: 229–267, 2003.[CrossRef][Medline]

Bazhenov M, Timofeev I, Steriade M, and Sejnowski TJ. Model of thalamocortical slow-wave sleep oscillations and transitions to activated states. J Neurosci 22: 8691–8704, 2002.[Abstract/Free Full Text]

Bazhenov M, Timofeev I, Steriade M, and Sejnowski TJ. Potassium model for slow (2–3 Hz) in vivo neocortical paroxysmal oscillations. J Neurophysiol 92: 1116–1132, 2004.[Abstract/Free Full Text]

Blatow M, Rozov A, Katona I, Hormuzdi SG, Meyer AH, Whittington MA, Caputi A, and Monyer H. A novel network of multipolar bursting interneurons generates theta frequency oscillations in neocortex. Neuron 38: 805–817, 2003.[CrossRef][Web of Science][Medline]

Bower JM and Beeman D. The Book of GENESIS. New York: Springer-Verlag, 1998.

Bragin A, Engel J, Wilson CL, Fried I, and Buzsáki G. High-frequency oscillations in human brain. Hippocampus 9: 137–142, 1999a.[CrossRef][Web of Science][Medline]

Bragin A, Engel J, Wilson CL, Fried I, and Mathern GW. Hippocampal and entorhinal cortex high-frequency oscillations (100–150 Hz) in human epileptic brain and in kainic acid-treated rats with chronic seizures. Epilepsia 40: 127–137, 1999b.[CrossRef][Web of Science][Medline]

Bragin A, Engel J, Wilson CL, Vizentin E, and Mathern GW. Electrophysiologic analysis of a chronic seizure model after unilateral hippocampal KA injection. Epilepsia 40: 1210–1221, 1999c.[CrossRef][Web of Science][Medline]

Bragin A, Modi I, Wilson CL, and Engel J. Local generation of fast ripples in epileptic brain. J Neurosci 22: 2012–2021, 2002.[Abstract/Free Full Text]

Bremer F. Considérations sur l’origine et la nature des "ondes" cérébrales. Electroencephalogr Clin Neurophysiol 1: 177–193, 1949.[Web of Science][Medline]

Bush PC and Sejnowski TJ. Reduced compartmental models of neocortical pyramidal cells. J Neurosci Methods 46: 159–166, 1993.[CrossRef][Web of Science][Medline]

Butera RJ, Rinzel J, and Smith JC. Models of respiratory rhythm generation in the pre-Bötzinger complex. I. Bursting pacemaker neurons. J Neurophysiol 82: 382–397, 1999.[Abstract/Free Full Text]

Connors BW and Gutnick MJ. Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci 13: 365–366, 1990.[Web of Science][Medline]

Crill WE. Functional implications of dendritic voltage-dependent conductances. J Physiol (Paris) 93: 17–21, 1999.

Cunningham MO, Whittington MA, Bibbig A, Roopun A, LeBeau FE, Vogt A, Monyer H, Buhl EH, and Traub RD. A role for fast rhythmic bursting neurons in cortical gamma oscillations in vitro. Proc Natl Acad Sci USA 101: 7152–7157, 2004.[Abstract/Free Full Text]

Darbon P, Yvon C, Legrand JC, and Streit J. INaP underlies intrinsic spiking and rhythm generation in networks of cultured rat spinal cord neurons. Eur J Neurosci 20: 976–988, 2004.[CrossRef][Web of Science][Medline]

DeFelipe J, Alonso-Nanclares L, and Arellano JI. Microstructure of the neocortex: comparative aspects. J Neurocytol 31: 299–316, 2002.[CrossRef][Web of Science][Medline]

Dégenètais E, Thierry A-M, Glowinski J, and Gioanni Y. Electrophysiological properties of pyramidal neurons in the rat prefrontal cortex: an in vivo intracellular recording study. J Neurophysiol 89: 1541–1566, 2003.[Abstract/Free Full Text]

Del Negro CA, Morgado-Valle C, and Feldman JL. Respiratory rhythm: an emergent network property? Neuron 34: 821–830, 2002.[CrossRef][Web of Science][Medline]

Del Negro CA, Morgado-Valle C, Hayes JA, Mackay DD, Pace RW, Crowder EA, and Feldman JL. Sodium and calcium current-mediated pacemaker neurons and respiratory rhythm generation. J Neurosci 25: 446–453, 2005.[Abstract/Free Full Text]

Destexhe A and Paré D. Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. J Neurophysiol 81: 1531–1547, 1999.[Abstract/Free Full Text]

Elsen FP and Ramirez JM. Calcium currents of rhythmic neurons recorded in the isolated respiratory network of neonatal mice. J Neurosci 18: 10652–10662, 1998.[Abstract/Free Full Text]

Feldmeyer D and Sakmann B. Synaptic efficacy and reliability of excitatory connections between principal neurones of the input (layer 4) and output (layer 5) of the neocortex. J Physiol 525: 31–39, 2000.[Abstract/Free Full Text]

Fisher RS, Webber WR, Lesser RP, Arroyo S, and Uematsu S. High-frequency EEG activity at the start of seizures. J Clin Neurophysiol 9: 441–448, 1992.[Web of Science][Medline]

Fleidervish IA and Gutnick MJ. Kinetics of slow inactivation of persistent sodium current in layer V neurons of mouse neocortical slices. J Neurophysiol 76: 2125–2130, 1996.[Abstract/Free Full Text]

Foehring RC, van Brederode JF, Kinney GA, and Spain WJ. Serotonergic modulation of supragranular neurons in rat sensorimotor cortex, J Neurosci 22: 8238–8250, 2002.[Abstract/Free Full Text]

Foehring RC and Wyler AR. Two patterns of firing in human neocortical neurons. Neurosci Lett 110: 279–285, 1990.[CrossRef][Web of Science][Medline]

Golomb D and Amitai Y. Propagating neuronal discharges in neocortical slices: computational and experimental study. J Neurophysiol 78: 1199–1211, 1997.[Abstract/Free Full Text]

Gray CM and McCormick DA. Chattering cells: superficial pyramidal neurons contributing to the generation of synchronous oscillations in the visual cortex. Science 274: 109–113, 1996.[Abstract/Free Full Text]

Grenier F, Timofeev I, and Steriade M. Focal synchronization of ripples (80–200 Hz) in neocortex and their neuronal correlates. J Neurophysiol 86: 1884–1898, 2001.[Abstract/Free Full Text]

Guatteo E, Franceschetti S, Bacci A, Avanzini G, and Wanke E. A TTX-sensitive conductance underlying burst firing in isolated pyramidal neurons from rat neocortex. Brain Res 741: 1–12, 1996.[CrossRef][Web of Science][Medline]

Hellwig B. A quantitative analysis of the local connectivity between pyramidal neurons in layer 2/3 of the rat visual cortex. Biol Cybern 82: 111–121, 2000.[CrossRef][Web of Science][Medline]

Hodgkin AL and Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 117: 500–544, 1952.[Free Full Text]

Huang CS, Song, JH, Nagata K, Yeh JZ, and Narahashi T. Effects of the neuroprotective agent riluzole on the high voltage-activated calcium channels of rat dorsal root ganglion neurons. J Pharmacol Exp Ther 282: 1280–1290, 1997.[Abstract/Free Full Text]

Istvan PJ and Zarzecki P. Intrinsic discharge patterns and somatosensory inputs for neurons in raccoon primary somatosensory cortex. J Neurophysiol 72: 2827–2839, 1994.[Abstract/Free Full Text]

Jacobs KM, Kharazia VN, and Prince DA. Mechanisms underlying epileptogenesis in cortical malformations. Epilepsy Res 36: 165–188, 1999.[CrossRef][Web of Science][Medline]

Jan CR, Lu YC, Jiann BP, Chang HT, and Huang JK. Effect of riluzole on cytosolic Ca2+ increase in human osteosarcoma cells. Pharmacology 66: 120–127, 2002.[CrossRef][Web of Science][Medline]

Köhling R. Voltage-gated sodium channels in epilepsy. Epilepsia 43: 1278–1295, 2002.[CrossRef][Web of Science][Medline]

Krimer LS and Goldman-Rakic PS. Prefrontal microcircuits: membrane properties and excitatory input of local, medium and wide arbor interneurons. J Neurosci 21: 3788–3796, 2001.[Abstract/Free Full Text]

Mantegazza M, Fanceschetti S, and Avanzi G. Anemone toxin (ATX II)-induced increase in persistent sodium current: effects on the firing properties of rat neocortical pyramidal neurons. J Physiol 507: 105–116, 1998.[Abstract/Free Full Text]

Marder E and Calabrese RL. Principles of rhythmic motor pattern generation. Physiol Rev 76: 687–717, 1996.[Abstract/Free Full Text]

Mountcastle VB. The columnar organization of the neocortex. Brain 120: 701–722, 1997.[Abstract/Free Full Text]

Nieuwenhuys R. The neocortex: an overview of its evolutionary development, structural organization and synaptology. Anat Embryol 190: 307–337, 1994.[Medline]

Nowak LG, Azouz R, Sanchez-Vives MV, Gray CM, and McCormick DA. Electrophysiological classes of cat primary visual cortical neurons in vivo as revealed by quantitative analyses. J Neurophysiol 89: 1541–1566, 2003.[Abstract/Free Full Text]

Nunez PL. Electric Fields of the Brain: The Neurophysics of EEG. Oxford, UK: Oxford Univ. Press, 1981.

O’Neill MJ, Bath CP, Dell CP, Hicks CA, Gilmore J, Ambler SJ, Ward MA, and Bleakman D. Effects of Ca2+ and Na+ channel inhibitors in vitro and in global cerebral ischaemia in vivo. Eur J Pharmacol 332: 121–131, 1997.[CrossRef][Web of Science][Medline]

Peña F and Ramirez JM. Endogenous activation of serotonin-2A receptors is required for respiratory rhythm generation in vitro. J Neurosci 22: 11055–11064, 2002.[Abstract/Free Full Text]

Peña F and Ramirez JM. Substance P-mediated modulation of pacemaker properties in the mammalian respiratory network. J Neurosci 24: 7549–7556, 2004.[Abstract/Free Full Text]

Peña F and Tapia R. Seizures and neurodegeneration induced by 4-aminopyridine in rat hippocampus in vivo: role of glutamate- and GABA-mediated neurotransmission and of ion channels. Neuroscience 101: 547–561, 2000.[CrossRef][Web of Science][Medline]

Ramirez JM, Tryba AK, and Peña F. Pacemaker neurons and neuronal networks: an integrative view. Curr Opin Neurobiol 14: 665–674, 2004.[CrossRef][Web of Science][Medline]

Sanabria ER, da Silva AV, Spreafico R, and Cavalheiro EA. Damage, reorganization, and abnormal neocortical hyperexcitability in the pilocarpine model of temporal lobe epilepsy. Epilepsia 43, Suppl. 5: 96–106, 2002.[Web of Science][Medline]

Sanabria ER, Su H, and Yaari Y. Initiation of network bursts by Ca2+-dependent intrinsic bursting in the rat pilocarpine model of temporal lobe epilepsy. J Physiol 532: 205–216, 2001.[Abstract/Free Full Text]

Sanchez-Vives MV and McCormick DA. Cellular and network mechanisms of rhythmic recurrent activity in neocortex. Nat Neurosci 3: 1027–1034, 2000.[CrossRef][Web of Science][Medline]

Schwindt P and Crill W. Mechanisms underlying burst and regular spiking evoked by dendritic depolarization in layer 5 cortical pyramidal neurons. J Neurophysiol 81: 1341–1354, 1999.[Abstract/Free Full Text]

Siniscalchi A, Bonci A, Mercuri NB, and Bernardi G. Effects of riluzole on rat cortical neurones: an in vitro electrophysiological study. Br J Pharmacol 120: 225–230, 1997.[CrossRef][Web of Science][Medline]

Somogyi P, Tamas G, Lujan R, and Buhl EH. Salient features of synaptic organization in the cerebral cortex. Brain Res Rev 26: 113–135, 1998.[CrossRef][Medline]

Spadoni F, Hainsworth AH, Mercuri NB, Caputi L, Martella G, Lavaroni F, Bernardi G, and Stefani A. Lamotrigine derivatives and riluzole inhibit I_Na, P in cortical neurons. Neuroreport 13: 1167–1170, 2002.[CrossRef][Web of Science][Medline]

Staiger JF, Schubert D, Zuschratter W, Kotter R, Luhmann HJ, and Zilles K. Innervation of interneurons immunoreactive for VIP by intrinsically bursting pyramidal cells and fast-spiking interneurons in infragranular layers of juvenile rat neocortex. Eur J Neurosci 16: 11–20, 2002.[CrossRef][Web of Science][Medline]

Staley KJ, Bains JS, Yee A, Hellier J, and Longacher JM. Statistical model relating CA3 burst probability to recovery from burst-induced depression at recurrent collateral synapses. J Neurophysiol 86: 2736–2747, 2001.[Abstract/Free Full Text]

Stefani A, Spadoni F, and Bernardi G. Differential inhibition by riluzole, lamotrigine, and phenytoin of sodium and calcium currents in cortical neurons: implications for neuroprotective strategies. Exp Neurol 147: 115–122, 1997.[CrossRef][Web of Science][Medline]

Steriade M. Neocortical cell classes are flexible entities. Nat Rev Neurosci 5: 121–134, 2004.[CrossRef][Web of Science][Medline]

Steriade M and Amzica F. Sleep oscillations developing into seizures in corticothalamic systems. Epilepsia 44, Suppl. 12: 9–20, 2003.

Steriade M, Amzica F, Neckelmann D, and Timofeev I. Spike-wave complexes and fast components of cortically generated seizures. II. Extra- and intracellular patterns. J Neurophysiol 80: 1456–1479, 1998.[Abstract/Free Full Text]

Steriade M, Timofeev I, and Grenier F. Natural waking and sleep states: a view from inside neocortical neurons. J Neurophysiol 85: 1969–1985, 2001.[Abstract/Free Full Text]

Steriade M, Timofeev I, Grenier F, and Durmuller N. Role of thalamic and cortical neurons in augmenting responses and self-sustained activity: dual intracellular recordings in vivo. J Neurosci 18: 6425–6443, 1998.[Abstract/Free Full Text]

Timofeev I, Grenier F, Bazhenov M, Sejnowski TJ, and Steriade M. Origin of slow cortical oscillations in deafferented cortical slabs. Cereb Cortex 10: 1185–1199, 2000a.[Abstract/Free Full Text]

Timofeev I, Grenier F, and Steriade M. Impact of intrinsic properties and synaptic factors on the activity of neocortical networks in vivo. J Physiol (Paris) 94: 343–355, 2000b.

Timofeev I, Grenier F, and Steriade M. Contribution of intrinsic neuronal factors in the generation of cortically driven electrographic seizures. J Neurophysiol 92: 1133–1143, 2004.[Abstract/Free Full Text]

Timofeev I and Steriade M. Neocortical seizures: initiation, development and cessation. Neuroscience 123: 299–336, 2004.[CrossRef][Web of Science][Medline]

Topolnik L, Steriade M, and Timofeev I. Hyperexcitability of intact neurons underlies acute development of trauma-related electrographic seizures in cats in vivo. Eur J Neurosci 18: 486–496, 2003.[CrossRef][Web of Science][Medline]

Traub RD, Buhl EH, Gloveli T, and Whittington MA. Fast rhythmic bursting can be induced in later 2/3 cortical neurons by enhancing persistent Na+ conductance or by blocking BK channels. J Neurophysiol 89: 909–921, 2003.[Abstract/Free Full Text]

Traub RD, Contreras D, Cunningham MO, Murray H, LeBeau FE, Roopun A, Bibbig A, Wilent WB, Higley MJ, and Whittington MA. Single-column thalamocortical network model exhibiting gamma oscillations, sleep spindles, and epileptogenic bursts. J Neurophysiol 93: 2194–2232, 2005.[Abstract/Free Full Text]

Traub RD, Jefferys JGR, Miles R, Whittington MA, and Toth K. A branching dendritic model of a rodent CA3 pyramidal neurone. J Physiol 481: 79–95, 1994.[Abstract/Free Full Text]

Traub RD, Whittington MA, Buhl EH, LeBeau FE, Bibbig A, Boyd S, Cross H, and Baldeweg T. A possible role for gap junctions in generation of very fast EEG oscillations preceding the onset of, or perhaps initiating, seizures. Epilepsia 42: 153–170, 2001.[CrossRef][Web of Science][Medline]

Urbani A and Belluzzi O. Riluzole inhibits the persistent sodium current in mammalian CNS neurons. Eur J Neurosci 12: 3567–3574, 2000.[CrossRef][Web of Science][Medline]

Van Drongelen W, Koch H, Marcuccilli C, Peña F, and Ramirez JM. Synchrony levels during evoked seizure-like bursts in mouse neocortical slices. J Neurophysiol 90: 1571–1580, 2003a.[Abstract/Free Full Text]

Van Drongelen W, Koch H, Peña F, Tryba A, Parkis M, Loweth J, Hecox KE, Kohrman M, Frim D, Chico MS, Ramirez JM, and Marcuccilli CJ. Differences in bursting properties between least and most abnormal tissue obtained from pediatric patients with intractable epilepsy. Program No. 99.15.2003. 2003 Abstract Viewer/Itinerary Planner. Washington, DC: Society for Neuroscience, 2003b, Online.

Van Drongelen W, Lee HC, Hereld M, Jones D, Cohoon M, Elsen F, Papka ME, and Stevens RL. Simulation of neocortical epileptiform activity using parallel computing. Neurocomputing 58–60: 1203–1209, 2004.

Van Drongelen W, Lee HC, Koch, Hereld M, Elsen F, Chen Z, and Stevens RL. Emergent epileptiform activity in neural networks with weak excitatory synapses. IEEE Trans Neural Syst Rehab 13: 236–241, 2005.[CrossRef]

Zona C, Siniscalchi A, Mercuri NB, and Bernardi G. Riluzole interacts with voltage-activated sodium and potassium currents in cultured rat cortical neurons. Neuroscience 85: 931–938, 1998.[CrossRef][Web of Science][Medline]




This article has been cited by other articles:


Home page
J. Neurophysiol.Home page
A. Czarnecki, V. Magloire, and J. Streit
Modulation of Intrinsic Spiking in Spinal Cord Neurons
J Neurophysiol, October 1, 2009; 102(4): 2441 - 2452.
[Abstract] [Full Text] [PDF]


Home page
J. Neurosci.Home page
S. Tazerart, L. Vinay, and F. Brocard
The Persistent Sodium Current Generates Pacemaker Activities in the Central Pattern Generator for Locomotion and Regulates the Locomotor Rhythm
J. Neurosci., August 20, 2008; 28(34): 8577 - 8589.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
A. Enomoto, J. M. Han, C.-F. Hsiao, and S. H. Chandler
Sodium Currents in Mesencephalic Trigeminal Neurons From Nav1.6 Null Mice
J Neurophysiol, August 1, 2007; 98(2): 710 - 719.
[Abstract] [Full Text] [PDF]


Home page
J. Neurosci.Home page
G. Zhong, M. A. Masino, and R. M. Harris-Warrick
Persistent Sodium Currents Participate in Fictive Locomotion Generation in Neonatal Mouse Spinal Cord
J. Neurosci., April 25, 2007; 27(17): 4507 - 4518.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
96/5/2564    most recent
00446.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by van Drongelen, W.
Right arrow Articles by Ramirez, J.-M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by van Drongelen, W.
Right arrow Articles by Ramirez, J.-M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the The American Physiological Society.