Anatomically based lower limb nerve model for electrical stimulation
© Kim et al; licensee BioMed Central Ltd. 2007
Received: 30 July 2007
Accepted: 17 December 2007
Published: 17 December 2007
Functional Electrical Stimulation (FES) is a technique that aims to rehabilitate or restore functionality of skeletal muscles using external electrical stimulation. Despite the success achieved within the field of FES, there are still a number of questions that remain unanswered. One way of providing input to the answers is through the use of computational models.
This paper describes the development of an anatomically based computer model of the motor neurons in the lower limb of the human leg and shows how it can be used to simulate electrical signal propagation from the beginning of the sciatic nerve to a skeletal muscle. One-dimensional cubic Hermite finite elements were used to represent the major portions of the lower limb nerves. These elements were fit to data that had been digitised using images from the Visible Man project. Nerves smaller than approximately 1 mm could not be seen in the images, and thus a tree-branching algorithm was used to connect the ends of the fitted nerve model to the respective skeletal muscle. To simulate electrical propagation, a previously published mammalian nerve model was implemented and solved on the anatomically based nerve mesh using a finite difference method. The grid points for the finite difference method were derived from the fitted finite element mesh. By adjusting the tree-branching algorithm, it is possible to represent different levels of motor-unit recruitment.
To illustrate the process of a propagating nerve stimulus to a muscle in detail, the above method was applied to the nerve tree that connects to the human semitendinosus muscle. A conduction velocity of 89.8 m/s was obtained for a 15 μm diameter nerve fibre. This signal was successfully propagated down the motor neurons to a selected group of motor units in the muscle.
An anatomically and physiologically based model of the posterior motor neurons in the human lower limb was developed. This model can be used to examine the effect of external stimulation on nerve and muscle activity, as may occur, for example, in the field of FES.
Functional Electrical Stimulation (FES) is a technique, using external electrical stimulation, which is capable of rehabilitating or restoring functionality of skeletal muscles that has been paralysed, for instance, by central nervous system lesions. The stimulation electrode can be located on the skin or implanted within the body to directly stimulate the muscles or nerves .
The direct stimulation of nerve trunks offers many advantages over transcutaneous or intramuscular stimulation. These advantages include the potential for controlling many muscles with a single implant, lower power requirements and the ability to place the electrode far from the contracting muscles . Also, further development on nerve electrodes, like FINE (flat interface nerve electrodes) , allow selective stimulation of fascicles in peripheral nerves, which could be important in selectively activating individual muscles as well as reducing muscle fatigue that often occurs with FES.
Despite the success of FES achieved to date, there are a number of questions that cannot yet be completely answered. For instance (i) how and where should the electrodes be arranged on the nerves to achieve a specific muscle response? (ii) what stimulation protocol should be used on a nerve to achieve a desired level of force while minimising the symptoms of muscle fatigue?
The answers to these questions strongly depend on how well one can predict the effects of external stimuli with respect to the initiation and propagation of nerve impulses within a nerve tree. An effective way of analysing such nerve responses is by appealing to a computational nerve model which is based on anatomical features and the electrophysiological behaviour of neurons.
Early nerve models used the Hodgkin-Huxley equations which are derived from experimental data on squid giant axons . A decade later, Fankenhaeuser and Huxley developed a new model, the so-called Fankenhaeuser-Huxley (or FH) model that was based on experimental data from myelinated nerve fibres of the toad Xenopus laevis . In 1976 McNeal developed the first nerve model for a myelinated fibre using stimulus electrodes outside a fibre. In 1985, the SENN model (Spatially Extended Nonlinear Node model), based on amphibian data, was presented by Reilly et al. . This model can accommodate arbitrary spatial and temporal distributions of the extracellular potentials and thus it makes it possible to simulate the behaviour of myelinated nerve fibres for different extracellular electrode configurations and stimulus protocols. In 1987, Sweeney et al.  published the first model for warm blooded nerves using the data of Chiu et al. , which stem from myelinated nerves of rabbits. This model has subsequently been called the CRRSS model, named after the investigators Chiu, Ritchie, Rogert, Stagg and Sweeney. Schwarz and Eikhof also presented their model, which was based on data from rat nerves . In 1992, a modified SENN model was presented by Frijns et al. . This model has adapted the SENN model to fit mammalian nerve fibre data by taking into account the influence of body temperature on nerve kinetics. In 2001 McIntyre et al.  described a model of mammalian nerve fibres using the concept of a double cable that had been earlier proposed by Halt and Clark .
While all the above mentioned models have previously been used within the field of FES [12–15], the CRRSS model is by far most popular one. This nerve model describes mammalian nerves more realisticly than earlier ones and it does not require the use of a double cable structure. For all these reasons this paper uses the CRRSS model in order to examine signal propagation in an anatomically based representation of the nerves in the human lower limb.
To the authors' knowledge, none of the previous works within the field of FES have used an anatomically based geometry of the nerve tree in the lower limb. The geometry of the nerve, however, plays a crucial role in analysing functional responses of muscles, e.g. the order of recruitment of muscles after an external source initiated a nerve stimulus. Therefore we have created an anatomically based geometric model of the posterior motor neurons in human lower limb to analyse the nerve responses more realistically and to allow us to look at the effect of external stimulation on nerve and muscle activity in the future.
Construction of the model of the lower limb nerves
The sciatic nerve is a major signal supplier to the legs and impulses conducted through its branches provide movement and feeling to the hip, thigh, knee, calf, ankle, foot and toes. The sciatic nerve starts from the lower spine and passes through the pelvic area and down through the hip and back of the leg branching into the tibial and common fibular nerves around the popliteal fissa area. The tibial nerve descends down to the feet along the back of the leg and the common fibular nerve is separated into the superficial and the deep fibular nerve descending along the front of the knee down to the feet [16, 17]. Our model of the nerve tree starts from the sciatic nerve and forms nerve branches (tibial, common fibular, superficial fibular, deep fibular nerves) and the branches to muscles – hamstring (semitendinosus, (part of) adductor magnus, semimembranosus, short biceps, long biceps), calf muscles (gastrocnemius, soleus), tibialis posterior, flexor digitorum longus, flexor hallucis longus, fibular longus, fibular bravis, tibialis anterior, extenxor digitorum longus, and the extenxor hallucis longus.
Digitising and fitting of main nerve trunk
Nerve entry points
The biggest challenge in constructing a complete nerve tree from the images of the VM was the fact that nerves smaller than 1 mm cannot be identified. Hence an algorithm needed to be developed that connects the main nerve trunk with the motor endplates of a muscle.
Motor endplates are the flattened ends of motor neurons that transmit neural impulses to a muscle. While the motor endplates cannot be determined directly from the VM images, their locations can be approximated using anatomical features specific to human skeletal muscle. In humans, muscle fibres run in most human leg muscles from tendon to tendon and the endplates are situated around the midpoint of the muscle fibre (the exceptions are the satorious and the gracillis muscles). Moreover, the innervating endplates form either a transverse band through the middle of the muscle or a concave endplate band depending on the shape of the specific muscle . This is also true for the semitendinosus muscle which is separated into two areas by an internal tendon. Thus the semitendinosus muscle has two bands of motor endplates and two separate nerve branches leading to the muscle. Using the geometrical models of the muscles in the lower limb the nerve entry points can hence be computed by determining the middle points of muscle fibres that run from tendon to tendon.
Each fibre in a nerve is capable of generating an action potential which can be propagated along its length without loss of amplitude, duration, or velocity. The individual fibres respond to stimulation according to the all-or-none principle: If the neuron does not reach the critical threshold level, then no action potential will be generated. On the other hand, if the threshold level is reached, an action potential of a fixed size occurs.
As an axon reaches the vicinity of the skeletal muscle fibre it loses its myelin sheath . The non-myelination and the small diameter of the nerve endings do both contribute to the reduction of the conduction velocity by a factor of between 1/50 and 1/100 when compared to the speed at which the impulse travels in the myelinated unbranched parent axon . In our model, however, the details of the nerve terminals were not modelled as the main focus of this study is on the conduction velocity of the unbranched axon. Hence, we assumed that the nerve fibre was myelinated along its entire length and the diameter of the fibre was constant along the whole nerve tree.
In the CRRSS model, the nodes of Ranvier are modeled by voltage-dependent sodium channels, voltage-independent leakage channels and nodal capacitance. There are no voltage-dependent potassium channels in the ionic current equations of the CRRSS model since these channel types are scarce at nodes of Ranvier in mammalian myelinated nerve fibres. Hence action potential repolarization in mammalian myelinated axons is attributed to rapid sodium inactivation and a large leakage current . A complete list of the CRRSS ionic currents equations together with the corresponding list of constants and variable names is given in APPENDIX A (See Additional file 1: appendixa).
With our proposed focus on external (extracellular) stimulation, the bidomain model was chosen to represent the electrical activity occurring along the nerve tree. The bidomain equations have the advantage that they explicitly model the extracellular potential and are given by
∇·((σ i + σ e ) ∇ φ e ) = -∇·(σ i ∇ V m ) + I stim , (1)
σ i : intracellular conductance
σ e : extracellular conductance
φ e : extracellular potential
V m : transmembrane potential
I ion : total ionic current, as given by the CRRSS model
C m : membrane capacitance
A m : the ratio of membrane surface area to the volume
I stim :extracelluar injected current
To use the governing two bidomain equations to model myelinated fibres, the left hand side of the equation (2) needed to be multiplied by l m /l n . Here, l m is a internodal length (myelinated part) and l n is a nodal gap length (node of Ranvier)(Figure 5(b)).
The modified bidomain equations were discretised using the same finite different scheme as proposed in . The resulting equations were then solved using an implicit numerical integration scheme (LSODA) [31, 32] with a uniform spatial discretisation and a fixed time step of 0.001 ms as implemented within CMISS . Solving the modified bidomain equations with values as given in  resulted in a slow conduction velocity. Therefore, the value of the axoplasm resistivity, the ratio of fibre diameter to axon diameter and the ratio of internodal length to axon diameter were changed in such a way that an appropriate action potential conduction velocity for humans was generated. In Waxman's review  of the conduction velocity in myelinated nerve fibres, he determined theoretically that to achieve maximal conduction velocity per axon diameter, the ratio of fibre diameter to axon diameter should be 0.6–0.7 and the ratio of the internodal length to axon diameter was found to be 100–200. These ratios change for different nerves and species.
Guided by the work in , we set the ratio of axon diameter to fibre diameter to be 0.67 and the ratio of internodal length to axon diameter to be 200.
Signal propagation to a muscle fibre
To illustrate the process of a propagating nerve stimulus to a muscle in detail, the above method was applied to the nerve tree that connected to the human semitendinosus muscle. The bidomain model was also used to simulate electrical propagation within the muscle fibres. Here though, the CRRSS model was replaced by the muscle cell model .
R i (1/σ i )
R e (1/σ e )
The ratio of membrane surface to volume(1/mm)
The conduction velocity in the main nerve trunk was calculated by running the simulation for 10 ms and measuring the distance which the stimulus has propagated from a point close to the beginning of the sciatic nerve. The calculated conduction velocity was 89.8 m/s. This result agrees well with previously published theoretical results of 90 m/s [10, 38]. It also agrees well with the experimental work of Boyd and Kalu , in which the relationship between conduction velocity and nerve fibre diameter was investigated experimentally using ten different muscle nerves in two cats.
In our work a conduction velocity of 85.35 m/s is predicted from their approximate relationship.
Simulation parameters for the 10 and 20 μm diameter fibres.
fibre diameter (μm)
axon diameter (μm)
internodal length (mm)
axoplasm resistivity (Ohmm)
We have created an anatomically based geometric model of the posterior motor neurons of the human lower limb. The geometry was mainly based on data from the Visible Man dataset and was supplemented, where necessary, using existing anatomical and physiological texts. Using this geometrical representation, action potentials were simulated using a modified representation of the bidomain equations that incorporated the CRRSS ionic conductance model. When an appropriate external stimulus was applied, an action potential was propagated along the nerve tree with a constant velocity of 89.8 m/s for a 15 μm diameter fibre. Most importantly the shape of an action potential was not altered by the simulation. The computed action potential conduction velocity agrees well with previously published theoretical and experimental estimates [10, 38, 39].
The axoplasm resistivity we used for the 15 μm diameter fibre was lower than that used in the original CRRSS model (0.24 Ohmm compared to 0.547 Ohmm ). However, in the original CRRSS paper, an intracellular stimulus was used with a 10 μm diameter fibre. It has also been noted that resistivity values are difficult to determine experimentally and agreed values are not yet well established .
The amplitude of the computed action potential (85 mV) was lower that reported in  which is about 105 – 110 mV. However, it should be noted that the amplitudes from  are obtained using experimental data from rat and cat nerves as a guide. One of the reasons for the different amplitude between the published values and the one presented in this paper could be the sodium (Na) transient since the membrane potential of the neuron is greatly affected by Na rushing into the cell when voltage gated channels open. If there is less Na available outside the cell, the inward driving force is decreased, creating a smaller depolarisation. To examine this in more detail, we compared the Na Nernst potential E Na between our model and the model of Schwarz et al. . In our case the value was 35.64 mV (which is the same as that used in the original CRRSS model) whereas the model of Schawrz et al.  used a value of 60 mV. Changing the Na Nernst potential to 55.64 mV in our model resulted in the action potential amplitude reaching approximately 100 mV together with an increased conduction velocity of 120 m/s. An additional increase of the axoplasm resistivity to 0.32 Ohmm decreased the conduction velocity to 90 m/s without any alteration to the action potential amplitude. Thus, the proposed model appears to be able to generate appropriate conduction velocities and action potential amplitudes using the axoplasm resistivity and Na Nernst potential. We note that currently good values for these parameters for the human motor nerves are not available.
One of the limitations of our model is the graded recruitment which was only mimicked by selecting a certain percentage of motor units and generating a nerve tree based on the motor entry points within these motor units. In other words, a graded recruitment was not achieved by the different strength of stimuli but by specifying the number of motor units to recruit a priori. This limitation can be overcome by modelling the sciatic nerve as a bundle of nerve fibres. With such a representation, each nerve fibre would be only connected to a specified set of motor units and the activation of each nerve fibre could then be defined according to the strength of stimuli. With our single nerve model replaced by a nerve bundle we can then investigate the relationship between surface electrode placement, stimulation strength, duration and activation of nerve fibres within the fibre bundle. The question of what external electrode configuration is required to achieve activation of a particular group of nerve fibres can then begin to be addressed. We intend to investigate this approach soon.
We have further shown that the nerve model can be connected to a skeletal muscle fibre model using the semitendinosus muscle as an illustration. The (neuron) action potentials propagated along the nerve to initiate a contractile response within the muscle. A novel approach of effectively bridging the scales of the electrophysiological cellular behaviour and the mechanical response of a skeletal muscle has been recently developed . This combination can be used as a powerful tool to investigate and answer open questions in the field of FES, e.g the use of specific stimulation protocols to achieve a desired movement or to delay the fatigue process.
This paper introduces an anatomically based geometric model of the posterior motor neurons in human lower limb based on data from the Visible Man dataset. Using this geometrical representation, action potentials were computed using a bidomain representation that incorporated the published ionic conductance model. The results indicate that this model can be used to look at the effect of external stimulation on nerve and muscle activity in the field of FES.
Funding for this work is provided by the Foundation for Research Science and Technology (FRST) under contract number C10X0402.
- Peripheral nerve and muscle stimulation [http://www.case.edu/groups/ANCL/pages/99/Ch4–2-Preprint.pdf]
- Leventhal D, Durand D: Subfascicle stimulation selectivity with the flat interface nerve electrode. Annals of Biomedical Enginering 2003, 31: 643–652. 10.1114/1.1569266View ArticleGoogle Scholar
- Hodgkin A, Huxley A: A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology 1952, 117: 500–544.View ArticleGoogle Scholar
- Frankenhaeuser B, Huxley A: The action potential in the myelinated nerve fibre of xenopus laevis as computed on the basis of voltage clamp data. The Journal of Physiology 1964, 171: 302–315.View ArticleGoogle Scholar
- Reilly J, Freeman V, Larkin W: Sensory effects of transient electrical stimulation evaluation with a neuroelectric model. IEEE Transactions on Biomedical Engineering 1985, 32: 1001–1011. 10.1109/TBME.1985.325509View ArticleGoogle Scholar
- Sweeney J, Mortimer J, Durand D: Modelling of mammalian myelinated nerve for functional neuromuscular stimulation. IEEE 9th Annual Conference of the Engineering in Medicine and Biology Society 1987, 1577–1578.Google Scholar
- Chiu S, Ritchie J, Rogart R, Stagg D: A quantitative description of membrane currents in rabbit myelinated nerve. The Journal of Physiology 1979, 292: 149–166.View ArticleGoogle Scholar
- Rattay F, Aberham M: Modeling Axon membranes for functional electrical stimulation. IEEE Transactions on Biomedical Engineering 1993, 40: 1201–1209. 10.1109/10.250575View ArticleGoogle Scholar
- Frijins J, Mooij J, ten Kate J: A quantitative approach to modeling mammalian myelinated nerve fibers for electrical prosthesis design. IEEE Transactions on Biomedical Engineering 1994, 41: 556–566. 10.1109/10.293243View ArticleGoogle Scholar
- McIntyre C, Richardson A, Grill W: Modeling the excitability of mammalian nerve fibers:influence of afterpotentials on the recovery cycle. J Neurophysiol 2001,87(2):995–1006.Google Scholar
- Halter J, Clark J: A distributed-parameter model of the myelinated nerve fiber. Journal of Theoretical Biology 1991, 148: 345–382. 10.1016/S0022-5193(05)80242-5View ArticleGoogle Scholar
- Litvak L, Delqutte B, Eddington D: Auditory nerve fiber responses to electric stimulation: modulated and unmodulated pulse trains. The Journal of the Acoustical Society of America 2001, 110: 368–379. 10.1121/1.1375140View ArticleGoogle Scholar
- Mahnam A, Firoozabadi S, Hashemi S, Janahmadi M: The use of ramp-shaped stimulus waveform for selective stimulation of neural fibers in peripheral nerves. The 2nd European Biomedical Engineering Conference EMBEC 05 2005
- Struijk J, Holsheimer J, van der Heide G, Boom H: Recruitment of dorsal column fibres in spinal cord stimulation: influence of collateral branching. IEEE Transactions on Biomedical Engineering 1992, 9: 903–912. 10.1109/10.256423View ArticleGoogle Scholar
- Colombo J, Parkins C: A model of electrical excitation of the mammalian auditory-nerve neuron. Hearing Research 1987, 31: 287–311. 10.1016/0378-5955(87)90197-3View ArticleGoogle Scholar
- Williams P, (Ed): Gray's anatomy. Churchill Livingstone; 1995.Google Scholar
- Agur A, Lee M: Grant's atlas of anatomy. Lippincott Williams & Wilkins; 1999.Google Scholar
- Spitzer V, Ackerman M, Scherzinger A, Whitlock D: The visible human male: a technical report. J of the American Med Informatics Ass 1996,3(2):118–130.View ArticleGoogle Scholar
- Eycleshymer A, Schoemaker D: A cross-section anatomy. London:Butterworths; 1970.Google Scholar
- Woodley S, Mercer S: Hamstring muscles: architecture and innervation. Cells Tissues Organs 2005, 179: 125–141. 10.1159/000085004View ArticleGoogle Scholar
- Sunderland S, Hughes E: Medical and non-medical features of the muscular branches of the sciatic nerve and its medial and lateral popliteal divisions. The Journal of Comparative Neurology 1946, 85: 205–222. 10.1002/cne.900850206View ArticleGoogle Scholar
- Blyth P: Modelling the nerves of the lower limb. Tech. rep., Department of bioengineering science, University of Auckland; 2001.Google Scholar
- Bradley C, Pullan A, Hunter P: Geometric modeling of the human torso. Annals of Biomedical Engineering 1997,27(1):96–111. 10.1007/BF02738542View ArticleGoogle Scholar
- Fernandez J, Mithraratne P, Thrupp S, Tawhai M, Hunter P: Anatomically based geometric modelling of the musculo-skeletal system and other organs. Journal Biomechanics and Modeling in Mechanobiology 2004, 2: 139–155. 10.1007/s10237-003-0036-1View ArticleGoogle Scholar
- Aquilonius S, Asknark H, Gillberg P, Nandedkar S, Olsson Y, Stalberg E: Topographial localization of motor endplates in cryosections of whole human muscles. Muscle & Nerve 1984, 7: 287–293. 10.1002/mus.880070406View ArticleGoogle Scholar
- Guyton A: Textbook of medical physiology. Philadelphia; London, Saunders; 1986.Google Scholar
- Neuron-no labels.png [http://en.wikipedia.org/wiki/Image:Neuron-no_labels.png]
- Gartner L, Hiatt J: Color atlas of histology. Lippincott Williams & Wilkins; 2006.Google Scholar
- Katz B, Miledi R: Propagation of electric activity in motor nerve terminals. Proceedings of the Royal Society of London, Series B, Biological Sciences 1965, 161: 453–482.View ArticleGoogle Scholar
- Rattay F: Analysis of models for external stimulation of axons. IEEE Transactions on Biomedical Engineering 1986, 33: 974–977. 10.1109/TBME.1986.325670View ArticleGoogle Scholar
- Hindmarsh A: ODEPACK, "A systematized collection of ODE solvers", in scientific computing. North-Holland, Amsterdam; 1982.Google Scholar
- Petzold L: Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential equations. SIAM Journal on Scientific Computing 1983, 4: 136–148. 10.1137/0904010View ArticleMATHGoogle Scholar
- CMISS [http://www.cmiss.org]
- Waxman S: Determinants of conduction velocity in myelinated nerve fibers. Muscle & Nerve 1980, 3: 141–150. 10.1002/mus.880030207View ArticleGoogle Scholar
- Reinard L, Bischhausen R: How are sheath dimensions affected by axon caliber and internode length? Brain Research 1982, 235: 335–350. 10.1016/0006-8993(82)91012-5View ArticleGoogle Scholar
- Shorten P, O'Callaghan P, Davidson J, Soboleva T: A mathematical model of fatigue in skeletal muscle force contraction. J Muscle Res Cell Motil 2007, in press.Google Scholar
- Schwarz J, Eikhof G: Na currents and action potentials in rat myelinated nerve fibers at 20 and 37° C . Pflgers Archiv European Journal of Physiology 1987, 409: 569–577. 10.1007/BF00584655View ArticleGoogle Scholar
- Rushton W: A theory of the effects of fiber size in medullated nerve. The Journal of Physiology 1951, 115: 101–122.View ArticleGoogle Scholar
- Boyd I, Kalu K: Scaling factor relating conduction velocity and diameter for myelinated afferent nerve fibers in the cat hind limb. The Journal of Physiology 1979, 289: 277–297.View ArticleGoogle Scholar
- Wesselink W, Holsheimer J, Sönmez Z, Boom H: A simulation model of the electrical characteristics of human myelinated sensory nerve fibers. Med Biol Eng Comput 1997,37(2):2028–2035.Google Scholar
- Röhrle O, Davidson J, Pullan A: Bridging scales: a three-dimensional electromechanical finite element model of skeletal muscle. SIAM Journal on Scientific Computing 2007, in press.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.