A distributed-parameter model of the myelinated nerve fiber

https://doi.org/10.1016/S0022-5193(05)80242-5Get rights and content

This paper presents a new model for the characterization of electrical activity in the nodal, paranodal and internodal regions of isolated amphibian and mammalian myelinated nerve fibers. It differs from previous models in the following ways: (1) in its ability to incorporate detailed anatomical and electrophysiological data; (2) in its approach to the myelinated nerve fiber as a multi-axial cable; and (3) in the numerical algorithm used to obtain distributed model equation solutions for potential and current. The morphometric properties are taken from detailed electron microscopic anatomical studies (Berthold & Rydmark, 1983a, Experientia 39, 964–976). The internodal axolemma is characterized as an excitable membrane and model-generated nodal and internodal membrane action potentials are presented. A system of describing equations for the equivalent network model is derived, based on the application of Kirchoff's Current Law, which take the form of multiple cross-coupled parabolic partial differential equations. An implicit numerical integration method is developed and the numerical solution implemented on a parallel processor. Non-uniform spatial step sizes are used, enabling detailed representation of the nodal region while minimizing the number of total segments necessary to represent the overall fiber. Conduction velocities of 20·2 m sec−1 at 20°C for a 15 μm diameter amphibian fiber and 57·6 m sec−1 at 37°C for a 17·5 μm diameter mammalian fiber are achieved, which agrees qualitatively with published experimental data at similar temperatures (Huxley & Stämpfli, 1949, J. Physiol., Lond. 108, 315–339; Rasminsky, 1973, Arch, Neurol. 28, 287–292). The simulation results demonstrate the ability of this model to produce detailed representations of the transaxonal, transmyelin and transfiber potentials and currents, as well as the longitudinal extra-axonal, periaxonal and intra-axonal currents. Also indicated is the potential contribution of the paranodal axolemma to nodal activity as well as the presence of significant longitudinal currents in the periaxonal space adjacent to the node of Ranvier.

