Staircase currents in motoneurons: insight into the spatial arrangement of calcium channels in the dendritic tree

Bifurcation analysis of motoneuronal excitability mechanisms under normal and ALS conditions

Introduction

Bifurcation analysis allows the examination of steady-state, non-linear dynamics of neurons and their effects on cell firing, yet its usage in neuroscience is limited to single-compartment models of highly reduced states. This is primarily due to the difficulty in developing high-fidelity neuronal models with 3D anatomy and multiple ion channels in XPPAUT, the primary bifurcation analysis software in neuroscience.

Methods

To facilitate bifurcation analysis of high-fidelity neuronal models under normal and disease conditions, we developed a multi-compartment model of a spinal motoneuron (MN) in XPPAUT and verified its firing accuracy against its original experimental data and against an anatomically detailed cell model that incorporates known MN non-linear firing mechanisms. We used the new model in XPPAUT to study the effects of somatic and dendritic ion channels on the MN bifurcation diagram under normal conditions and after amyotrophic lateral sclerosis (ALS) cellular changes.

Results

Our results show that somatic small-conductance Ca2+-activated K (SK) channels and dendritic L-type Ca2+ channels have the strongest effects on the bifurcation diagram of MNs under normal conditions. Specifically, somatic SK channels extend the limit cycles and generate a subcritical Hopf bifurcation node in the V-I bifurcation diagram of the MN to replace a supercritical node Hopf node, whereas L-type Ca2+ channels shift the limit cycles to negative currents. In ALS, our results show that dendritic enlargement has opposing effects on MN excitability, has a greater overall impact than somatic enlargement, and dendritic overbranching offsets the dendritic enlargement hyperexcitability effects.

Discussion

Together, the new multi-compartment model developed in XPPAUT facilitates studying neuronal excitability in health and disease using bifurcation analysis.

Maturation of persistent and hyperpolarization-activated inward currents shapes the differential activation of motoneuron subtypes during postnatal development

The size principle underlies the orderly recruitment of motor units; however, motoneuron size is a poor predictor of recruitment amongst functionally defined motoneuron subtypes. Whilst intrinsic properties are key regulators of motoneuron recruitment, the underlying currents involved are not well defined. Whole-cell patch-clamp electrophysiology was deployed to study intrinsic properties, and the underlying currents, that contribute to the differential activation of delayed and immediate firing motoneuron subtypes. Motoneurons were studied during the first three postnatal weeks in mice to identify key properties that contribute to rheobase and may be important to establish orderly recruitment. We find that delayed and immediate firing motoneurons are functionally homogeneous during the first postnatal week and are activated based on size, irrespective of subtype. The rheobase of motoneuron subtypes becomes staggered during the second postnatal week, which coincides with the differential maturation of passive and active properties, particularly persistent inward currents. Rheobase of delayed firing motoneurons increases further in the third postnatal week due to the development of a prominent resting hyperpolarization-activated inward current. Our results suggest that motoneuron recruitment is multifactorial, with recruitment order established during postnatal development through the differential maturation of passive properties and sequential integration of persistent and hyperpolarization-activated inward currents.

Linking Motoneuron PIC Location to Motor Function in Closed-Loop Motor Unit System Including Afferent Feedback: A Computational Investigation

The goal of this study is to investigate how the activation location of persistent inward current (PIC) over motoneuron dendrites is linked to motor output in the closed-loop motor unit. Here, a physiologically realistic model of a motor unit including afferent inputs from muscle spindles was comprehensively analyzed under intracellular stimulation at the soma and synaptic inputs over the dendrites during isometric contractions over a full physiological range of muscle lengths. The motor output of the motor unit model was operationally assessed by evaluating the rate of force development, the degree of force potentiation and the capability of self-sustaining force production. Simulations of the model motor unit demonstrated a tendency for a faster rate of force development, a greater degree of force potentiation, and greater capacity for self-sustaining force production under both somatic and dendritic stimulation of the motoneuron as the PIC channels were positioned farther from the soma along the path of motoneuron dendrites. Interestingly, these effects of PIC activation location on force generation significantly differed among different states of muscle length. The rate of force development and the degree of force potentiation were systematically modulated by the variation of PIC channel location for shorter-than-optimal muscles but not for optimal and longer-than-optimal muscles. Similarly, the warm-up behavior of the motor unit depended on the interplay between PIC channel location and muscle length variation. These results suggest that the location of PIC activation over motoneuron dendrites may be distinctively reflected in the motor performance during shortening muscle contractions.

Calpain fosters the hyperexcitability of motoneurons after spinal cord injury and leads to spasticity

Up-regulation of the persistent sodium current (INaP) and down-regulation of the potassium/chloride extruder KCC2 lead to spasticity after spinal cord injury (SCI). We here identified calpain as the driver of the up- and down-regulation of INaP and KCC2, respectively, in neonatal rat lumbar motoneurons. Few days after SCI, neonatal rats developed behavioral signs of spasticity with the emergence of both hyperreflexia and abnormal involuntary muscle contractions on hindlimbs. At the same time, in vitro isolated lumbar spinal cords became hyperreflexive and displayed numerous spontaneous motor outputs. Calpain-I expression paralleled with a proteolysis of voltage-gated sodium (Nav) channels and KCC2. Acute inhibition of calpains reduced this proteolysis, restored the motoneuronal expression of Nav and KCC2, normalized INaP and KCC2 function, and curtailed spasticity. In sum, by up- and down-regulating INaP and KCC2, the calpain-mediated proteolysis of Nav and KCC2 drives the hyperexcitability of motoneurons which leads to spasticity after SCI.

