×

Reduced order models of myelinated axonal compartments. (English) Zbl 1427.92027

Summary: The paper presents a hierarchical series of computational models for myelinated axonal compartments. Three classes of models are considered, either with distributed parameters (2.5D EQS-ElectroQuasi Static, 1D TL-Transmission Lines) or with lumped parameters (0D). They are systematically analyzed with both analytical and numerical approaches, the main goal being to identify the best procedure for order reduction of each case. An appropriate error estimator is proposed in order to assess the accuracy of the models. This is the foundation of a procedure able to find the simplest reduced model having an imposed precision. The most computationally efficient model from the three geometries proved to be the analytical 1D one, which is able to have accuracy less than 0.1%. By order reduction with vector fitting, a finite model is generated with a relative difference of \(10^{- 4}\) for order 5. The dynamical models thus extracted allow an efficient simulation of neurons and, consequently, of neuronal circuits. In such situations, the linear models of the myelinated compartments coupled with the dynamical, nonlinear models of the Ranvier nodes, neuronal body (soma) and dendritic tree give global reduced models. In order to ease the simulation of large-scale neuronal systems, the sub-models at each level, including those of myelinated compartments should have the lowest possible order. The presented procedure is a first step in achieving simulations of neural systems with accuracy control.

MSC:

92C20 Neural biology
92B20 Neural networks for/in biological studies, artificial life and related topics
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Abramowitz, Milton; Stegun, Irene A.; Romer, Robert H., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, American Journal of Physics, 56, 958-958 (1988)
[2] Antoulas, AC; Sorensen, DC; Gugercin, S., A survey of model reduction methods for large-scale systems, Contemporary Mathematics, 280, 193-220 (2001) · Zbl 1048.93014
[3] Bărbulescu, R., Ioan, D., Ciuprina, G. (2018). Modeling the saltatory conduction in myelinated axons by order reduction. In 2018 20th international conference on neuroinformatics and computational neuroscience (ICNCN) (pp. 2059-2062).
[4] Bărbulescu, R., Ioan, D., Ciuprina, G. (2019). Coupled macromodels for the simulation of the saltatory conduction. Scientific Bulletin Series C (under review).
[5] Bărbulescu, R., Ioan, D., Ciurea, J. (2016). Simple 1d models for neuro-signals transmission along axons. In 2016 international conference and exposition on electrical and power engineering (EPE) (pp. 313-319): IEEE.
[6] Berljafa, M., & Güttel, S. (2014). A rational Krylov toolbox for matlab. · Zbl 1319.65028
[7] Bower, J., & Beeman, D. (1998). The Book of GENESIS: exploring realistic neural models with the GEneral NEural SImulation System. Interpretation of Music; 55. TELOS. https://books.google.ro/books?id=WapQAAAAMAAJ. Accessed 14 Dec 2017. · Zbl 0820.68136
[8] Brain Facts and Figures. (2017). https://faculty.washington.edu/chudler/facts.html#brain. Accessed 14 Dec 2017.
[9] Brown, A.M., & Hamann, M. (2014). Computational modeling of the effects of auditory nerve dysmyelination. Frontiers in Neuroanatomy 8.
[10] Burger, J.R. (2009). Human memory modeled with standard analog and digital circuits: inspiration for man-made computers. New Jersey: Wiley.
[11] Carnevale, N.T., & Hines, M.L. (2006). The NEURON book. Cambridge: Cambridge University Press.
[12] Ciuprina, G., Villena, J.F., Ioan, D., Ilievski, Z., Kula, S., ter Maten, E.J.W., Mohaghegh, K., Pulch, R., Schilders, W.H., Silveira, L.M., et al. (2015). Parameterized model order reduction. In Coupled multiscale simulation and optimization in nanoelectronics (pp. 267-359): Springer.
[13] Clark, J.; Plonsey, R., A mathematical evaluation of the core conductor model, Biophysical Journal, 6, 95-112 (1966)
[14] De Sterck, H., & Ullrich, P. (2009). Introduction to computational pdes. Course Notes for Amath 442.
[15] Digicortex. (2017). http://www.digicortex.net/. Accessed: 2017-12-14.
[16] Elmore, WC, The transient response of damped linear networks with particular regard to wideband amplifiers, Journal of Applied Physics, 19, 55-63 (1948)
[17] Fitzhugh, R., Computation of impulse initiation and saltatory conduction in a myelinated nerve fiber, Biophysical Journal, 2, 11-21 (1962)
[18] FitzHugh, R. (1966). Mathematical models of excitation and propagation in nerve. Publisher Unknown.
[19] 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, 171, 302-315 (1964)
[20] Ganapathy, N.; Clark, J., Extracellular currents and potentials of the active myelinated nerve fiber, Biophysical Journal, 52, 749-761 (1987)
[21] Goldman, L.; Albus, JS, Computation of impulse conduction in myelinated fibers; theoretical basis of the velocity-diameter relation, Biophysical Journal, 8, 596-607 (1968)
[22] Gow, A.; Devaux, J., A model of tight junction function in central nervous system myelinated axons, Neuron Glia Biology, 4, 307-317 (2008)
[23] Graham, B., Gillies, A., Willshaw, D. (2011). Principles of computational modelling in neuroscience.
[24] Graupner, M., & Brunel, N. (2010). Mechanisms of induction and maintenance of spike-timing dependent plasticity in biophysical synapse models. Frontiers in Computational Neuroscience 4.
[25] Gurney, K. (1997). An introduction to neural networks. Boca Raton: CRC Press.
[26] Gustavsen, B.; Semlyen, A., Rational approximation of frequency domain responses by vector fitting, IEEE Transactions on Power Delivery, 14, 1052-1061 (1999)
[27] Halter, JA; Clark, JWJr, A distributed-parameter model of the myelinated nerve fiber, Journal of Theoretical Biology, 148, 345-382 (1991)
[28] Haykin, S. (1999). Neural networks: a comprehensive foundation. International edition. Prentice Hall. https://books.google.ro/books?id=M5abQgAACAAJ. · Zbl 0934.68076
[29] Heres, P.P. (2005). Robust and efficient Krylov subspace methods for model order reduction. Ph.D. thesis, Technische Universiteit Eindhoven.
[30] Hodgkin, AL; Huxley, AF, A quantitative description of membrane current and its application to conduction and excitation in nerve, The Journal of Physiology, 117, 500-544 (1952)
[31] Huxley, A.; Stämpeli, R., Evidence for saltatory conduction in peripheral myelinated nerve fibres, The Journal of Physiology, 108, 315-339 (1949)
[32] Ioan, D. (1988). Metode pentru calculul campului electromagnetic. Separarea variabilelor. IPB.
[33] Ioan, D.; Ciuprina, G.; Radulescu, M.; Seebacher, E., Compact modeling and fast simulation of on-chip interconnect lines, IEEE Transactions on Magnetics, 42, 547-550 (2006)
[34] Ioan, D.; Munteanu, I., Missing link rediscovered: the electromagnetic circuit element concept, JSAEM Studies in Applied Electromagnetics and Mechanics, 8, 302-320 (1999)
[35] Izhikevich, EM, Simple model of spiking neurons, IEEE Transactions on Neural Networks, 14, 1569-1572 (2003)
[36] Izhikevich, E.M. (2007). Dynamical systems in neuroscience. Cambridge: MIT Press.
[37] Jin, J.M. (2015). The finite element method in electromagnetics. New Jersey: Wiley.
[38] Joucla, S., Glière, A., Yvert, B. (2014). Current approaches to model extracellular electrical neural microstimulation. Frontiers in Computational Neuroscience 8.
[39] Keener, J., & Sneyd, J. (2010). Mathematical physiology: I: cellular physiology. Springer Science & Business Media. · Zbl 1273.92017
[40] Kuokkanen, P. (2012). On the origin of the extracellular potential in the nucleus laminaris of the barn owl. Ph.D. thesis, Humboldt-Universität zu Berlin, Mathematisch-Naturwissenschaftliche Fakultät I.
[41] Lindsay, K., Ogden, J., Halliday, D., Rosenberg, J. (1999). An introduction to the principles of neuronal modelling. In Modern techniques in neuroscience research (pp. 213-306): Springer.
[42] Lozier, DW, Nist digital library of mathematical functions, Annals of Mathematics and Artificial Intelligence, 38, 105-119 (2003) · Zbl 1019.65001
[43] Maass, W., Networks of spiking neurons: the third generation of neural network models, Neural Networks, 10, 1659-1671 (1997)
[44] MATLAB. (2015). Balanced model truncation via square root method. The MathWorks Inc., Natick, Massachusetts.
[45] MATLAB. (2015). Singular value decomposition. The MathWorks Inc., Natick, Massachusetts.
[46] Moore, J.; Ramon, F.; Joyner, R., Axon voltage-clamp simulations. i. methods and tests, Biophysical Journal, 15, 11-24 (1975)
[47] Moore, JW; Joyner, RW; Brill, MH; Waxman, SD; Najar-Joa, M., Simulations of conduction in uniform myelinated fibers. relative sensitivity to changes in nodal and internodal parameters, Biophysical Journal, 21, 147-160 (1978)
[48] Morrison, A.; Diesmann, M.; Gerstner, W., Phenomenological models of synaptic plasticity based on spike timing, Biological Cybernetics, 98, 459-478 (2008) · Zbl 1145.92306
[49] Niebur, E., Neuronal cable theory, Scholarpedia, 3, 2674 (2008)
[50] Panzer, H. (2014). Model order reduction by Krylov subspace methods with global error bounds and automatic choice of parameters. Verlag Dr. Hut.
[51] Parasuram, H., Nair, B., D’Angelo, E., Hines, M., Naldi, G., Diwakar, S. (2016). Computational modeling of single neuron extracellular electric potentials and network local field potentials using lfpsim. Frontiers in Computational Neuroscience 10.
[52] Paugam-Moisy, H., & Bohte, S. (2012). Computing with spiking neuron networks. In Handbook of natural computing (pp. 335-376): Springer.
[53] Rapetti, F.; Rousseaux, G., On quasi-static models hidden in Maxwell’s equations, Applied Numerical Mathematics, 79, 92-106 (2014) · Zbl 1291.78019
[54] Rattay, F.; Potrusil, T.; Wenger, C.; Wise, AK; Glueckert, R.; Schrott-Fischer, A., Impact of morphometry, myelinization and synaptic current strength on spike conduction in human and cat spiral ganglion neurons, PLoS One, 8, e79,256 (2013)
[55] Robinson, P.; Rennie, C.; Rowe, D.; O’Connor, S.; Gordon, E., Multiscale brain modelling, Philosophical Transactions of the Royal Society of London B: Biological Sciences, 360, 1043-1050 (2005)
[56] Salimbahrami, B., & Lohmann, B. (2002). Krylov subspace methods in linear model order reduction: introduction and invariance properties. In Sci. Rep. Inst. of Automation. Citeseer. · Zbl 1354.65076
[57] Simulating action potential with the Hodgkin-Huxley model. (2018). https://www.comsol.com/model/simulating-action-potential-with-the-hodgkin-huxley-model-47121https://www.comsol.com/model/simulating-action-potential-with-the-hodgkin-huxley-model-47121. Accessed: 2018-06-10.
[58] Stephanova, D., Myelin as longitudinal conductor: a multi-layered model of the myelinated human motor nerve fibre, Biological Cybernetics, 84, 301-308 (2001)
[59] Stoica, P., Moses, R.L., et al. (2005). Spectral analysis of signals Vol. 452. Upper Saddle River: Pearson Prentice Hall.
[60] Struijk, JJ; Holsheimer, J.; Heide, GG; Boom, HB, Recruitment of dorsal column fibers in spinal cord stimulation: influence of collateral branching, IEEE Transactions on Biomedical Engineering, 39, 903-912 (1992)
[61] Timotin, A. (2004). La structure de la fibre nerveuse: Un projet optimal. Proceedings of the Romanian Academy Series A 5(1).
[62] The Nervous System. (2017). (Structure and function) (Nursing) part 1. http://what-when-how.com/nursing/the-nervous-system-structure-and-function-nursing-part-1/http://what-when-how.com/nursing/the-nervous-system-structure-and-function-nursing-part-1/ . Accessed: 2017-12-14.
[63] Geit, W.; Schutter, E.; Achard, P., Automated neuron model optimization techniques: a review, Biological Cybernetics, 99, 241-251 (2008) · Zbl 1154.92013
[64] Villapecellín-Cid, MM; Medina, F.; Roa, LM, Internodal myelinated segments: delay and rgc time-domain green function model, IEEE Transactions on Biomedical Engineering, 51, 389-391 (2004)
[65] Young, R. (2015). Mathematical modeling of the evolution and development of myelin. Ph.D. thesis, University of Hawaii at Manoa.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.