Motoneuron model of self-sustained firing after spinal cord injury

Effect of Non-invasive Spinal Stimulation on Self-sustained Firing Motoneuron Model: In-Silico Study Using Human Body Model

Transcutaneous spinal electrical stimulation (tSCS) is a non-invasive neuromodulation approach using a low intensity direct current. Recent developments in the technique have opened the possibility that tSCS can help restore motor function after spinal cord injury (SCI). However, the exact mechanism of action tSCS has on the spinal circuits is still unknown. Due to the complexity of experimental synthesis in a human model to delineate the mechanisms, models that link the stimulation paradigm and circuit behaviors are advantageous. Thus, this study aims to simulate the underlying changes in motor circuit firing rates in response to external stimuli induced by tSCS. Serial stimulations combining a high-fidelity finite element model with the human torso and spinal cord with a lumped motor neuron model is constructed. The parameters for both components of the model were derived from previous studies. We focused our analysis on a lumped motor neuron model that describes sustained firing behavior of the motor neuron driven primarily by persistent inward current (PIC), a signature behavior of the motor neuron after SCI. Modulation of the PIC behaviors was achieved by stimulating voltage-dependent calcium and sodium channels in the dendrite using a tSCS-induced electric field (E-field) expressed at different a spatial locations of the motor neuron in the gray matter. The PIC behaviors of spinal motor neurons in the left ventral horn were suppressed, while for the most part invariant in the right ventral horn. These initial simulations will provide a steppingstone for future examinations that incorporate additional neuronal models of inhibitory and excitatory interneurons to access the circuit-level effect of spinal stimulation.

Trans-Spinal Direct Current Stimulation Targets Ca2+ Channels to Induce Persistent Motor Unit Responses

Trans-spinal direct current stimulation (tsDCS) is a neuromodulatory approach to augment spinal cord activity to improve function after neurological disease and injury. Little is known about the mechanisms underlying tsDCS actions on the motor system. The purpose of this study is to determine the role for a persistent inward current (PIC)-like response in motoneurons in mediating tsDCS actions. We recorded single motor units from the extensor and flexor carpi radialis muscles in healthy sedated rats and measured unit activity changes produced by cervical enlargement cathodal and anodal tsDCS (c-tsDCS; a-tsDCS). Both c-tsDCS and a-tsDCS immediately increased spontaneous motor unit firing during stimulation. After c-tsDCS was stopped, spontaneous firing persisted for a substantial period (165 ± 5s), yet after a-tsDCS activity shortly returned to baseline (27 ± 7s). Administration of the L-type calcium channel blocker Nimodipine reduced spontaneous motor unit firing during c-tsDCS and blocked the persistent response. By contrast, Nimodipine did not change unit firing during a-tsDCS but the short persistent response was blocked. Computer simulation using a two-compartment neuronal model replicated the main experimental observations: larger and more persistent responses during and after c-tsDCS than a-tsDCS. Using reduced Ca²⁺ conductance to model Nimodipine action, a reduced response during c-tsDCS and elimination of the persistent response was observed. Our experimental findings, supported by computer simulation, show that c-tsDCS can target Ca²⁺ conductances to augment motoneuron activity. As tsDCS is well-tolerated in humans, this knowledge informs therapeutic treatment strategies to achieve rehabilitation goals after injury; in particular, to increase muscle force.

Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings

Thoracic paravertebral sympathetic postganglionic neurons (tSPNs) comprise the final integrative output of the distributed sympathetic nervous system controlling vascular and thermoregulatory systems. Considered a non-integrating relay, what little is known of tSPN intrinsic excitability has been determined by sharp microelectrodes with presumed impalement injury. We thus undertook the first electrophysiological characterization of tSPN cellular properties using whole-cell recordings and coupled results with a conductance-based model to explore the principles governing their excitability in adult mice of both sexes. Recorded membrane resistance and time constant values were an order of magnitude greater than values previously obtained, leading to a demonstrable capacity for synaptic integration in driving recruitment. Variation in membrane resistivity was the primary determinant controlling cell excitability with vastly lower currents required for tSPN recruitment. Unlike previous microelectrode recordings in mouse which observed inability to sustain firing, all tSPNs were capable of repetitive firing. Computational modeling demonstrated that observed differences are explained by introduction of a microelectrode impalement injury conductance. Overall, tSPNs largely linearly encoded injected current magnitudes over a broad frequency range with distinct subpopulations differentiable based on repetitive firing signatures. Thus, whole-cell recordings reveal tSPNs have more dramatically amplified excitability than previously thought, with greater intrinsic capacity for synaptic integration and with the ability for maintained firing to support sustained actions on vasomotor tone and thermoregulatory function. Rather than acting as a relay, these studies support a more responsive role and possible intrinsic capacity for tSPNs to drive sympathetic autonomic function.

PyMUS: Python-Based Simulation Software for Virtual Experiments on Motor Unit System

