Search

Pain-related sensory input is processed in the spinal dorsal horn (SDH) before being relayed to the brain. That processing profoundly influences whether stimuli are correctly or incorrectly perceived as painful. Significant advances have been made in identifying the types of excitatory and inhibitory neurons that comprise the SDH, and there is some information about how neuron types are connected, but it remains unclear how the overall circuit processes sensory input or how that processing is disrupted under chronic pain conditions. To explore SDH function, we developed a computational model of the circuit that is tightly constrained by experimental data. Our model comprises conductance-based neuron models that reproduce the characteristic firing patterns of spinal neurons. Excitatory and inhibitory neuron populations, defined by their expression of genetic markers, spiking pattern, or morphology, were synaptically connected according to available qualitative data. Using a genetic algorithm, synaptic weights were tuned to reproduce projection neuron firing rates (model output) based on primary afferent firing rates (model input) across a range of mechanical stimulus intensities. Disparate synaptic weight combinations could produce equivalent circuit function, revealing degeneracy that may underlie heterogeneous responses of different circuits to perturbations or pathological insults. To validate our model, we verified that it responded to reduction of inhibition (i.e. disinhibition) and ablation of specific neuron types in a manner consistent with experiments. Thus validated, our model offers a valuable resource for interpreting experimental results and testing hypotheses in silico to plan experiments for examining normal and pathological SDH circuit function.Significance StatementWe developed a multiscale computer model of the posterior part of spinal cord gray matter (spinal dorsal horn), involved in perception of touch and pain. The model reproduces several experimental observations and makes predictions about how specific types of spinal neurons and synapses influence projection neurons that send information to the brain. Misfiring of these projection neurons can produce anomalous sensations associated with chronic pain. Our computer model will not only assist in planning future experiments, but will also be useful for developing new pharmacotherapy for chronic pain disorders, connecting the effect of drugs acting at the molecular scale with emergent properties of neurons and circuits that shape the pain experience.

The lionfish (Pterois volitans) is a venomous invasive species found in the Caribbean and Northwestern Atlantic. It poses a growing health problem because of the increase in frequency of painful stings, for which no treatment or antidote exists, and the long-term disability caused by the pain. Understanding the venom's algogenic properties can help identify better treatment for these envenomations. In this study, we provide the first characterization of the pain and inflammation caused by lionfish venom and examine the mechanisms through which it causes pain using a combination of in vivo and in vitro approaches including behavioral, physiological, calcium imaging, and electrophysiological testing. Intraplantar injections of the venom produce a significant increase in pain behavior, as well as a marked increase in mechanical sensitivity for up to 24 hours after injection. The algogenic substance(s) are heat-labile peptides that cause neurogenic inflammation at the site of injection and induction of Fos and microglia activation in the superficial layers of the dorsal horn. Finally, calcium imaging and electrophysiology experiments show that the venom acts predominantly on nonpeptidergic, TRPV1-negative, nociceptors, a subset of neurons implicated in sensing mechanical pain. These data provide the first characterization of the pain and inflammation caused by lionfish venom, as well as the first insight into its possible cellular mechanism of action.

KEY POINTS Primary somatosensory neurons normally respond to somatic depolarization with transient spiking but can switch to repetitive spiking under pathological conditions. This switch in spiking pattern reflects a qualitative change in spike initiation dynamics and contributes to the hyperexcitability associated with chronic pain. Neurons can be converted to repetitive spiking by adding a virtual conductance using dynamic clamp. By titrating the conductance to determine how much must be added to cause repetitive spiking, we found that small cells are more susceptible to switching (i.e. required less added conductance) than medium-large cells. By measuring how much less conductance is required to cause repetitive spiking when dynamic clamp was combined with other pathomimetic manipulations (e.g. application of inflammatory mediators), we measured how much each manipulation facilitated repetitive spiking. Our results suggest that many pathological factors facilitate repetitive spiking but that the switch to repetitive spiking requires the cumulative effect of many co-occurring factors. ABSTRACT Primary somatosensory neurons become hyperexcitable in many chronic pain conditions. Hyperexcitability can include a switch from transient to repetitive spiking during sustained somatic depolarization. This switch results from diverse pathological processes that impact ion channel expression or function. Because multiple pathological processes co-occur, isolating how much each contributes to switching the spiking pattern is difficult. Our approach to this challenge involves adding a virtual sodium conductance via dynamic clamp. The magnitude of that conductance was titrated to determine the minimum required to enable rheobasic stimulation to evoke repetitive spiking. The minimum required conductance, termed g¯ Na ∗, was re-measured before and during manipulations designed to model various pathological processes in vitro. The reduction in g¯ Na ∗ caused by each pathomimetic manipulation reflects how much the modelled process contributes to switching the spiking pattern. We found that elevating extracellular potassium or applying inflammatory mediators reduced g¯ Na ∗ whereas direct hyperpolarization had no effect. Inflammatory mediators reduced g¯ Na ∗ more in medium-large (>30 μm diameter) neurons than in small (⩽30 μm diameter) neurons, but had equivalent effects in cutaneous and muscle afferents. The repetitive spiking induced by dynamic clamp was also found to differ between small and medium-large neurons, thus revealing latent differences in adaptation. Our study demonstrates a novel way to determine to what extent individual pathological factors facilitate repetitive spiking. Our results suggest that most factors facilitate but do not cause repetitive spiking on their own, and, therefore, that a switch to repetitive spiking results from the cumulative effect of many co-occurring factors.

