Analytical solutions of the Frankenhaeuser-Huxley equations I: minimal model for backpropagation of action potentials in sparsely excitable dendrites

Exact solutions to cable equations in branching neurons with tapering dendrites

Neurons are biological cells with uniquely complex dendritic morphologies that are not present in other cell types. Electrical signals in a neuron with branching dendrites can be studied by cable theory which provides a general mathematical modelling framework of spatio-temporal voltage dynamics. Typically such models need to be solved numerically unless the cell membrane is modelled either by passive or quasi-active dynamics, in which cases analytical solutions can be reduced to calculation of the Green's function describing the fundamental input-output relationship in a given morphology. Such analytically tractable models often assume individual dendritic segments to be cylinders. However, it is known that dendritic segments in many types of neurons taper, i.e. their radii decline from proximal to distal ends. Here we consider a generalised form of cable theory which takes into account both branching and tapering structures of dendritic trees. We demonstrate that analytical solutions can be found in compact algebraic forms in an arbitrary branching neuron with a class of tapering dendrites studied earlier in the context of single neuronal cables by Poznanski (Bull. Math. Biol. 53(3):457-467, 1991). We apply this extended framework to a number of simplified neuronal models and contrast their output dynamics in the presence of tapering versus cylindrical segments.

Action potential propagation and synchronisation in myelinated axons

With the advent of advanced MRI techniques it has become possible to study axonal white matter non-invasively and in great detail. Measuring the various parameters of the long-range connections of the brain opens up the possibility to build and refine detailed models of large-scale neuronal activity. One particular challenge is to find a mathematical description of action potential propagation that is sufficiently simple, yet still biologically plausible to model signal transmission across entire axonal fibre bundles. We develop a mathematical framework in which we replace the Hodgkin-Huxley dynamics by a spike-diffuse-spike model with passive sub-threshold dynamics and explicit, threshold-activated ion channel currents. This allows us to study in detail the influence of the various model parameters on the action potential velocity and on the entrainment of action potentials between ephaptically coupled fibres without having to recur to numerical simulations. Specifically, we recover known results regarding the influence of axon diameter, node of Ranvier length and internode length on the velocity of action potentials. Additionally, we find that the velocity depends more strongly on the thickness of the myelin sheath than was suggested by previous theoretical studies. We further explain the slowing down and synchronisation of action potentials in ephaptically coupled fibres by their dynamic interaction. In summary, this study presents a solution to incorporate detailed axonal parameters into a whole-brain modelling framework.

With more and more data becoming available on white-matter tracts, the need arises to develop modelling frameworks that incorporate these data at the whole-brain level. This requires the development of efficient mathematical schemes to study parameter dependencies that can then be matched with data, in particular the speed of action potentials that cause delays between brain regions. Here, we develop a method that describes the formation of action potentials by threshold activated currents, often referred to as spike-diffuse-spike modelling. A particular focus of our study is the dependence of the speed of action potentials on structural parameters. We find that the diameter of axons and the thickness of the myelin sheath have a strong influence on the speed, whereas the length of myelinated segments and node of Ranvier length have a lesser effect. In addition to examining single axons, we demonstrate that action potentials between nearby axons can synchronise and slow down their propagation speed.

Induced mitochondrial membrane potential for modeling solitonic conduction of electrotonic signals

A cable model that includes polarization-induced capacitive current is derived for modeling the solitonic conduction of electrotonic potentials in neuronal branchlets with microstructure containing endoplasmic membranes. A solution of the nonlinear cable equation modified for fissured intracellular medium with a source term representing charge ‘soakage’ is used to show how intracellular capacitive effects of bound electrical charges within mitochondrial membranes can influence electrotonic signals expressed as solitary waves. The elastic collision resulting from a head-on collision of two solitary waves results in localized and non-dispersing electrical solitons created by the nonlinearity of the source term. It has been shown that solitons in neurons with mitochondrial membrane and quasi-electrostatic interactions of charges held by the microstructure (i.e., charge ‘soakage’) have a slower velocity of propagation compared with solitons in neurons with microstructure, but without endoplasmic membranes. When the equilibrium potential is a small deviation from rest, the nonohmic conductance acts as a leaky channel and the solitons are small compared when the equilibrium potential is large and the outer mitochondrial membrane acts as an amplifier, boosting the amplitude of the endogenously generated solitons. These findings demonstrate a functional role of quasi-electrostatic interactions of bound electrical charges held by microstructure for sustaining solitons with robust self-regulation in their amplitude through changes in the mitochondrial membrane equilibrium potential. The implication of our results indicate that a phenomenological description of ionic current can be successfully modeled with displacement current in Maxwell’s equations as a conduction process involving quasi-electrostatic interactions without the inclusion of diffusive current. This is the first study in which solitonic conduction of electrotonic potentials are generated by polarization-induced capacitive current in microstructure and nonohmic mitochondrial membrane current.