Sequential activation of multiple persistent inward currents induces staircase currents in serotonergic neurons of medulla in ePet-EYFP mice

Persistent inward currents (PICs) are widely reported in rodent spinal neurons. A distinctive pattern observed recently is staircase-like PICs induced by voltage ramp in serotonergic neurons of mouse medulla. The mechanism underlying this pattern of PICs is unclear. Combining electrophysiological, pharmacological, and computational approaches, we investigated the staircase PICs in serotonergic neurons of medulla in ePet-EYFP transgenic mice (postnatal days 1-7). Staircase PICs induced by 10-s voltage biramps were observed in 70% of serotonergic neurons (n = 73). Staircase PICs activated at -48.8 ± 5 mV and consisted of two components, with the first PIC of 45.8 ± 51 pA and the second PIC of 197.3 ± 126 pA (n = 51). Staircase PICs were also composed of low-voltage-activated sodium PIC (Na-PIC; onset -46.2 ± 5 mV, n = 34), high-voltage-activated calcium PIC (Ca-PIC; onset -29.3 ± 6 mV, n = 23), and high-voltage-activated tetrodotoxin (TTX)- and dihydropyridine-resistant sodium PIC (TDR-PIC; onset -16.8 ± 4 mV, n = 28). Serotonergic neurons expressing Na-PIC, Ca-PIC, and TDR-PIC were evenly distributed in medulla. Bath application of 1-2 μM TTX blocked the first PIC and decreased the second PIC by 36% (n = 23, P < 0.05). Nimodipine (25 μM) reduced the second PIC by 38% (n = 34, P < 0.001) without altering the first PIC. TTX and nimodipine removed the first PIC and reduced the second PIC by 59% (n = 28, P < 0.01). A modeling study mimicked the staircase PICs and verified experimental conclusions that sequential activation of Na-PIC, Ca-PIC, and TDR-PIC in order of voltage thresholds induced staircase PICs in serotonergic neurons. Further experimental results suggested that the multiple components of staircase PICs play functional roles in regulating excitability of serotonergic neurons in medulla.NEW & NOTEWORTHY Staircase persistent inward currents (PICs) are mediated by activation of L-type calcium channels in dendrites of mouse spinal motoneurons. A novel mechanism is explored in this study. Here we report that the staircase PICs are mediated by sequentially activating sodium and calcium PICs in serotonergic neurons of mouse medulla.

Motoneuron output regulated by ionic channels: a modeling study of motoneuron frequency-current relationships during fictive locomotion

Cat lumbar motoneurons display changes in membrane properties during fictive locomotion. These changes include reduction of input resistance and afterhyperpolarization, hyperpolarization of voltage threshold, and voltage-dependent excitation of the motoneurons. The state-dependent alteration of membrane properties leads to dramatic changes in frequency-current (F-I) relationship. The mechanism underlying these changes remains unknown. Using a motoneuron model combined with electrophysiological data, we investigated the channel mechanisms underlying the regulation of motoneuronal excitability and motor output. Simulation results showed that upregulation of transient sodium, persistent sodium, or Ca_v1.3 calcium conductances or downregulation of calcium-activated potassium or KCNQ/K_v7 potassium conductances could increase motoneuronal excitability and motor output through hyperpolarizing (left shifting) the F-I relationships or increasing the F-I slopes, whereas downregulation of input resistance or upregulation of potassium-mediated leak conductance produced the opposite effects. The excitatory phase of locomotor drive potentials (LDPs) also substantially hyperpolarized the F-I relationships and increased the F-I slopes, whereas the inhibitory phase of the LDPs had opposite effects to a similar extent. The simulation results also showed that none of the individual channel modulations could produce all the changes in the F-I relationships. The effects of modulation of Ca_v1.3 and KCNQ/K_v7 on F-I relationships were supported by slice experiments with the Ca_v1.3 agonist Bay K8644 and the KCNQ/K_v7 antagonist XE-991. The conclusion is that the varying changes in F-I relationships during fictive locomotion could be regulated by multichannel modulations. This study provides insight into the ionic basis for control of motor output in walking. NEW & NOTEWORTHY Mammalian spinal motoneurons have their excitability adapted to facilitate recruitment and firing during locomotion. Cat lumbar motoneurons display dramatic changes in membrane properties during fictive locomotion. These changes lead to a varying alteration of frequency-current relationship. The mechanisms underlying the changes remain unknown. In particular, little is known about the ionic basis for regulation of motoneuronal excitability and thus control of the motor output for walking by the spinal motor system.

Persistent Sodium Current Drives Excitability of Immature Renshaw Cells in Early Embryonic Spinal Networks