Primary afferents transduce environmental stimuli into electrical activity that is transmitted centrally to be decoded into corresponding sensations. However, it remains unknown how afferent populations encode different somatosensory inputs. To address this, we performed two-photon Ca2+ imaging from thousands of dorsal root ganglion (DRG) neurons in anesthetized mice while applying mechanical and thermal stimuli to hind paws. We found that approximately half of all neurons are polymodal and that heat and cold are encoded very differently. As temperature increases, more heating-sensitive neurons are activated, and most individual neurons respond more strongly, consistent with graded coding at population and single-neuron levels, respectively. In contrast, most cooling-sensitive neurons respond in an ungraded fashion, inconsistent with graded coding and suggesting combinatorial coding, based on which neurons are co-activated. Although individual neurons may respond to multiple stimuli, our results show that different stimuli activate distinct combinations of diversely tuned neurons, enabling rich population-level coding.

Background Kilohertz-frequency electric field stimulation (kEFS) applied to the spinal cord can reduce chronic pain without causing the buzzing sensation (paresthesia) associated with activation of dorsal column fibers. This suggests that high-rate spinal cord stimulation (SCS) has a mode of action distinct from conventional, parasthesia-based SCS. A recent study reported that kEFS hyperpolarizes spinal neurons, yet this potentially transformative mode of action contradicts previous evidence that kEFS induces depolarization and was based on patch clamp recordings whose accuracy in the presence of kEFS has not been verified. Objectives We sought to elucidate the basis for kEFS-induced hyperpolarization and to validate the effects of kEFS observed in patch clamp recordings by comparing with independent optical methods. Methods Using patch clamp electrophysiology and voltage-sensitive dye (VSD) imaging, we measured the response to kEFS applied in vitro to hippocampal and spinal neurons. Results The kEFS-induced hyperpolarization observed with current clamp recordings was corroborated by VSD imaging and rheobase measurements in patched neurons. However, no hyperpolarization was observed when imaging unpatched neurons or when recording with a voltage-follower amplifier rather than with a patch clamp amplifier (PCA). We found that EFS induced an artifactual current in PCAs that was injected back into current clamped neurons. The artifactual current induced by single, charge-balanced EFS pulses caused modest hyperpolarization, but these unitary hyperpolarizations accumulated when EFS pulses were repeated at kilohertz frequencies. Conclusion Our results rule out hyperpolarization as the mechanism underlying kEFS-mediated analgesia and highlight the risk of recording artifacts caused by extracellular electrical stimulation.

