Synaptic control of the shape of the motoneuron pool input-output function

NeuroMechanics: Electrophysiological and computational methods to accurately estimate the neural drive to muscles in humans in vivo

The ultimate neural signal for muscle control is the neural drive sent from the spinal cord to muscles. This neural signal comprises the ensemble of action potentials discharged by the active spinal motoneurons, which is transmitted to the innervated muscle fibres to generate forces. Accurately estimating the neural drive to muscles in humans in vivo is challenging since it requires the identification of the activity of a sample of motor units (MUs) that is representative of the active MU population. Current electrophysiological recordings usually fail in this task by identifying small MU samples with over-representation of higher-threshold with respect to lower-threshold MUs. Here, we describe recent advances in electrophysiological methods that allow the identification of more representative samples of greater numbers of MUs than previously possible. This is obtained with large and very dense arrays of electromyographic electrodes. Moreover, recently developed computational methods of data augmentation further extend experimental MU samples to infer the activity of the full MU pool. In conclusion, the combination of new electrode technologies and computational modelling allows for an accurate estimate of the neural drive to muscles and opens new perspectives in the study of the neural control of movement and in neural interfacing.

Intrinsic motoneuron properties in typical human development

Motoneuron properties and their firing patterns undergo significant changes throughout development and in response to neuromodulators such as serotonin. Here, we examined the age-related development of self-sustained firing and general excitability of tibialis anterior motoneurons in a young development (7-17 years), young adult (18-28 years) and adult (32-53 years) group, as well as in a separate group of participants taking selective serotonin reuptake inhibitors (SSRIs, aged 11-28 years). Self-sustained firing, as measured by ΔF, was larger in the young development (∼5.8 Hz, n = 20) compared to the young adult (∼4.9 Hz, n = 13) and adult (∼4.8 Hz, n = 8) groups, consistent with a developmental decrease in self-sustained firing mediated by persistent inward currents (PIC). ΔF was also larger in participants taking SSRIs (∼6.5 Hz, n = 9) compared to their age-matched controls (∼5.3 Hz, n = 26), consistent with increased levels of spinal serotonin facilitating the motoneuron PIC. Participants in the young development and SSRI groups also had higher firing rates and a steeper acceleration in initial firing rates (secondary ranges), consistent with the PIC producing a steeper acceleration in membrane depolarization at the onset of motoneuron firing. In summary, both the young development and SSRI groups exhibited increased intrinsic motoneuron excitability compared to the adults, which, in the young development group, was also associated with a larger unsteadiness in the dorsiflexion torque profiles. We propose several intrinsic and extrinsic factors that affect both motoneuron PICs and cell discharge which vary during development, with a time course similar to the changes in motoneuron firing behaviour observed in the present study. KEY POINTS: Neurons in the spinal cord that activate muscles in the limbs (motoneurons) undergo increases in excitability shortly after birth to help animals stand and walk. We examined whether the excitability of human ankle flexor motoneurons also continues to change from child to adulthood by recording the activity of the muscle fibres they innervate. Motoneurons in children and adolescents aged 7-17 years (young development group) had higher signatures of excitability that included faster firing rates and more self-sustained activity compared to adults aged ≥18 years. Participants aged 11-28 years of age taking serotonin reuptake inhibitors had the highest measures of motoneuron excitability compared to their age-matched controls. The young development group also had more unstable contractions, which might partly be related to the high excitability of the motoneurons.

Motor unit firing rate modulation is more impaired during flexion synergy-driven contractions of the biceps brachii in chronic stroke

Following a hemiparetic stroke, individuals exhibit altered motor unit firing patterns during voluntary muscle contractions, including impairments in firing rate modulation and recruitment. These individuals also exhibit abnormal muscle coactivation through multi-joint synergies (e.g., flexion synergy). Here, we investigate whether motor unit firing activity during flexion synergy-driven contractions of the paretic biceps brachii differs from that of voluntary contractions and use these differences to predict changes in descending motor commands. To accomplish this, we characterized motor unit firing patterns of the biceps brachii in individuals with chronic hemiparetic stroke during voluntary isometric elbow flexion contractions in the paretic and non-paretic limbs, as well as during contractions driven by voluntary effort and by flexion synergy expression in the paretic limb. We observed significant reductions in motor unit firing rate modulation from the non-paretic to paretic limb (non-paretic - paretic: 0.14 pps/%MVT, 95% CI: [0.09 0.19]) that were further reduced during synergy-driven contractions (voluntary paretic - synergy driven: 0.19 pps/%MVT, 95% CI: [0.14 0.25]). Moreover, using recently developed metrics, we evaluated how a stroke-induced reliance on indirect motor pathways alters the inputs that motor units receive and revealed progressive increases in neuromodulatory and inhibitory drive to the motor pool in the paretic limb, with the changes greatest during synergy-driven contractions. These findings suggest that an interplay between heightened neuromodulatory drive and alterations in inhibitory command structure may account for the observed motor unit impairments, further illuminating underlying neural mechanisms involved in the flexion synergy and its impact on motor unit firing patterns post-stroke.

Multivariable closed-loop control of deep brain stimulation for Parkinson's disease

Objective. Closed-loop deep brain stimulation (DBS) methods for Parkinson's disease (PD) to-date modulate either stimulation amplitude or frequency to control a single biomarker. While good performance has been demonstrated for symptoms that are correlated with the chosen biomarker, suboptimal regulation can occur for uncorrelated symptoms or when the relationship between biomarker and symptom varies. Control of stimulation-induced side-effects is typically not considered.Approach.A multivariable control architecture is presented to selectively target suppression of either tremor or subthalamic nucleus beta band oscillations. DBS pulse amplitude and duration are modulated to maintain amplitude below a threshold and avoid stimulation of distal large diameter axons associated with stimulation-induced side effects. A supervisor selects between a bank of controllers which modulate DBS pulse amplitude to control rest tremor or beta activity depending on the level of muscle electromyographic (EMG) activity detected. A secondary controller limits pulse amplitude and modulates pulse duration to target smaller diameter axons lying close to the electrode. The control architecture was investigated in a computational model of the PD motor network which simulated the cortico-basal ganglia network, motoneuron pool, EMG and muscle force signals.Main results.Good control of both rest tremor and beta activity was observed with reduced power delivered when compared with conventional open loop stimulation, The supervisor avoided over- or under-stimulation which occurred when using a single controller tuned to one biomarker. When DBS amplitude was constrained, the secondary controller maintained the efficacy of stimulation by increasing pulse duration to compensate for reduced amplitude. Dual parameter control delivered effective control of the target biomarkers, with additional savings in the power delivered.Significance.Non-linear multivariable control can enable targeted suppression of motor symptoms for PD patients. Moreover, dual parameter control facilitates automatic regulation of the stimulation therapeutic dosage to prevent overstimulation, whilst providing additional power savings.

