Role of persistent sodium current in bursting activity of mouse neocortical networks in vitro

Interdependence of cellular and network properties in respiratory rhythmogenesis

How breathing is generated by the preBötzinger Complex (preBötC) remains divided between two ideological frameworks, and the persistent sodium current (INaP) lies at the heart of this debate. Although INaP is widely expressed, the pacemaker hypothesis considers it essential because it endows a small subset of neurons with intrinsic bursting or "pacemaker" activity. In contrast, burstlet theory considers INaP dispensable because rhythm emerges from "pre-inspiratory" spiking activity driven by feed-forward network interactions. Using computational modeling, we discover that changes in spike shape can dissociate INaP from intrinsic bursting. Consistent with many experimental benchmarks, conditional effects on spike shape during simulated changes in oxygenation, development, extracellular potassium, and temperature alter the prevalence of intrinsic bursting and pre-inspiratory spiking without altering the role of INaP. Our results support a unifying hypothesis where INaP and excitatory network interactions, but not intrinsic bursting or pre-inspiratory spiking, are critical interdependent features of preBötC rhythmogenesis.

Guanfacine inhibits interictal epileptiform events and sodium currents in prefrontal cortex pyramidal neurons

BACKGROUND

Guanfacine (an alpha-2A receptor agonist) is a commonly used drug with recognized efficacy in the treatment of attention deficit hyperactivity disorder (ADHD). This study aimed to assess the effects of guanfacine on short-lasting (interictal) epileptiform discharges in cortical neurons. Moreover, we assessed the effects of guanfacine on voltage-gated sodium currents.

METHODS

We conducted patch-clamp recordings in prefrontal cortex pyramidal neurons obtained from young rats. Interictal epileptiform events were evoked in cortical slices in a zero magnesium proepileptic extracellular solution with an elevated concentration of potassium ions.

RESULTS

Interictal epileptiform discharges were spontaneous depolarisations, which triggered action potentials. Guanfacine (10 and 100 µM) inhibited the frequency of epileptiform discharges. The effect of guanfacine on interictal events persisted in the presence of alpha-2 adrenergic receptor antagonist idazoxan. The tested drug inhibited neuronal excitability. Tonic NMDA currents were not influenced by guanfacine. Recordings from dispersed neurons showed that the tested drug (10 and 100 µM) inhibited persistent and fast inactivating voltage-gated sodium currents.

CONCLUSIONS

This study shows that guanfacine inhibits interictal discharges in cortical neurons independently of alpha-2A adrenergic receptors. This effect may be mediated by voltage-gated sodium currents. Inhibition of interictal activity by guanfacine may be of clinical importance because interictal events often occur in patients with ADHD and may contribute to symptoms of this disease.

Selective inhibition of excitatory synaptic transmission alters the emergent bursting dynamics of in vitro neural networks

Neurons in vitro connect to each other and form neural networks that display emergent electrophysiological activity. This activity begins as spontaneous uncorrelated firing in the early phase of development, and as functional excitatory and inhibitory synapses mature, the activity typically emerges as spontaneous network bursts. Network bursts are events of coordinated global activation among many neurons interspersed with periods of silencing and are important for synaptic plasticity, neural information processing, and network computation. While bursting is the consequence of balanced excitatory-inhibitory (E/I) interactions, the functional mechanisms underlying their evolution from physiological to potentially pathophysiological states, such as decreasing or increasing in synchrony, are still poorly understood. Synaptic activity, especially that related to maturity of E/I synaptic transmission, is known to strongly influence these processes. In this study, we used selective chemogenetic inhibition to target and disrupt excitatory synaptic transmission in in vitro neural networks to study functional response and recovery of spontaneous network bursts over time. We found that over time, inhibition resulted in increases in both network burstiness and synchrony. Our results indicate that the disruption in excitatory synaptic transmission during early network development likely affected inhibitory synaptic maturity which resulted in an overall decrease in network inhibition at later stages. These findings lend support to the importance of E/I balance in maintaining physiological bursting dynamics and, conceivably, information processing capacity in neural networks.

Autocrine Neuromodulation and Network Activity Patterns in the Locus Coeruleus of Newborn Rat Slices

Already in newborns, the locus coeruleus (LC) controls multiple brain functions and may have a complex organization as in adults. Our findings in newborn rat brain slices indicate that LC neurons (i) generate at ~1 Hz a ~0.3 s-lasting local field potential (LFP) comprising summated phase-locked single spike discharge, (ii) express intrinsic ‘pacemaker’ or ‘burster’ properties and (iii) receive solely excitatory or initially excitatory−secondary inhibitory inputs. μ-opioid or ɑ2 noradrenaline receptor agonists block LFP rhythm at 100−250 nM whereas slightly lower doses transform its bell-shaped pattern into slower crescendo-shaped multipeak bursts. GABAA and glycine receptors hyperpolarize LC neurons to abolish rhythm which remains though unaffected by blocking them. Rhythm persists also during ionotropic glutamate receptor (iGluR) inhibition whereas <10 mV depolarization during iGluR agonists accelerates spiking to cause subtype-specific fast (spindle-shaped) LFP oscillations. Similar modest neuronal depolarization causing a cytosolic Ca2+ rise occurs (without effect on neighboring astrocytes) during LFP acceleration by CNQX activating a TARP-AMPA-type iGluR complex. In contrast, noradrenaline lowers neuronal Ca2+ baseline via ɑ2 receptors, but evokes an ɑ1 receptor-mediated ‘concentric’ astrocytic Ca2+ wave. In summary, the neonatal LC has a complex (possibly modular) organization to enable discharge pattern transformations that might facilitate discrete actions on target circuits.

Ion dynamics at the energy-deprived tripartite synapse

The anatomical and functional organization of neurons and astrocytes at 'tripartite synapses' is essential for reliable neurotransmission, which critically depends on ATP. In low energy conditions, synaptic transmission fails, accompanied by a breakdown of ion gradients, changes in membrane potentials and cell swelling. The resulting cellular damage and cell death are causal to the often devastating consequences of an ischemic stroke. The severity of ischemic damage depends on the age and the brain region in which a stroke occurs, but the reasons for this differential vulnerability are far from understood. In the present study, we address this question by developing a comprehensive biophysical model of a glutamatergic synapse to identify key determinants of synaptic failure during energy deprivation. Our model is based on fundamental biophysical principles, includes dynamics of the most relevant ions, i.e., Na+, K+, Ca2+, Cl- and glutamate, and is calibrated with experimental data. It confirms the critical role of the Na+/K+-ATPase in maintaining ion gradients, membrane potentials and cell volumes. Our simulations demonstrate that the system exhibits two stable states, one physiological and one pathological. During energy deprivation, the physiological state may disappear, forcing a transit to the pathological state, which can be reverted when blocking voltage-gated Na+ and K+ channels. Our model predicts that the transition to the pathological state is favoured if the extracellular space fraction is small. A reduction in the extracellular space volume fraction, as, e.g. observed with ageing, will thus promote the brain's susceptibility to ischemic damage. Our work provides new insights into the brain's ability to recover from energy deprivation, with translational relevance for diagnosis and treatment of ischemic strokes.

Differential Inhibition of Human Nav1.2 Resurgent and Persistent Sodium Currents by Cannabidiol and GS967