We constructed a physiologically plausible computationally efficient model of a motor unit and developed simulation software that allows for integrative investigations of the input-output processing in the motor unit system. The model motor unit was first built by coupling the motoneuron model and muscle unit model to a simplified axon model. To build the motoneuron model, we used a recently reported two-compartment modeling approach that accurately captures the key cell-type-related electrical properties under both passive conditions (somatic input resistance, membrane time constant, and signal attenuation properties between the soma and the dendrites) and active conditions (rheobase current and afterhyperpolarization duration at the soma and plateau behavior at the dendrites). To construct the muscle unit, we used a recently developed muscle modeling approach that reflects the experimentally identified dependencies of muscle activation dynamics on isometric, isokinetic and dynamic variation in muscle length over a full range of stimulation frequencies. Then, we designed the simulation software based on the object-oriented programing paradigm and developed the software using open-source Python language to be fully operational using graphical user interfaces. Using the developed software, separate simulations could be performed for a single motoneuron, muscle unit and motor unit under a wide range of experimental input protocols, and a hierarchical analysis could be performed from a single channel to the entire system behavior. Our model motor unit and simulation software may represent efficient tools not only for researchers studying the neural control of force production from a cellular perspective but also for instructors and students in motor physiology classroom settings.

Intensity sensitive modulation effect of theta burst form of median nerve stimulation on the monosynaptic spinal reflex

The effects of electrical stimulation of median nerve with a continuous theta burst pattern (EcTBS) on the spinal H-reflex were studied. Different intensities and durations of EcTBS were given to the median nerve to 11 healthy individuals. The amplitude ratio of the H-reflex to maximum M wave (H/M ratio), corticospinal excitability and inhibition measured using motor evoked potentials (MEPs), short-interval intracortical inhibition and facilitation (SICI/ICF), spinal reciprocal inhibition (RI), and postactivation depression (PAD) were measured before and after EcTBS. In result, the H/M ratio was reduced followed by EcTBS at 90% H-reflex threshold, and the effect lasted longer after 1200 pulses than after 600 pulses of EcTBS. In contrast, EcTBS at 110% threshold facilitated the H/M ratio, while at 80% threshold it had no effect. Maximum M wave, MEPs, SICI/ICF, RI, and PAD all remained unchanged after EcTBS. In conclusion, EcTBS produced lasting effects purely on the H-reflex, probably, through effects on postsynaptic plasticity. The effect of EcTBS depends on the intensity and duration of stimulation. EcTBS is beneficial to research on mechanisms of human plasticity. Moreover, its ability to modulate spinal excitability is expected to have therapeutic benefits on neurological disorders involving spinal cord dysfunction.

Differential contributions of somatic and dendritic calcium-dependent potassium currents to the control of motoneuron excitability following spinal cord injury

The hyperexcitability of alpha-motoneurons and accompanying spasticity following spinal cord injury (SCI) have been attributed to enhanced persistent inward currents (PICs), including L-type calcium and persistent sodium currents. Factors controlling PICs may offer new therapies for managing spasticity. Such factors include calcium-activated potassium (KCa) currents, comprising in motoneurons an after-hyperpolarization-producing current (I KCaN) activated by N/P-type calcium currents, and a second current (I KCaL) activated by L-type calcium currents (Li and Bennett in J neurophysiol 97:767-783, 2007). We hypothesize that these two currents offer differential control of PICs and motoneuron excitability based on their probable somatic and dendritic locations, respectively. We reproduced SCI-induced PIC enhancement in a two-compartment motoneuron model that resulted in persistent dendritic plateau potentials. Removing dendritic I KCaL eliminated primary frequency range discharge and produced an abrupt transition into tertiary range firing without significant changes in the overall frequency gain. However, I KCaN removal mainly increased the gain. Steady-state analyses of dendritic membrane potential showed that I KCaL limits plateau potential magnitude and strongly modulates the somatic injected current thresholds for plateau onset and offset. In contrast, I KCaN had no effect on the plateau magnitude and thresholds. These results suggest that impaired function of I KCaL may be an important intrinsic mechanism underlying PIC-induced motoneuron hyperexcitability following SCI.

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.

Modulation of inhibitory strength and kinetics facilitates regulation of persistent inward currents and motoneuron excitability following spinal cord injury

Spasticity is commonly observed after chronic spinal cord injury (SCI) and many other central nervous system disorders (e.g., multiple sclerosis, stroke). SCI-induced spasticity has been associated with motoneuron hyperexcitability partly due to enhanced activation of intrinsic persistent inward currents (PICs). Disrupted spinal inhibitory mechanisms also have been implicated. Altered inhibition can result from complex changes in the strength, kinetics, and reversal potential (E(Cl(-))) of γ-aminobutyric acid A (GABA(A)) and glycine receptor currents. Development of optimal therapeutic strategies requires an understanding of the impact of these interacting factors on motoneuron excitability. We employed computational methods to study the effects of conductance, kinetics, and E(Cl(-)) of a dendritic inhibition on PIC activation and motoneuron discharge. A two-compartment motoneuron with enhanced PICs characteristic of SCI and receiving recurrent inhibition from Renshaw cells was utilized in these simulations. This dendritic inhibition regulated PIC onset and offset and exerted its strongest effects at motoneuron recruitment and in the secondary range of the current-frequency relationship during PIC activation. Increasing inhibitory conductance compensated for moderate depolarizing shifts in E(Cl(-)) by limiting PIC activation and self-sustained firing. Furthermore, GABA(A) currents exerted greater control on PIC activation than glycinergic currents, an effect attributable to their slower kinetics. These results suggest that modulation of the strength and kinetics of GABA(A) currents could provide treatment strategies for uncontrollable spasms.