Motor Unit Discharge Characteristics and Conduction Velocity of the Vastii Muscles in Long-Term Resistance-Trained Men

PURPOSE

Adjustments in motor unit (MU) discharge properties have been shown after short-term resistance training; however, MU adaptations in long-term resistance-trained (RT) individuals are less clear. Here, we concurrently assessed MU discharge characteristics and MU conduction velocity in long-term RT and untrained (UT) men.

METHODS

Motor unit discharge characteristics (discharge rate, recruitment, and derecruitment threshold) and MU conduction velocity were assessed after the decomposition of high-density electromyograms recorded from vastus lateralis (VL) and vastus medialis (VM) of RT (>3 yr; n = 14) and UT ( n = 13) during submaximal and maximal isometric knee extension.

RESULTS

Resistance-trained men were on average 42% stronger (maximal voluntary force [MVF], 976.7 ± 85.4 N vs 685.5 ± 123.1 N; P < 0.0001), but exhibited similar relative MU recruitment (VL, 21.3% ± 4.3% vs 21.0% ± 2.3% MVF; VM, 24.5% ± 4.2% vs 22.7% ± 5.3% MVF) and derecruitment thresholds (VL, 20.3% ± 4.3% vs 19.8% ± 2.9% MVF; VM, 24.2% ± 4.8% vs 22.9% ± 3.7% MVF; P ≥ 0.4543). There were also no differences between groups in MU discharge rate at recruitment and derecruitment or at the plateau phase of submaximal contractions (VL, 10.6 ± 1.2 pps vs 10.3 ± 1.5 pps; VM, 10.7 ± 1.6 pps vs 10.8 ± 1.7 pps; P ≥ 0.3028). During maximal contractions of a subsample population (10 RT, 9 UT), MU discharge rate was also similar in RT compared with UT (VL, 21.1 ± 4.1 pps vs 14.0 ± 4.5 pps; VM, 19.5 ± 5.0 pps vs 17.0 ± 6.3 pps; P = 0.7173). Motor unit conduction velocity was greater in RT compared with UT individuals in both VL (4.9 ± 0.5 m·s -1 vs 4.5 ± 0.3 m·s -1 ; P < 0.0013) and VM (4.8 ± 0.5 m·s -1 vs 4.4 ± 0.3 m·s -1 ; P < 0.0073).

CONCLUSIONS

Resistance-trained and UT men display similar MU discharge characteristics in the knee extensor muscles during maximal and submaximal contractions. The between-group strength difference is likely explained by superior muscle morphology of RT as suggested by greater MU conduction velocity.

Inhibitory Synaptic Influences on Developmental Motor Disorders

During development, GABA and glycine play major trophic and synaptic roles in the establishment of the neuromotor system. In this review, we summarise the formation, function and maturation of GABAergic and glycinergic synapses within neuromotor circuits during development. We take special care to discuss the differences in limb and respiratory neuromotor control. We then investigate the influences that GABAergic and glycinergic neurotransmission has on two major developmental neuromotor disorders: Rett syndrome and spastic cerebral palsy. We present these two syndromes in order to contrast the approaches to disease mechanism and therapy. While both conditions have motor dysfunctions at their core, one condition Rett syndrome, despite having myriad symptoms, has scientists focused on the breathing abnormalities and their alleviation-to great clinical advances. By contrast, cerebral palsy remains a scientific quagmire or poor definitions, no widely adopted model and a lack of therapeutic focus. We conclude that the sheer abundance of diversity of inhibitory neurotransmitter targets should provide hope for intractable conditions, particularly those that exhibit broad spectra of dysfunction-such as spastic cerebral palsy and Rett syndrome.

A geometric approach to quantifying the neuromodulatory effects of persistent inward currents on individual motor unit discharge patterns

Objective.All motor commands flow through motoneurons, which entrain control of their innervated muscle fibers, forming a motor unit (MU). Owing to the high fidelity of action potentials within MUs, their discharge profiles detail the organization of ionotropic excitatory/inhibitory as well as metabotropic neuromodulatory commands to motoneurons. Neuromodulatory inputs (e.g. norepinephrine, serotonin) enhance motoneuron excitability and facilitate persistent inward currents (PICs). PICs introduce quantifiable properties in MU discharge profiles by augmenting depolarizing currents upon activation (i.e. PIC amplification) and facilitating discharge at lower levels of excitatory input than required for recruitment (i.e. PIC prolongation).Approach. Here, we introduce a novel geometric approach to estimate neuromodulatory and inhibitory contributions to MU discharge by exploiting discharge non-linearities introduced by PIC amplification during time-varying linear tasks. In specific, we quantify the deviation from linear discharge ('brace height') and the rate of change in discharge (i.e. acceleration slope, attenuation slope, angle). We further characterize these metrics on a simulated motoneuron pool with known excitatory, inhibitory, and neuromodulatory inputs and on human MUs (number of MUs; Tibialis Anterior: 1448, Medial Gastrocnemius: 2100, Soleus: 1062, First Dorsal Interosseus: 2296).Main results. In the simulated motor pool, we found brace height and attenuation slope to consistently indicate changes in neuromodulation and the pattern of inhibition (excitation-inhibition coupling), respectively, whereas the paired MU analysis (ΔF) was dependent on both neuromodulation and inhibition pattern. Furthermore, we provide estimates of these metrics in human MUs and show comparable variability in ΔFand brace height measures for MUs matched across multiple trials.Significance. Spanning both datasets, we found brace height quantification to provide an intuitive method for achieving graded estimates of neuromodulatory and inhibitory drive to individual MUs. This complements common techniques and provides an avenue for decoupling changes in the level of neuromodulatory and pattern of inhibitory motor commands.

Estimation of the firing behaviour of a complete motoneuron pool by combining electromyography signal decomposition and realistic motoneuron modelling

Our understanding of the firing behaviour of motoneuron (MN) pools during human voluntary muscle contractions is currently limited to electrophysiological findings from animal experiments extrapolated to humans, mathematical models of MN pools not validated for human data, and experimental results obtained from decomposition of electromyographical (EMG) signals. These approaches are limited in accuracy or provide information on only small partitions of the MN population. Here, we propose a method based on the combination of high-density EMG (HDEMG) data and realistic modelling for predicting the behaviour of entire pools of motoneurons in humans. The method builds on a physiologically realistic model of a MN pool which predicts, from the experimental spike trains of a smaller number of individual MNs identified from decomposed HDEMG signals, the unknown recruitment and firing activity of the remaining unidentified MNs in the complete MN pool. The MN pool model is described as a cohort of single-compartment leaky fire-and-integrate (LIF) models of MNs scaled by a physiologically realistic distribution of MN electrophysiological properties and driven by a spinal synaptic input, both derived from decomposed HDEMG data. The MN spike trains and effective neural drive to muscle, predicted with this method, have been successfully validated experimentally. A representative application of the method in MN-driven neuromuscular modelling is also presented. The proposed approach provides a validated tool for neuroscientists, experimentalists, and modelers to infer the firing activity of MNs that cannot be observed experimentally, investigate the neuromechanics of human MN pools, support future experimental investigations, and advance neuromuscular modelling for investigating the neural strategies controlling human voluntary contractions.