A nonlinear cable framework for bidirectional synaptic plasticity

Finding the rules underlying how axons of cortical neurons form neural circuits and modify their corresponding synaptic strength is the still subject of intense research. Experiments have shown that internal calcium concentration, and both the precise timing and temporal order of pre and postsynaptic action potentials, are important constituents governing whether the strength of a synapse located on the dendrite is increased or decreased. In particular, previous investigations focusing on spike timing-dependent plasticity (STDP) have typically observed an asymmetric temporal window governing changes in synaptic efficacy. Such a temporal window emphasizes that if a presynaptic spike, arriving at the synaptic terminal, precedes the generation of a postsynaptic action potential, then the synapse is potentiated; however if the temporal order is reversed, then depression occurs. Furthermore, recent experimental studies have now demonstrated that the temporal window also depends on the dendritic location of the synapse. Specifically, it was shown that in distal regions of the apical dendrite, the magnitude of potentiation was smaller and the window for depression was broader, when compared to observations from the proximal region of the dendrite. To date, the underlying mechanism(s) for such a distance-dependent effect is (are) currently unknown. Here, using the ionic cable theory framework in conjunction with the standard calcium based plasticity model, we show for the first time that such distance-dependent inhomogeneities in the temporal learning window for STDP can be largely explained by both the spatial and active properties of the dendrite.

Thermal noise due to surface-charge effects within the Debye layer of endogenous structures in dendrites

An assumption commonly used in cable theory is revised by taking into account electrical amplification due to intracellular capacitive effects in passive dendritic cables. A generalized cable equation for a cylindrical volume representation of a dendritic segment is derived from Maxwell's equations under assumptions: (i) the electric-field polarization is restricted longitudinally along the cable length; (ii) extracellular isopotentiality; (iii) quasielectrostatic conditions; and (iv) homogeneous medium with constant conductivity and permittivity. The generalized cable equation is identical to Barenblatt's equation arising in the theory of infiltration in fissured strata with a known analytical solution expressed in terms of a definite integral involving a modified Bessel function and the solution to a linear one-dimensional classical cable equation. Its solution is used to determine the impact of thermal noise on voltage attenuation with distance at any particular time. A regular perturbation expansion for the membrane potential about the linear one-dimensional classical cable equation solution is derived in terms of a Green's function in order to describe the dynamics of free charge within the Debye layer of endogenous structures in passive dendritic cables. The asymptotic value of the first perturbative term is explicitly evaluated for small values of time to predict how the slowly fluctuating (in submillisecond range) electric field attributed to intracellular capacitive effects alters the amplitude of the membrane potential. It was found that capacitive effects are almost negligible for cables with electrotonic lengths L>0.5 , contributes up to 10% of the signal for cables with electrotonic lengths in the range between 0.25<L<0.5 , and dominates the membrane potential for electrotonically short cables (L<0.2) . These results show that electrotonically short dendritic cables with both ends sealed are prone to significant neurobiological thermal noise due to intracellular capacitive effects. The presence of significant thermal noise weakens the assumption of intracellular isopotentiality when approximating dendrites with compartments.

Analytical solutions for nonlinear cable equations with calcium dynamics. II. Saltatory transmission in a sparsely excitable cable model

In order to gain a better theoretical understanding of the interaction between voltage and calcium influx, we present the simulation results for saltatory transmission in a sparsely excitable model of a continuous cylindrical segment of nerve fiber, where calcium diffuses internally and various ion channels are distributed as hotspots along the cable. A standard set of ion channel descriptions is used to illustrate how different numbers and distributions of ion channel hotspots affect the propagation and transmission of a single action potential and/or a spike train and how such hotspots affect calcium influx and diffusion within continuous cylindrical segment of nerve fiber.

Analytical solutions for nonlinear cable equations with calcium dynamics. I: Derivations

