Model-based prediction of fusimotor activity and its effect on muscle spindle activity during voluntary wrist movements

Linking cortex and contraction-Integrating models along the corticomuscular pathway

Computational models of the neuromusculoskeletal system provide a deterministic approach to investigate input-output relationships in the human motor system. Neuromusculoskeletal models are typically used to estimate muscle activations and forces that are consistent with observed motion under healthy and pathological conditions. However, many movement pathologies originate in the brain, including stroke, cerebral palsy, and Parkinson's disease, while most neuromusculoskeletal models deal exclusively with the peripheral nervous system and do not incorporate models of the motor cortex, cerebellum, or spinal cord. An integrated understanding of motor control is necessary to reveal underlying neural-input and motor-output relationships. To facilitate the development of integrated corticomuscular motor pathway models, we provide an overview of the neuromusculoskeletal modelling landscape with a focus on integrating computational models of the motor cortex, spinal cord circuitry, α-motoneurons  and skeletal muscle in regard to their role in generating voluntary muscle contraction. Further, we highlight the challenges and opportunities associated with an integrated corticomuscular pathway model, such as challenges in defining neuron connectivities, modelling standardisation, and opportunities in applying models to study emergent behaviour. Integrated corticomuscular pathway models have applications in brain-machine-interaction, education, and our understanding of neurological disease.

Feedforward modulation of gamma motor neuron activity can improve motor command accuracy

Objective. Coactivation of gamma and alpha motor neuron activity ensures that muscle spindle responsiveness is maintained during muscle contractions. However, some evidence suggests that the activity of gamma motor neurons is phase-advanced with respect to that of alpha motor neurons during manual control tasks. We hypothesized that this might be a deliberate control strategy to maximize movement accuracy.Approach. Using a computational model of the neural activation of a muscle and its type Ia sensory feedback to the motor neurons, we systematically investigated the impact of the phase difference between oscillatory descending input to alpha and dynamic gamma motor neurons. Specifically, the amplification of the alpha motor neuron drive to the muscle was investigated as a function of the frequency of the synaptic input (1-9 Hz individually or superimposed) and the alpha-gamma phase difference (0-2π).Main results. Simulation results showed that when the phase advance of the dynamic gamma drive resulted in delays between muscle velocity and type Ia afferent feedback similar to those previously observed experimentally, low-frequency components (1 and 2 Hz) of the motor neuron synaptic input were amplified (gain up to 1.7). On the other hand, synaptic input at higher frequencies was unaffected.Significance. This finding suggests that by imposing a phase advance of the input to dynamic gamma motor neurons, components of the neural drive usually associated with voluntary control are amplified. In this way, our study suggests that this neural strategy increases the control-to-neural-noise ratio of the motor output during movement.

Proprioceptive Recognition with Artificial Neural Networks Based on Organizations of Spinocerebellar Tract and Cerebellum

Muscle kinematics and kinetics are nonlinearly encoded by proprioceptors, and the changes in muscle length and velocity are integrated into Ia afferent. Besides, proprioceptive signals from multiple muscles are probably mixed in afferent pathways, which all lead to difficulties in proprioceptive recognition for the cerebellum. In this study, the artificial neural networks, whose organizations are biologically based on the spinocerebellar tract and cerebellum, are utilized to decode the proprioceptive signals. Consistent with the controversy of the proprioceptive division in the dorsal spinocerebellar tract, the spinocerebellar tract networks incorporated two distinct inferences, (1) the centralized networks, which mixed Ia, II, and Ib and processed them together; (2) the decentralized networks, which processed Ia, II, and Ib afferents separately. The cerebellar networks were based on the Marr-Albus model to recognize the kinematic states. The networks were trained by a specific movement, and the trained networks were subsequently required to predict kinematic states of six movements. The results demonstrated that the centralized networks, which were more consistent with the physiological findings in recent years, had better recognition accuracy than the decentralized networks, and the networks were still effective even when proprioceptive afferents from multiple muscles were integrated.

Physiological tremor increases when skeletal muscle is shortened: implications for fusimotor control

KEY POINTS

In tonic, isometric, plantarflexion contractions, physiological tremor increases as the ankle joint becomes plantarflexed. Modulation of physiological tremor as a function of muscle stretch differs from that of the stretch reflex amplitude. Amplitude of physiological tremor may be altered as a function of reflex pathway gains. Healthy humans likely increase their γ-static fusimotor drive when muscles shorten. Quantification of physiological tremor by manipulation of joint angle may be a useful experimental probe of afferent gains and/or the integrity of automatic fusimotor control.

ABSTRACT