Our experimental understanding of the firing behaviour of motoneuron (MN) pools during human voluntary muscle contractions is currently limited to the observation of small samples of active MNs obtained from EMG decomposition. EMG decomposition therefore provides an important but incomplete description of the role of individual MNs in the firing activity of the complete MN pool, which limits our understanding of the neural strategies of the whole MN pool and of how the firing activity of each MN contributes to the neural drive to muscle. Here, we combine decomposed high-density EMG (HDEMG) data and a physiologically realistic model of MN population to predict the unknown recruitment and firing activity of the remaining unidentified MNs in the complete MN pool. In brief, an experimental estimation of the synaptic current is input to a cohort of MN models, which are calibrated using the available decomposed HDEMG data, and predict the MN spike trains fired by the entire MN population. This novel approach is experimentally validated and applied to muscle force prediction from neuromuscular modelling.

Chronic exercise increases excitability of lamina X neurons through enhancement of persistent inward currents and dendritic development in mice

KEY POINTS

Chronic exercise alters adaptability of spinal motor system in rodents. Multiple mechanisms are responsible for the adaptation, including regulation of neuronal excitability and change in dendritic morphology. Spinal interneurons in lamina X are a cluster of heterogeneous neurons playing multifunctional roles in the spinal cord, especially in regulating locomotor activity. Chronic exercise in juvenile mice increased excitability of these interneurons and facilitated dendritic development. Mechanisms underlying these changes remain unknown. Lamina X neurons expressed persistent inward currents (PICs) composed of calcium (Ca-PIC) and sodium (Na-PIC) components. The exercise-increased excitability of lamina X neurons was mediated by enhancing Ca-PIC and Na-PIC components and facilitating dendritic length. Na-PIC contributed more to lowering of PIC onset and Ca-PIC to increase of PIC amplitude. This study unveiled novel morphological and ionic mechanisms underlying adaptation of lamina X neurons in rodents during chronic exercise.

ABSTRACT

Chronic exercise has been shown to enhance excitability of spinal interneurons in rodents. However, the mechanisms underlying this enhancement remain unclear. In this study we investigated adaptability of lamina X neurons with three-week treadmill exercise in mice of P21-P24. Whole-cell path-clamp recording was performed on the interneurons from slices of T12-L4. The experimental results included: (1) Treadmill exercise reduced rheobase by 7.4±2.2 pA (control: 11.3±6.1 pA, n = 12; exercise: 3.8±4.6 pA, n = 13; P = 0.002) and hyperpolarized voltage threshold by 7.1±1.5 mV (control: -36.6±4.6 mV, exercise: -43.7±2.7 mV; P = 0.001). (2) Exercise enhanced persistent inward currents (PICs) with increase of amplitude (control: 140.6±56.3 pA, n = 25; exercise: 225.9±62.5 pA, n = 17; P = 0.001) and hyperpolarization of onset (control: -50.3±3.6 mV, exercise: -56.5±5.5 mV; P = 0.001). (3) PICs consisted of dihydropyridine-sensitive calcium (Ca-PIC) and tetrodotoxin-sensitive sodium (Na-PIC) components. Exercise increased amplitude of both components but hyperpolarized onset of Na-PIC only. (4) Exercise reduced derecruitment current of repetitive firing evoked by current bi-ramp and prolonged firing in falling phase of the bi-ramp. The derecruitment reduction was eliminated by bath application of 3 μM riluzole or 25 μM nimodipine, suggesting that both Na-PIC and Ca-PIC contributed to the exercise-prolonged hysteresis of firing. (5) Exercise facilitated dendritic development with significant increase in dendritic length by 285.1±113 μm (control: 457.8±171.8 μm, n = 12; exercise: 742.9±357 μm, n = 14; P = 0.019). We concluded that three-week treadmill exercise increased excitability of lamina X interneurons through enhancement of PICs and increase of dendritic length. This study provided insight into cellular and channel mechanisms underlying adaptation of the spinal motor system in exercise. Abstract figure legend A. B6 mice were randomly divided into control group and exercise group. Control group mice remained sedentary in the cage; exercise group mice completed 60 min treadmill runs each day (6 days/week) for a period of 3 weeks. B. Whole-cell patch clamp recordings were made from lumbar lamina X neurons after three-weeks exercise. C. Exercise facilitated development of dendrites of lamina X neurons. D. Exercise enhanced persistent inward currents. E. Exercise increased excitability of lamina X neurons by hyperpolarizing voltage threshold for action potential generation. This article is protected by copyright. All rights reserved.

Locomotor Pattern and Force Generation Modulated by Ionic Channels: A Computational Study of Spinal Networks Underlying Locomotion

Locomotion is a fundamental movement in vertebrates produced by spinal networks known as central pattern generators (CPG). During fictive locomotion cat lumbar motoneurons (MNs) exhibit changes in membrane properties, including hyperpolarization of voltage threshold, reduction of afterhyperpolarization and input resistance, and amplification of nonlinear membrane properties. Both modeling and electrophysiological studies suggest that these changes can be produced by upregulating voltage-gated sodium channel (VGSC), persistent sodium (NaP), or L-type calcium channel (LTCC) or downregulating delayed-rectifier potassium (K(DR)) or calcium-dependent potassium channel (KCa) in spinal MNs. Further studies implicate that these channel modulations increase motor output and facilitate MN recruitment. However, it remains unknown how the channel modulation of CPG networks or MN pools affects the rhythmic generation of locomotion and force production of skeletal muscle during locomotion. In order to investigate this issue, we built a two-level CPG model composed of excitatory interneuron pools (Exc-INs), coupled reciprocally with inhibitory interneuron pools (Inh-INs), and projected to the flexor-extensor MN pools innervating skeletal muscles. Each pool consisted of 100 neurons with membrane properties based on cat spinal neurons. VGSC, K(DR), NaP, KCa, LTCC, and H-current channels were included in the model. Simulation results showed that (1) upregulating VGSC, NaP, or LTCC or downregulating KCa in MNs increased discharge rate and recruitment of MNs, thus facilitating locomotor pattern formation, increased amplitude of electroneurogram (ENG) bursting, and enhanced force generation of skeletal muscles. (2) The same channel modulation in Exc-INs increased the firing frequency of the Exc-INs, facilitated rhythmic generation, and increased flexor-extensor durations of step cycles. (3) Contrarily, downregulation of NaP or LTCC in MNs or Exc-INs or both CPG (Exc-INs and Inh-INs) and MNs disrupted locomotor pattern and reduced or even blocked the ENG bursting of MNs and force generation of skeletal muscles. (4) Pharmacological experiments showed that bath application of 25 μM nimodipine or 2 μM riluzole completely blocked fictive locomotion in isolated rat spinal cord, consistent with simulation results. We concluded that upregulation of VGSC, NaP, or LTCC or downregulation of KCa facilitated rhythmic generation and force production during walking, with NaP and LTCC playing an essential role.

