Predicting electromyographic signals under realistic conditions using a multiscale chemo–electro–mechanical finite element model

This paper presents a novel multiscale finite element-based framework for modelling electromyographic (EMG) signals. The framework combines (i) a biophysical description of the excitation–contraction coupling at the half-sarcomere level, (ii) a model of the action potential (AP) propagation along muscle fibres, (iii) a continuum-mechanical formulation of force generation and deformation of the muscle, and (iv) a model for predicting the intramuscular and surface EMG. Owing to the biophysical description of the half-sarcomere, the model inherently accounts for physiological properties of skeletal muscle. To demonstrate this, the influence of membrane fatigue on the EMG signal during sustained contractions is investigated. During a stimulation period of 500 ms at 100 Hz, the predicted EMG amplitude decreases by 40% and the AP propagation velocity decreases by 15%. Further, the model can take into account contraction-induced deformations of the muscle. This is demonstrated by simulating fixed-length contractions of an idealized geometry and a model of the human tibialis anterior muscle (TA). The model of the TA furthermore demonstrates that the proposed finite element model is capable of simulating realistic geometries, complex fibre architectures, and can include different types of heterogeneities. In addition, the TA model accounts for a distributed innervation zone, different fibre types and appeals to motor unit discharge times that are based on a biophysical description of the α motor neurons.

[1]  H. Milner-Brown,et al.  Muscle membrane excitation and impulse propagation velocity are reduced during muscle fatigue , 1986, Muscle & nerve.

[2]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[3]  Todd A. Kuiken,et al.  Volume conduction in an anatomically based surface EMG model , 2004, IEEE Transactions on Biomedical Engineering.

[4]  Hiroshi Nishiura,et al.  Age-Dependent Estimates of the Epidemiological Impact of Pandemic Influenza (H1N1-2009) in Japan , 2013, Comput. Math. Methods Medicine.

[5]  Bernd Markert,et al.  A general polyconvex strain‐energy function for fiber‐reinforced materials , 2005 .

[6]  S H Bryant,et al.  Cable properties of external intercostal muscle fibres from myotonic and nonmyotonic goats , 1969, The Journal of physiology.

[7]  L. Mesin,et al.  Advances in surface electromyographic signal simulation with analytical and numerical descriptions of the volume conductor , 2004, Medical and Biological Engineering and Computing.

[8]  D. Farina,et al.  Decorrelation of cortical inputs and motoneuron output. , 2011, Journal of neurophysiology.

[9]  K B Campbell,et al.  Stiffness-distortion sarcomere model for muscle simulation. , 1999, Journal of applied physiology.

[10]  Oliver Röhrle,et al.  Modeling the Chemoelectromechanical Behavior of Skeletal Muscle Using the Parallel Open-Source Software Library OpenCMISS , 2013, Comput. Math. Methods Medicine.

[11]  Zhaohua Ding,et al.  Quantitative diffusion tensor MRI-based fiber tracking of human skeletal muscle. , 2007, Journal of applied physiology.

[12]  D. Winter,et al.  Models of recruitment and rate coding organization in motor-unit pools. , 1993, Journal of neurophysiology.

[13]  Hisashi Kawai,et al.  Investigation of Innervation Zone Shift with Continuous Dynamic Muscle Contraction , 2013, Comput. Math. Methods Medicine.

[14]  Luca Mesin Analytical Generation Model of Surface Electromyogram for Multi Layer Volume Conductors , 2005 .

[15]  G V Dimitrov,et al.  Calculation of extracellular potentials produced by an inclined muscle fibre at a rectangular plate electrode. , 1999, Medical engineering & physics.

[16]  J. B. Davidson,et al.  A Physiologically Based, Multi-Scale Model of Skeletal Muscle Structure and Function , 2012, Front. Physio..

[17]  G Rau,et al.  The presence of unknown layer of skin and fat is an obstacle to a correct estimation of the motor unit size from surface detected potentials. , 2002, Electromyography and clinical neurophysiology.

[18]  V. Spitzer,et al.  The visible human dataset: The anatomical platform for human simulation , 1998, The Anatomical record.

[19]  Luca Mesin,et al.  Volume conductor models in surface electromyography: Computational techniques , 2013, Comput. Biol. Medicine.

[20]  Todd A. Kuiken,et al.  A multiple-layer finite-element model of the surface EMG signal , 2002, IEEE Transactions on Biomedical Engineering.

[21]  Andrew J. Pullan,et al.  Solving the cardiac bidomain equations for discontinuous conductivities , 2006, IEEE Transactions on Biomedical Engineering.

[22]  Roger M. Enoka,et al.  Physiology of the motor neuron and the motor unit , 2004 .

[23]  Antonio I Cuesta-Vargas,et al.  Electromyography and sonomyography analysis of the tibialis anterior: a cross sectional study , 2014, Journal of Foot and Ankle Research.

[24]  Dario Farina,et al.  A surface EMG generation model with multilayer cylindrical description of the volume conductor , 2004, IEEE Transactions on Biomedical Engineering.

[25]  Roberto Merletti,et al.  The extraction of neural strategies from the surface EMG. , 2004, Journal of applied physiology.

