Motoneurons receive synaptic inputs from tens of thousands of connections that cause membrane potential to fluctuate continuously (synaptic noise), which introduces variability in discharge times of action potentials. We hypothesized that the influence of synaptic noise on force steadiness during voluntary contractions is limited to low muscle forces. The hypothesis was examined with an analytical description of transduction of motor unit spike trains into muscle force, a computational model of motor unit recruitment and rate coding, and experimental analysis of interspike interval variability during steady contractions with the abductor digiti minimi muscle. Simulations varied contraction force, level of synaptic noise, size of motor unit population, recruitment range, twitch contraction times, and level of motor unit short-term synchronization. Consistent with the analytical derivations, simulations and experimental data showed that force variability at target forces above a threshold was primarily due to low-frequency oscillations in neural drive, whereas the influence of synaptic noise was almost completely attenuated by two low-pass filters, one related to convolution of motoneuron spike trains with motor unit twitches (temporal summation) and the other attributable to summation of single motor unit forces (spatial summation). The threshold force above which synaptic noise ceased to influence force steadiness depended on recruitment range, size of motor unit population, and muscle contractile properties. This threshold was low (<10% of maximal force) for typical values of these parameters. Results indicate that motor unit recruitment and muscle properties of a typical muscle are tuned to limit the influence of synaptic noise on force steadiness to low forces and that the inability to produce a constant force during stronger contractions is mainly attributable to the common low-frequency oscillations in motoneuron discharge rates.

Tremor is one of the most prevalent movement disorders. There is a large proportion of patients (around 25%) in whom current treatments do not attain a significant tremor reduction. This paper proposes a tremor suppression strategy that detects tremor from the electromyographic signals of the muscles from which tremor originates and counteracts it by delivering electrical stimulation to the antagonist muscles in an out of phase manner. The detection was based on the iterative Hilbert transform and stimulation was delivered above the motor threshold (motor stimulation) and below the motor threshold (sensory stimulation). The system was tested on six patients with predominant wrist flexion/extension tremor (four with Parkinson disease and two with Essential tremor) and led to an average tremor reduction in the range of 46%-81% and 35%-48% across five patients when using the motor and sensory stimulation, respectively. In one patient, the system did not attenuate tremor. These results demonstrate that tremor attenuation might be achieved by delivering electrical stimulation below the motor threshold, preventing muscle fatigue and discomfort for the patients, which sets the basis for the development of an alternative treatment for tremor.

The amplitude of the surface EMG does not reach the level achieved during a maximal voluntary contraction force at the end of a sustained, submaximal contraction, despite near-maximal levels of voluntary effort. The depression of EMG amplitude may be explained by several neural and muscular adjustments during fatiguing contractions, including decreased net neural drive to the muscle, changes in the shape of the motor unit action potentials, and EMG amplitude cancellation. The changes in these parameters for the entire motor unit pool, however, cannot be measured experimentally. The present study used a computational model to simulate the adjustments during sustained isometric contractions and thereby determine the relative importance of these factors in explaining the submaximal levels of EMG amplitude at task failure. The simulation results indicated that the amount of amplitude cancellation in the simulated EMG (∼40%) exhibited a negligible change during the fatiguing contractions. Instead, the main determinant of the submaximal EMG amplitude at task failure was a decrease in muscle activation (number of muscle fiber action potentials), due to a reduction in the net synaptic input to motor neurons, with a lesser contribution from changes in the shape of the motor unit action potentials. Despite the association between the submaximal EMG amplitude and reduced muscle activation, the deficit in EMG amplitude at task failure was not consistently associated with the decrease in neural drive (number of motor unit action potentials) to the muscle. This indicates that the EMG amplitude cannot be used as an index of neural drive.

This paper presents the fully automatic identification of motor unit spike trains from high-density surface electromyograms (EMG) in pathological tremor. First, a mathematical derivation is provided to theoretically prove the possibility of decomposing noise-free high-density surface EMG signals into motor unit spike trains with high correlation, which are typical of tremor contractions. Further, the proposed decomposition method is tested on simulated signals with different levels of noise and on experimental signals from 14 tremor-affected patients. In the case of simulated tremor with central frequency ranging from 5 Hz to 11 Hz and signal-to-noise ratio of 20 dB, the method identified ∼8 motor units per contraction with sensitivity in spike timing identification ≥ 95% and false alarm and miss rates ≤ 5%. In experimental signals, the number of identified motor units varied substantially (range 0-21) across patients and contraction types, as expected. The behaviour of the identified motor units was consistent with previous data obtained by intramuscular EMG decomposition. These results demonstrate for the first time the possibility of a fully non-invasive investigation of motor unit behaviour in tremor-affected patients. The method provides a new means for physiological investigations of pathological tremor.

The ability to produce rapid forces requires quick motor unit recruitment, high motor unit discharge rates, and fast motor unit force twitches. The relative importance of these parameters for maximum rate of force development (RFD), however, is poorly understood. In this study, we systematically investigated these relationships using a computational model of motor unit pool activity and force. Across simulations, neural and muscular properties were systematically varied in experimentally observed ranges. Motor units were recruited over an interval starting from contraction onset (range: 22–233 ms). Upon recruitment, discharge rates declined from an initial rate (range: 89–212 pulses per second), with varying likelihood of doublet (interspike interval of 3 ms; range: 0–50%). Finally, muscular adaptations were modeled by changing average twitch contraction time (range: 42–78 ms). Spectral analysis showed that the effective neural drive to the simulated muscle had smaller bandwidths than the average motor unit twitch, indicating that the bandwidth of the motor output, and thus the capacity for explosive force, was limited mainly by neural properties. The simulated RFD increased by 1,050 ± 281% maximal voluntary contraction force per second from the longest to the shortest recruitment interval. This effect was more than fourfold higher than the effect of increasing the initial discharge rate, more than fivefold higher than the effect of increasing the chance of doublets, and more than sixfold higher than the effect of decreasing twitch contraction times. The simulated results suggest that the physiological variation of the rate by which motor units are recruited during ballistic contractions is the main determinant for the variability in RFD across individuals. NEW & NOTEWORTHY An important limitation of human performance is the ability to generate explosive movements by means of rapid development of muscle force. The physiological determinants of this ability, however, are poorly understood. In this study, we show using extensive simulations that the rate by which motor units are recruited is the main limiting factor for maximum rate of force development.