Spontaneous network activity (SNA) emerges in the spinal cord (SC) before the formation of peripheral sensory inputs and central descending inputs. SNA is characterized by recurrent giant depolarizing potentials (GDPs). Because GDPs in motoneurons (MNs) are mainly evoked by prolonged release of GABA, they likely necessitate sustained firing of interneurons. To address this issue we analyzed, as a model, embryonic Renshaw cell (V1^R) activity at the onset of SNA (E12.5) in the embryonic mouse SC (both sexes). V1^R are one of the interneurons known to contact MNs, which are generated early in the embryonic SC. Here, we show that V1^R already produce GABA in E12.5 embryo, and that V1^R make synaptic-like contacts with MNs and have putative extrasynaptic release sites, while paracrine release of GABA occurs at this developmental stage. In addition, we discovered that V1^R are spontaneously active during SNA and can already generate several intrinsic activity patterns including repetitive-spiking and sodium-dependent plateau potential that rely on the presence of persistent sodium currents (I_Nap). This is the first demonstration that I_Nap is present in the embryonic SC and that this current can control intrinsic activation properties of newborn interneurons in the SC of mammalian embryos. Finally, we found that 5 μm riluzole, which is known to block I_NaP, altered SNA by reducing episode duration and increasing inter-episode interval. Because SNA is essential for neuronal maturation, axon pathfinding, and synaptogenesis, the presence of I_NaP in embryonic SC neurons may play a role in the early development of mammalian locomotor networks.SIGNIFICANCE STATEMENT The developing spinal cord (SC) exhibits spontaneous network activity (SNA) involved in the building of nascent locomotor circuits in the embryo. Many studies suggest that SNA depends on the rhythmic release of GABA, yet intracellular recordings of GABAergic neurons have never been performed at the onset of SNA in the SC. We first discovered that embryonic Renshaw cells (V1^R) are GABAergic at E12.5 and spontaneously active during SNA. We uncover a new role for persistent sodium currents (I_NaP) in driving plateau potential in V1^R and in SNA patterning in the embryonic SC. Our study thus sheds light on a role for INaP in the excitability of V1^R and the developing SC.

Impact of the localization of dendritic calcium persistent inward current on the input-output properties of spinal motoneuron pool: a computational study

The goal of this study is to investigate how the dendritic Ca-PIC location influences nonlinear input-output properties and depends on the type of motoneurons across the motoneuron pool. A model motoneuron pool consisting of 10 motoneurons was constructed using a recently developed two-compartment modeling approach that reflected key cell type-associated properties experimentally identified. The dendritic excitability and firing output depended systematically on both the PIC location and the motoneuron type. The PIC onset and offset in the current-voltage (I-V) relationship tended to occur at more hyperpolarized voltages as the path length to the PIC channels from the soma increased and as the cell type shifted from high- to low-threshold motoneurons. At the same time, the firing acceleration and frequency hysteresis in the frequency-current (F-I) relationship became faster and larger, respectively. However, the PIC onset-offset hysteresis increased as the path length and the recruitment threshold increased. Furthermore, the gain of frequency-current function before full PIC activation was larger for PIC channels located over distal dendritic regions in low- compared with high-threshold motoneurons. When compared with previously published experimental observations, the modeling concurred when Ca-PIC channels were placed closer to the soma in high- than low-threshold motoneurons in the model motoneuron pool. All of these results suggest that the negative relationship of Ca-PIC location and cell recruitment threshold may underlie the systematic variation in I-V and F-I transformation across the motoneuron pool.NEW & NOTEWORTHY How does the dendritic location of calcium persistent inward current (Ca-PIC) influence dendritic excitability and firing behavior across the spinal motoneuron pool? This issue was investigated developing a model motoneuron pool that reflected key motoneuron type-specific properties experimentally identified. The simulation results point out the negative relationship between the distance of Ca-PIC source from the soma and cell recruitment threshold as a basis underlying the systematic variation in input-output properties of motoneurons over the motoneuron pool.

Motor Neurons

Motor neurons translate synaptic input from widely distributed premotor networks into patterns of action potentials that orchestrate motor unit force and motor behavior. Intercalated between the CNS and muscles, motor neurons add to and adjust the final motor command. The identity and functional properties of this facility in the path from synaptic sites to the motor axon is reviewed with emphasis on voltage sensitive ion channels and regulatory metabotropic transmitter pathways. The catalog of the intrinsic response properties, their underlying mechanisms, and regulation obtained from motoneurons in in vitro preparations is far from complete. Nevertheless, a foundation has been provided for pursuing functional significance of intrinsic response properties in motoneurons in vivo during motor behavior at levels from molecules to systems. © 2017 American Physiological Society. Compr Physiol 7:463-484, 2017.

Voltage-gated calcium channels are abnormal in cultured spinal motoneurons in the G93A-SOD1 transgenic mouse model of ALS

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease characterized by progressive loss of motoneurons. Hyperexcitability and excitotoxicity have been implicated in the early pathogenesis of ALS. Studies addressing excitotoxic motoneuron death and intracellular Ca(2+) overload have mostly focused on Ca(2+) influx through AMPA glutamate receptors. However, intrinsic excitability of motoneurons through voltage-gated ion channels may also have a role in the neurodegeneration. In this study we examined the function and localization of voltage-gated Ca(2+) channels in cultured spinal cord motoneurons from mice expressing a mutant form of human superoxide dismutase-1 with a Gly93→Ala substitution (G93A-SOD1). Using whole-cell patch-clamp recordings, we showed that high voltage activated (HVA) Ca(2+) currents are increased in G93A-SOD1 motoneurons, but low voltage activated Ca(2+) currents are not affected. G93A-SOD1 motoneurons also have altered persistent Ca(2+) current mediated by L-type Ca(2+) channels. Quantitative single-cell RT-PCR revealed higher levels of Ca1a, Ca1b, Ca1c, and Ca1e subunit mRNA expression in G93A-SOD1 motoneurons, indicating that the increase of HVA Ca(2+) currents may result from upregulation of Ca(2+) channel mRNA expression in motoneurons. The localizations of the Ca1B N-type and Ca1D L-type Ca(2+) channels in motoneurons were examined by immunocytochemistry and confocal microscopy. G93A-SOD1 motoneurons had increased Ca1B channels on the plasma membrane of soma and dendrites. Ca1D channels are similar on the plasma membrane of soma and lower on the plasma membrane of dendrites of G93A-SOD1 motoneurons. Our study demonstrates that voltage-gated Ca(2+) channels have aberrant functions and localizations in ALS mouse motoneurons. The increased HVA Ca(2+) currents and PCCa current could contribute to early pathogenesis of ALS.