The Cellular Basis for the Generation of Firing Patterns in Human Motor Units

Motor units, which comprise a motoneuron and the set of muscle fibers it innervates, are the fundamental neuromuscular transducers for all motor commands. The one to one relationship between a motoneuron and its innervated muscle fibers allow motoneuron firing patterns to be readily measured in humans. In this chapter, we summarize the current understanding of the cellular basis for the generation of firing patterns in human motor units. We provide a brief review of landmark insights from classic studies and then proceed to consider the features of motor unit firing patterns that are most likely to be sensitive estimators of motoneuron inputs and properties. In addition, we discuss recent advances in technology for recording human motor unit firing patterns and highly realistic computer simulations of motoneurons. The final section presents our recent efforts to use the power of supercomputers for implementation of the motoneuron models, with a goal of achieving a true "reverse engineering" approach that maximizes the insights from motor unit firing patterns into the synaptic structure of motor commands.

Firing rate trajectories of human occipitofrontalis motor units in response to triangular voluntary contraction intensity

During voluntary contractions, limb muscle motor unit (MU) firing rates accelerate over a small force range and saturate in response to increasing contraction intensity. In comparison, facial muscles are cranially innervated, and some function without crossing joints. Therefore, the MU firing rate behaviour and characteristics of saturation were explored in a facial muscle that moves skin and facia during voluntary contractions. We evaluated the firing rate trajectory in response to triangular voluntary contraction ramps in the occipitofrontalis muscle of 11 adult participants. Intramuscular electromyography of the frontalis aspect was used to record single MU trains followed up to maximal voluntary contraction intensities. Firing rates were measured from each MU sample, with the firing rate trajectory fit as both exponential (i.e., saturation) and linear models that were compared statistically. The rate coding behaviour of frontalis MUs was broad, as the peak firing rate (mean 76 Hz) was ninefold greater than the firing rate at recruitment threshold (mean 8 Hz). Across 20 MU samples, only 40% (8 MU samples) were determined to have a firing rate trajectory that saturated and had slow acceleration in response to increasing voluntary drive until maximum. The exponential curve of the firing rate trajectory had ~ tenfold lower acceleration as compared to prior reports in limb muscles. These results across all MU samples indicated that voluntary control of the frontalis muscle requires relatively slower accelerating or linear MU firing rate trajectories, suggesting that movements of facial muscles may be directly representative of extrinsic synaptic inputs.

Neural control of the healthy pectoralis major from low-to-moderate isometric contractions

The pectoralis major critically enables arm movement in several directions. However, its neural control remains unknown. High-density electromyography (HD-sEMG) was acquired from the pectoralis major in two sets of experiments in healthy young adults. Participants performed ramp-and-hold isometric contractions in: adduction, internal rotation, flexion, and horizontal adduction at three force levels: 15%, 25%, and 50% scaled to task-specific maximal voluntary force (MVF). HD-sEMG signals were decomposed into motor unit spike trains using a convolutive blind source separation algorithm and matched across force levels using a motor unit matching algorithm. The mean discharge rate and coefficient of variation were quantified across the hold and compared between 15% and 25% MVF across all tasks, whereas comparisons between 25% and 50% MVF were made where available. Mean motor unit discharge rate was not significantly different between 15% and 25% MVF (all P > 0.05) across all tasks or between 25% and 50% MVF in horizontal adduction (P = 0.11), indicating an apparent saturation across force levels and the absence of rate coding. These findings suggest that the pectoralis major likely relies on motor unit recruitment to increase force, providing first-line evidence of motor unit recruitment in this muscle and paving the way for more deliberate investigations of the pectoralis major involvement in shoulder function.NEW & NOTEWORTHY This work is the first to investigate the relative contribution of rate coding and motor unit recruitment in the pectoralis major muscle in several functionally relevant tasks and across varying force levels in healthy adults. Our results demonstrate the absence of motor unit rate coding with an increase in EMG amplitude with increases in force level in all tasks examined, indicating that the pectoralis major relies on motor unit recruitment to increase force.

In-silico development and assessment of a Kalman filter motor decoder for prosthetic hand control

Up to 50% of amputees abandon their prostheses, partly due to rapid degradation of the control systems, which require frequent recalibration. The goal of this study was to develop a Kalman filter-based approach to decoding motoneuron activity to identify movement kinematics and thereby provide stable, long-term, accurate, real-time decoding. The Kalman filter-based decoder was examined via biologically varied datasets generated from a high-fidelity computational model of the spinal motoneuron pool. The estimated movement kinematics controlled a simulated MuJoCo prosthetic hand. This clear-box approach showed successful estimation of hand movements under eight varied physiological conditions with no retraining. The mean correlation coefficient of 0.98 and mean normalized root mean square error of 0.06 over these eight datasets provide proof of concept that this decoder would improve long-term integrity of performance while performing new, untrained movements. Additionally, the decoder operated in real-time (~0.3 ms). Further results include robust performance of the Kalman filter when re-trained to more severe post-amputation limitations in the type and number of motoneurons remaining. An additional analysis shows that the decoder achieves better accuracy when using the firing of individual motoneurons as input, compared to using aggregate pool firing. Moreover, the decoder demonstrated robustness to noise affecting both the trained decoder parameters and the decoded motoneuron activity. These results demonstrate the utility of a proof of concept Kalman filter decoder that can support prosthetics' control systems to maintain accurate and stable real-time movement performance.

GABAergic Mechanisms Can Redress the Tilted Balance between Excitation and Inhibition in Damaged Spinal Networks

Correct operation of neuronal networks depends on the interplay between synaptic excitation and inhibition processes leading to a dynamic state termed balanced network. In the spinal cord, balanced network activity is fundamental for the expression of locomotor patterns necessary for rhythmic activation of limb extensor and flexor muscles. After spinal cord lesion, paralysis ensues often followed by spasticity. These conditions imply that, below the damaged site, the state of balanced networks has been disrupted and that restoration might be attempted by modulating the excitability of sublesional spinal neurons. Because of the widespread expression of inhibitory GABAergic neurons in the spinal cord, their role in the early and late phases of spinal cord injury deserves full attention. Thus, an early surge in extracellular GABA might be involved in the onset of spinal shock while a relative deficit of GABAergic mechanisms may be a contributor to spasticity. We discuss the role of GABA A receptors at synaptic and extrasynaptic level to modulate network excitability and to offer a pharmacological target for symptom control. In particular, it is proposed that activation of GABA A receptors with synthetic GABA agonists may downregulate motoneuron hyperexcitability (due to enhanced persistent ionic currents) and, therefore, diminish spasticity. This approach might constitute a complementary strategy to regulate network excitability after injury so that reconstruction of damaged spinal networks with new materials or cell transplants might proceed more successfully.