Key points Distinct spiking patterns may arise from qualitative differences in ion channel expression (i.e. when different neurons express distinct ion channels) and/or when quantitative differences in expression levels qualitatively alter the spike generation process. We hypothesized that spiking patterns in neurons of the superficial dorsal horn (SDH) of spinal cord reflect both mechanisms. We reproduced SDH neuron spiking patterns by varying densities of KV 1- and A-type potassium conductances. Plotting the spiking patterns that emerge from different density combinations revealed spiking-pattern regions separated by boundaries (bifurcations). This map suggests that certain spiking pattern combinations occur when the distribution of potassium channel densities straddle boundaries, whereas other spiking patterns reflect distinct patterns of ion channel expression. The former mechanism may explain why certain spiking patterns co-occur in genetically identified neuron types. We also present algorithms to predict spiking pattern proportions from ion channel density distributions, and vice versa. Abstract Neurons are often classified by spiking pattern. Yet, some neurons exhibit distinct patterns under subtly different test conditions, which suggests that they operate near an abrupt transition, or bifurcation. A set of such neurons may exhibit heterogeneous spiking patterns not because of qualitative differences in which ion channels they express, but rather because quantitative differences in expression levels cause neurons to operate on opposite sides of a bifurcation. Neurons in the spinal dorsal horn, for example, respond to somatic current injection with patterns that include tonic, single, gap, delayed and reluctant spiking. It is unclear whether these patterns reflect five cell populations (defined by distinct ion channel expression patterns), heterogeneity within a single population, or some combination thereof. We reproduced all five spiking patterns in a computational model by varying the densities of a low-threshold (KV 1-type) potassium conductance and an inactivating (A-type) potassium conductance and found that single, gap, delayed and reluctant spiking arise when the joint probability distribution of those channel densities spans two intersecting bifurcations that divide the parameter space into quadrants, each associated with a different spiking pattern. Tonic spiking likely arises from a separate distribution of potassium channel densities. These results argue in favour of two cell populations, one characterized by tonic spiking and the other by heterogeneous spiking patterns. We present algorithms to predict spiking pattern proportions based on ion channel density distributions and, conversely, to estimate ion channel density distributions based on spiking pattern proportions. The implications for classifying cells based on spiking pattern are discussed.

Neuropathic pain is a debilitating condition caused by the abnormal processing of somatosensory input. Synaptic inhibition in the spinal dorsal horn plays a key role in that processing. Mechanical allodynia – the misperception of light touch as painful – occurs when inhibition is compromised. Disinhibition is due primarily to chloride dysregulation caused by hypofunction of the potassium-chloride co-transporter KCC2. Here we show, in rats, that excitatory neurons are disproportionately affected. This is not because chloride is differentially dysregulated in excitatory and inhibitory neurons, but, rather, because excitatory neurons rely more heavily on inhibition to counterbalance strong excitation. Receptive fields in both cell types have a center-surround organization but disinhibition unmasks more excitatory input to excitatory neurons. Differences in intrinsic excitability also affect how chloride dysregulation affects spiking. These results deepen understanding of how excitation and inhibition are normally balanced in the spinal dorsal horn, and how their imbalance disrupts somatosensory processing.

A new therapeutic option for treating arthritis pain., Chronic joint pain such as mechanical allodynia is the most debilitating symptom of arthritis, yet effective therapies are lacking. We identify the pannexin-1 (Panx1) channel as a therapeutic target for alleviating mechanical allodynia, a cardinal sign of arthritis. In rats, joint pain caused by intra-articular injection of monosodium iodoacetate (MIA) was associated with spinal adenosine 5′-triphosphate (ATP) release and a microglia-specific up-regulation of P2X7 receptors (P2X7Rs). Blockade of P2X7R or ablation of spinal microglia prevented and reversed mechanical allodynia. P2X7Rs drive Panx1 channel activation, and in rats with mechanical allodynia, Panx1 function was increased in spinal microglia. Specifically, microglial Panx1-mediated release of the proinflammatory cytokine interleukin-1β (IL-1β) induced mechanical allodynia in the MIA-injected hindlimb. Intrathecal administration of the Panx1-blocking peptide 10panx suppressed the aberrant discharge of spinal laminae I-II neurons evoked by innocuous mechanical hindpaw stimulation in arthritic rats. Furthermore, mice with a microglia-specific genetic deletion of Panx1 were protected from developing mechanical allodynia. Treatment with probenecid, a clinically used broad-spectrum Panx1 blocker, resulted in a striking attenuation of MIA-induced mechanical allodynia and normalized responses in the dynamic weight-bearing test, without affecting acute nociception. Probenecid reversal of mechanical allodynia was also observed in rats 13 weeks after anterior cruciate ligament transection, a model of posttraumatic osteoarthritis. Thus, Panx1-targeted therapy is a new mechanistic approach for alleviating joint pain.