Contribution of intrinsic motoneuron properties to discharge hysteresis and its estimation based on paired motor unit recordings: a simulation study

Motoneuron activity is strongly influenced by the activation of persistent inward currents (PICs) mediated by voltage-gated sodium and calcium channels. However, the amount of PIC contribution to the activation of human motoneurons can only be estimated indirectly. Simultaneous recordings of pairs of motor units have been used to provide an estimate of the PIC contribution by using the firing rate of the lower threshold unit to provide an estimate of the common synaptic drive to both units, and the difference in firing rate (ΔF) of this lower threshold unit at recruitment and de-recruitment of the higher threshold unit to estimate the PIC contribution to activation of the higher threshold unit. It has recently been suggested that a number of factors other than PIC can contribute to ΔF values, including mechanisms underlying spike frequency adaptation and spike threshold accommodation. In the present study, we used a set of compartmental models representing a sample of 20 motoneurons with a range of thresholds to investigate how several different intrinsic motoneuron properties can potentially contribute to variations in ΔF values. We drove the models with linearly increasing and decreasing noisy conductance commands of different rate of rise and duration and determined the influence of different intrinsic mechanisms on discharge hysteresis (the difference in excitatory drive at recruitment and de-recruitment) and ΔF. Our results indicate that, although other factors can contribute, variations in discharge hysteresis and ΔF values primarily reflect the contribution of dendritic PICs to motoneuron activation.

Comparison of dendritic calcium transients in juvenile wild type and SOD1(G93A) mouse lumbar motoneurons

Previous studies of spinal motoneurons in the SOD1 mouse model of amyotrophic lateral sclerosis have shown alterations long before disease onset, including increased dendritic branching, increased persistent Na⁺ and Ca²⁺ currents, and impaired axonal transport. In this study dendritic Ca²⁺ entry was investigated using two photon excitation fluorescence microscopy and whole-cell patch-clamp of juvenile (P4-11) motoneurons. Neurons were filled with both Ca²⁺ Green-1 and Texas Red dextrans, and line scans performed throughout. Steps were taken to account for different sources of variability, including (1) dye filling and laser penetration, (2) dendritic anatomy, and (3) the time elapsed from the start of recording. First, Ca²⁺ Green-1 fluorescence was normalized by Texas Red; next, neurons were reconstructed so anatomy could be evaluated; finally, time was recorded. Customized software detected the largest Ca²⁺ transients (area under the curve) from each line scan and matched it with parameters above. Overall, larger dendritic diameter and shorter path distance from the soma were significant predictors of larger transients, while time was not significant up to 2 h (data thereafter was dropped). However, Ca²⁺ transients showed additional variability. Controlling for previous factors, significant variation was found between Ca²⁺ signals from different processes of the same neuron in 3/7 neurons. This could reflect differential expression of Ca²⁺ channels, local neuromodulation or other variations. Finally, Ca²⁺ transients in SOD1^G93A motoneurons were significantly smaller than in non-transgenic motoneurons. In conclusion, motoneuron processes show highly variable Ca²⁺ transients, but these transients are smaller overall in SOD1^G93A motoneurons.

Development of modified cable models to simulate accurate neuronal active behaviors

In large network and single three-dimensional (3-D) neuron simulations, high computing speed dictates using reduced cable models to simulate neuronal firing behaviors. However, these models are unwarranted under active conditions and lack accurate representation of dendritic active conductances that greatly shape neuronal firing. Here, realistic 3-D (R3D) models (which contain full anatomical details of dendrites) of spinal motoneurons were systematically compared with their reduced single unbranched cable (SUC, which reduces the dendrites to a single electrically equivalent cable) counterpart under passive and active conditions. The SUC models matched the R3D model's passive properties but failed to match key active properties, especially active behaviors originating from dendrites. For instance, persistent inward currents (PIC) hysteresis, frequency-current (FI) relationship secondary range slope, firing hysteresis, plateau potential partial deactivation, staircase currents, synaptic current transfer ratio, and regional FI relationships were not accurately reproduced by the SUC models. The dendritic morphology oversimplification and lack of dendritic active conductances spatial segregation in the SUC models caused significant underestimation of those behaviors. Next, SUC models were modified by adding key branching features in an attempt to restore their active behaviors. The addition of primary dendritic branching only partially restored some active behaviors, whereas the addition of secondary dendritic branching restored most behaviors. Importantly, the proposed modified models successfully replicated the active properties without sacrificing model simplicity, making them attractive candidates for running R3D single neuron and network simulations with accurate firing behaviors. The present results indicate that using reduced models to examine PIC behaviors in spinal motoneurons is unwarranted.

Neuromodulation impact on nonlinear firing behavior of a reduced model motoneuron with the active dendrite

Neuromodulatory inputs from brainstem systems modulate the normal function of spinal motoneurons by altering the activation properties of persistent inward currents (PICs) in their dendrites. However, the effect of the PIC on firing outputs also depends on its location in the dendritic tree. To investigate the interaction between PIC neuromodulation and PIC location dependence, we used a two-compartment model that was biologically realistic in that it retains directional and frequency-dependent electrical coupling between the soma and the dendrites, as seen in multi-compartment models based on full anatomical reconstructions of motoneurons. Our two-compartment approach allowed us to systematically vary the coupling parameters between the soma and the dendrite to accurately reproduce the effect of location of the dendritic PIC on the generation of nonlinear (hysteretic) motoneuron firing patterns. Our results show that as a single parameter value for PIC activation was either increased or decreased by 20% from its default value, the solution space of the coupling parameter values for nonlinear firing outputs was drastically reduced by approximately 80%. As a result, the model tended to fire only in a linear mode at the majority of dendritic PIC sites. The same results were obtained when all parameters for the PIC activation simultaneously changed only by approximately ±10%. Our results suggest the democratization effect of neuromodulation: the neuromodulation by the brainstem systems may play a role in switching the motoneurons with PICs at different dendritic locations to a similar mode of firing by reducing the effect of the dendritic location of PICs on the firing behavior.