Serotonergic Modulation of Persistent Inward Currents in Serotonergic Neurons of Medulla in ePet-EYFP Mice

Serotonergic (5-HT) neurons in the medulla play multiple functional roles associated with many symptoms and motor activities. The descending serotonergic pathway from medulla is essential for initiating locomotion. However, the ionic properties of 5-HT neurons in the medulla remain unclear. Using whole-cell patch-clamp technique, we studied the biophysical and modulatory properties of persistent inward currents (PICs) in 5-HT neurons of medulla in ePet-EYFP transgenic mice (P3–P6). PICs were recorded by a family of voltage bi-ramps (10-s duration, 40-mV peak step), and the ascending and descending PICs were mirrored to analyze the PIC hysteresis. PICs were found in 77% of 5-HT neurons (198/258) with no significant difference between parapyramidal region (n = 107) and midline raphe nuclei (MRN) (n = 91) in either PIC onset (−47.4 ± 10 mV and −48.7 ± 7 mV; P = 0.44) or PIC amplitude (226.9 ± 138 pA and 259.2 ± 141 pA; P = 0.29). Ninety-six percentage (191/198) of the 5-HT neurons displayed counterclockwise hysteresis and four percentage (7/198) exhibited the clockwise hysteresis. The composite PICs could be differentiated as calcium component (Ca_PIC) by bath application of nimodipine (25 μM), sodium component (Na_PIC) by tetrodotoxin (TTX, 2 μM), and TTX- and dihydropyridine-resistance component (TDR_PIC) by TTX and nimodipine. Ca_PIC, Na_PIC and TDR_PIC all contributed to upregulation of excitability of 5-HT neurons. 5-HT (15 μM) enhanced the PICs, including a 26% increase in amplitude of the compound currents of Ca_PIC and TDR_PIC (P < 0.001, n = 9), 3.6 ± 5 mV hyperpolarization of Na_PIC and TDR_PIC onset (P < 0.05, n = 12), 30% increase in amplitude of TDR_PIC (P < 0.01), and 2.0 ± 3 mV hyperpolarization of TDR_PIC onset (P < 0.05, n = 18). 5-HT also facilitated repetitive firing of 5-HT neurons through modulation of composite PIC, Na_PIC and TDR_PIC, and Ca_PIC and TDR_PIC, respectively. In particular, the high voltage-activated TDR_PIC facilitated the repetitive firing in higher membrane potential, and this facilitation could be amplified by 5-HT. Morphological data analysis indicated that the dendrites of 5-HT neurons possessed dense spherical varicosities intensively crossing 5-HT neurons in medulla. We characterized the PICs in 5-HT neurons and unveiled the mechanism underlying upregulation of excitability of 5-HT neurons through serotonergic modulation of PICs. This study provided insight into channel mechanisms responsible for the serotonergic modulation of serotonergic neurons in brainstem.

Three-Week Treadmill Exercise Enhances Persistent Inward Currents, Facilitates Dendritic Plasticity, and Upregulates the Excitability of Dorsal Raphe Serotonin Neurons in ePet-EYFP Mice

Exercise plays a key role in preventing or treating mental or motor disorders caused by dysfunction of the serotonergic system. However, the electrophysiological and ionic channel mechanisms underlying these effects remain unclear. In this study, we investigated the effects of 3-week treadmill exercise on the electrophysiological and channel properties of dorsal raphe nucleus (DRN). Serotonin (5-HT) neurons in ePet-EYFP mice, using whole-cell patch clamp recording. Treadmill exercise was induced in ePet-EYFP mice of P21–24 for 3 weeks, and whole-cell patch clamp recording was performed on EYFP-positive 5-HT neurons from DRN slices of P42–45 mice. Experiment data showed that 5-HT neurons in the DRN were a heterogeneous population with multiple firing patterns (single firing, phasic firing, and tonic firing). Persistent inward currents (PICs) with multiple patterns were expressed in 5-HT neurons and composed of Cav1.3 (Ca-PIC) and sodium (Na-PIC) components. Exercise hyperpolarized the voltage threshold for action potential (AP) by 3.1 ± 1.0 mV (control: n = 14, exercise: n = 18, p = 0.005) and increased the AP amplitude by 6.7 ± 3.0 mV (p = 0.031) and firing frequency by more than 22% especially within a range of current stimulation stronger than 70 pA. A 3-week treadmill exercise was sufficient to hyperpolarize PIC onset by 2.6 ± 1.3 mV (control: −53.4 ± 4.7 mV, n = 28; exercise: −56.0 ± 4.7 mV, n = 25, p = 0.050) and increase the PIC amplitude by 28% (control: 193.6 ± 81.8 pA; exercise: 248.5 ± 105.4 pA, p = 0.038). Furthermore, exercise hyperpolarized Na-PIC onset by 3.8 ± 1.8 mV (control: n = 8, exercise: n = 9, p = 0.049) and increased the Ca-PIC amplitude by 23% (p = 0.013). The exercise-induced enhancement of the PIC amplitude was mainly mediated by Ca-PIC and hyperpolarization of PIC onset by Na-PIC. Moreover, exercise facilitated dendritic plasticity, which was shown as the increased number of branch points by 1.5 ± 0.5 (p = 0.009) and dendritic branches by 2.1 ± 0.6 (n = 20, p = 0.001) and length by 732.0 ± 100.1 μm (p < 0.001) especially within the range of 50–200 μm from the soma. Functional analysis suggested that treadmill exercise enhanced Na-PIC for facilitation of spike initiation and Ca-PIC for regulation of repetitive firing. We concluded that PICs broadly existed in DRN 5-HT neurons and could influence serotonergic neurotransmission in juvenile mice and that 3-week treadmill exercise induced synaptic adaptations, enhanced PICs, and thus upregulated the excitability of the 5-HT neurons.

A modeling study of spinal motoneuron recruitment regulated by ionic channels during fictive locomotion