Many epilepsy patients are refractory to conventional antiepileptic drugs. Resurgent and persistent currents can be enhanced by epilepsy mutations in the Nav1.2 channel, but conventional antiepileptic drugs inhibit normal transient currents through these channels, along with aberrant resurgent and persistent currents that are enhanced by Nav1.2 epilepsy mutations. Pharmacotherapies that specifically target aberrant resurgent and/or persistent currents would likely have fewer unwanted side effects and be effective in many patients with refractory epilepsy. This study investigated the effects of cannbidiol (CBD) and GS967 (each at 1 μM) on transient, resurgent, and persistent currents in human embryonic kidney (HEK) cells stably expressing wild-type hNav1.2 channels. We found that CBD preferentially inhibits resurgent currents over transient currents in this paradigm; and that GS967 preferentially inhibits persistent currents over transient currents. Therefore, CBD and GS967 may represent a new class of more targeted and effective antiepileptic drugs.

Resurgent and Gating Pore Currents Induced by De Novo SCN2A Epilepsy Mutations

Over 150 mutations in the SCN2A gene, which encodes the neuronal Nav1.2 protein, have been implicated in human epilepsy cases. Of these, R1882Q and R853Q are two of the most commonly reported mutations. This study utilized voltage-clamp electrophysiology to characterize the biophysical effects of the R1882Q and R853Q mutations on the hNav1.2 channel, including their effects on resurgent current and gating pore current, which are not typically investigated in the study of Nav1.2 channel mutations. HEK cells transiently transfected with DNA encoding either wild-type (WT) or mutant hNav1.2 revealed that the R1882Q mutation induced a gain-of-function phenotype, including slowed fast inactivation, depolarization of the voltage dependence of inactivation, and increased persistent current. In this model system, the R853Q mutation primarily produced loss-of-function effects, including reduced transient current amplitude and density, hyperpolarization of the voltage dependence of inactivation, and decreased persistent current. The presence of a Navβ4 peptide (KKLITFILKKTREK-OH) in the pipette solution induced resurgent currents, which were increased by the R1882Q mutation and decreased by the R853Q mutation. Further study of the R853Q mutation in Xenopus oocytes indicated a reduced surface expression and revealed a robust gating pore current at negative membrane potentials, a function absent in the WT channel. This not only shows that different epileptogenic point mutations in hNav1.2 have distinct biophysical effects on the channel, but also illustrates that individual mutations can have complex consequences that are difficult to identify using conventional analyses. Distinct mutations may, therefore, require tailored pharmacotherapies in order to eliminate seizures.

Effects of persistent sodium current blockade in respiratory circuits depend on the pharmacological mechanism of action and network dynamics

The mechanism(s) of action of most commonly used pharmacological blockers of voltage-gated ion channels are well understood; however, this knowledge is rarely considered when interpreting experimental data. Effects of blockade are often assumed to be equivalent, regardless of the mechanism of the blocker involved. Using computer simulations, we demonstrate that this assumption may not always be correct. We simulate the blockade of a persistent sodium current (I_NaP), proposed to underlie rhythm generation in pre-Bötzinger complex (pre-BötC) respiratory neurons, via two distinct pharmacological mechanisms: (1) pore obstruction mediated by tetrodotoxin and (2) altered inactivation dynamics mediated by riluzole. The reported effects of experimental application of tetrodotoxin and riluzole in respiratory circuits are diverse and seemingly contradictory and have led to considerable debate within the field as to the specific role of I_NaP in respiratory circuits. The results of our simulations match a wide array of experimental data spanning from the level of isolated pre-BötC neurons to the level of the intact respiratory network and also generate a series of experimentally testable predictions. Specifically, in this study we: (1) provide a mechanistic explanation for seemingly contradictory experimental results from in vitro studies of I_NaP block, (2) show that the effects of I_NaP block in in vitro preparations are not necessarily equivalent to those in more intact preparations, (3) demonstrate and explain why riluzole application may fail to effectively block I_NaP in the intact respiratory network, and (4) derive the prediction that effective block of I_NaP by low concentration tetrodotoxin will stop respiratory rhythm generation in the intact respiratory network. These simulations support a critical role for I_NaP in respiratory rhythmogenesis in vivo and illustrate the importance of considering mechanism when interpreting and simulating data relating to pharmacological blockade.

The application of pharmacological agents that affect transmembrane ionic currents in neurons is a commonly used experimental technique. A simplistic interpretation of experiments involving these agents suggests that antagonist application removes the impacted current and that subsequently observed changes in activity are attributable to the loss of that current’s effects. The more complex reality, however, is that different drugs may have distinct mechanisms of action, some corresponding not to a removal of a current but rather to a changing of its properties. We use computational modeling to explore the implications of the distinct mechanisms associated with two drugs, riluzole and tetrodotoxin, that are often characterized as sodium channel blockers. Through this approach, we offer potential explanations for disparate findings observed in experiments on neural respiratory circuits and show that the experimental results are consistent with a key role for the persistent sodium current in respiratory rhythm generation.

Role of paroxysmal depolarization in focal seizure activity

We analyze the role of inhibition in sustaining focal epileptic seizure activity. We review ongoing seizure activity at the mesoscopic scale that can be observed with microelectrode arrays as well as at the macroscale of standard clinical EEG. We provide clinical, experimental, and modeling data to support the hypothesis that paroxysmal depolarization (PD) is a critical component of the ictal machinery. We present dual-patch recordings in cortical cultures showing reduced synaptic transmission associated with presynaptic occurrence of PD, and we find that the PD threshold is cell size related. We further find evidence that optically evoked PD activity in parvalbumin neurons can promote propagation of neuronal excitation in neocortical networks in vitro. Spike sorting results from microelectrode array measurements around ictal wave propagation in human focal seizures demonstrate a strong increase in putative inhibitory firing with an approaching excitatory wave, followed by a sudden reduction of firing at passage. At the macroscopic level, we summarize evidence that this excitatory ictal wave activity is strongly correlated with oscillatory activity across a centimeter-sized cortical network. We summarize Wilson-Cowan-type modeling showing how inhibitory function is crucial for this behavior. Our findings motivated us to develop a network motif of neurons in silico, governed by a reduced version of the Hodgkin-Huxley formalism, to show how feedforward, feedback, PD, and local failure of inhibition contribute to observed dynamics across network scales. The presented multidisciplinary evidence suggests that the PD not only is a cellular marker or epiphenomenon but actively contributes to seizure activity.NEW & NOTEWORTHY We present mechanisms of ongoing focal seizures across meso- and macroscales of microelectrode array and standard clinical recordings, respectively. We find modeling, experimental, and clinical evidence for a dual role of inhibition across these scales: local failure of inhibition allows propagation of a mesoscopic ictal wave, whereas inhibition elsewhere remains intact and sustains macroscopic oscillatory activity. We present evidence for paroxysmal depolarization as a mechanism behind this dual role of inhibition in shaping ictal activity.

Kinetic properties of persistent Na+ current orchestrate oscillatory bursting in respiratory neurons

