Soliton model in neuroscience

History

Ichiji Tasaki pioneered a thermodynamic approach to the phenomenon of nerve pulse propagation which identified several phenomena that were not included in the Hodgkin–Huxley model. Along with measuring various non electrical components of a nerve impulse, Tasaki investigated the physical chemistry of phase transitions in nerve fibers and its importance for nerve pulse propagation. Based on Tasaki's work, Konrad Kaufman proposed sound waves as a physical basis for nerve pulse propagation in an unpublished manuscript. The basic idea at the core of the soliton model is the balancing of intrinsic dispersion of the two dimensional sound waves in the membrane by nonlinear elastic properties near a phase transition. The initial impulse can acquire a stable shape under such circumstances, in general known as a solitary wave. Solitons are the simplest solution of the set of nonlinear wave equations governing such phenomenon and were applied to model nerve impulse in 2005 by Thomas Heimburg and Andrew D. Jackson, both at the Niels Bohr Institute of the University of Copenhagen. Heimburg heads the institute's Membrane Biophysics Group. The biological physics group of Matthias Schneider has studied propagation of two-dimensional sound waves in lipid interfaces and their possible role in biological signalling

Justification

The model starts with the observation that cell membranes always have a freezing point (the temperature below which the consistency changes from fluid to gel-like) only slightly below the organism's body temperature, and this allows for the propagation of solitons. An action potential traveling along a mixed nerve results in a slight increase in temperature followed by a decrease in temperature. Soliton model proponents claim that no net heat is released during the overall pulse and that the observed temperature changes are inconsistent with the Hodgkin-Huxley model. However, this is untrue: the Hodgkin Huxley model predicts a biphasic release and absorption of heat. In addition, the AP causes a slight local thickening of the membrane and a force acting outwards; this effect is not predicted by the Hodgkin–Huxley model but does not contradict it, either.

The soliton model attempts to explain the electrical currents associated with the action potential as follows: the traveling soliton locally changes density and thickness of the membrane, and since the membrane contains many charged and polar substances, this will result in an electrical effect, akin to piezoelectricity. Indeed, such nonlinear sound waves have now been shown to exist at lipid interfaces that show superficial similarity to action potentials (electro-opto-mechanical coupling, velocities, biphasic pulse shape, threshold for excitation etc.). Furthermore, the waves remain localized in the membrane and do not spread out in the surrounding due to an impedance mismatch. However these mechanical waves differ in fundamental ways from the properties described by the Hodgkin-Huxley Model in that they do not exhibit the very specific ionic selectivity, voltage sensitivity, time-course, or sensitivity to toxins that are integral to the Hodgkin-Huxley model.

Formalism

The soliton representing the action potential of nerves is the solution of the partial differential equation

∂ 2 Δ ρ ∂ t 2 = ∂ ∂ x [ ( c 0 2 + p Δ ρ + q Δ ρ 2 ) ∂ Δ ρ ∂ x ] − h ∂ 4 Δ ρ ∂ x 4 ,

where t is time and x is the position along the nerve axon. Δρ is the change in membrane density under the influence of the action potential, c₀ is the sound velocity of the nerve membrane, p and q describe the nature of the phase transition and thereby the nonlinearity of the elastic constants of the nerve membrane. The parameters c₀, p and q are dictated by the thermodynamic properties of the nerve membrane and cannot be adjusted freely. They have to be determined experimentally. The parameter h describes the frequency dependence of the sound velocity of the membrane (dispersion relation). The above equation does not contain any fit parameters. It is formally related to the Boussinesq approximation for solitons in water canals. The solutions of the above equation possess a limiting maximum amplitude and a minimum propagation velocity that is similar to the pulse velocity in myelinated nerves. Under restrictive assumptions, there exist periodic solutions that display hyperpolarization and refractory periods.

Role of ion channels

Advocates of the soliton model claim that it explains several aspects of the action potential, which are not explained by the Hodgkin–Huxley model. Since it is of thermodynamic nature it does not address the properties of single macromolecules like ion channel proteins on a molecular scale. It is rather assumed that their properties are implicitly contained in the macroscopic thermodynamic properties of the nerve membranes. The soliton model predicts membrane current fluctuations during the action potential. These currents are of similar appearance as those reported for ion channel proteins. They are thought to be caused by lipid membrane pores spontaneously generated by the thermal fluctuations. Unfortunately for the model, however, the similarities are superficial. Such thermal fluctuations cannot explain the specific ionic selectivity, the specific time-course of the response to voltage changes, or the pharmacology of neuron channels.