During fictive locomotion cat lumbar motoneurons exhibit changes in membrane proprieties including a decrease in voltage threshold (V_th), afterhyperpolarization (AHP) and input resistance (R_in) and an increase in non-linear membrane property. The impact of these changes on the motoneuron recruitment remains unknown. Using modeling approach we investigated the channel mechanism regulating the motoneuron recruitment. Three types of motoneuron pools including slow (S), fatigue-resistant (FR) and fast-fatigable (FF) motoneurons were constructed based on the membrane proprieties of cat lumbar motoneurons. The transient sodium (NaT), persistent sodium (NaP), delayed-rectifier potassium [K(DR)], Ca²⁺-dependent K⁺ [K(AHP)] and L-type calcium (Ca_L) channels were included in the models. Simulation results showed that (1) Strengthening synaptic inputs increased the number of recruitments in all three types of motoneurons following the size principle. (2) Increasing NaT or NaP or decreasing K(DR) or K(AHP) lowered rheobase of spike generation thus increased recruitment of motoneuron pools. (3) Decreasing R_in reduced recruitment in all three types of motoneurons. (4) The FF-type motoneuron pool, followed by FR- and S-type, were the most sensitive to increase of synaptic inputs, reduction of R_in, upregulation of NaT and NaP, and downregulation of K(DR) and K(AHP). (5) Increasing Ca_L enhanced overall discharge rate of motoneuron pools with little effect on the recruitment. Simulation results suggested that modulation of ionic channels altered the output of motoneuron pools with either modulating the number of recruited motoneurons or regulating the overall discharge rate of motoneuron pools. Multiple channels contributed to the recruitment of motoneurons with interaction of excitatory and inhibitory synaptic inputs during walking.

Estimation of self-sustained activity produced by persistent inward currents using firing rate profiles of multiple motor units in humans

Persistent inward calcium and sodium currents (I_P) activated during motoneuron recruitment help synaptic inputs maintain self-sustained firing until derecruitment. Here, we estimate the contribution of the I_P to self-sustained firing in human motoneurons of varying recruitment threshold by measuring the difference in synaptic input needed to maintain minimal firing once the I_P is fully activated compared with the larger synaptic input required to initiate firing before full I_P activation. Synaptic input to ≈20 dorsiflexor motoneurons simultaneously recorded during ramp contractions was estimated from firing profiles of motor units decomposed from high-density surface electromyography (EMG). To avoid errors introduced when using high-threshold units firing in their nonlinear range, we developed methods where the lowest threshold units firing linearly with force were used to construct a composite (control) unit firing rate profile to estimate synaptic input to higher threshold (test) units. The difference in the composite firing rate (synaptic input) at the time of test unit recruitment and derecruitment (ΔF = F_recruit - F_derecruit) was used to measure I_P amplitude that sustained firing. Test units with recruitment thresholds 1-30% of maximum had similar ΔF values, which likely included both slow and fast motor units activated by small and large motoneurons, respectively. This suggests that the portion of the I_P that sustains firing is similar across a wide range of motoneuron sizes.NEW & NOTEWORTHY A new method of estimating synaptic drive to multiple, simultaneously recorded motor units provides evidence that the portion of the depolarizing drive from persistent inward currents that contributes to self-sustained firing is similar across motoneurons of different sizes.

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.

Changes in motoneuron excitability during voluntary muscle activity in humans with spinal cord injury

The excitability of resting motoneurons increases following spinal cord injury (SCI). The extent to which motoneuron excitability changes during voluntary muscle activity in humans with SCI, however, remains poorly understood. To address this question, we measured F waves by using supramaximal electrical stimulation of the ulnar nerve at the wrist and cervicomedullary motor-evoked potentials (CMEPs) by using high-current electrical stimulation over the cervicomedullary junction in the first dorsal interosseous muscle at rest and during 5 and 30% of maximal voluntary contraction into index finger abduction in individuals with chronic cervical incomplete SCI and aged-matched control participants. We found higher persistence (number of F waves present in each set) and amplitude of F waves at rest in SCI compared with control participants. With increasing levels of voluntary contraction, the amplitude, but not the persistence, of F waves increased in both groups but to a lesser extent in SCI compared with control participants. Similarly, the CMEP amplitude increased in both groups but to a lesser extent in SCI compared with controls. These results were also found at matched absolutely levels of electromyographic activity, suggesting that these changes were not related to decreases in voluntary motor output after SCI. F-wave and CMEP amplitudes were positively correlated across conditions in both groups. These results support the hypothesis that the responsiveness of the motoneuron pool during voluntary activity decreases following SCI, which could alter the generation and strength of voluntary muscle contractions.NEW & NOTEWORTHY How the excitability of motoneurons changes during voluntary muscle activity in humans with spinal cord injury (SCI) remains poorly understood. We found that F-wave and cervicomedullary motor-evoked potential amplitude, outcomes reflecting motoneuronal excitability, increased during voluntary activity compared with rest in SCI participants but to a lesser extent that in controls. These results suggest that the responsiveness of motoneurons during voluntary activity decreases following SCI, which might affect functionally relevant plasticity after the injury.

Investigating the Effect of Persistent Inward Currents on Motor Unit Firing Rates and Beta-Band Coherence in a Model of the First Dorsal Interosseous Muscle

Neuromodulatory drive resulting in the generation of persistent inward currents (PICs) within motoneuron dendrites has been demonstrated to introduce nonlinearities into the motoneuron input-output function for a given motor command. It is less understood, however, as to what role PICs play during voluntary contractions or on the correlation between motoneuron firings arising as a result of common synaptic inputs to the motoneuron pool. To examine this, a biophysical model of the motoneuron pool representing the first dorsal interosseous (FDI) muscle was used to simulate the effects of PICs on motor unit firing patterns and beta-band (15-30 Hz) motor unit coherence at 20, 30, and 40 percent of maximum voluntary contraction (MVC). The contribution of PICs at each MVC was quantified by calculating the difference in the mean firing rate of each motoneuron within the pool and assessing changes in the mean firing rate distribution and motor unit coherence with and without PICs present. The results of the current study demonstrated that increased activation of PICs progressively reduced motor unit coherence, however, changes in coherence were modest when investigating activation levels consistent with experimentally observed mean motor unit firing rates in the FDI muscle during isometric voluntary contraction.

Nonlinear Input-Output Functions of Motoneurons

All movements are generated by the activation of motoneurons, and hence their input-output properties define the final step in processing of all motor commands. A major challenge to understanding this transformation has been the striking nonlinear behavior of motoneurons conferred by the activation of persistent inward currents (PICs) mediated by their voltage-gated Na+ and Ca2+ channels. In this review, we focus on the contribution that these PICs make to motoneuronal discharge and how the nonlinearities they engender impede the construction of a comprehensive model of motor control.

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.

Understanding exercise-dependent plasticity of motoneurons using intracellular and intramuscular approaches

Spinal motoneurons (MN) exhibit exercise-dependent adaptations to increased activity, such as exercise and locomotion, as well as decreased activity associated with disuse, spinal cord injury, and aging. The development of several experimental approaches, in both human and animal models, has contributed significantly to our understanding of this plasticity. The purpose of this review is to summarize how intracellular recordings in an animal model and motor unit recordings in a human model have, together, contributed to our current understanding of exercise-dependent MN plasticity. These approaches and techniques will allow neuroscientists to continue to advance our understanding of MN physiology and the plasticity of the "final common path" of the motor system, and to design experiments to answer the critical questions that are emerging in this field.