The rhythmic pattern of breathing depends on the pre-Bötzinger complex (preBötC) in the brainstem, a vital circuit that contains a population of neurons with intrinsic oscillatory bursting behavior. Here, we investigate the specific kinetic properties that enable voltage-gated sodium channels to establish oscillatory bursting in preBötC inspiratory neurons, which exhibit an unusually large persistent Na⁺ current (I_NaP). We first characterize the kinetics of I_NaP in neonatal rat brainstem slices in vitro, using whole-cell patch-clamp and computational modeling, and then test the contribution of I_NaP to rhythmic bursting in live neurons, using the dynamic clamp technique. We provide evidence that subthreshold activation, persistence at suprathreshold potentials, slow inactivation, and slow recovery from inactivation are kinetic features of I_NaP that regulate all aspects of intrinsic rhythmic bursting in preBötC neurons. The slow and cumulative inactivation of I_NaP during the burst active phase controls burst duration and termination, while the slow recovery from inactivation controls the duration of the interburst interval. To demonstrate this mechanism, we develop a Markov state model of I_NaP that explains a comprehensive set of voltage clamp data. By adding or subtracting a computer-generated I_NaP from a live neuron via dynamic clamp, we are able to convert nonbursters into intrinsic bursters, and vice versa. As a control, we test a model with inactivation features removed. Adding noninactivating I_NaP into nonbursters results in a pattern of random transitions between sustained firing and quiescence. The relative amplitude of I_NaP is the key factor that separates intrinsic bursters from nonbursters and can change the fraction of intrinsic bursters in the preBötC. I_NaP could thus be an important target for regulating network rhythmogenic properties.

mGluR5-dependent modulation of dendritic excitability in CA1 pyramidal neurons mediated by enhancement of persistent Na+ currents

KEY POINTS

High-frequency stimulation (HFS) of the Schaffer collateral pathway activates metabotropic glutamate receptor 5 (mGluR5) signalling in the proximal apical dendrites of CA1 pyramidal neurons. The synaptic activation of mGluR5-mediated calcium signalling causes a significant increase in persistent sodium current (INa,P ) in the dendrites. Increased INa,P by HFS underlies potentiation of synaptic inputs at both the proximal and distal dendrite, leading to an enhanced probability of action potential firing associated with decreased action potential thresholds. Therefore, HFS-induced activation of intracellular mGluR5 serves an important role as an instructive signal for potentiation of upcoming inputs by increasing dendritic excitability.

ABSTRACT

Dendritic Na+ channels in pyramidal neurons are known to amplify synaptic signals, thereby facilitating action potential (AP) generation. However, the mechanisms that modulate dendritic Na+ channels have remained largely uncharacterized. Here, we report a new form of short-term plasticity in which proximal excitatory synaptic inputs to hippocampal CA1 pyramidal neurons transiently elevate dendritic excitability. High-frequency stimulations (HFS) to the Schaffer collateral (SC) pathway activate mGluR5-dependent Ca2+ signalling in the apical dendrites, which, with calmodulin, upregulates specifically Nav1.6 channel-mediated persistent Na+ currents (INa,P ) in the dendrites. This HFS-induced increase in dendritic INa,P results in transient increases in the amplitude of excitatory postsynaptic potentials induced by both proximal SC and distal perforant path stimulation, leading to the enhanced probability of AP firing associated with decreased AP thresholds. Taken together, our study identifies dendritic INa,P as a novel target for mediating activity-dependent modulation of dendritic integration and neuronal output.

The Dynamic Basis of Respiratory Rhythm Generation: One Breath at a Time

Rhythmicity is a universal timing mechanism in the brain, and the rhythmogenic mechanisms are generally dynamic. This is illustrated for the neuronal control of breathing, a behavior that occurs as a one-, two-, or three-phase rhythm. Each breath is assembled stochastically, and increasing evidence suggests that each phase can be generated independently by a dedicated excitatory microcircuit. Within each microcircuit, rhythmicity emerges through three entangled mechanisms: ( a) glutamatergic transmission, which is amplified by ( b) intrinsic bursting and opposed by ( c) concurrent inhibition. This rhythmogenic triangle is dynamically tuned by neuromodulators and other network interactions. The ability of coupled oscillators to reconfigure and recombine may allow breathing to remain robust yet plastic enough to conform to nonventilatory behaviors such as vocalization, swallowing, and coughing. Lessons learned from the respiratory network may translate to other highly dynamic and integrated rhythmic systems, if approached one breath at a time.

Cross-scale effects of neural interactions during human neocortical seizure activity

Small-scale neuronal networks may impose widespread effects on large network dynamics. To unravel this relationship, we analyzed eight multiscale recordings of spontaneous seizures from four patients with epilepsy. During seizures, multiunit spike activity organizes into a submillimeter-sized wavefront, and this activity correlates significantly with low-frequency rhythms from electrocorticographic recordings across a 10-cm-sized neocortical network. Notably, this correlation effect is specific to the ictal wavefront and is absent interictally or from action potential activity outside the wavefront territory. To examine the multiscale interactions, we created a model using a multiscale, nonlinear system and found evidence for a dual role for feedforward inhibition in seizures: while inhibition at the wavefront fails, allowing seizure propagation, feedforward inhibition of the surrounding centimeter-scale networks is activated via long-range excitatory connections. Bifurcation analysis revealed that distinct dynamical pathways for seizure termination depend on the surrounding inhibition strength. Using our model, we found that the mesoscopic, local wavefront acts as the forcing term of the ictal process, while the macroscopic, centimeter-sized network modulates the oscillatory seizure activity.

Pro-excitatory alterations in sodium channel activity facilitate subiculum neuron hyperexcitability in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is a common form of adult epilepsy involving the limbic structures of the temporal lobe. Subiculum neurons act to provide a major output from the hippocampus and consist of a large population of endogenously bursting excitatory neurons. In TLE, subiculum neurons are largely spared, become hyperexcitable and show spontaneous epileptiform activity. The basis for this hyperexcitability is unclear, but is likely to involve alterations in the expression levels and function of various ion channels. In this study, we sought to determine the importance of sodium channel currents in facilitating neuronal hyperexcitability of subiculum neurons in the continuous hippocampal stimulation (CHS) rat model of TLE. Subiculum neurons from TLE rats were hyperexcitable, firing a higher frequency of action potentials after somatic current injection and action potential (AP) bursts after synaptic stimulation. Voltage clamp recordings revealed increases in resurgent (I_NaR) and persistent (I_NaP) sodium channel currents and pro-excitatory shifts in sodium channel activation and inactivation parameters that would facilitate increases in AP generation. Attenuation of I_NaR and I_NaP currents with 4,9-anhydro-tetrodotoxin (4,9-ah TTX; 100nM), a toxin with increased potency against Na_v1.6 channels, suppressed neuronal firing frequency and inhibited AP bursting induced by synaptic stimulation in TLE neurons. These findings support an important role of sodium channels, particularly Na_v1.6, in facilitating subiculum neuron hyperexcitability in TLE and provide further support for the importance of I_NaR and I_NaP currents in establishing epileptiform activity of subiculum neurons.

Ionic mechanisms maintaining action potential conduction velocity at high firing frequencies in an unmyelinated axon

