Abstract:A partial differential equation model for the three-dimensional current flow in an excitable, unmyelinated axon is considered. Where the axon radius is significantly below a critical value R(crit) (that depends upon intra- and extra-cellular conductivity and ion channel conductance) the resistance of the intracellular space is significantly higher than that of the extracellular space, such that the potential outside the axon is uniformly small whilst the intracellular potential is approximated by the transmembrane potential. In turn, since the current flow is predominantly axial, it can be shown that the transmembrane potential is approximated by a solution to the one-dimensional cable equation. It is noted that the radius of the squid giant axon, investigated by (Hodgkin and Huxley 1952e), lies close to R(crit). This motivates us to apply the three-dimensional model to the squid giant axon and compare the results thus found to those obtained using the cable equation. In the context of the in vitro experiments conducted in (Hodgkin and Huxley 1952e) we find only a small difference between the wave profiles determined using these two different approaches and little difference between the speeds of action potential propagation predicted. This suggests that the cable equation approximation is accurate in this scenario. However when applied to the it in vivo setting, in which the conductivity of the surrounding tissue is considerably lower than that of the axoplasm, there are marked differences in both wave profile and speed of action potential propagation calculated using the two approaches. In particular, the cable equation significantly over predicts the increase in the velocity of propagation as axon radius increases. The consequences of these results are discussed in terms of the evolutionary costs associated with increasing the speed of action potential propagation by increasing axon radius.

Solving Hodgkin-Huxley equations using the compact difference scheme -tapering dendrite

Solving the Cable Equation Using a Compact Difference Scheme -- Passive Soma Dendrite

Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

Efficient Numerical Methods for Solving the Membrane Equations of Human Ventricular Cell Model

A Computational Approach for the inverse problem of neuronal conductances determination

Simulation of Three-Dimensional Dendritic Growth Using H-adaptive Finite Element Method

Electro-diffusive modeling and the role of spine geometry on action potential propagation in neurons

An Implicit Parallel Multigrid Computing Scheme to Solve Coupled Thermal-Solute Phase-Field Equations for Dendrite Evolution

Modelling in vivo action potential propagation along a giant axon

Well-posedness and regularity of solutions to neural field problems with dendritic processing

A novel central compact finite-difference scheme for third derivatives with high spectral resolution

Dendrites and conformal symmetry

Fractional Cable Model for Signal Conduction in Spiny Neuronal Dendrites

Dirac assisted tree method for 1D heterogeneous Helmholtz equations with arbitrary variable wave numbers

Computational Modeling of Three-Dimensional Electrodiffusion in Biological Systems: Application to the Node of Ranvier

Convergence analysis of exponential time differencing scheme for the nonlocal Cahn-Hilliard equation

Spatio-temporal modeling of saltatory conduction in neurons using Poisson-Nernst-Planck treatment and estimation of conduction velocity

Electrodiffusion Models of Axon and Extracellular Space Using the Poisson-Nernst-Planck Equations

Dynamics from seconds to hours in Hodgkin-Huxley model with time-dependent ion concentrations and buffer reservoirs

Dendritic growth at very low undercoolings