References (67)

  • BakerM. et al.

    Accommodation in rat myelinated axons depends on two pharmacologically distinct types of potassium channels

    J. Physiol., Lond.

    (1985)
  • BarrettE.F. et al.

    Intracellular recording from vertebrate myelinated axons: mechanism of the depolarizing afterpotential

    J. Physiol., Lond.

    (1982)
  • BertholdC.-H. et al.

    Axon diameter and myelin sheath thickness in nerve fibers of the ventral spinal root of the seventh lumbar nerve of the adult and developing cat

    J. Anat.

    (1983)
  • BertholdC.-H. et al.

    Anatomy of the paranode-node-paranode in the cat

    Experientia

    (1983)
  • BertholdC.-H. et al.

    Electron microscopic serial section analysis of node of Ranvier in lumbosacral spinal roots of the cat: ultrastructural organization of nodal compartments in fibres of different sizes

    J. Neurocytol.

    (1983)
  • BostockH. et al.

    The internodal axon membrane: electrical excitability and continuous conduction in segmental demyelination

    J. Physiol., Lond.

    (1978)
  • BostockH. et al.

    The effect of 4-aminopyridine and tetraethylammonium ions on normal and demyelinated mammalian nerve fibres

    J. Physiol., Lond.

    (1981)
  • BrillM.H. et al.

    Conduction velocity and spike configuration in myelinated fibres: computed dependence on internode distance

    J. Neurol. Neurosurg. Psychiat.

    (1977)
  • BrismarT.

    Electrical properties of isolated demyelinated rat nerve fibres

    Acta Physiol. Scand.

    (1981)
  • BrismarT.

    Potential clamp of demyelinated rat nerve fibres

    Acta Physiol. Scand.

    (1981)
  • BrismarT. et al.

    Potassium permeability in rat myelinated nerve fibres

    Acta Physiol. Scand.

    (1985)
  • ChiuS.Y.

    Asymmetry currents in the mammalian myelinated nerve fiber

    J. Physiol., Lond.

    (1980)
  • ChiuS.Y.

    Sodium currents in axon-associated Schwann cells from adult rabbits

    J. Physiol., Lond.

    (1987)
  • ChiuS.Y. et al.

    Evidence for the presence of potassium channels in the paranodal region of acutely demyelinated mammalian single nerve fibres

    J. Physiol., Lond.

    (1981)
  • ChiuS.Y. et al.

    Evidence for the presence of potassium channels in the internode of frog myelinated nerve fibers

    J. Physiol., Lond.

    (1982)
  • ChiuS.Y. et al.

    A quantitative description of membrane currents in rabbit myelinated nerve

    J. Physiol., Lond.

    (1979)
  • CrankJ. et al.

    A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type

  • DodgeF.A.

    A study of ionic permeability changes underlying excitation in myelinated nerve fibers of the frog

    Ph.D. thesis. The Rockerfeller Institute, New York

    (1963)
  • DongarraJ.J. et al.

    LINPACK User's Guide

    (1979)
  • DuboisJ.M.

    Evidence for the existence of three types of potassium channels in the frog Ranvier node membrane

    J. Physiol., Lond.

    (1981)
  • DuboisJ.M.

    Simultaneous changes in the equilibrium potential and potassium conductance in voltage clamped Ranvier node in the frog

    J. Physiol., Lond.

    (1981)
  • FrankenhaeuserB. et al.

    The action potentials in the myelinated nerve fibre of Xenopus laevis as computed on the basis of voltage clamp data

    J. Physiol., Lond.

    (1964)
  • FrankenhaeuserB. et al.

    The effect of temperature on the sodium and potassium permeability changes in myelinated fibres of Xenopus laevis

    J. Physiol., Lond.

    (1963)
  • Cited by (89)

    • Saltatory Conduction along Myelinated Axons Involves a Periaxonal Nanocircuit

      2020, Cell
      Citation Excerpt :

      The emergent DC properties were consistent with the low capacitance produced by ∼15 myelin lamellae, a conductive nanoscale periaxonal space (∼50 Ω cm) and electrical sealing of the paranodes (∼550 Ω cm; Figures 2 and 3). These findings are in agreement with earlier predicted incomplete insulation by the myelin sheath on the basis of experimental recordings (Barrett and Barrett, 1982; Blight and Someya, 1985; Funch and Faber, 1984) and provide a quantitative framework for axoplasmic, periaxonal, and paranodal resistivity, substantially advancing our understanding of the parameters to model submyelin conduction (Arancibia-Carcámo et al., 2017; Blight, 1985; Dimitrov, 2005; Gow and Devaux, 2008; Halter and Clark, 1991; McIntyre et al., 2002; Richardson et al., 2000; Stephanova and Bostock, 1995; Young et al., 2013). The spatial and temporal evolution of the transaxonal potentials beneath myelin is strikingly complex with gradually attenuating waves toward the middle of the internode (Figures 5 and 6).

    • Modeling conduction delays in the corpus callosum using MRI-measured g-ratio

      2019, NeuroImage
      Citation Excerpt :

      Due to the myelin sheath, the voltage progression in the internode can be described as passive. The progression of voltage along both segments of the axons have already been modelled and numerically simulated (e.g. Rushton, 1951; Moore et al., 1978; Haltert and Clark, 1991; Basser, 1993). This understanding of the axon signal conduction, together with knowledge on the axon's electrical properties (Debanne et al., 2011; Bakiri et al., 2011) and the ability to measure its geometry, will allow to implement a model that links white matter structure with function in biophysical terms.

    • Evolution of rapid nerve conduction

      2016, Brain Research
      Citation Excerpt :

      This was confirmed in another study, which found that the optimal internode length (i.e. the length that corresponds to the maximum conduction velocity) depended on the thickness of the myelin sheath, with the optimal length shifting to longer internodes as the thickness increased (and capacitance dropped) (Castelfranco and Hartline, 2015). Another condition that affects the theoretical results is whether the models describe the myelin sheath as distinct from the axolemma, including current flow in the submyelin space as in “double cable” models (Blight, 1985; Halter and Clark, 1991; Hines and Shrager, 1991; McIntyre et al., 2002). Overall these detailed models show a relationship between conduction velocity and internode length similar to that described by Brill et al. (1977), but the predictions of the effect on conduction velocity of a specific change in internode length can differ (Lasiene et al., 2008; Young et al., 2013a; Young et al., 2013b).

    • Computer Models of Peripheral Nerves

      2015, Nerves and Nerve Injuries: Pain, Treatment, Injury, Disease and Future Directions: Vol 2
    View all citing articles on Scopus
    View full text