Measuring the excitability of individual axons is complicated by the prohibitive difficulty in obtaining intracellular recordings. Here, we present an innovative methodology that enables local excitability to be measured anywhere in a channelrhodopsin (ChR2)-expressing neuron. The approach hinges on activating ChR2 in a spatially and temporally precise manner while recording the resulting spike train from a remote site. We validated this approach in primary afferent neurons (PANs). Initial encoding of somatosensory stimuli relies on transduction of the physical stimulus into a receptor potential and transformation of the receptor potential into a spike train; the transformation process depends on the excitability of the most distal PAN endings but, as explained above, is extraordinarily difficult to study in situ using traditional methods. Using ChR2-based photoactivation, we show 1) that excitability differs between the distal endings and more proximal portions of PAN axons, 2) that the transformation process differs between PANs, and 3) that the transformation process is directly affected by inflammation. Beyond presenting an innovative method by which to study axonal excitability, this study has validated its utility in helping to decipher the earliest stages of somatosensory encoding.

Pain, itch, heat, cold, and touch represent different percepts arising from somatosensory input. How stimuli give rise to these percepts has been debated for over a century. Recent work supports the view that primary afferents are highly specialized to transduce and encode specific stimulus modalities. However, cross-modal interactions (e.g. inhibition or exacerbation of pain by touch) support convergence rather than specificity in central circuits. We outline how peripheral specialization together with central convergence could enable spinal microcircuits to combine inputs from distinctly specialized, co-activated afferents and to modulate the output signals thus formed through computations like normalization. These issues will be discussed alongside recent advances in our understanding of microcircuitry in the superficial dorsal horn.

Neuropathic pain, which arises from damage to the nervous system, is a major unmet clinical challenge. Reversing the neuronal hyperexcitability induced by nerve damage is a logical treatment strategy but has proven frustratingly difficult. Here, we propose a novel explanation for that difficulty. Changes in several different ion channels are individually sufficient to cause hyperexcitability in primary somatosensory neurons. Despite offering multiple drug targets, this scenario is problematic: if multiple sufficient changes are triggered by nerve injury, then no single change is necessary for hyperexcitability. This so-called degeneracy compromises therapeutic interventions because drug effects on any one ion channel can be circumvented by changes occurring in other ion channels. Overcoming degeneracy demands a more integrative approach to drug discovery.

Nociceptive signals originating in the periphery must be transmitted to the brain to evoke pain. Rather than being conveyed unchanged, those signals undergo extensive processing in the spinal dorsal horn. Synaptic inhibition plays a crucial role in that processing. On the one hand, neuropathy and inflammation are associated with reduced spinal inhibition; on the other hand, the hypersensitivity associated with inflammatory and neuropathic pain can be reproduced by blocking inhibition at the spinal level. To understand the consequences of disinhibition and how to therapeutically reverse it, one must understand how synaptic inhibition normally operates. To that end, this chapter will discuss the structure and function of GABAA and glycine receptors together with the role of associated molecules involved in transmitter handling and chloride regulation. Mechanisms by which inhibition modulates cellular excitability will be described. The chapter will end with discussion of how inhibition goes awry under pathological conditions and what the implications are for the treatment of resulting pain.

Lamina I of the spinal dorsal horn plays an important role in processing and relaying nociceptive information to the brain. It comprises physiologically distinct cell types that process information in fundamentally different ways: tonic neurons fire repetitively during stimulation and display prolonged EPSPs, suggesting operation as integrators, whereas single-spike neurons act like coincidence detectors. Using whole-cell recordings from a rat spinal slice preparation, we set out to determine the basis for prolonged EPSPs in tonic cells and the implications for signal processing. Kinetics of synaptic currents could not explain differences in EPSP kinetics. Instead, tonic neurons were found to express a persistent sodium current,INa,P, that amplified and prolonged depolarization in response to brief stimulation. Tonic neurons also expressed a persistent calcium current,ICa,P, that contributed to prolongation but not to amplification. Simulations using NEURON software demonstrated thatINa,Pwas necessary and sufficient to explain amplification, whereasINa,PandICa,Pacted synergistically to prolong depolarization: initial activation of the slower current (ICa,P) depended on the faster current (INa,P) but maintained activation of the faster current likewise depended on the slower current. Additional investigation revealed thatINa,PandICa,Pcould dramatically increase integration time (>30×) and thereby encourage temporal summation but at the expense of spike time precision. Thus, by prolonging subthreshold depolarization, intrinsic inward currents allow tonic neurons in spinal lamina I to specialize as integrators that are optimally suited to encode stimulus intensity.