Application to anesthesia

The authors claim that their model explains the previously obscure mode of action of numerous anesthetics. The Meyer–Overton observation holds that the strength of a wide variety of chemically diverse anesthetics is proportional to their lipid solubility, suggesting that they do not act by binding to specific proteins such as ion channels but instead by dissolving in and changing the properties of the lipid membrane. Dissolving substances in the membrane lowers the membrane's freezing point, and the resulting larger difference between body temperature and freezing point inhibits the propagation of solitons. By increasing pressure, lowering pH or lowering temperature, this difference can be restored back to normal, which should cancel the action of anesthetics: this is indeed observed. The amount of pressure needed to cancel the action of an anesthetic of a given lipid solubility can be computed from the soliton model and agrees reasonably well with experimental observations.

Differences between model predictions and experimental observations

The following is a list of some of the disagreements between experimental observations and the "Soliton model":

Antidromic invasion of soma from axon

An action potential initiated anywhere on an axon will travel in an antidromic (backward) direction to the neuron soma (cell body) without loss of amplitude and produce a full-amplitude action potential in the soma. As the membrane area of the soma is orders of magnitude larger than the area of the axon, conservation of energy requires that an adiabatic mechanical wave decrease in amplitude. Since the absence of heat production is one of the claimed justifications of the 'soliton model', this is particularly difficult to explain within that model.

Persistence of action potential over wide temperature range

An important assumption of the soliton model is the presence of a phase transition near the ambient temperature of the axon ("Formalism", above). Then, rapid change of temperature away from the phase transition temperature would necessarily cause large changes in the action potential. Below the phase transition temperature, the solition wave would not be possible. Yet, action potentials are present at 0 °C. The time course is slowed in a manner predicted by the measured opening and closing kinetics of the Hodgkin-Huxley ion channels.

Collisions

Nerve impulses traveling in opposite directions annihilate each other on collision. On the other hand, mechanical waves do not annihilate but pass through each other. Solition model proponents have attempted to show that action potentials can pass through a collision; however, collision annihilation of orthodromic and antidromic action potentials is a routinely observed phenomenon in neuroscience laboratories and are the basis of a standard technique for identification of neurons. Solitons pass each other on collision (Figure--"Collision of Solitions"), solitary waves in general can pass, annihilate or bounce of each other and solitons are only a special case of such solitary waves.

Ionic currents under voltage clamp

The voltage clamp, used by Hodgkin and Huxley (1952) (Hodgkin-Huxley Model) to experimentally dissect the action potential in the squid giant axon, uses electronic feedback to measure the current necessary to hold membrane voltage constant at a commanded value. A silver wire, inserted into the interior of the axon, forces a constant membrane voltage along the length of the axon. Under these circumstances, there is no possibility of a traveling 'soliton'. Any thermodynamic changes are very different from those resulting from an action potential. Yet, the measured currents accurately reproduce the action potential. (citation needed)

Single channel currents

The patch clamp technique isolates a microscopic patch of membrane on the tip of a glass pipette. It is then possible to record currents from single ionic channels. There is no possibility of propagating solitons or thermodynamic changes. Yet, the properties of these channels (temporal response to voltage jumps, ionic selectivity) accurately predict the properties of the macroscopic currents measured under conventional voltage clamp. (citation needed)

Saltatory transmission at nodes of Ranvier

Rapid transmission of the action potential in myelinated axons depends on electrical currents across inactive regions of membrane. The 'soliton model' cannot account for this phenomenon. (citation needed)

Selective ionic conductivity

The current underlying the action potential depolarization is selective for sodium. Repolarization depends on a selective potassium current. These currents have very specific responses to voltage changes which quantitatively explain the action potential. Substitution of non-permeant ions for sodium abolishes the action potential. The 'solition model' cannot explain either the ionic selectivity or the responses to voltage changes. ( Hodgkin-Huxley Model )(citation needed)

Pharmacology

The drug tetrodotoxin (TTX) blocks action potentials at extremely low concentrations. The site of action of TTX on the sodium channel has been identified. Dendrotoxins block the potassium channels. These drugs produce quantitatively predictable changes in the action potential. The 'soliton model' provides no explanation for these pharmacological effects.

Action waves

A recent theoretical model, proposed by Ahmed El Hady and Benjamin Machta, proposes that there is a mechanical surface wave which co-propagates with the electrical action potential. These surface waves are called "action waves". In the El Hady–Machta's model, these co-propagating waves are driven by voltage changes across the membrane caused by the action potential.