Asymmetry in signal propagation between the soma and dendrites plays a key role in determining dendritic excitability in motoneurons

It is widely recognized that propagation of electrophysiological signals between the soma and dendrites of neurons differs depending on direction, i.e. it is asymmetric. How this asymmetry influences the activation of voltage-gated dendritic channels, and consequent neuronal behavior, remains unclear. Based on the analysis of asymmetry in several types of motoneurons, we extended our previous methodology for reducing a fully reconstructed motoneuron model to a two-compartment representation that preserved asymmetric signal propagation. The reduced models accurately replicated the dendritic excitability and the dynamics of the anatomical model involving a persistent inward current (PIC) dispersed over the dendrites. The relationship between asymmetric signal propagation and dendritic excitability was investigated using the reduced models while varying the asymmetry in signal propagation between the soma and the dendrite with PIC density constant. We found that increases in signal attenuation from soma to dendrites increased the activation threshold of a PIC (hypo-excitability), whereas increases in signal attenuation from dendrites to soma decreased the activation threshold of a PIC (hyper-excitability). These effects were so strong that reversing the asymmetry in the soma-to-dendrite vs. dendrite-to-soma attenuation, reversed the correlation between PIC threshold and distance of this current source from the soma. We propose the tight relation of the asymmetric signal propagation to the input resistance in the dendrites as a mechanism underlying the influence of the asymmetric signal propagation on the dendritic excitability. All these results emphasize the importance of maintaining the physiological asymmetry in dendritic signaling not only for normal function of the cells but also for biophysically realistic simulations of dendritic excitability.

Regulation of persistent activity in hippocampal mossy cells by inhibitory synaptic potentials

The hippocampal formation receives strong cholinergic input from the septal/diagonal band complex. Although the functional effects of cholinergic activation have been extensively studied in pyramidal neurons within the hippocampus and entorhinal cortex, less is known about the role of cholinergic receptors on dentate gyrus neurons. Using intracellular recordings from rat dentate hilar neurons, we find that activation of m1-type muscarinic receptors selectively increases the excitability of glutamatergic mossy cells but not of hilar interneurons. Following brief stimuli, cholinergic modulation reveals a latent afterdepolarization response in mossy cells that can extend the duration of stimulus-evoked depolarization by >100 msec. Depolarizing stimuli also could trigger persistent firing in mossy cells exposed to carbachol or an m1 receptor agonist. Evoked IPSPs attenuated the ADP response in mossy cells. The functional effect of IPSPs was amplified during ADP responses triggered in the presence of cholinergic receptor agonists but not during slowly decaying simulated ADPs, suggesting that modulation of ADP responses by IPSPs arises from destabilization of the intrinsic currents underlying the ADP. Evoked IPSPs also could halt persistent firing triggered by depolarizing stimuli. These results show that through intrinsic properties modulated by muscarinic receptors, mossy cells can prolong depolarizing responses to excitatory input and extend the time window where multiple synaptic inputs can summate. By actively regulating the intrinsic response to synaptic input, inhibitory synaptic input can dynamically control the integration window that enables detection of coincident inputs and shape the spatial pattern of hilar cell activity.

Early-life cocaine interferes with BDNF-mediated behavioral plasticity

An important aspect of goal-directed action selection is differentiating between actions that are more or less likely to be reinforced. With repeated performance or psychostimulant exposure, however, actions can assume stimulus-elicited-or "habitual"-qualities that are resistant to change. We show that selective knockdown of prelimbic prefrontal cortical Brain-derived neurotrophic factor (Bdnf) increases sensitivity to response-outcome associations, blocking habit-like behavioral inflexibility. A history of adolescent cocaine exposure, however, occludes the "beneficial" effects of Bdnf knockdown. This finding highlights a challenge in treating addiction-that drugs of abuse may bias decision-making toward habit systems even in individuals with putative neurobiological resiliencies.

Sequential acquisition of cacophony calcium currents, sodium channels and voltage-dependent potassium currents affects spike shape and dendrite growth during postembryonic maturation of an identified Drosophila motoneuron

During metamorphosis the CNS undergoes profound changes to accommodate the switch from larval to adult behaviors. In Drosophila and other holometabolous insects, adult neurons differentiate either from respecified larval neurons, newly born neurons, or are born embryonically but remain developmentally arrested until differentiation during pupal life. This study addresses the latter in the identified Drosophila flight motoneuron 5. In situ patch-clamp recordings, intracellular dye fills and immunocytochemistry address the interplay between dendritic shape, excitability and ionic current development. During pupal life, changes in excitability and spike shape correspond to a stereotyped, progressive appearance of voltage-gated ion channels. High-voltage-activated calcium current is the first current to appear at pupal stage P4, prior to the onset of dendrite growth. This is followed by voltage-gated sodium as well as transient potassium channel expression, when first dendrites grow, and sodium-dependent action potentials can be evoked by somatic current injection. Sustained potassium current appears later than transient potassium current. During the early stages of rapid dendritic growth, sodium-dependent action potentials are broadened by a calcium component. Narrowing of spike shape coincides with sequential increases in transient and sustained potassium currents during stages when dendritic growth ceases. Targeted RNAi knockdown of pupal calcium current significantly reduces dendritic growth. These data indicate that the stereotyped sequential acquisition of different voltage-gated ion channels affects spike shape and excitability such that activity-dependent calcium influx serves as a partner of genetic programs during critical stages of motoneuron dendrite growth.