Gain control in the sensorimotor system

Coordinated movement depends on constant interaction between neural circuits that produce motor output and those that report sensory consequences. Fundamental to this process are mechanisms for controlling the influence that sensory signals have on motor pathways - for example, reducing feedback gains when they are disruptive and increasing gains when advantageous. Sensory gain control comes in many forms and serves diverse purposes - in some cases sensory input is attenuated to maintain movement stability and filter out irrelevant or self-generated signals, or enhanced to facilitate salient signals for improved movement execution and adaptation. The ubiquitous presence of sensory gain control across species at multiple levels of the nervous system reflects the importance of tuning the impact that feedback information has on behavioral output.

Variability in common synaptic input to motor neurons modulates both force steadiness and pegboard time in young and older adults

KEY POINTS

The fluctuations in force during a steady isometric contraction (force steadiness) are associated with oscillations in common synaptic input to the involved motor neurons. Decreases in force steadiness are associated with increases in pegboard times in older adults, although a mechanism for this link has not been established. We used a state-space model to estimate the variability in common synaptic input to motor neurons during steady, isometric contractions. The estimate of common synaptic input was derived from the discharge times of motor units as recorded with high-density surface electrodes. We found that the variability in common synaptic input to motor neurons modulates force steadiness for young and older adults, as well as pegboard time for older adults.

ABSTRACT

We investigated the associations between grooved pegboard times, force steadiness (coefficient of variation for force) and variability in an estimate of the common synaptic input to motor neurons innervating the wrist extensor muscles during steady contractions performed by young and older adults. The discharge times of motor units were derived from recordings obtained with high-density surface electrodes when participants performed steady isometric contractions at 10% and 20% of maximal voluntary contraction force. The steady contractions were performed with a pinch grip and wrist extension, both independently (single action) and concurrently (double action). The variance in common synaptic input to motor neurons was estimated with a state-space model of the latent common input dynamics. There was a statistically significant association between the coefficient of variation for force during the steady contractions and the estimated variance in common synaptic input in young (r2  = 0.31) and older (r2  = 0.39) adults, although not between either the mean or the coefficient of variation for interspike interval of single motor units with the coefficient of variation for force. Moreover, the estimated variance in common synaptic input during the double-action task with the wrist extensors at the 20% target was significantly associated with grooved pegboard time (r2  = 0.47) for older adults but not young adults. These findings indicate that longer pegboard times of older adults were associated with worse force steadiness and greater fluctuations in the estimated common synaptic input to motor neurons during steady contractions.

Robust and accurate decoding of motoneuron behaviour and prediction of the resulting force output

KEY POINTS

The spinal alpha motoneuron is the only cell in the human CNS whose discharge can be routinely recorded in humans. We have reengineered motor unit collection and decomposition approaches, originally developed in humans, to measure the neural drive to muscle and estimate muscle force generation in the in vivo cat model. Experimental, computational, and predictive approaches are used to demonstrate the validity of this approach across a wide range of modes to activate the motor pool. The utility of this approach is shown through the ability to track individual motor units across trials, allowing for better predictions of muscle force than the electromyography signal, and providing insights in to the stereotypical discharge characteristics in response to synaptic activation of the motor pool. This approach now allows for a direct link between the intracellular data of single motoneurons, the discharge properties of motoneuron populations, and muscle force generation in the same preparation.

ABSTRACT

The discharge of a spinal alpha motoneuron and the resulting contraction of its muscle fibres represents the functional quantum of the motor system. Recent advances in the recording and decomposition of the electromyographic signal allow for the identification of several tens of concurrently active motor units. These detailed population data provide the potential to achieve deep insights into the synaptic organization of motor commands. Yet most of our understanding of the synaptic input to motoneurons is derived from intracellular recordings in animal preparations. Thus, it is necessary to extend the new electrode and decomposition methods to recording of motor unit populations in these same preparations. To achieve this goal, we use high-density electrode arrays and decomposition techniques, analogous to those developed for humans, to record and decompose the activity of tens of concurrently active motor units in a hindlimb muscle in the in vivo cat. Our results showed that the decomposition method in this animal preparation was highly accurate, with conventional two-source validation providing rates of agreement equal to or superior to those found in humans. Multidimensional reconstruction of the motor unit action potential provides the ability to accurately track the same motor unit across multiple contractions. Additionally, correlational analyses demonstrate that the composite spike train provides better estimates of whole muscle force than conventional estimates obtained from the electromyographic signal. Lastly, stark differences are observed between the modes of activation, in particular tendon vibration produced quantal interspike intervals at integer multiples of the vibration period.

The effects of model composition design choices on high-fidelity simulations of motoneuron recruitment and firing behaviors

OBJECTIVE

Computational models often require tradeoffs, such as balancing detail with efficiency; yet optimal balance should incorporate sound design features that do not bias the results of the specific scientific question under investigation. The present study examines how model design choices impact simulation results.

APPROACH

We developed a rigorously-validated high-fidelity computational model of the spinal motoneuron pool to study three long-standing model design practices which have yet to be examined for their impact on motoneuron recruitment, firing rate, and force simulations. The practices examined were the use of: (1) generic cell models to simulate different motoneuron types, (2) discrete property ranges for different motoneuron types, and (3) biological homogeneity of cell properties within motoneuron types.

MAIN RESULTS

Our results show that each of these practices accentuates conditions of motoneuron recruitment based on the size principle, and minimizes conditions of mixed and reversed recruitment orders, which have been observed in animal and human recordings. Specifically, strict motoneuron orderly size recruitment occurs, but in a compressed range, after which mixed and reverse motoneuron recruitment occurs due to the overlap in electrical properties of different motoneuron types. Additionally, these practices underestimate the motoneuron firing rates and force data simulated by existing models.

SIGNIFICANCE

Our results indicate that current modeling practices increase conditions of motoneuron recruitment based on the size principle, and decrease conditions of mixed and reversed recruitment order, which, in turn, impacts the predictions made by existing models on motoneuron recruitment, firing rate, and force. Additionally, mixed and reverse motoneuron recruitment generated higher muscle force than orderly size motoneuron recruitment in these simulations and represents one potential scheme to increase muscle efficiency. The examined model design practices, as well as the present results, are applicable to neuronal modeling throughout the nervous system.

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.

Altered Neuromodulatory Drive May Contribute to Exaggerated Tonic Vibration Reflexes in Chronic Hemiparetic Stroke