The descending contralateral movement detector (DCMD) is a high‐performance interneuron in locusts with an axon capable of transmitting action potentials (AP) at more than 500 Hz. We investigated biophysical mechanisms for fidelity of high‐frequency transmission in this axon. We measured conduction velocities (CVs) at room temperature during exposure to 10 mmol/L cadmium, a calcium current antagonist, and found significant reduction in CV with reduction at frequencies >200 Hz of ~10%. Higher temperatures induced greater CV reductions during exposure to cadmium across all frequencies of ~20–30%. Intracellular recordings during 15 min of exposure to cadmium or nickel, also a calcium current antagonist, revealed an increase in the magnitude of the afterhyperpolarization potential (AHP) and the time to recover to baseline after the AHP (Medians for Control: −19.8%; Nickel: 167.2%; Cadmium: 387.2%), that could be due to a T‐type calcium current. However, the removal of extracellular calcium did not mimic divalent cation exposure suggesting calcium currents are not the cause of the AHP increase. Computational modeling showed that the effects of the divalent cations could be modeled with a persistent sodium current which could be blocked by high concentrations of divalent cations. Persistent sodium current shortened the AHP duration in our models and increased CV for high‐frequency APs. We suggest that faithful, high‐frequency axonal conduction in the DCMD is enabled by a mechanism that shortens the AHP duration like a persistent or resurgent sodium current.

Development of pacemaker properties and rhythmogenic mechanisms in the mouse embryonic respiratory network

Breathing is a vital rhythmic behavior generated by hindbrain neuronal circuitry, including the preBötzinger complex network (preBötC) that controls inspiration. The emergence of preBötC network activity during prenatal development has been described, but little is known regarding inspiratory neurons expressing pacemaker properties at embryonic stages. Here, we combined calcium imaging and electrophysiological recordings in mouse embryo brainstem slices together with computational modeling to reveal the existence of heterogeneous pacemaker oscillatory properties relying on distinct combinations of burst-generating INaP and ICAN conductances. The respective proportion of the different inspiratory pacemaker subtypes changes during prenatal development. Concomitantly, network rhythmogenesis switches from a purely INaP/ICAN-dependent mechanism at E16.5 to a combined pacemaker/network-driven process at E18.5. Our results provide the first description of pacemaker bursting properties in embryonic preBötC neurons and indicate that network rhythmogenesis undergoes important changes during prenatal development through alterations in both circuit properties and the biophysical characteristics of pacemaker neurons.

Multiscale Aspects of Generation of High-Gamma Activity during Seizures in Human Neocortex

High-gamma (HG; 80-150 Hz) activity in macroscopic clinical records is considered a marker for critical brain regions involved in seizure initiation; it is correlated with pathological multiunit firing during neocortical seizures in the seizure core, an area identified by correlated multiunit spiking and low frequency seizure activity.

High-gamma (HG; 80-150 Hz) activity in macroscopic clinical records is considered a marker for critical brain regions involved in seizure initiation; it is correlated with pathological multiunit firing during neocortical seizures in the seizure core, an area identified by correlated multiunit spiking and low frequency seizure activity. However, the effects of the spatiotemporal dynamics of seizure on HG power generation are not well understood. Here, we studied HG generation and propagation, using a three-step, multiscale signal analysis and modeling approach. First, we analyzed concurrent neuronal and microscopic network HG activity in neocortical slices from seven intractable epilepsy patients. We found HG activity in these networks, especially when neurons displayed paroxysmal depolarization shifts and network activity was highly synchronized. Second, we examined HG activity acquired with microelectrode arrays recorded during human seizures (n = 8). We confirmed the presence of synchronized HG power across microelectrode records and the macroscale, both specifically associated with the core region of the seizure. Third, we used volume conduction-based modeling to relate HG activity and network synchrony at different network scales. We showed that local HG oscillations require high levels of synchrony to cross scales, and that this requirement is met at the microscopic scale, but not within macroscopic networks. Instead, we present evidence that HG power at the macroscale may result from harmonics of ongoing seizure activity. Ictal HG power marks the seizure core, but the generating mechanism can differ across spatial scales.

Network burst activity in hippocampal neuronal cultures: the role of synaptic and intrinsic currents

The goal of this work was to define the contributions of intrinsic and synaptic mechanisms toward spontaneous network-wide bursting activity, observed in dissociated rat hippocampal cell cultures. This network behavior is typically characterized by short-duration bursts, separated by order of magnitude longer interburst intervals. We hypothesize that while short-timescale synaptic processes modulate spectro-temporal intraburst properties and network-wide burst propagation, much longer timescales of intrinsic membrane properties such as persistent sodium (Nap) currents govern burst onset during interburst intervals. To test this, we used synaptic receptor antagonists picrotoxin, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), and 3-(2-carboxypiperazine-4-yl)propyl-1-phosphonate (CPP) to selectively block GABAA, AMPA, and NMDA receptors and riluzole to selectively block Nap channels. We systematically compared intracellular activity (recorded with patch clamp) and network activity (recorded with multielectrode arrays) in eight different synaptic connectivity conditions: GABAA + NMDA + AMPA, NMDA + AMPA, GABAA + AMPA, GABAA + NMDA, AMPA, NMDA, GABAA, and all receptors blocked. Furthermore, we used mixed-effects modeling to quantify the aforementioned independent and interactive synaptic receptor contributions toward spectro-temporal burst properties including intraburst spike rate, burst activity index, burst duration, power in the local field potential, network connectivity, and transmission delays. We found that blocking intrinsic Nap currents completely abolished bursting activity, demonstrating their critical role in burst onset within the network. On the other hand, blocking different combinations of synaptic receptors revealed that spectro-temporal burst properties are uniquely associated with synaptic functionality and that excitatory connectivity is necessary for the presence of network-wide bursting. In addition to confirming the critical contribution of direct excitatory effects, mixed-effects modeling also revealed distinct combined (nonlinear) contributions of excitatory and inhibitory synaptic activity to network bursting properties.

Modeling focal epileptic activity in the Wilson-cowan model with depolarization block

Measurements of neuronal signals during human seizure activity and evoked epileptic activity in experimental models suggest that, in these pathological states, the individual nerve cells experience an activity driven depolarization block, i.e. they saturate. We examined the effect of such a saturation in the Wilson-Cowan formalism by adapting the nonlinear activation function; we substituted the commonly applied sigmoid for a Gaussian function. We discuss experimental recordings during a seizure that support this substitution. Next we perform a bifurcation analysis on the Wilson-Cowan model with a Gaussian activation function. The main effect is an additional stable equilibrium with high excitatory and low inhibitory activity. Analysis of coupled local networks then shows that such high activity can stay localized or spread. Specifically, in a spatial continuum we show a wavefront with inhibition leading followed by excitatory activity. We relate our model simulations to observations of spreading activity during seizures.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1186/s13408-015-0019-4) contains supplementary material 1.

Large-scale modeling of epileptic seizures: scaling properties of two parallel neuronal network simulation algorithms

Our limited understanding of the relationship between the behavior of individual neurons and large neuronal networks is an important limitation in current epilepsy research and may be one of the main causes of our inadequate ability to treat it. Addressing this problem directly via experiments is impossibly complex; thus, we have been developing and studying medium-large-scale simulations of detailed neuronal networks to guide us. Flexibility in the connection schemas and a complete description of the cortical tissue seem necessary for this purpose. In this paper we examine some of the basic issues encountered in these multiscale simulations. We have determined the detailed behavior of two such simulators on parallel computer systems. The observed memory and computation-time scaling behavior for a distributed memory implementation were very good over the range studied, both in terms of network sizes (2,000 to 400,000 neurons) and processor pool sizes (1 to 256 processors). Our simulations required between a few megabytes and about 150 gigabytes of RAM and lasted between a few minutes and about a week, well within the capability of most multinode clusters. Therefore, simulations of epileptic seizures on networks with millions of cells should be feasible on current supercomputers.