Modern pain-control theory predicts that a loss of inhibition (disinhibition) in the dorsal horn of the spinal cord is a crucial substrate for chronic pain syndromes. However, the nature of the mechanisms that underlie such disinhibition has remained controversial. Here we present evidence for a novel mechanism of disinhibition following peripheral nerve injury. It involves a trans-synaptic reduction in the expression of the potassium-chloride exporter KCC2, and the consequent disruption of anion homeostasis in neurons of lamina I of the superficial dorsal horn, one of the main spinal nociceptive output pathways. In our experiments, the resulting shift in the transmembrane anion gradient caused normally inhibitory anionic synaptic currents to be excitatory, substantially driving up the net excitability of lamina I neurons. Local blockade or knock-down of the spinal KCC2 exporter in intact rats markedly reduced the nociceptive threshold, confirming that the reported disruption of anion homeostasis in lamina I neurons was sufficient to cause neuropathic pain.

Neurons rely on action potentials, or spikes, to encode information. But spikes can encode different stimulus features in different neurons. We show here through simulations and experiments how neurons encode the integral or derivative of their input based on the distinct tuning properties conferred upon them by subthreshold currents. Slow-activating subthreshold inward (depolarizing) current mediates positive feedback control of subthreshold voltage, sustaining depolarization and allowing the neuron to spike on the basis of its integrated stimulus waveform. Slow-activating subthreshold outward (hyperpolarizing) current mediates negative feedback control of subthreshold voltage, truncating depolarization and forcing the neuron to spike on the basis of its differentiated stimulus waveform. Depending on its direction, slow-activating subthreshold current cooperates or competes with fast-activating inward current during spike initiation. This explanation predicts that sensitivity to the rate of change of stimulus intensity differs qualitatively between integrators and differentiators. This was confirmed experimentally in spinal sensory neurons that naturally behave as specialized integrators or differentiators. Predicted sensitivity to different stimulus features was confirmed by covariance analysis. Integration and differentiation, which are themselves inverse operations, are thus shown to be implemented by the slow feedback mediated by oppositely directed subthreshold currents expressed in different neurons.

The central terminals of primary afferent fibers experience depolarization upon activation of GABAA receptors (GABAAR) because their intracellular chloride concentration is maintained above electrochemical equilibrium. Primary afferent depolarization (PAD) normally mediates inhibition via sodium channel inactivation and shunting but can evoke spikes under certain conditions. Antidromic (centrifugal) conduction of these spikes may contribute to neurogenic inflammation while orthodromic (centripetal) conduction could contribute to pain in the case of nociceptive fibers. PAD-induced spiking is assumed to override presynaptic inhibition. Using computer simulations and dynamic clamp experiments, we sought to identify which biophysical changes are required to enable PAD-induced spiking and whether those changes necessarily compromise PAD-mediated inhibition. According to computational modeling, a depolarizing shift in GABA reversal potential (EGABA) and increased intrinsic excitability (manifest as altered spike initiation properties) were necessary for PAD-induced spiking, whereas increased GABAAR conductance density (ḡGABA) had mixed effects. We tested our predictions experimentally by using dynamic clamp to insert virtual GABAAR conductances with different EGABA and kinetics into acutely dissociated dorsal root ganglion (DRG) neuron somata. Comparable experiments in central axon terminals are prohibitively difficult but the biophysical requirements for PAD-induced spiking are arguably similar in soma and axon. Neurons from naïve (i.e. uninjured) rats were compared before and after pharmacological manipulation of intrinsic excitability, and against neurons from nerve-injured rats. Experimental data confirmed that, in most neurons, both predicted changes were necessary to yield PAD-induced spiking. Importantly, such changes did not prevent PAD from inhibiting other spiking or from blocking spike propagation. In fact, since the high value of ḡGABA required for PAD-induced spiking still mediates strong inhibition, we conclude that PAD-induced spiking does not represent failure of presynaptic inhibition. Instead, diminished PAD caused by reduction of ḡGABA poses a greater risk to presynaptic inhibition and the sensory processing that relies upon it., Author Summary Postsynaptic GABAAR mediate inhibition by causing hyperpolarization or by preventing (shunting) the depolarization caused by concurrent excitatory input. Presynaptic GABAAR work differently, in the spinal cord at least. Because of their higher-than-equilibrium intracellular chloride concentration, the central terminals of primary afferent fibers are depolarized by activation of GABAAR. This so-called primary afferent depolarization, or PAD, nonetheless reduces spike propagation and synaptic release from those fibers because of shunting effects and sodium channel inactivation. But those inhibitory effects can be diminished under certain pathological conditions; in fact, the emergence of dorsal root reflexes suggests that PAD can become paradoxically excitatory. The biophysical basis for this paradoxical excitation has been hinted at by experiments, but here, for the first time, we use computational modeling and dynamic clamp experiments to decipher how distinct contributing factors interact to enable PAD-induced spiking. Our results suggest that PAD-induced spiking requires a shift in GABA reversal potential plus changes in intrinsic excitability that allow for repetitive spiking during sustained depolarization. Inhibitory effects of PAD are retained under conditions in which GABAAR activation causes transient spiking and are only lost if GABAAR activation can evoke repetitive spiking.