Exaggerated stretch-sensitive reflexes are a common finding in elbow flexors of the contralesional arm in chronic hemiparetic stroke, particularly when muscles are not voluntarily activated prior to stretch. Previous investigations have suggested that this exaggeration could arise either from an abnormal tonic ionotropic drive to motoneuron pools innervating the paretic limbs, which could bring additional motor units near firing threshold, or from an increased influence of descending monoaminergic neuromodulatory pathways, which could depolarize motoneurons and amplify their responses to synaptic inputs. However, previous investigations have been unable to differentiate between these explanations, leaving the source(s) of this excitability increase unclear. Here, we used tonic vibration reflexes (TVRs) during voluntary muscle contractions of increasing magnitude to infer the sources of spinal motor excitability in individuals with chronic hemiparetic stroke. We show that when the paretic and non-paretic elbow flexors are preactivated to the same percentage of maximum prior to vibration, TVRs remain significantly elevated in the paretic arm. We also show that the rate of vibration-induced torque development increases as a function of increasing preactivation in the paretic limb, even though the amplitude of vibration-induced torque remains conspicuously unchanged as preactivation increases. It is highly unlikely that these findings could be explained by a source that is either purely ionotropic or purely neuromodulatory, because matching preactivation should control for the effects of a potential ionotropic drive (and lead to comparable tonic vibration reflex responses between limbs), while a purely monoaminergic mechanism would increase reflex magnitude as a function of preactivation. Thus, our results suggest that increased excitability of motor pools innervating the paretic limb post-stroke is likely to arise from both ionotropic and neuromodulatory mechanisms.

A latent low-dimensional common input drives a pool of motor neurons: a probabilistic latent state-space model

Motor neurons appear to be activated with a common input signal that modulates the discharge activity of all neurons in the motor nucleus. It has proven difficult for neurophysiologists to quantify the variability in a common input signal, but characterization of such a signal may improve our understanding of how the activation signal varies across motor tasks. Contemporary methods of quantifying the common input to motor neurons rely on compiling discrete action potentials into continuous time series, assuming the motor pool acts as a linear filter, and requiring signals to be of sufficient duration for frequency analysis. We introduce a space-state model in which the discharge activity of motor neurons is modeled as inhomogeneous Poisson processes and propose a method to quantify an abstract latent trajectory that represents the common input received by motor neurons. The approach also approximates the variation in synaptic noise in the common input signal. The model is validated with four data sets: a simulation of 120 motor units, a pair of integrate-and-fire neurons with a Renshaw cell providing inhibitory feedback, the discharge activity of 10 integrate-and-fire neurons, and the discharge times of concurrently active motor units during an isometric voluntary contraction. The simulations revealed that a latent state-space model is able to quantify the trajectory and variability of the common input signal across all four conditions. When compared with the cumulative spike train method of characterizing common input, the state-space approach was more sensitive to the details of the common input current and was less influenced by the duration of the signal. The state-space approach appears to be capable of detecting rather modest changes in common input signals across conditions.NEW & NOTEWORTHY We propose a state-space model that explicitly delineates a common input signal sent to motor neurons and the physiological noise inherent in synaptic signal transmission. This is the first application of a deterministic state-space model to represent the discharge characteristics of motor units during voluntary contractions.

Modification of Spastic Stretch Reflexes at the Elbow by Flexion Synergy Expression in Individuals With Chronic Hemiparetic Stroke

OBJECTIVE

To systematically characterize the effect of flexion synergy expression on the manifestation of elbow flexor stretch reflexes poststroke, and to relate these findings to elbow flexor stretch reflexes in individuals without neurologic injury.

DESIGN

Controlled cohort study.

SETTING

Academic medical center.

PARTICIPANTS

Participants (N=20) included individuals with chronic hemiparetic stroke (n=10) and a convenience sample of individuals without neurologic or musculoskeletal injury (n=10).

INTERVENTIONS

Participants with stroke were interfaced with a robotic device that precisely manipulated flexion synergy expression (by regulating shoulder abduction loading) while delivering controlled elbow extension perturbations over a wide range of velocities. This device was also used to elicit elbow flexor stretch reflexes during volitional elbow flexor activation, both in the cohort of individuals with stroke and in a control cohort. In both cases, the amplitude of volitional elbow flexor preactivation was matched to that generated involuntarily during flexion synergy expression.

MAIN OUTCOME MEASURES

The amplitude of short- and long-latency stretch reflexes in the biceps brachii, assessed by electromyography, and expressed as a function of background muscle activation and stretch velocity.

RESULTS

Increased shoulder abduction loading potentiated elbow flexor stretch reflexes via flexion synergy expression in the paretic arm. Compared with stretch reflexes in individuals without neurologic injury, paretic reflexes were larger at rest but were approximately equal to control muscles at matched levels of preactivation.

CONCLUSIONS

Because flexion synergy expression modifies stretch reflexes in involved muscles, interventions that reduce flexion synergy expression may confer the added benefit of reducing spasticity during functional use of the arm.

A motor unit-based model of muscle fatigue

Muscle fatigue is a temporary decline in the force and power capacity of skeletal muscle resulting from muscle activity. Because control of muscle is realized at the level of the motor unit (MU), it seems important to consider the physiological properties of motor units when attempting to understand and predict muscle fatigue. Therefore, we developed a phenomenological model of motor unit fatigue as a tractable means to predict muscle fatigue for a variety of tasks and to illustrate the individual contractile responses of MUs whose collective action determines the trajectory of changes in muscle force capacity during prolonged activity. An existing MU population model was used to simulate MU firing rates and isometric muscle forces and, to that model, we added fatigue-related changes in MU force, contraction time, and firing rate associated with sustained voluntary contractions. The model accurately estimated endurance times for sustained isometric contractions across a wide range of target levels. In addition, simulations were run for situations that have little experimental precedent to demonstrate the potential utility of the model to predict motor unit fatigue for more complicated, real-world applications. Moreover, the model provided insight into the complex orchestration of MU force contributions during fatigue, that would be unattainable with current experimental approaches.

The potential for understanding the synaptic organization of human motor commands via the firing patterns of motoneurons

Motoneurons are unique in being the only neurons in the CNS whose firing patterns can be easily recorded in human subjects. This is because of the one-to-one relationship between the motoneuron and muscle cell behavior. It has long been appreciated that the connection of motoneurons to their muscle fibers allows their action potentials to be amplified and recorded, but only recently has it become possible to simultaneously record the firing pattern of many motoneurons via array electrodes placed on the skin. These firing patterns contain detailed information about the synaptic organization of motor commands to the motoneurons. This review focuses on parameters in these firing patterns that are directly linked to specific features of this organization. It is now well established that motor commands consist of three components, excitation, inhibition, and neuromodulation; the importance of the third component has become increasingly evident. Firing parameters linked to each of the three components are discussed, along with consideration of potential limitations in their utility for understanding the underlying organization of motor commands. Future work based on realistic computer simulations of motoneurons may allow quantitative "reverse engineering" of human motoneuron firing patterns to provide good estimates of the relative amplitudes and temporal patterns of all three components of motor commands.