The cellular building blocks of breathing

Respiratory brainstem neurons fulfill critical roles in controlling breathing: they generate the activity patterns for breathing and contribute to various sensory responses including changes in O2 and CO2. These complex sensorimotor tasks depend on the dynamic interplay between numerous cellular building blocks that consist of voltage-, calcium-, and ATP-dependent ionic conductances, various ionotropic and metabotropic synaptic mechanisms, as well as neuromodulators acting on G-protein coupled receptors and second messenger systems. As described in this review, the sensorimotor responses of the respiratory network emerge through the state-dependent integration of all these building blocks. There is no known respiratory function that involves only a small number of intrinsic, synaptic, or modulatory properties. Because of the complex integration of numerous intrinsic, synaptic, and modulatory mechanisms, the respiratory network is capable of continuously adapting to changes in the external and internal environment, which makes breathing one of the most integrated behaviors. Not surprisingly, inspiration is critical not only in the control of ventilation, but also in the context of "inspiring behaviors" such as arousal of the mind and even creativity. Far-reaching implications apply also to the underlying network mechanisms, as lessons learned from the respiratory network apply to network functions in general.

The effect of morphology upon electrophysiological responses of retinal ganglion cells: simulation results

Retinal ganglion cells (RGCs) display differences in their morphology and intrinsic electrophysiology. The goal of this study is to characterize the ionic currents that explain the behavior of ON and OFF RGCs and to explore if all morphological types of RGCs exhibit the phenomena described in electrophysiological data. We extend our previous single compartment cell models of ON and OFF RGCs to more biophysically realistic multicompartment cell models and investigate the effect of cell morphology on intrinsic electrophysiological properties. The membrane dynamics are described using the Hodgkin - Huxley type formalism. A subset of published patch-clamp data from isolated intact mouse retina is used to constrain the model and another subset is used to validate the model. Two hundred morphologically distinct ON and OFF RGCs are simulated with various densities of ionic currents in different morphological neuron compartments. Our model predicts that the differences between ON and OFF cells are explained by the presence of the low voltage activated calcium current in OFF cells and absence of such in ON cells. Our study shows through simulation that particular morphological types of RGCs are capable of exhibiting the full range of phenomena described in recent experiments. Comparisons of outputs from different cells indicate that the RGC morphologies that best describe recent experimental results are ones that have a larger ratio of soma to total surface area.

Modeling Neural Activity

This paper provides an overview of different types of models for studying activity of nerve cells and their networks with a special emphasis on neural oscillations. One part describes the neuronal models based on the Hodgkin and Huxley formalism first described in the 1950s. It is discussed how further simplifications of this formalism enable mathematical analysis of the process of neural excitability. The focus of the paper’s second component is on network activity. Understanding network function is one of the important frontiers remaining in neuroscience. At present, experimental techniques can only provide global recordings or samples of the activity of the huge networks that form the nervous system. Models in neuroscience can therefore play a critical role by providing a framework for integration of necessarily incomplete datasets, thereby providing insight into the mechanisms of neural function. Network models can either explicitly contain individual network nodes that model the neurons, or they can be based on representations of compound population activity. The latter approach was pioneered by Wilson and Cowan in the 1970s. Finally I provide an overview and discuss how network models are employed in the study of neuronal network pathology such as epilepsy.

The role of voltage dependence of the NMDA receptor in cellular and network oscillation

Unraveling the mechanisms underlying oscillatory behavior is critical for understanding normal and pathological brain processes. Here we used electrophysiology in mouse neocortical slices and principles of nonlinear dynamics to demonstrate how an increase in the N-methyl-d-aspartic acid receptor (NMDAR) conductance can create a nonlinear whole-cell current-voltage (I-V) relationship which leads to changes in cellular stability. We discovered two behaviorally and morphologically distinct pyramidal cell populations. Under control conditions, both cell types responded to depolarizing current injection with regular spiking patterns. However, upon NMDAR activation, an intrinsic oscillatory (IO) cell type (n = 44) showed a nonlinear whole-cell I-V relationship, intrinsic voltage-dependent oscillations plus amplification of alternating input current, and these properties persisted after disabling action potential generation with tetrodotoxin (TTX). The other non-oscillatory (NO) neuronal population (n = 24) demonstrated none of these behaviors. Simultaneous intra- and extracellular recordings demonstrated the NMDAR's capacity to promote low-frequency seizure-like network oscillations via its effects on intrinsic neuronal properties. The two pyramidal cell types demonstrated different relationships with network oscillation--the IO cells were leaders that were activated early in the population activity cycle while the activation of the NO cell type was distributed across network bursts. The properties of IO neurons disappeared in a low-magnesium environment where the voltage dependence of the receptor is abolished; concurrently, the cellular contribution to network oscillation switched to synchronous firing. Thus, depending upon the efficacy of NMDAR in altering the linearity of the whole-cell I-V relationship, the two cell populations played different roles in sustaining network oscillation.

Cycle-by-cycle assembly of respiratory network activity is dynamic and stochastic

Rhythmically active networks are typically composed of neurons that can be classified as silent, tonic spiking, or rhythmic bursting based on their intrinsic activity patterns. Within these networks, neurons are thought to discharge in distinct phase relationships with their overall network output, and it has been hypothesized that bursting pacemaker neurons may lead and potentially trigger cycle onsets. We used multielectrode recording from 72 experiments to test these ideas in rhythmically active slices containing the pre-Bötzinger complex, a region critical for breathing. Following synaptic blockade, respiratory neurons exhibited a gradient of intrinsic spiking to rhythmic bursting activities and thus defied an easy classification into bursting pacemaker and nonbursting categories. Features of their firing activity within the functional network were analyzed for correlation with subsequent rhythmic bursting in synaptic isolation. Higher firing rates through all phases of fictive respiration statistically predicted bursting pacemaker behavior. However, a cycle-by-cycle analysis indicated that respiratory neurons were stochastically activated with each burst. Intrinsically bursting pacemakers led some population bursts and followed others. This variability was not reproduced in traditional fully interconnected computational models, while sparsely connected network models reproduced these results both qualitatively and quantitatively. We hypothesize that pacemaker neurons do not act as clock-like drivers of the respiratory rhythm but rather play a flexible and dynamic role in the initiation and stabilization of each burst. Thus, at the behavioral level, each breath can be thought of as de novo assembly of a stochastic collaboration of network topology and intrinsic properties.

Metabotropic glutamate receptors (mGluR5) activate transient receptor potential canonical channels to improve the regularity of the respiratory rhythm generated by the pre-Bötzinger complex in mice