Correlated spiking has been widely observed, but its impact on neural coding remains controversial. Correlation arising from comodulation of rates across neurons has been shown to vary with the firing rates of individual neurons. This translates into rate and correlation being equivalently tuned to the stimulus; under those conditions, correlated spiking does not provide information beyond that already available from individual neuron firing rates. Such correlations are irrelevant and can reduce coding efficiency by introducing redundancy. Using simulations and experiments in rat hippocampal neurons, we show here that pairs of neurons receiving correlated input also exhibit correlations arising from precise spike-time synchronization. Contrary to rate comodulation, spike-time synchronization is unaffected by firing rate, thus enabling synchrony- and rate-based coding to operate independently. The type of output correlation depends on whether intrinsic neuron properties promote integration or coincidence detection: “ideal” integrators (with spike generation sensitive to stimulus mean) exhibit rate comodulation, whereas ideal coincidence detectors (with spike generation sensitive to stimulus variance) exhibit precise spike-time synchronization. Pyramidal neurons are sensitive to both stimulus mean and variance, and thus exhibit both types of output correlation proportioned according to which operating mode is dominant. Our results explain how different types of correlations arise based on how individual neurons generate spikes, and why spike-time synchronization and rate comodulation can encode different stimulus properties. Our results also highlight the importance of neuronal properties for population-level coding insofar as neural networks can employ different coding schemes depending on the dominant operating mode of their constituent neurons.

Fast axonal conduction of action potentials in mammals relies on myelin insulation. Demyelination can cause slowed, blocked, desynchronized, or paradoxically excessive spiking that underlies the symptoms observed in demyelination diseases. The diversity and timing of such symptoms are poorly understood, often intermittent, and uncorrelated with disease progress. We modeled the effects of demyelination (and secondary remodeling) on intrinsic axonal excitability using Hodgkin–Huxley and reduced Morris–Lecar models. Simulations and analysis suggested a simple explanation for the breadth of symptoms and revealed that the ratio of sodium to leak conductance, g Na /g L , acted as a four-way switch controlling excitability patterns that included spike failure, single spike transmission, afterdischarge, and spontaneous spiking. Failure occurred when this ratio fell below a threshold value. Afterdischarge occurred at g Na /g L just below the threshold for spontaneous spiking and required a slow inward current that allowed for two stable attractor states, one corresponding to quiescence and the other to repetitive spiking. A neuron prone to afterdischarge could function normally unless it was switched to its “pathological” attractor state; thus, although the underlying pathology may develop slowly by continuous changes in membrane conductances, a discontinuous change in axonal excitability can occur and lead to paroxysmal symptoms. We conclude that tonic and paroxysmal positive symptoms as well as negative symptoms may be a consequence of varying degrees of imbalance between g Na and g L after demyelination. The KCNK family of g L potassium channels may be an important target for new drugs to treat the symptoms of demyelination.

Lamina I of the spinal dorsal horn plays an important role in the processing and relay of nociceptive information. Signal processing depends, in part, on neuronal membrane properties. Intrinsic membrane properties of lamina I neurons were therefore investigated using whole cell patch clamp recordings in a slice preparation of adult rat spinal cord. Based on responses to somatic current injection, four cell types were identified: tonic, which fire comparatively slowly but continuously throughout stimulation; phasic, which fire a high frequency burst of variable duration; delayed onset, which fire irregularly and with a marked delay to the first spike; and single spike, which typically fire only one action potential even when strongly depolarised. Classification by spiking pattern was further refined by identification of characteristic stimulus-response curves and quantification of several response parameters. Objectivity of the classification was confirmed by cluster analysis. Responses to stimulus trains and synaptic input as well as the kinetics of spontaneous synaptic events revealed differences in the signal processing characteristics of the cell types: tonic and delayed onset cells appeared to act predominantly as integrators whereas phasic and single spike cells acted as coincidence detectors. Intracellular labelling revealed a significant correlation between morphological and physiological cell types: tonic cells were typically fusiform, phasic cells were pyramidal, and delayed onset and single spike cells were multipolar. Thus, there are multiple physiological cells types in lamina I with specific morphological correlates and distinctive signal processing characteristics that confer significant differences in the transduction of input into spike trains.