The interaction between membrane potential and internal calcium concentration plays many important roles in regulating synaptic integration and neuronal firing. In order to gain a better theoretical understanding between the voltage-calcium interaction, a nonlinear cable equation with calcium dynamics is solved analytically. This general reaction-diffusion system represents a model of a cylindrical dendritic segment in which calcium diffuses internally in the presence of buffers, pumps and exchangers, and where ion channels are sparsely distributed over the membrane,in the form of hotspots, acting as point current sources along the dendritic membrane. In order to proceed, the reaction-diffusion system is recast into a system of coupled nonlinear integral equations, with which a perturbative expansion in dimensionless voltage and calcium concentration are used to find analytical solutions to this general system. The resulting solutions can be used to investigate, the interaction between the membrane potential and underlying calcium dynamics in a natural (non-discretized) setting.

fMRI MODELS OF DENDRITIC AND ASTROCYTIC NETWORKS

In order to elucidate the relationships between hierarchical structures within the neocortical neuropil and the information carried by an ensemble of neurons encompassing a single voxel, it is essential to predict through volume conductor modeling LFPs representing average extracellular potentials, which are expressed in terms of interstitial potentials of individual cells in networks of gap-junctionally connected astrocytes and synaptically connected neurons. These relationships have been provided and can then be used to investigate how the underlying neuronal population activity can be inferred from the measurement of the BOLD signal through electrovascular coupling mechanisms across the blood-brain barrier. The importance of both synaptic and extrasynaptic transmission as the basis of electrophysiological indices triggering vascular responses between dendritic and astrocytic networks, and sequential configurations of firing patterns in composite neural networks is emphasized. The purpose of this review is to show how fMRI data may be used to draw conclusions about the information transmitted by individual neurons in populations generating the BOLD signal.

Analytical solution of the cable equation with synaptic reversal potentials for passive neurons with tip-to-tip dendrodendritic coupling

A passive cable model is presented for a pair of electrotonically coupled neurons in order to investigate the effects of tip-to-tip dendrodendritic gap junctions on the interaction between excitation and either pre or postsynaptic inhibition. The model represents each dendritic tree by a tapered equivalent cylinder attached to an isopotential soma. Analytical solution of the cable equation with synaptic reversal potentials is considered for each neuron to yield a system of Volterra integral equations for the voltage. The solution to the system of linear integral equations (expressed as a Neumann series) is used to determine the current spread within the two coupled neurons, and to re-examine the sensitivity of the soma potentials (in particular) to the coupling resistance for various loci of synaptic inputs. The model is actually posed generally, so that active as well as passive properties could be considered. In the active case, a system of non-linear integral equations is derived for the voltage.

Rallian "equivalent" cylinders reconsidered: comparisons with literal compartments

In Rall's "equivalent" cylinder morphological-to-electrical transformation, neuronal arborizations are reduced to single unbranched core-conductors. The conventional assumption that such an "equivalent" reconstructs the electrical properties of the fibers it represents was tested directly; electrical properties and responses of "equivalent" cylinders were compared with those of their literal branch constituents for fibers with a single symmetrical bifurcation. The numerical solution methods were validated independently by their accurate reconstruction of the responses of an analog circuit configured with compartmental architecture to solve the cable equation for passive fibers with a symmetrical bifurcation. In passive fibers, "equivalent" cylinders misestimated the spatial distribution of voltage amplitudes and steady-state input resistance, partly due to the lack of axial current bifurcation. In active fibers with a single propagating action potential, the spatial distributions of point-to-point conduction velocity values (measured in meters/second) for a literal branch point differed significantly from those of their "equivalent" cylinders. "Equivalent" cylinders also underestimated the diameter-dependent delay in propagation through the branch point and branches, due to the larger "equivalent" diameter. Corrections to the "equivalent" cylinder did not reconcile differences between "equivalent" and literal models. However, "equivalent" and literal branch fibers had the same (a) steady-state resistance "looking into" an isolated symmetrical branch point and (b) geometry-independent point-to-point propagation velocity when measured in space constants per millisecond except within +/-1 space constant from the geometrical inhomogeneity. In summary, Rall's "equivalent" cylinders did not accurately reconstruct all passive or active electrophysiological properties and responses of their literal compartments. For the modeling of individual neurons, the requirement of single-branch resolution is discussed.