Metabotropic glutamate receptors (mGluRs) are hypothesized to play a key role in generating the central respiratory rhythm and other rhythmic activities driven by central pattern generators (e.g. locomotion). However, the functional role of mGluRs in rhythmic respiratory activity and many motor patterns is very poorly understood. Here, we used mouse respiratory brain-slice preparations containing the pre-Bötzinger complex (pre-BötC) to identify the role of group I mGluRs (mGluR1 and mGluR5) in respiratory rhythm generation. We found that activation of mGluR1/5 is not required for the pre-BötC to generate a respiratory rhythm. However, our data suggest that mGluR1 and mGluR5 differentially modulate the respiratory rhythm. Blocking endogenous mGluR5 activity with 2-Methyl-6-(phenylethynyl)pyridine (MPEP) decreases the inspiratory burst duration, burst area and frequency, whereas it increases the irregularity of the fictive eupneic inspiratory rhythm generated by the pre-BötC. In contrast, blocking mGluR1 reduces the frequency. Moreover, the mGluR1/5 agonist 3,5-dihydroxyphenylglycine increases the frequency and decreases the irregularity of the respiratory rhythm. Based on previous studies, we hypothesized that mGluR signaling decreases the irregularity of the respiratory rhythm by activating transient receptor potential canonical (TRPC) channels, which carry a non-specific cation current (ICAN). Indeed, 3,5-dihydroxyphenylglycine (DHPG) application reduces cycle-by-cycle variability and subsequent application of the TRPC channel blocker 1-[2-(4-methoxyphenyl)-2-[3-(4-methoxyphenyl)propoxy]ethyl]imidazole (SKF-96365) hydrochloride reverses this effect. Our data suggest that mGluR5 activation of ICAN-carrying TRPC channels plays an important role in governing the cycle-by-cycle variability of the respiratory rhythm. These data suggest that modulation of TRPC channels may correct irregular respiratory rhythms in some central neuronal diseases.

Analysis of stability and bifurcations of fixed points and periodic solutions of a lumped model of neocortex with two delays

A lumped model of neural activity in neocortex is studied to identify regions of multi-stability of both steady states and periodic solutions. Presence of both steady states and periodic solutions is considered to correspond with epileptogenesis. The model, which consists of two delay differential equations with two fixed time lags is mainly studied for its dependency on varying connection strength between populations. Equilibria are identified, and using linear stability analysis, all transitions are determined under which both trivial and non-trivial fixed points lose stability. Periodic solutions arising at some of these bifurcations are numerically studied with a two-parameter bifurcation analysis.

Rhythmic intrinsic bursting neurons in human neocortex obtained from pediatric patients with epilepsy

Neocortical oscillations result from synchronized activity of a synaptically coupled network and can be strongly influenced by the intrinsic firing properties of individual neurons. As such, the intrinsic electroresponsive properties of individual neurons may have important implications for overall network function. Rhythmic intrinsic bursting (rIB) neurons are of particular interest, as they are poised to initiate and/or strongly influence network oscillations. Although neocortical rIB neurons have been recognized in multiple species, the current study is the first to identify and characterize rIB neurons in the human neocortex. Using whole-cell current-clamp recordings, rIB neurons (n = 12) are identified in human neocortical tissue resected from pediatric patients with intractable epilepsy. In contrast to human regular spiking neurons (n = 12), human rIB neurons exhibit rhythmic bursts of action potentials at frequencies of 0.1-4 Hz. These bursts persist after blockade of fast excitatory neurotransmission and voltage-gated calcium channels. However, bursting is eliminated by subsequent application of the persistent sodium current (I(NaP)) blocker, riluzole. In the presence of riluzole (either 10 or 20 μm), human rIB neurons no longer burst, but fire tonically like regular spiking neurons. These data demonstrate that I(NaP) plays a critical role in intrinsic oscillatory activity observed in rIB neurons in the human neocortex. It is hypothesized that aberrant changes in I(NaP) expression and/or function may ultimately contribute to neurological diseases that are linked to abnormal network activity, such as epilepsy.

Network reconfiguration and neuronal plasticity in rhythm-generating networks

Neuronal networks are highly plastic and reconfigure in a state-dependent manner. The plasticity at the network level emerges through multiple intrinsic and synaptic membrane properties that imbue neurons and their interactions with numerous nonlinear properties. These properties are continuously regulated by neuromodulators and homeostatic mechanisms that are critical to maintain not only network stability and also adapt networks in a short- and long-term manner to changes in behavioral, developmental, metabolic, and environmental conditions. This review provides concrete examples from neuronal networks in invertebrates and vertebrates, and illustrates that the concepts and rules that govern neuronal networks and behaviors are universal.

From neuron to behavior: dynamic equation-based prediction of biological processes in motor control

This article presents the use of continuous dynamic models in the form of differential equations to describe and predict temporal changes in biological processes and discusses several of its important advantages over discontinuous bistable ones, exemplified on the stick insect walking system. In this system, coordinated locomotion is produced by concerted joint dynamics and interactions on different dynamical scales, which is therefore difficult to understand. Modeling using differential equations possesses, in general, the potential for the inclusion of biological detail, the suitability for simulation, and most importantly, parameter manipulation to make predictions about the system behavior. We will show in this review article how, in case of the stick insect walking system, continuous dynamical system models can help to understand coordinated locomotion.

N-methyl-D-aspartate-induced oscillatory properties in neocortical pyramidal neurons from patients with epilepsy

N-Methyl-D-aspartate (NMDA) receptors have been implicated in epileptogenesis, but how these receptors contribute to epilepsy remains unknown. In particular, their role is likely to be complicated because of their voltage-dependent behavior. Here, the authors investigate how activation of NMDA receptors can affect the intrinsic production of oscillation and the resonance properties of neocortical pyramidal neurons from children with intractable epilepsy. Intracellular whole-cell patch clamp recordings in cortical slices from these patients revealed that pyramidal neurons do not produce spontaneous oscillation under control conditions. However, they did exhibit resonance around 1.5 Hz. On NMDA receptor activation, with bath-applied NMDA (10 μM), the majority of neurons produced voltage-dependent intrinsic oscillation associated with a change in the stability of the neuronal system as reflected by the whole-cell I-V curve. Furthermore, the degree of resonance was amplified while the frequency of resonance was shifted to lower frequencies (∼1 Hz) in NMDA. These results suggest that NMDA receptors may both promote the production of low-frequency oscillation and sharpen the response of the cell to lower frequencies. Both these behaviors may be amplified in tissue from patients with epilepsy, resulting in an increased propensity to generate seizures.

Neuronal bursting properties in focal and parafocal regions in pediatric neocortical epilepsy stratified by histology

To test the hypothesis that focal and parafocal neocortical tissue from pediatric patients with intractable epilepsy exhibits cellular and synaptic differences, the authors characterized the propensity of these neurons to generate (a) voltage-dependent bursting and (b) synaptically driven paroxysmal depolarization shifts. Neocortical slices were prepared from tissue resected from patients with intractable epilepsy. Multiunit network activity and simultaneous whole-cell patch recordings were made from neurons from three patient groups: (1) those with normal histology; (2) those with mild and severe cortical dysplasia; and (3) those with abnormal pathology but without cortical dysplasia. Seizure-like activity was characterized by population bursting with concomitant bursting in intracellularly recorded cortical neurons (n = 59). The authors found significantly more N-methyl-D-aspartic acid-driven voltage-dependent bursting neurons in focal versus parafocal tissue in patients with severe cortical dysplasia (P < 0.01). Occurrence of paroxysmal depolarization shifts and burst amplitude and burst duration were significantly related to tissue type: focal or parafocal (P < 0.05). The authors show that functional differences between focal and parafocal tissue in patients with severe cortical dysplasia exist. There are functional differences between patient groups with different histology, and bursting properties can be significantly associated with the distinction between focal and parafocal tissue.

The role of spiking and bursting pacemakers in the neuronal control of breathing