Mechanisms underlying subunit independence in pyramidal neuron dendrites

Pyramidal neuron (PN) dendrites compartmentalize voltage signals and can generate local spikes, which has led to the proposal that their dendrites act as independent computational subunits within a multilayered processing scheme. However, when a PN is strongly activated, back-propagating action potentials (bAPs) sweeping outward from the soma synchronize dendritic membrane potentials many times per second. How PN dendrites maintain the independence of their voltage-dependent computations, despite these repeated voltage resets, remains unknown. Using a detailed compartmental model of a layer 5 PN, and an improved method for quantifying subunit independence that incorporates a more accurate model of dendritic integration, we first established that the output of each dendrite can be almost perfectly predicted by the intensity and spatial configuration of its own synaptic inputs, and is nearly invariant to the rate of bAP-mediated "cross-talk" from other dendrites over a 100-fold range. Then, through an analysis of conductance, voltage, and current waveforms within the model cell, we identify three biophysical mechanisms that together help make independent dendritic computation possible in a firing neuron, suggesting that a major subtype of neocortical neuron has been optimized for layered, compartmentalized processing under in-vivo-like spiking conditions.

Sodium-mediated plateau potentials in lumbar motoneurons of neonatal rats

The development and the ionic nature of bistable behavior in lumbar motoneurons were investigated in rats. One week after birth, almost all (∼80%) ankle extensor motoneurons recorded in whole-cell configuration displayed self-sustained spiking in response to a brief depolarization that emerged when the temperature was raised >30°C. The effect of L-type Ca(2+) channel blockers on self-sustained spiking was variable, whereas blockade of the persistent sodium current (I(NaP)) abolished them. When hyperpolarized, bistable motoneurons displayed a characteristic slow afterdepolarization (sADP). The sADPs generated by repeated depolarizing pulses summed to promote a plateau potential. The sADP was tightly associated with the emergence of Ca(2+) spikes. Substitution of extracellular Na(+) or chelation of intracellular Ca(2+) abolished both sADP and the plateau potential without affecting Ca(2+) spikes. These data suggest a key role of a Ca(2+)-activated nonselective cation conductance ((CaN)) in generating the plateau potential. In line with this, the blockade of (CaN) by flufenamate abolished both sADP and plateau potentials. Furthermore, 2-aminoethoxydiphenyl borate (2-APB), a common activator of thermo-sensitive vanilloid transient receptor potential (TRPV) cation channels, promoted the sADP. Among TRPV channels, only the selective activation of TRPV2 channels by probenecid promoted the sADP to generate a plateau potential. To conclude, bistable behaviors are, to a large extent, determined by the interplay between three currents: L-type I(Ca), I(NaP), and a Na(+)-mediated I(CaN) flowing through putative TRPV2 channels.

Motor unit

Movement is accomplished by the controlled activation of motor unit populations. Our understanding of motor unit physiology has been derived from experimental work on the properties of single motor units and from computational studies that have integrated the experimental observations into the function of motor unit populations. The article provides brief descriptions of motor unit anatomy and muscle unit properties, with more substantial reviews of motoneuron properties, motor unit recruitment and rate modulation when humans perform voluntary contractions, and the function of an entire motor unit pool. The article emphasizes the advances in knowledge on the cellular and molecular mechanisms underlying the neuromodulation of motoneuron activity and attempts to explain the discharge characteristics of human motor units in terms of these principles. A major finding from this work has been the critical role of descending pathways from the brainstem in modulating the properties and activity of spinal motoneurons. Progress has been substantial, but significant gaps in knowledge remain.

Statistical computer model analysis of the reciprocal and recurrent inhibitions of the Ia-EPSP in α-motoneurons

We simulate the inhibition of Ia-glutamatergic excitatory postsynaptic potential (EPSP) by preceding it with glycinergic recurrent (REN) and reciprocal (REC) inhibitory postsynaptic potentials (IPSPs). The inhibition is evaluated in the presence of voltage-dependent conductances of sodium, delayed rectifier potassium, and slow potassium in five α-motoneurons (MNs). We distribute the channels along the neuronal dendrites using, alternatively, a density function of exponential rise (ER), exponential decay (ED), or a step function (ST). We examine the change in EPSP amplitude, the rate of rise (RR), and the time integral (TI) due to inhibition. The results yield six major conclusions. First, the EPSP peak and the kinetics depending on the time interval are either amplified or depressed by the REC and REN shunting inhibitions. Second, the mean EPSP peak, its TI, and RR inhibition of ST, ER, and ED distributions turn out to be similar for analogous ranges of G. Third, for identical G, the large variations in the parameters' values can be attributed to the sodium conductance step (g(Na_step)) and the active dendritic area. We find that small g(Na_step) on a few dendrites maintains the EPSP peak, its TI, and RR inhibition similar to the passive state, but high g(Na_step) on many dendrites decrease the inhibition and sometimes generates even an excitatory effect. Fourth, the MN's input resistance does not alter the efficacy of EPSP inhibition. Fifth, the REC and REN inhibitions slightly change the EPSP peak and its RR. However, EPSP TI is depressed by the REN inhibition more than the REC inhibition. Finally, only an inhibitory effect shows up during the EPSP TI inhibition, while there are both inhibitory and excitatory impacts on the EPSP peak and its RR.

Ca(v)2 channels mediate low and high voltage-activated calcium currents in Drosophila motoneurons