The involuntary force fluctuations associated with physiological (as distinct from pathological) tremor are an unavoidable component of human motor control. While the origins of physiological tremor are known to depend on muscle afferentation, it is possible that the mechanical properties of muscle-tendon systems also affect its generation, amplification and maintenance. In this paper, we investigated the dependence of physiological tremor on muscle length in healthy individuals. We measured physiological tremor during tonic, isometric plantarflexion torque at 30% of maximum at three ankle angles. The amplitude of physiological tremor increased as calf muscles shortened in contrast to the stretch reflex whose amplitude decreases as muscle shortens. We used a published closed-loop simulation model of afferented muscle to explore the mechanisms responsible for this behaviour. We demonstrate that changing muscle lengths does not suffice to explain our experimental findings. Rather, the model consistently required the modulation of  γ-static fusimotor drive to produce increases in physiological tremor with muscle shortening - while successfully replicating the concomitant reduction in stretch reflex amplitude. This need to control γ-static fusimotor drive explicitly as a function of muscle length has important implications. First, it permits the amplitudes of physiological tremor and stretch reflex to be decoupled. Second, it postulates neuromechanical interactions that require length-dependent γ drive modulation to be independent from α drive to the parent muscle. Lastly, it suggests that physiological tremor can be used as a simple, non-invasive measure of the afferent mechanisms underlying healthy motor function, and their disruption in neurological conditions.

Proprioceptive Feedback through a Neuromorphic Muscle Spindle Model

Connecting biologically inspired neural simulations to physical or simulated embodiments can be useful both in robotics, for the development of a new kind of bio-inspired controllers, and in neuroscience, to test detailed brain models in complete action-perception loops. The aim of this work is to develop a fully spike-based, biologically inspired mechanism for the translation of proprioceptive feedback. The translation is achieved by implementing a computational model of neural activity of type Ia and type II afferent fibers of muscle spindles, the primary source of proprioceptive information, which, in mammals is regulated through fusimotor activation and provides necessary adjustments during voluntary muscle contractions. As such, both static and dynamic γ-motoneurons activities are taken into account in the proposed model. Information from the actual proprioceptive sensors (i.e., motor encoders) is then used to simulate the spindle contraction and relaxation, and therefore drive the neural activity. To assess the feasibility of this approach, the model is implemented on the NEST spiking neural network simulator and on the SpiNNaker neuromorphic hardware platform and tested on simulated and physical robotic platforms. The results demonstrate that the model can be used in both simulated and real-time robotic applications to translate encoder values into a biologically plausible neural activity. Thus, this model provides a completely spike-based building block, suitable for neuromorphic platforms, that will enable the development of sensory-motor closed loops which could include neural simulations of areas of the central nervous system or of low-level reflexes.

Emergence of gamma motor activity in an artificial neural network model of the corticospinal system

Muscle spindle discharge during active movement is a function of mechanical and neural parameters. Muscle length changes (and their derivatives) represent its primary mechanical, fusimotor drive its neural component. However, neither the action nor the function of fusimotor and in particular of γ-drive, have been clearly established, since γ-motor activity during voluntary, non-locomotor movements remains largely unknown. Here, using a computational approach, we explored whether γ-drive emerges in an artificial neural network model of the corticospinal system linked to a biomechanical antagonist wrist simulator. The wrist simulator included length-sensitive and γ-drive-dependent type Ia and type II muscle spindle activity. Network activity and connectivity were derived by a gradient descent algorithm to generate reciprocal, known target α-motor unit activity during wrist flexion-extension (F/E) movements. Two tasks were simulated: an alternating F/E task and a slow F/E tracking task. Emergence of γ-motor activity in the alternating F/E network was a function of α-motor unit drive: if muscle afferent (together with supraspinal) input was required for driving α-motor units, then γ-drive emerged in the form of α-γ coactivation, as predicted by empirical studies. In the slow F/E tracking network, γ-drive emerged in the form of α-γ dissociation and provided critical, bidirectional muscle afferent activity to the cortical network, containing known bidirectional target units. The model thus demonstrates the complementary aspects of spindle output and hence γ-drive: i) muscle spindle activity as a driving force of α-motor unit activity, and ii) afferent activity providing continuous sensory information, both of which crucially depend on γ-drive.

Human muscle spindle sensitivity reflects the balance of activity between antagonistic muscles

Muscle spindles are commonly considered as stretch receptors encoding movement, but the functional consequence of their efferent control has remained unclear. The "α-γ coactivation" hypothesis states that activity in a muscle is positively related to the output of its spindle afferents. However, in addition to the above, possible reciprocal inhibition of spindle controllers entails a negative relationship between contractile activity in one muscle and spindle afferent output from its antagonist. By recording spindle afferent responses from alert humans using microneurography, I show that spindle output does reflect antagonistic muscle balance. Specifically, regardless of identical kinematic profiles across active finger movements, stretch of the loaded antagonist muscle (i.e., extensor) was accompanied by increased afferent firing rates from this muscle compared with the baseline case of no constant external load. In contrast, spindle firing rates from the stretching antagonist were lowest when the agonist muscle powering movement (i.e., flexor) acted against an additional resistive load. Stepwise regressions confirmed that instantaneous velocity, extensor, and flexor muscle activity had a significant effect on spindle afferent responses, with flexor activity having a negative effect. Therefore, the results indicate that, as consequence of their efferent control, spindle sensitivity (gain) to muscle stretch reflects the balance of activity between antagonistic muscles rather than only the activity of the spindle-bearing muscle.