3 Dendrites and Synapses

3.3 Spatial Structure: Axons

Any given neuron has a single axon that leaves the soma to make synaptic contacts. Like dendrites, axons have a range of different morphologies. Some axons project mainly to neurons close by. This is the case for neurons in the layer 2-3 of cortex, their axons branch out in all directions from the soma forming a star-shaped axonal arbor called a ‘daisy’. Other neurons such as pyramidal neurons situated deeper in the cortex have axons that plunge in the white matter and may cross the whole brain to reach another brain area. There are even longer axons that leave the central nervous system and travel down the spinal cord to reach muscles at the tip of the foot.

In terms of propagation dynamics, we distinguish two types of axons: the myelinated and the unmyelinated axons. We will see that myelin is useful to increase propagation speed in far-reaching projections. This is the case for cortical projections passing through the white matter, or for axons crossing the spinal cord. Short projections on the other hand use axons devoid of myelin.

3.3.1 Unmyelinated Axons

Mathematical description of the membrane potential in the axon is identical to that of dendrites with active ion channels. Unmyelinated axons contain sodium and potassium channels uniformly distributed over their entire length. The classical example is the squid giant axon investigated by Hodgkin and Huxley. The Hodgkin-Huxley model described in Chapter 2 was developed for a small axon segment. The general equation for ion channels imbedded on a passive membrane is

crLtu(t,x)=2x2u(t,x)-rL(u(t,x)-El)-rLiion[u](t,x).c\,r_{L}\frac{\partial}{\partial t}\,u(t,x)=\frac{\partial^{2}}{\partial x^{2}% }\,u(t,x)-r_{L}(u(t,x)-E_{l})-r_{L}i_{\text{ion}}[u](t,x)\,. (3.37)

where we have reverted to a variable uu in units of mV from the equation of active dendrites seen in Sect. 3.2.3. For the giant squid axon, the ionic current are described by the Hodgkin-Huxley model

iion[u](t,x)=gNam3(t,x)h(t,x)(u(t,x)-ENa)+gKn4(t,x)(u(t,x)-EK).i_{\text{ion}}[u](t,x)=g_{\rm Na}\,m^{3}(t,x)h(t,x)\,(u(t,x)-E_{\rm Na})+g_{% \rm K}\,n^{4}(t,x)\,(u(t,x)-E_{\rm K}). (3.38)

In other systems, the axon may be covered with other types of sodium or potassium ion channels.

When an action potential is fired in the axon initial segment, the elevated membrane potential will depolarize the adjacent axonal segments. Sodium channels farther down the axon which were previously closed will start to open, thereby depolarizing the membrane further. The action potential propagates by activating sodium channels along the cable rather than by spreading the charges such as in a passive dendrite. The properties of the ion channels strongly influence conduction velocity. In the unmyelinated axons of the hippocampus, the conduction velocity axons is 0.25 m/s.

The dynamics described by Eqs. (3.37)-(3.38) reproduce many properties of real axons. In particular, two spikes traveling in opposite direction will collide and annihilate each other. This is unlike wave propagating on water. Another property is reflection at branch points. When the impedance mismatch at the point where a single axon splits into two can is significant, the action potential can reflect and start travelling in the direction it came from.

The solution of Eq. (3.37) with sodium and postassium ion channels such as in Eq. (3.38) can not be written in a closed form. Properties of axonal propagation are either studied numerically (see Sect. 3.4) or with reduced models of ion channels.

Example: Speed of Propagation with Simplified Action Potential Dynamics

For the sake of studying propagation properties, we can replace the active properties of a small axonal segment by a bistable switch (152; 356) . We can write the time- and space-dependent membrane potential as

crLtu(t,x)=2x2u(t,x)-rLg1-au(t,x)(u(t,x)-1)(u(t,x)-a).c\,r_{L}\frac{\partial}{\partial t}\,u(t,x)=\frac{\partial^{2}}{\partial x^{2}% }\,u(t,x)-\frac{r_{L}g}{1-a}u(t,x)(u(t,x)-1)(u(t,x)-a)\,. (3.39)

where a<1/2a<1/2 and gg are parameters. The membrane potential is scaled such that it rests at zero but may be activated to u=1u=1. The reduced model can switch between u=0u=0 and u=1u=1 if it is pushed above u=au=a, but does not reproduce the full upswing followed by downswing of action potentials.

It turns out that Eq. 3.39 can also be interpreted as a model of flame front propagation. The solution of this equation follows (564)