Breathing is controlled by a distributed network involving areas in the neocortex, cerebellum, pons, medulla, spinal cord, and various other subcortical regions. However, only one area seems to be essential and sufficient for generating the respiratory rhythm: the preBötzinger complex (preBötC). Lesioning this area abolishes breathing and following isolation in a brain slice the preBötC continues to generate different forms of respiratory activities. The use of slice preparations led to a thorough understanding of the cellular mechanisms that underlie the generation of inspiratory activity within this network. Two types of inward currents, the persistent sodium current (I(NaP)) and the calcium-activated non-specific cation current (I(CAN)), play important roles in respiratory rhythm generation. These currents give rise to autonomous pacemaker activity within respiratory neurons, leading to the generation of intrinsic spiking and bursting activity. These membrane properties amplify as well as activate synaptic mechanisms that are critical for the initiation and maintenance of inspiratory activity. In this review, we describe the dynamic interplay between synaptic and intrinsic membrane properties in the generation of the respiratory rhythm and we relate these mechanisms to rhythm generating networks involved in other behaviors.

Chapter 3--networks within networks: the neuronal control of breathing

Breathing emerges through complex network interactions involving neurons distributed throughout the nervous system. The respiratory rhythm generating network is composed of micro networks functioning within larger networks to generate distinct rhythms and patterns that characterize breathing. The pre-Bötzinger complex, a rhythm generating network located within the ventrolateral medulla assumes a core function without which respiratory rhythm generation and breathing cease altogether. It contains subnetworks with distinct synaptic and intrinsic membrane properties that give rise to different types of respiratory rhythmic activities including eupneic, sigh, and gasping activities. While critical aspects of these rhythmic activities are preserved when isolated in in vitro preparations, the pre-Bötzinger complex functions in the behaving animal as part of a larger network that receives important inputs from areas such as the pons and parafacial nucleus. The respiratory network is also an integrator of modulatory and sensory inputs that imbue the network with the important ability to adapt to changes in the behavioral, metabolic, and developmental conditions of the organism. This review summarizes our current understanding of these interactions and relates the emerging concepts to insights gained in other rhythm generating networks.

Slow intrinsic oscillations in thick neocortical slices of hypoxia tolerant deep diving seals

Direct evidence that the mammalian neocortex is an important generator of intrinsic activity comes from isolated neocortical slices that spontaneously generate multiple rhythms including those in the beta, delta and gamma range. These oscillations are also seen in intact animals where they interact with other areas including the hippocampus, thalamus and basal ganglia. Here we show that thick isolated neocortical slices from hooded seals intrinsically generate persistent spontaneous activities, both repetitive non-rhythmic activity with activity states lasting for several minutes, and oscillating activity with rhythms that are much slower (<0.1 Hz) than the rhythms previously described in vitro. These intrinsic activities were very robust and persisted for up to 1 h even in severely hypoxic conditions. We hypothesize that the remarkable hypoxia tolerance of the hooded seal nervous system made it possible to maintain functional integrity in slices thick enough to preserve intact neuronal networks capable of generating these slow oscillations. The observed activities in seal neocortical slices support the notion that mammalian cortical networks intrinsically generate multiple states of activity that include oscillatory activity all the way down to <0.1 Hz. This intrinsic neocortical excitability is an important contributor not only to sleep but also to the default awake state of the neocortex.

Systematic comparison of the behaviors produced by computational models of epileptic neocortex

Two existing models of brain dynamics in epilepsy, one detailed (i.e., realistic) and one abstract (i.e., simplified) are compared in terms of behavioral range and match to in vitro mouse recordings. A new method is introduced for comparing across computational models that may have very different forms. First, high-level metrics were extracted from model and in vitro output time series. A principal components analysis was then performed over these metrics to obtain a reduced set of derived features. These features define a low-dimensional behavior space in which quantitative measures of behavioral range and degree of match to real data can be obtained. The detailed and abstract models and the mouse recordings overlapped considerably in behavior space. Both the range of behaviors and similarity to mouse data were similar between the detailed and abstract models. When no high-level metrics were used and principal components analysis was computed over raw time series, the models overlapped minimally with the mouse recordings. The method introduced here is suitable for comparing across different kinds of model data and across real brain recordings. It appears that, despite differences in form and computational expense, detailed and abstract models do not necessarily differ in their behaviors.

Prostaglandin E2-induced synaptic plasticity in neocortical networks of organotypic slice cultures

Traumatic brain injury (TBI) is a major cause of epilepsy, yet the mechanisms underlying the progression from TBI to epilepsy are unknown. TBI induces the expression of COX-2 (cyclooxygenase-2) and increases levels of prostaglandin E2 (PGE2). Here, we demonstrate that acutely applied PGE2 (2 mum) decreases neocortical network activity by postsynaptically reducing excitatory synaptic transmission in acute and organotypic neocortical slices of mice. In contrast, long-term exposure to PGE2 (2 mum; 48 h) presynaptically increases excitatory synaptic transmission, leading to a hyperexcitable network state that is characterized by the generation of paroxysmal depolarization shifts (PDSs). PDSs were also evoked as a result of depriving organotypic slices of activity by treating them with tetrodotoxin (TTX, 1 mum; 48 h). This treatment predominantly increased postsynaptically excitatory synaptic transmission. The network and cellular effects of PGE2 and TTX treatments reversed within 1 week. Differences in the underlying mechanisms (presynaptic vs postsynaptic) as well as occlusion experiments in which slices were exposed to TTX plus PGE2 suggest that the two substances evoke distinct forms of homeostatic plasticity, both of which result in a hyperexcitable network state. PGE2 and TTX (alone or together with PGE2) also increased levels of apoptotic cell death in organotypic slices. Thus, we hypothesize that the increase in excitability and apoptosis may constitute the first steps in a cascade of events that eventually lead to epileptogenesis triggered by TBI.

Persistent sodium current and its role in epilepsy

Sodium currents are essential for the initiation and propagation of neuronal firing. Alterations of sodium currents can lead to abnormal neuronal activity, such as occurs in epilepsy. The transient voltage-gated sodium current mediates the upstroke of the action potential. A small fraction of sodium current, termed the persistent sodium current (I(NaP)), fails to inactivate significantly, even with prolonged depolarization. I(NaP) is activated in the subthreshold voltage range and is capable of amplifying a neuron's response to synaptic input and enhancing its repetitive firing capability. A burgeoning literature is documenting mutations in sodium channels that underlie human disease, including epilepsy. Some of these mutations lead to altered neuronal excitability by increasing I(NaP). This review focuses on the pathophysiological effects of I(NaP) in epilepsy.

Modulation of intrinsic spiking in spinal cord neurons