Different blends of membrane currents underlie distinct functions of neurons in the brain. A major step towards understanding neuronal function, therefore, is to identify the genes that encode different ionic currents. This study combined in situ patch clamp recordings of somatodendritic calcium currents in an identified adult Drosophila motoneuron with targeted genetic manipulation. Voltage clamp recordings revealed transient low voltage-activated (LVA) currents with activation between –60 mV and –70 mV as well as high voltage-activated (HVA) current with an activation voltage around –30 mV. LVA could be fully inactivated by prepulses to –50 mV and was partially amiloride sensitive. Recordings from newly generated mutant flies demonstrated that DmαG (Ca(v)3 homolog) encoded the amiloride-sensitive portion of the transient LVA calcium current. We further demonstrated that the Ca(v)2 homolog, Dmca1A, mediated the amiloride-insensitive component of LVA current. This novel role of Ca(v)2 channels was substantiated by patch clamp recordings from conditional mutants, RNAi knock-downs, and following Dmca1A overexpression. In addition, we show that Dmca1A underlies the HVA somatodendritic calcium currents in vivo. Therefore, the Drosophila Ca(v)2 homolog, Dmca1A, underlies HVA and LVA somatodendritic calcium currents in the same neuron. Interestingly, DmαG is required for regulating LVA and HVA derived from Dmca1A in vivo. In summary, each vertebrate gene family for voltage-gated calcium channels is represented by a single gene in Drosophila, namely Dmca1D (Ca(v)1), Dmca1A (Ca(v)2) and DmαG (Ca(v)3), but the commonly held view that LVA calcium currents are usually mediated by Ca(v)3 rather than Ca(v)2 channels may require reconsideration.

Contribution of intrinsic properties and synaptic inputs to motoneuron discharge patterns: a simulation study

Motoneuron discharge patterns reflect the interaction of synaptic inputs with intrinsic conductances. Recent work has focused on the contribution of conductances mediating persistent inward currents (PICs), which amplify and prolong the effects of synaptic inputs on motoneuron discharge. Certain features of human motor unit discharge are thought to reflect a relatively stereotyped activation of PICs by excitatory synaptic inputs; these features include rate saturation and de-recruitment at a lower level of net excitation than that required for recruitment. However, PIC activation is also influenced by the pattern and spatial distribution of inhibitory inputs that are activated concurrently with excitatory inputs. To estimate the potential contributions of PIC activation and synaptic input patterns to motor unit discharge patterns, we examined the responses of a set of cable motoneuron models to different patterns of excitatory and inhibitory inputs. The models were first tuned to approximate the current- and voltage-clamp responses of low- and medium-threshold spinal motoneurons studied in decerebrate cats and then driven with different patterns of excitatory and inhibitory inputs. The responses of the models to excitatory inputs reproduced a number of features of human motor unit discharge. However, the pattern of rate modulation was strongly influenced by the temporal and spatial pattern of concurrent inhibitory inputs. Thus, even though PIC activation is likely to exert a strong influence on firing rate modulation, PIC activation in combination with different patterns of excitatory and inhibitory synaptic inputs can produce a wide variety of motor unit discharge patterns.

The beginning of intracellular recording in spinal neurons: facts, reflections, and speculations

Intracellular (IC) recording of action potentials in neurons of the vertebrate central nervous system (CNS) was first reported by John Eccles and two colleagues, Walter Brock and John Coombs, in Dunedin, NZL in 1951/1952 and by Walter Woodbury and Harry Patton in Seattle, WA, USA in 1952. Both groups studied spinal cord neurons of the adult cat. In this review, we discuss the precedents to their notable achievement and reflect and speculate on some of the scientific and personal nuances of their work and its immediate and later impact. We then briefly discuss early achievements in IC recording in the study of CNS neurobiology in other laboratories around the world, and some of the methods that led to enhancement of CNS IC-recording techniques. Our modern understanding of CNS neurophysiology directly emanates from the pioneering endeavors of the five who wrote the seminal 1951/1952 articles.

Whither motoneurons?

In the preceding series of articles, the history of vertebrate motoneuron and motor unit neurobiological studies has been discussed. In this article, we select a few examples of recent advances in neuroscience and discuss their application or potential application to the study of motoneurons and the control of movement. We conclude, like Sherrington, that in order to understand normal, traumatized, and diseased human behavior, it is critical to continue to study motoneuron biology using all available and emerging tools. This article is part of a Special Issue entitled Historical Review.

Tetrodotoxin-, dihydropyridine-, and riluzole-resistant persistent inward current: novel sodium channels in rodent spinal neurons

Recently, we reported the tetrodotoxin (TTX)- and dihydropyridine (DHP)-resistant (TDR) inward currents in neonatal mouse spinal neurons. In this study, we further characterized these currents in the presence of 1-5 μM TTX and 20-30 μM DHP (nifedipine, nimodipine, or isradipine). TDR inward currents were recorded by voltage ramp (persistent inward current, TDR-PIC) and step (TDR-I(p)) protocols. TDR-PIC and TDR-I(p) were found in 80.2% of recorded neurons (101/126) crossing laminae I to X from T12 to L6. TDR-PIC activated at -28.6 ± 13 mV with an amplitude of 80.6 ± 75 pA and time constant of 470.6 ± 240 ms (n = 75). TDR-I(p) had an amplitude of 151.2 ± 151 pA and a voltage threshold of -17.0 ± 9 mV (n = 54) with a wide range of kinetics parameters. The half-maximal activation was -21.5 ± 8 mV (-37 to -12 mV, n = 29) with a time constant of 5.2 ± 2 ms (1.2-11.2 ms, n = 19), whereas the half-maximal inactivation was -26.9 ± 9 mV (-39 to -18 mV, n = 14) with a time constant of 1.4 ± 0.4 s (0.5-2.2 s, n = 19). TDR-PIC and TDR-I(p) could be reduced by 60% in zero calcium and completely removed in zero sodium solutions, suggesting that they were mediated by sodium ions. Furthermore, the reversal potential of TDR-I(p) was estimated as 56.6 ± 3 mV (n = 10). TDR-PIC and TDR-I(p) persisted in 1-205 μM TTX, 20-100 μM DHP, 3-30 μM riluzole, 50-300 μM flufenamic acid, and 2-30 mM intracellular BAPTA. They also persisted with T-, N-, P/Q-, and R-type calcium channel blockers. In conclusion, we demonstrated novel TTX-, DHP-, and riluzole-resistant sodium channels in neonatal rodent spinal neurons. The unique pharmacological and electrophysiological properties would allow these channels to play a functional role in spinal motor system.