The tentacle withdrawal reflex of the snail Helix aspersa exhibits a complex combination of habituation and sensitization consistent with the dual-process theory of plasticity. Habituation, sensitization, or a combination of both were elicited by varying stimulation parameters and lesion condition. Analysis of response plasticity shows that the late phase of the response is selectively enhanced by sensitization, whereas all phases are decreased by habituation. Previous data have shown that tentacle withdrawal is mediated conjointly by parallel monosynaptic and polysynaptic pathways. The former mediates the early phase, whereas the latter mediates the late phase of the response. Plastic loci were identified by stimulating and recording at different points within the neural circuit, in combination with selective lesions. Results indicate that depression occurs at an upstream locus, before circuit divergence, and is therefore expressed in all pathways, whereas facilitation requires downstream facilitatory neurons and is selectively expressed in polysynaptic pathways. Differential expression of plasticity between pathways helps explain the behavioral manifestation of depression and facilitation. A simple mathematical model is used to show how serial positioning of depression and facilitation can explain the kinetics of dual-process learning. These results illustrate how the position of cellular plasticity in the network affects behavioral change and how forms of plasticity can interact to determine the kinetics of the net changes.

Prescott, Steven A., Nishi Gill, and Ronald Chase. Neural circuit mediating tentacle withdrawal in Helix aspersa, with specific reference to the competence of the motor neuron C3. J. Neurophysiol. 78: 2951–2965, 1997. The tentacle withdrawal reflex in the terrestrial snail Helix aspersa involves bending and retraction of the tentacles. When elicited by mechanical stimulation of the tentacle, the reflex is mediated by the conjoint action of the central and peripheral nervous systems. The neural circuit underlying the stimulus-response pathways was studied in vitro using a combination of morphological and physiological techniques. Sensory input caused by stimulation of the nose (situated at the superior tentacle's tip) first passes into the tentacle ganglion. Motor fibers are likely excited in the tentacle ganglion to form a peripheral stimulus-response pathway. While still in the tentacle ganglion, the excitation caused by a brief stimulus is transformed into a prolonged neuronal discharge. This modified signal travels, via the olfactory nerve, to the cerebral ganglion where it excites the giant motor neuron C3 along with numerous smaller motor neurons. Afferent input to C3 also arrives from several other sources. The afferent convergence is followed by a marked divergence of C3's output. C3 innervates the muscles mediating both tentacle retraction and tentacle bending through multiple cerebral nerves. Thus C3's pattern of effector innervation allows this single cell to elicit and coordinate both components of the tentacle withdrawal reflex. Lesion experiments indicate that C3 is responsible for 85% of the central contribution to tentacle retraction, though C3 is actually sufficient to mediate maximal muscle contraction as evidenced by intracellular stimulation. In addition to C3, three groups of putative central motor neurons were identified through nerve backfills and nerve recordings. The additional motor neurons mediating tentacle retraction are important for maximizing the rate of muscle contraction, whereas those mediating tentacle bending are likely more important for nondefensive behaviors. These neurons are arranged in parallel with C3, but unlike C3, each of these neurons innervates only a single effector or portion thereof. Given C3's direct innervation of multiple effectors and its sufficiency to evoke strong responses in those effectors, we conclude that C3 is paramount in eliciting and coordinating tentacle withdrawal.

Chloride homeostasis is a critical determinant of the strength and robustness of inhibition mediated by GABAA receptors (GABAARs). The impact of changes in steady state Cl− gradient is relatively straightforward to understand, but how dynamic interplay between Cl− influx, diffusion, extrusion and interaction with other ion species affects synaptic signaling remains uncertain. Here we used electrodiffusion modeling to investigate the nonlinear interactions between these processes. Results demonstrate that diffusion is crucial for redistributing intracellular Cl− load on a fast time scale, whereas Cl−extrusion controls steady state levels. Interaction between diffusion and extrusion can result in a somato-dendritic Cl− gradient even when KCC2 is distributed uniformly across the cell. Reducing KCC2 activity led to decreased efficacy of GABAAR-mediated inhibition, but increasing GABAAR input failed to fully compensate for this form of disinhibition because of activity-dependent accumulation of Cl−. Furthermore, if spiking persisted despite the presence of GABAAR input, Cl− accumulation became accelerated because of the large Cl− driving force that occurs during spikes. The resulting positive feedback loop caused catastrophic failure of inhibition. Simulations also revealed other feedback loops, such as competition between Cl− and pH regulation. Several model predictions were tested and confirmed by [Cl−]i imaging experiments. Our study has thus uncovered how Cl− regulation depends on a multiplicity of dynamically interacting mechanisms. Furthermore, the model revealed that enhancing KCC2 activity beyond normal levels did not negatively impact firing frequency or cause overt extracellular K− accumulation, demonstrating that enhancing KCC2 activity is a valid strategy for therapeutic intervention., Author Summary Fast synaptic inhibition relies on chloride current to hyperpolarize the neuron or to prevent depolarization caused by concurrent excitatory input. Both scenarios necessarily involve chloride flux into the cell and, thus, a change in intracellular chloride concentration. The vast majority of models neglect changes in ion concentration despite experimental evidence that such changes occur and are not inconsequential. The importance of considering chloride homeostasis mechanisms is heightened by evidence that several neurological diseases are associated with deficient chloride extrusion capacity. Steady state chloride levels are altered in those disease states. Fast chloride dynamics are also likely affected, but those changes have yet to be explored. To this end, we built an electrodiffusion model that tracks changes in the concentration of chloride plus multiple other ion species. Simulations in this model revealed a multitude of complex, nonlinear interactions that have important consequences for the efficacy of synaptic inhibition. Several predictions from the model were tested and confirmed with chloride imaging experiments.

Multiple sclerosis and other demyelination diseases are associated with positive (gain of function), negative (loss of function), and paroxysmal (sudden and short duration) symptoms. These symptoms occur in a remarkable diversity of modalities including prolonged or, in the case of paroxysmal symptoms, brief episodes of pain, paresthesia (numbness), dysarthria (difficulty speaking), weakness, muscle spasms, etc. [1]. Using computational modeling, we have recently established the mechanisms underlying four distinct phases of action potential (AP) firing associated with demyelination and secondary compensatory changes [2]. Transitions in axonal excitability caused by varying the ratio of gNa and gL indicate that a single mechanism controlling conductance balance is sufficient to explain normal spiking, failed conduction, afterdischarge (AD) and spontaneous activity, and can thus account for the full range of negative, paroxysmal and tonic positive symptoms, respectively. Positive symptoms are thought to arise from ectopic spike discharge, but the biophysical mechanisms responsible for the abrupt onset and self limiting duration of paroxysmal AD are poorly understood and were thus a particular area of focus. We used two types of complementary models. The first was a Hodgkin-Huxley (HH), multi-compartment model of a focally demyelinated axon created using the NEURON Simulator (http://www.neuron.yale.edu/neuron, [3]) based on a geometrically simplified model to first understand the behaviors of electrical conductances. The second was a modified Morris-Lecar, single compartment model [4] that has been described in detail [5]. We explored this model with phase-plane and bifurcation analysis using XPPAUT in order to rigorously characterize the dynamical mechanisms underlying phenomena identified in the HH model. In both models, the paroxysmal AD can only occur when the system is bi-stable and initiation occurs when a perturbation abruptly switches the system between attractor states. Termination of paroxysms occurs when internal dynamics shift the basins of attraction until the system falls back to its original stable state. Temporal summation and refractoriness can also be simply explained within this framework. Disruption of the oligodendrocytes that form the myelin sheaths may lead to abnormal electrolyte homeostasis in the vicinity of the denuded axon. Our models indicate that the pathological accumulation of Na + or K + can influence the duration of paroxysmal firing by acting as an additional, ultra-slow feedback mechanism in the dynamical system. As such, a naked patch of previously myelinated axonal membrane becomes capable of competing with the same neuron’s soma for the initiation of burst firing patterns and interferes with normal firing patterns. Although not biologically detailed, our modeling nonetheless provides significant insight into how biologically realistic mechanisms can interact to produce neuronal discharges as well as other abnormal firing characteristics in the course of demyelination.