[26]  D. Winter,et al.  Crosstalk in surface electromyography: Theoretical and practical estimates. , 1994, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[27]  P. Komi,et al.  Innervation zone shift at different levels of isometric contraction in the biceps brachii muscle. , 2009, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[28]  Jason M. DeFreitas,et al.  An examination of innervation zone movement with increases in isometric torque production , 2008, Clinical Neurophysiology.

[29]  John B. Davidson,et al.  A mathematical model of fatigue in skeletal muscle force contraction , 2007, Journal of Muscle Research and Cell Motility.

[30]  R. Enoka Eccentric contractions require unique activation strategies by the nervous system. , 1996, Journal of applied physiology.

[31]  John B. Davidson,et al.  Bridging Scales: A Three-Dimensional Electromechanical Finite Element Model of Skeletal Muscle , 2008, SIAM J. Sci. Comput..

[32]  W. Wallinga,et al.  Modelling action potentials and membrane currents of mammalian skeletal muscle fibres in coherence with potassium concentration changes in the T-tubular system , 1999, European Biophysics Journal.

[33]  Dario Farina,et al.  A finite element model for describing the effect of muscle shortening on surface EMG , 2006, IEEE Transactions on Biomedical Engineering.

[34]  Oliver Röhrle,et al.  A multiscale chemo-electro-mechanical skeletal muscle model to analyze muscle contraction and force generation for different muscle fiber arrangements , 2014, Front. Physiol..

[35]  G. Somjen,et al.  FUNCTIONAL SIGNIFICANCE OF CELL SIZE IN SPINAL MOTONEURONS. , 1965, Journal of neurophysiology.

[36]  Dario Farina,et al.  Decoding the neural drive to muscles from the surface electromyogram , 2010, Clinical Neurophysiology.

[37]  Heye Zhang,et al.  OpenCMISS: a multi-physics & multi-scale computational infrastructure for the VPH/Physiome project. , 2011, Progress in biophysics and molecular biology.

[38]  Aslak Tveito,et al.  Optimal monodomain approximations of the bidomain equations , 2007, Appl. Math. Comput..

[39]  Roberto Merletti,et al.  Atlas of Muscle Innervation Zones , 2012 .

[40]  G. Somjen,et al.  Excitability and inhibitability of motoneurons of different sizes. , 1965, Journal of neurophysiology.

[41]  Winnie Jensen,et al.  Relationship between grasping force and features of single-channel intramuscular EMG signals , 2009, Journal of Neuroscience Methods.

[42]  O Röhrle Simulating the Electro-Mechanical Behavior of Skeletal Muscles , 2010, Computing in Science & Engineering.

[43]  O. A. Nikitin,et al.  Neither high-pass filtering nor mathematical differentiation of the EMG signals can considerably reduce cross-talk. , 2002, Journal of electromyography and kinesiology : official journal of the International Society of Electrophysiological Kinesiology.

[44]  R. W. Lau,et al.  The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. , 1996, Physics in medicine and biology.

[45]  Rogério Rodrigues Lima Cisi,et al.  Simulation system of spinal cord motor nuclei and associated nerves and muscles, in a Web-based architecture , 2008, Journal of Computational Neuroscience.

[46]  Andrew J. Pullan,et al.  Mathematically Modelling the Electrical Activity of the Heart: From Cell to Body Surface and Back Again , 2005 .

[47]  M. Karhanek,et al.  An allosteric model of the molecular interactions of excitation- contraction coupling in skeletal muscle , 1993, The Journal of general physiology.

[48]  Lothar A. Blatter,et al.  Role of β-hydroxybutyrate, its polymer poly-β-hydroxybutyrate and inorganic polyphosphate in mammalian health and disease , 2014, Front. Physiol..

[49]  C. Gleissner Macrophage Phenotype Modulation by CXCL4 in Atherosclerosis , 2011, Front. Physio..

[50]  Xing Cai,et al.  On the Computational Complexity of the Bidomain and the Monodomain Models of Electrophysiology , 2006, Annals of Biomedical Engineering.

[51]  C. Juel,et al.  Muscle action potential propagation velocity changes during activity , 1988, Muscle & nerve.

[52]  D. Farina,et al.  Estimating motor unit discharge patterns from high-density surface electromyogram , 2009, Clinical Neurophysiology.

[53]  S. Gandevia,et al.  Measurement of muscle contraction with ultrasound imaging , 2003, Muscle & nerve.

[54]  Luca Mesin,et al.  Simulation of Surface EMG Signals for a Multilayer Volume Conductor With a Superficial Bone or Blood Vessel , 2008, IEEE Transactions on Biomedical Engineering.

[55]  Natalia A. Trayanova,et al.  Computational techniques for solving the bidomain equations in three dimensions , 2002, IEEE Transactions on Biomedical Engineering.

[56]  R. H. Adrian,et al.  Reconstruction of the action potential of frog sartorius muscle , 1973, The Journal of physiology.

[57]  Roberto Merletti,et al.  Electromyography. Physiology, engineering and non invasive applications , 2005 .

[58]  S. Baylor,et al.  Model of Sarcomeric Ca2+ Movements, Including ATP Ca2+ Binding and Diffusion, during Activation of Frog Skeletal Muscle , 1998, The Journal of general physiology.

[59]  S. Rush,et al.  Resistivity of Body Tissues at Low Frequencies , 1963, Circulation research.

[60]  Dario Farina,et al.  A novel approach for precise simulation of the EMG signal detected by surface electrodes , 2001, IEEE Trans. Biomed. Eng..