Asymmetric electrotonic coupling between the soma and dendrites alters the bistable firing behaviour of reduced models

The goal of the study was to investigate the influence of asymmetric coupling, between the soma and dendrites, on the nonlinear dynamic behaviour of a two-compartment model. We used a recently published method for generating reduced two-compartment models that retain the asymmetric coupling of anatomically reconstructed motor neurons. The passive input-output relationship of the asymmetrically coupled model was analytically compared to the symmetrically coupled case. Predictions based on the analytic comparison were tested using numerical simulations. The simulations evaluated the nonlinear dynamics of the models as a function of coupling parameters. Analytical results showed that the input resistance at the dendrite of the asymmetric model was directly related to the degree of coupling asymmetry. In contrast, a comparable symmetric model had identical input resistances at both the soma and dendrite regardless of coupling strength. These findings lead to predictions that variations in dendritic excitability, subsequent to changes in input resistance, might change the current threshold and onset timing of the plateau potential generated in the dendrite. Since the plateau potential underlies bistable firing, these results further predicted that asymmetric coupling might alter nonlinear (i.e. bistable) firing patterns. The numerical simulations supported analytical predictions, showing that the fully bistable firing pattern of the asymmetric model depended on the degree of coupling asymmetry and its correlated dendritic excitability. The physiological property of asymmetric coupling plays an important role in generating and stabilizing the bistability of motor neurons by interacting with the excitability of dendritic branches.

Evidence from computer simulations for alterations in the membrane biophysical properties and dendritic processing of synaptic inputs in mutant superoxide dismutase-1 motoneurons

A critical step in improving our understanding of the development of amyotrophic lateral sclerosis (ALS) is to identify the factors contributing to the alterations in the excitability of motoneurons and assess their individual contributions. Here we investigated the early alterations in the passive electrical and morphological properties of neonatal spinal motoneurons that occur by 10 d after birth, long before disease onset. We identified some of the factors contributing to these alterations, and estimated their individual contributions. To achieve this goal, we undertook a computer simulation analysis using realistic morphologies of reconstructed wild-type (WT) and mutant superoxide dismutase-1 (mSOD1) motoneurons. Ion channel parameters of these models were then tuned to match the experimental data on electrical properties obtained from these same motoneurons. We found that the reduced excitability of mSOD1 models was accompanied with decreased specific membrane resistance by approximately 25% and efficacy of synaptic inputs (slow and fast) by 12-22%. Linearity of summation of synaptic currents was similar to WT. We also assessed the contribution of the alteration in dendritic morphology alone to this decreased excitability and found that it reduced the input resistance by 10% and the efficacy of synaptic inputs by 7-15%. Our results were also confirmed in models with dendritic active conductances. Our simulations indicated that the alteration in passive electrical properties of mSOD1 models resulted from concurrent alterations in their morphology and membrane biophysical properties, and consequently altered the motoneuronal dendritic processing of synaptic inputs. These results clarify new aspects of spinal motoneurons malfunction in ALS.

Distribution of vestibulospinal contacts on the dendrites of ipsilateral splenius motoneurons: an anatomical substrate for push-pull interactions during vestibulocollic reflexes

Excitatory and inhibitory synapses may control neuronal output through a push-pull mechanism--that is, increases in excitation are coupled to simultaneous decreases in inhibition or vice versa. This pattern of activity is characteristic of excitatory and inhibitory vestibulospinal axons that mediate vestibulocollic reflexes. Previously, we showed that medial vestibulospinal tract (MVST) neurons in the rostral descending vestibular nucleus (DVN), an excitatory pathway, primarily innervate the medial dendrites of contralateral splenius motoneurons. In the present study, we tested the hypothesis that the counterparts of the push-pull mechanism, the ipsilateral inhibitory MVST synapses, are distributed on the dendritic tree such that the interactions with excitatory MVST synapses are enhanced. We combined anterograde tracing and intracellular staining in adult felines and show that most contacts (approximately 70%) between inhibitory MVST neurons in the rostral DVN and ipsilateral splenius motoneurons are also located on medial dendrites. There was a weak bias towards proximal dendrites. Using computational methods, we further show that the organization of excitatory and inhibitory MVST synapses on splenius motoneurons increases their likelihood for interaction. We found that if either excitatory or inhibitory MVST synapses were uniformly distributed throughout the dendritic tree, the proportion of inhibitory contacts in close proximity to excitatory contacts decreased. Thus, the compartmentalized distribution of excitatory and inhibitory MVST synapses on splenius motoneurons may be specifically designed to enhance their interactions during vestibulocollic reflexes. This suggests that the push-pull modulation of motoneuron output is based, in part, on the spatial arrangement of synapses on the dendritic tree.