u(x,t)=11+exp(x-𝗏t2λ)u(x,t)=\frac{1}{1+\exp\left(\frac{x-{\sf v}t}{\sqrt{2}\lambda^{\ast}}\right)} (3.40)

with traveling speed

𝗏=c(1-2a)2(1-a)rL/g.{\sf v}=\frac{c(1-2a)}{\sqrt{2(1-a)r_{L}/g}}. (3.41)

The propagation velocity depends on the capacitance per unit length cc, the longitudinal resistance per unit length rLr_{L}, and the excitability parameters gg and aa.

How does the conduction velocity scale with axon size? Since rLr_{L}, cc and gg themselves depend on the diameter of the axon, we expect the velocity to reflect that relationship. The parameters cc and gg scale with the circumference of the cellular membrane and therefore scale linearly with the radius ρ\rho. The cytoplasmic resistance per unit length, however, scales with the cross-sectional area, rLρ2r_{L}\propto\rho^{2}. With these relations in mind, eq. (3.41) shows that the conduction velocity is proportional to the square root of the diameter 𝗏ρ{\sf v}\propto\sqrt{\rho}. Therefore, increasing the diameter improves propagation velocity. This is thought to be the reason why the unmyelinated axons that Hodgkin and Huxley studied were so large (up to ρ=500\rho=500 μ\mum).

3.3.2 Myelinated Axons

Myelinated axons have sodium and potassium channels only in restricted segments called nodes of Ranvier. These nodes form only 0.2% of the axonal length, the rest is considered a passive membrane that is wrapped into a myelin sheath. Myelin mainly decreases the membrane capacitance CC and increase the resistance RTR_{T} by a factor of up to 300 (118). Ions are trapped by myelin since it prevents them to either flow outside the axon or accumulate on the membrane. Instead, ions flow in and out of the the nodes such that an ion leaving a node of Ranvier forces another to enter the following node. Assuming that the nodes are equally separated by a myelinated segment of length LL, we can model the evolution of the membrane potential at each node unu_{n}. The dynamics of idealized myelinated axons follow Kirchoff’s equation with a resistance RL=LrLR_{L}=L\,r_{L} replacing the myelinated segment

Cdundt=1LrL(un+1(t)-2un(t)+un-1(t))-ionIion,n(t)C\frac{{\text{d}}u_{n}}{{\text{d}}t}=\frac{1}{L\,r_{L}}(u_{n+1}(t)-2u_{n}(t)+u% _{n-1}(t))-\sum_{ion}I_{ion,n}(t) (3.42)

where CC is the total capacitance of the node. This equation was encountered in the derivation of the Cable Equation (Section 3.17). The conduction velocity is greatly increased by myelin such that some nerves reach 70-80 m/s(118).

Example: Propagation Speed with Simplified Action Potential Dynamics

Using the simplification of the ion channel dynamics in Eq. (3.39) for each node

C   d   un   d   t=1LrL(un+1(t)-2un(t)+un-1(t))-g1-aun(t)(un(t)-1)(un(t)-a)C\frac{{\text{d}}u_{n}}{{\text{d}}t}=\frac{1}{L\,r_{L}}(u_{n+1}(t)-2u_{n}(t)+u% _{n-1}(t))-\frac{g}{1-a}u_{n}(t)(u_{n}(t)-1)(u_{n}(t)-a) (3.43)

where gg and a<1/2a<1/2 are parameters regulating the excitability of the node. Unlike the Eq. (3.39), the parameter gg has units of conductance per node since we nodes of Ranvier discrete segments. An activated node may fail to excite the adjacent nodes if the membrane potential does not reach u=au=a. In this model, the internodal distance must satisfy (142)

L<L=1-a4a2rLgL<L^{\ast}=\frac{1-a}{4a^{2}r_{L}g} (3.44)

for propagation to be sustained. When the internodal distance LL is smaller than LL^{\ast}, propagation will fail. When the internodal distance is larger, propagation will succeed.

The propagation velocity for small L-LL^{\ast}-L follows (59)

𝗏πga(1-a)CL(L-L){\sf v}\approx\frac{\pi g\,a}{(1-a)C}\sqrt{L(L^{\ast}-L)} (3.45)

which is maximum at L=L/2L=L^{\ast}/2. Since in most myelinated axons, internodal distance scales linearly with their radius (545), the velocity of myelinated axons also scales linearly with radius, 𝗏Lρ{\sf v}\propto L\propto\rho.