Evaluation of linear and non-linear activation dynamics models for insect muscle
Authors: Nalin Harischandra, Anthony J Clare, Jure Zakotnik, Laura M L Blackburn, Tom Matheson, Volker Dürr
Journal: PLOS Computational Biology
Accepted: 2019-09-25T00:00:00Z


In computational modelling of sensory-motor control, the dynamics of muscle contraction is an important determinant of movement timing and joint stiffness. This is particularly so in animals with many slow muscles, as is the case in insects—many of which are important models for sensory-motor control. A muscle model is generally used to transform motoneuronal input into muscle force. Although standard models exist for vertebrate muscle innervated by many motoneurons, there is no agreement on a parametric model for single motoneuron stimulation of invertebrate muscle. Although several different models have been proposed, they have never been evaluated using a common experimental data set. We evaluate five models for isometric force production of a well-studied model system: the locust hind leg tibial extensor muscle. The response of this muscle to motoneuron spikes is best modelled as a non-linear low-pass system. Linear first-order models can approximate isometric force time courses well at high spike rates, but they cannot account for appropriate force time courses at low spike rates. A linear third-order model performs better, but only non-linear models can account for frequency-dependent change of decay time and force potentiation at intermediate stimulus frequencies. Some of the differences among published models are due to differences among experimental data sets. We developed a comprehensive toolbox for modelling muscle activation dynamics, and optimised model parameters using one data set. The “Hatze-Zakotnik model” that emphasizes an accurate single-twitch time course and uses frequency-dependent modulation of the twitch for force potentiation performs best for the slow motoneuron. Frequency-dependent modulation of a single twitch works less well for the fast motoneuron. The non-linear “Wilson” model that optimises parameters to all data set parts simultaneously performs better here. Our open-access toolbox provides powerful tools for researchers to fit appropriate models to a range of insect muscles.
Author summary: Insects are important study organisms in the neuroscience of sensory-motor systems. Since the dynamics of muscle contraction and associated changes in force, torque or stiffness are central to our understanding of sensory-motor systems in general, the choice of the most appropriate model for insect muscle matters. Computational modelling of muscle properties typically separates activation dynamics from contraction dynamics. The former models the development of muscle force in response to motoneuron activity, whereas the latter describes how this force is affected by the current physical state of the muscle: its length and contraction velocity. We evaluate five published activation dynamics models for insect muscle. We explain differences between them, suggest how to decide which one to use, and provide an open-source toolbox for activation dynamics modelling. We further show that non-linear models are the best choice if: (i) the time course of a single muscle twitch is slow, (ii) the spike frequency ranges between one and thirty spikes per second, or (iii) the sensory-motor system tends to execute movements in a similar manner even if the demand on joint torque or stiffness changes.