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.

Motoneuron-driven computational muscle modelling with motor unit resolution and subject-specific musculoskeletal anatomy.

The computational simulation of human voluntary muscle contraction is possible with EMG-driven Hill-type models of whole muscles. Despite impactful applications in numerous fields, the neuromechanical information and the physiological accuracy such models provide remain limited because of multiscale simplifications that limit comprehensive description of muscle internal dynamics during contraction. We addressed this limitation by developing a novel motoneuron-driven neuromuscular model, that describes the force-generating dynamics of a population of individual motor units, each of which was described with a Hill-type actuator and controlled by a dedicated experimentally derived motoneuronal control. In forward simulation of human voluntary muscle contraction, the model transforms a vector of motoneuron spike trains decoded from high-density EMG signals into a vector of motor unit forces that sum into the predicted whole muscle force. The motoneuronal control provides comprehensive and separate descriptions of the dynamics of motor unit recruitment and discharge and decodes the subject's intention. The neuromuscular model is subject-specific, muscle-specific, includes an advanced and physiological description of motor unit activation dynamics, and is validated against an experimental muscle force. Accurate force predictions were obtained when the vector of experimental neural controls was representative of the discharge activity of the complete motor unit pool. This was achieved with large and dense grids of EMG electrodes during medium-force contractions or with computational methods that physiologically estimate the discharge activity of the motor units that were not identified experimentally. This neuromuscular model advances the state-of-the-art of neuromuscular modelling, bringing together the fields of motor control and musculoskeletal modelling, and finding applications in neuromuscular control and human-machine interfacing research.

Larger and Denser: An Optimal Design for Surface Grids of EMG Electrodes to Identify Greater and More Representative Samples of Motor Units.

The spinal motor neurons are the only neural cells whose individual activity can be noninvasively identified. This is usually done using grids of surface electromyographic (EMG) electrodes and source separation algorithms; an approach called EMG decomposition. In this study, we combined computational and experimental analyses to assess how the design parameters of grids of electrodes influence the number and the properties of the identified motor units. We first computed the percentage of motor units that could be theoretically discriminated within a pool of 200 simulated motor units when decomposing EMG signals recorded with grids of various sizes and interelectrode distances (IEDs). Increasing the density, the number of electrodes, and the size of the grids, increased the number of motor units that our decomposition algorithm could theoretically discriminate, i.e., up to 83.5% of the simulated pool (range across conditions: 30.5-83.5%). We then identified motor units from experimental EMG signals recorded in six participants with grids of various sizes (range: 2-36 cm2) and IED (range: 4-16 mm). The configuration with the largest number of electrodes and the shortest IED maximized the number of identified motor units (56 ± 14; range: 39-79) and the percentage of early recruited motor units within these samples (29 ± 14%). Finally, the number of identified motor units further increased with a prototyped grid of 256 electrodes and an IED of 2 mm. Taken together, our results showed that larger and denser surface grids of electrodes allow to identify a more representative pool of motor units than currently reported in experimental studies.

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.

Mathematical relationships between spinal motoneuron properties.

Our understanding of the behaviour of spinal alpha-motoneurons (MNs) in mammals partly relies on our knowledge of the relationships between MN membrane properties, such as MN size, resistance, rheobase, capacitance, time constant, axonal conduction velocity, and afterhyperpolarization duration. We reprocessed the data from 40 experimental studies in adult cat, rat, and mouse MN preparations to empirically derive a set of quantitative mathematical relationships between these MN electrophysiological and anatomical properties. This validated mathematical framework, which supports past findings that the MN membrane properties are all related to each other and clarifies the nature of their associations, is besides consistent with the Henneman's size principle and Rall's cable theory. The derived mathematical relationships provide a convenient tool for neuroscientists and experimenters to complete experimental datasets, explore the relationships between pairs of MN properties never concurrently observed in previous experiments, or investigate inter-mammalian-species variations in MN membrane properties. Using this mathematical framework, modellers can build profiles of inter-consistent MN-specific properties to scale pools of MN models, with consequences on the accuracy and the interpretability of the simulations.

Muscles receive their instructions through electrical signals carried by tens or hundreds of cells connected to the command centers of the body. These 'alpha-motoneurons' have various sizes and electrical characteristics which affect how they transmit signals. Previous experiments have shown that these properties are linked; for instance, larger motoneurons transfer electrical signals more quickly. The exact nature of the mathematical relationships between these characteristics, however, remains unclear. This limits our understanding of the behaviour of motoneurons from experimental data. To identify the equations linking eight motoneuron properties, Caillet et al. analysed published datasets from experimental studies on cat motoneurons. This approach uncovered simple mathematical associations: in fact, only one characteristic needs to be measured experimentally to calculate all the other properties. The relationships identified were also consistent with previously accepted approaches for modelling motoneuron activity. Caillet et al. then validated this mathematical framework with data from studies on rodents, showing that some of the equations hold true for different mammals. This work offers a quick and easy way for researchers to calculate the characteristics of a motoneuron based on a single observation. This will allow non-measured properties to be added to experimental datasets, and it could help to uncover the diversity of motoneurons at work within a population.