The vertebrate spinal cord is equipped with a number of neuronal networks that underlie repetitive patterns of behavior as locomotion. Activity in such networks is mediated not only by intrinsic cellular properties but also by synaptic coupling. In this study, we focused on the modulation of the intrinsic activity by 5-hydroxytryptamine (5-HT, serotonin) and the cholinergic agonist muscarine in spinal cord cultures (embryonic age 14 rats). We investigated theses cultures (slices and dissociated cells) at the network level using multielectrode arrays (MEAs) and at the cellular level using whole cell patch clamp. All cultures showed bursting network activity and intrinsic activity when gamma-aminobutyric acid, glycine, and glutamate transmission was blocked. Using MEAs, we observed an increase of the intrinsic activity in the ventral part of the slices with 5-HT and muscarine. In single-cell recordings we found that 43 and 35% of the cells that were silent in the absence of fast synaptic activity were transformed into intrinsically spiking cells by 5-HT and muscarine, respectively. We tested the hypothesis that these neuromodulators act via modulation of the persistent sodium currents (I(NaP)) in these neurons. We found that 5-HT increased threefold the amplitude of I(NaP), specifically in the nonintrinsically spiking cells, and thus switched these cells into intrinsically spiking cells via activation of 5-HT(2) receptor and the phospholipase C pathway. In contrast, the effect of muscarine on nonintrinsically spiking neurons seems to be independent of I(NaP). We conclude from these findings that serotoninergic and cholinergic modulation can turn silent into spontaneously spiking neurons and thus initiate new sources of activity for rhythm generation in spinal networks.

A generalized linear integrate-and-fire neural model produces diverse spiking behaviors

For simulations of neural networks, there is a trade-off between the size of the network that can be simulated and the complexity of the model used for individual neurons. In this study, we describe a generalization of the leaky integrate-and-fire model that produces a wide variety of spiking behaviors while still being analytically solvable between firings. For different parameter values, the model produces spiking or bursting, tonic, phasic or adapting responses, depolarizing or hyperpolarizing after potentials and so forth. The model consists of a diagonalizable set of linear differential equations describing the time evolution of membrane potential, a variable threshold, and an arbitrary number of firing-induced currents. Each of these variables is modified by an update rule when the potential reaches threshold. The variables used are intuitive and have biological significance. The model's rich behavior does not come from the differential equations, which are linear, but rather from complex update rules. This single-neuron model can be implemented using algorithms similar to the standard integrate-and-fire model. It is a natural match with event-driven algorithms for which the firing times are obtained as a solution of a polynomial equation.

Processing noxious information at the subnucleus reticularis dorsalis (SRD) of anesthetized cats: wind-up mechanisms

With the exception of one monkey's study, where wind-up was not reported, electrophysiological data from SRD neurons were obtained in rodents where they show wind-up. This work was designed to examine the response properties of SRD neurons in anesthetized cats to study how general the data from rats may be. Since cat's SRD cells showed wind-up, its underlying mechanisms were approached, an issue not previously addressed at supraspinal level. Electrical stimulation, extracellular (combined with microiontophoresis) and intracellular techniques revealed that A delta information reaches the SRD via the ventrolateral cord, whereas C information preferentially follows a dorsal route. Wind-up was usually generated by spinal and peripheral stimulation, but it was also evoked either by stimulating the nucleus reticularis gigantocellularis (NRGc), even after spinal cord section and bilateral full thickness removal of the cerebral cortex, or by applying microiontophoretic pulses of l-glutamate at 0.3-1 Hz. Wind-up relied on afferent repetitive activity gradually depolarizing the SRD neurons leading 3-4.5 Hz subthreshold membrane rhythmic activity to threshold. Riluzole retarded wind-up generation and decreased the number of spikes per stimulus during wind-up. GABA or glycine abolished spontaneous and sensory-evoked activity and bicuculline, but not strychnine, increased spontaneous and stimulus-evoked activity. These results demonstrate that wind-up at the SRD is not merely the reflection of spinal wind-up, but (i) can be locally generated, (ii) is partially dependent upon persistent sodium currents, and (iii) is under the modulation of a tonic GABAa-dependent inhibition decreasing SRD excitability. Injury and/or inflammation producing tonic C-fiber activation will surpass tonic inhibition generating wind-up.

The persistent sodium current generates pacemaker activities in the central pattern generator for locomotion and regulates the locomotor rhythm

Rhythm generation in neuronal networks relies on synaptic interactions and pacemaker properties. Little is known about the contribution of the latter mechanisms to the integrated network activity underlying locomotion in mammals. We tested the hypothesis that the persistent sodium current (I(NaP)) is critical in generating locomotion in neonatal rodents using both slice and isolated spinal cord preparations. After removing extracellular calcium, 75% of interneurons in the area of the central pattern generator (CPG) for locomotion exhibited bursting properties and I(NaP) was concomitantly upregulated. Putative CPG interneurons such as commissural and Hb9 interneurons also expressed I(NaP)-dependent (riluzole-sensitive) bursting properties. Most bursting cells exhibited a pacemaker-like behavior (i.e., burst frequency increased with depolarizing currents). Veratridine upregulated I(NaP), induced riluzole-sensitive bursting properties, and slowed down the locomotor rhythm. This study provides evidence that I(NaP) generates pacemaker activities in CPG interneurons and contributes to the regulation of the locomotor activity.

propagation of seizure-like activity in a model of neocortex

SUMMARY

Seizures in pediatric epilepsy are often associated with spreading, repetitive bursting activity in neocortex. The authors examined onset and propagation of seizure-like activity using a computational model of cortical circuitry. The model includes two pyramidal cell types and four types of inhibitory interneurons. Each neuron is represented by a multicompartmental model with biophysically realistic ion channels. The authors determined the role of bursting neurons and found that their capability of driving network oscillations is most prominent in networks with either weak or relatively strong excitatory synaptic coupling. Synaptic coupling strength was varied in a separate set of simulations to examine its role in network bursting. Oscillations both between cortical layers (vertical oscillations) and between cortical areas (horizontal oscillations) emerge at moderate excitatory coupling strengths. For horizontal propagation, existence of a fast-conducting fiber system and its properties are critical. Seizure-like oscillatory activity may originate from single neurons or small networks, and that activity may propagate in two principal fashions: one that can be represented by a unidirectional (pacemaker)-type process and the other as multi- or bidirectional propagating waves. The frequency of the bursting patterns relates to underlying propagating activity that can either sustain or disrupt the ongoing oscillation.

Sodium currents in mesencephalic trigeminal neurons from Nav1.6 null mice

Previous studies using pharmacological methods suggest that subthreshold sodium currents are critical for rhythmical burst generation in mesencephalic trigeminal neurons (Mes V). In this study, we characterized transient (I(NaT)), persistent (I(N)(aP)), and resurgent (I(res)) sodium currents in Na(v)1.6-null mice (med mouse, Na(v)1.6(-/-)) lacking expression of the sodium channel gene Scn8a. We found that peak transient, persistent, and resurgent sodium currents from med (Na(v)1.6(-/-)) mice were reduced by 18, 39, and 76% relative to their wild-type (Na(v)1.6(+/+)) littermates, respectively. Current clamp recordings indicated that, in response to sinusoidal constant amplitude current (ZAP function), all neurons exhibited membrane resonance. However, Mes V neurons from med mice had reduced peak amplitudes in the impedance-frequency relationship (resonant Q-value) and attenuated subthreshold oscillations despite the similar passive membrane properties compared with wild-type littermates. The spike frequency-current relationship exhibited reduced instantaneous discharge frequencies and spike block at low stimulus currents and seldom showed maintained spike discharge throughout the stimulus in the majority of med neurons compared with wild-type neurons. Importantly, med neurons never exhibited maintained stimulus-induced rhythmical burst discharge unlike those of wild-type littermates. The data showed that subthreshold sodium currents are critical determinants of Mes V electrogenesis and burst generation and suggest a role for resurgent sodium currents in control of spike discharge.