# Soliton model in neuroscience

Nonlinear electro-mechanical wave measured in an artificial lipid system

The soliton hypothesis in neuroscience is a model that claims to explain how action potentials are initiated and conducted along axons based on a thermodynamic theory of nerve pulse propagation.[1] It proposes that the signals travel along the cell's membrane in the form of certain kinds of solitary sound (or density) pulses that can be modeled as solitons. The model is proposed as an alternative to the Hodgkin–Huxley model[2] in which action potentials: voltage-gated ion channels in the membrane open and allow sodium ions to enter the cell (inward current). The resulting decrease in membrane potential opens nearby voltage-gated sodium channels, thus propagating the action potential. The transmembrane potential is restored by delayed opening of potassium channels. Soliton hypothesis proponents assert that energy is mainly conserved during propagation except dissipation losses; however, measured temperature changes are also consistent with the Hodgkin-Huxley model.[3][4]

The soliton model ( and sound waves in general) depends on adiabatic propagation in which the energy provided at the source of excitation is carried adiabatically through the medium, i.e. plasma membrane. The measurement of a temperature pulse and the claimed absence of heat release during an action potential[5][6] were the basis of the proposal that nerve impulses are an adiabatic phenomenon much like sound waves. Synaptically evoked action potentials in the electric organ of the electric eel are associated with substantial positive (only) heat production.[7] In the garfish olfactory nerve, the action potential is associated with a biphasic temperature change; however, there is a net production of heat.[8] These published results are consistent with the Hodgkin-Huxley Model and the authors interpret their work in terms of that model: The initial sodium current releases heat as the membrane capacitance is discharged; heat is absorbed during recharge of the membrane capacitance as potassium ions move with their concentration gradient but against the membrane potential. This mechanism is called the "Condenser Theory". Additional heat may be generated by membrane configuration changes driven by the changes in membrane potential. An increase in entropy during depolarization would release heat; entropy increase during repolarization would absorb heat.[9]

## 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.[10] 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.[11] 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.[12] 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,[13][14][15] 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 [16][17][18][19]

## 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.[20] 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.[9] In addition, the AP causes a slight local thickening of the membrane and a force acting outwards;[21] 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.).[17] Furthermore, the waves remain localized in the membrane and do not spread out in the surrounding due to an impedance mismatch.[22] 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

${\displaystyle {\frac {\partial ^{2}\Delta \rho }{\partial t^{2}}}={\frac {\partial }{\partial x}}\left[\left(c_{0}^{2}+p\Delta \rho +q\Delta \rho ^{2}\right){\frac {\partial \Delta \rho }{\partial x}}\right]-h{\frac {\partial ^{4}\Delta \rho }{\partial 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, c0 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 c0, 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.[23]

## 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.[24] They are thought to be caused by lipid membrane pores spontaneously generated by the thermal fluctuations. Such thermal fluctuations cannot explain the specific ionic selectivity or the specific time-course of the response to voltage changes.

## 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.[25] 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

Collision of solitons

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.[26][citation needed]
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 soliton 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.[27]
Collisions
Nerve impulses traveling in opposite directions annihilate each other on collision.[28] On the other hand, mechanical waves do not annihilate but pass through each other. Soliton model proponents have attempted to show that action potentials can pass through a collision;[29] 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.[30] Solitons pass each other on collision (Figure--"Collision of Solitons"), solitary waves in general can pass, annihilate or bounce of each other [31] and solitons are only a special case of such solitary waves.[32]
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.[33] (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 'soliton 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.[34] Dendrotoxins block the potassium channels. These drugs produce quantitatively predictable changes in the action potential.[33] 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".[35] In the El Hady–Machta's model, these co-propagating waves are driven by voltage changes across the membrane caused by the action potential.

## References

1. ^ Andersen, S; Jackson, A; Heimburg, T (2009). "Towards a thermodynamic theory of nerve pulse propagation" (PDF). Progress in Neurobiology. 88: 104–113. doi:10.1016/j.pneurobio.2009.03.002.
2. ^ Hodgkin AL, Huxley AF, Katz B (1952). "Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo". Journal of Physiology. 116 (4): 424&ndash, 448. doi:10.1113/jphysiol.1952.sp004717. PMC 1392213. PMID 14946713.
Hodgkin AL, Huxley AF (1952). "Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo". Journal of Physiology. 116 (4): 449&ndash, 472. doi:10.1113/jphysiol.1952.sp004717. PMC 1392213. PMID 14946713.
Hodgkin AL, Huxley AF (1952). "The components of membrane conductance in the giant axon of Loligo". J Physiol. 116 (4): 473&ndash, 496. doi:10.1113/jphysiol.1952.sp004718. PMC 1392209. PMID 14946714.
Hodgkin AL, Huxley AF (1952). "The dual effect of membrane potential on sodium conductance in the giant axon of Loligo". J Physiol. 116 (4): 497&ndash, 506. doi:10.1113/jphysiol.1952.sp004719. PMC 1392212. PMID 14946715.
Hodgkin AL, Huxley AF (1952). "A quantitative description of membrane current and its application to conduction and excitation in nerve". J Physiol. 117 (4): 500&ndash, 544. doi:10.1113/jphysiol.1952.sp004764. PMC 1392413. PMID 12991237.
3. ^ Margineanu, D.-G; Schoffeniels, E. (1977). "Molecular events and energy changes during the action potential" (PDF). PNAS. 74: 3810–3813. Bibcode:1977PNAS...74.3810M. doi:10.1073/pnas.74.9.3810. PMC 431740.
4. ^ Hasenstaub, A; Callaway, E; Otte, S; Sejnowski, T (2010). "Metabolic cost as a unifying principle governing neuronal biophysics". Proceedings of the National Academy of Sciences of the USA. 107 (27): 12329–12334. Bibcode:2010PNAS..10712329H. doi:10.1073/pnas.0914886107. PMC 2901447.
5. ^ Tasaki, Ichiji (13 October 1995). "Mechanical and Thermal Changes in the Torpedo Electric Organ Associated with Its Postsynaptic Potentials". Biochemical and Biophysical Research Communications. 215 (2): 654–658. doi:10.1006/bbrc.1995.2514. PMID 7488005.
6. ^ Howarth, J V; Keynes, R D; Ritchie, J M; Muralt, A von (1 Jul 1975). "The heat production associated with the passage of a single impulse in pike olfactory nerve fibres" (PDF). The Journal of Physiology. 249 (2): 349–368. doi:10.1113/jphysiol.1975.sp011019. PMC 1309578.
7. ^ Tasaki, I; Byrne, P. M. (1993). "Rapid heat production associated with electrical excitation of the electric organs of the electric eel". Biochem Biophys Res Commun. 197: 910–915. doi:10.1006/bbrc.1993.2565.
8. ^ Tasaki, K; Kusano, K; Byrne, PM (1989). "Rapid thermal and mechanical changes in garfish olfactory nerve associated with a propagated impulse". Biophys J. 55: 1033–1040. Bibcode:1989BpJ....55.1033T. doi:10.1016/s0006-3495(89)82902-9. PMC 1330571. PMID 2765644.
9. ^ a b Howarth, J. V. (1975). "Heat Production in Non-Myelinated Nerves". Philosophical Transactions of the Royal Society. 270: 425–432. Bibcode:1975RSPTB.270..425H. doi:10.1098/rstb.1975.0020. JSTOR 2417341.
10. ^ Tasaki, Ichiji (1982). Physiology and Electrochemistry of Nerve Fibers. Bethesda, Maryland: Academic Press Inc. (London). ISBN 0-12-683780-5.
11. ^ Kaufmann, Konrad (1989). Action Potentials and Electrochemical Coupling in the Macroscopic Chiral Phospholipid Membrane. Caruaru, Brazil.
12. ^ Xin-Yi, Wang (1985). "Solitary wave and nonequilibrium phase transition in liquid crystals". Physical Review A. 32: 3126–3129. Bibcode:1985PhRvA..32.3126X. doi:10.1103/PhysRevA.32.3126.
13. ^ Heimburg, T., Jackson, A.D. (12 July 2005). "On soliton propagation in biomembranes and nerves". Proc. Natl. Acad. Sci. U.S.A. 102 (2): 9790–9795. Bibcode:2005PNAS..102.9790H. doi:10.1073/pnas.0503823102. PMC 1175000. PMID 15994235.
14. ^ Heimburg, T., Jackson, A.D. (2007). "On the action potential as a propagating density pulse and the role of anesthetics". Biophys. Rev. Lett. 2: 57–78. arXiv:physics/0610117. Bibcode:2006physics..10117H. doi:10.1142/S179304800700043X.
15. ^ Andersen, S.S.L., Jackson, A.D., Heimburg, T. (2009). "Towards a thermodynamic theory of nerve pulse propagation". Prog. Neurobiol. 88 (2): 104–113. doi:10.1016/j.pneurobio.2009.03.002.
16. ^ Griesbauer, J; Bossinger, S; Wixforth, A; Schneider, M (9 May 2012). "Propagation of 2D Pressure Pulses in Lipid Monolayers and Its Possible Implications for Biology". Physical Review Letters. 108: 198103. arXiv:1211.4104. Bibcode:2012PhRvL.108s8103G. doi:10.1103/PhysRevLett.108.198103. PMID 23003093.
17. ^ a b Shrivastava, Shamit; Schneider, Matthias (18 June 2014). "Evidence for two-dimensional solitary sound waves in a lipid controlled interface and its implications for biological signalling". Royal Society Interface. 11: 20140098. doi:10.1098/rsif.2014.0098. PMC 4078894.
18. ^ Griesbauer, J; Bossinger, S; Wixforth, A; Schneider, M (19 Dec 2012). "Simultaneously propagating voltage and pressure pulses in lipid monolayers of pork brain and synthetic lipids". Physical Review E. 86. arXiv:1211.4105. Bibcode:2012PhRvE..86f1909G. doi:10.1103/PhysRevE.86.061909.
19. ^ Shrivastava, Shamit (Jan 2014). NON-LINEAR SOLITARY SOUND WAVES IN LIPID MEMBRANES AND THEIR POSSIBLE ROLE IN BIOLOGICAL SIGNALING (1st ed.). Boston, MA 02215 US: Thesis, Boston University.
20. ^ Abbott, B.C., Hill, A.V., Howarth, J.V. (1958). "The positive and negative heat associated with a nerve impulse". Proceedings of the Royal Society B. 148 (931): 149–187. Bibcode:1958RSPSB.148..149A. doi:10.1098/rspb.1958.0012.
21. ^ Iwasa, K., Tasaki I., Gibbons, R. (1980). "Swelling of nerve fibres associated with action potentials". Science. 210 (4467): 338–9. Bibcode:1980Sci...210..338I. doi:10.1126/science.7423196. PMID 7423196.
22. ^ Griesbauer, J; Wixforth, A; Schneider, M F (15 Nov 2009). "Wave Propagation in Lipid Monolayers" (PDF). Biophysical Journal. 97 (10): 2710–2716. Bibcode:2009BpJ....97.2710G. doi:10.1016/j.bpj.2009.07.049. PMC 2776282. PMID 19917224.
23. ^ Villagran Vargas, E., Ludu, A., Hustert, R., Gumrich, P., Jackson, A.D., Heimburg, T. (2010). "Periodic solutions and refractory periods in the soliton theory for nerves and the locust femoral nerve". Biophysical Chemistry. 153 (2–3): 159–167. arXiv:1006.3281. doi:10.1016/j.bpc.2010.11.001. PMID 21177017.
24. ^ Heimburg, T. (2010). "Lipid Ion Channels". Biophys. Chem. 150 (1–3): 2–22. arXiv:1001.2524. doi:10.1016/j.bpc.2010.02.018. PMID 20385440.
25. ^ Heimburg, T., Jackson, A.D. (2007). "The thermodynamics of general anesthesia". Biophys. J. 92 (9): 3159–65. arXiv:physics/0610147. Bibcode:2007BpJ....92.3159H. doi:10.1529/biophysj.106.099754. PMC 1852341. PMID 17293400.
26. ^ Rall, W and Shepherd, GM (1968) Theoretical reconstructions of dendrodendritic synaptic interactions in the olfactory bulb. J Neurophysiol 31, 884-915.http://jn.physiology.org/content/jn/31/6/884.full.pdf
27. ^ Hodgkin; Katz (1949). "The Effect of Temperature on the Electrical Activity of the Giant Axon of the Squid" (PDF). J. Physiol. 109: 240–249. doi:10.1113/jphysiol.1949.sp004388. PMC 1392577. PMID 15394322.
28. ^ Tasaki, Ichiji (1949). "Collision of two nerve impulses in the nerve fiber". Biochim Biophys Acta. 3: 494–497. doi:10.1016/0006-3002(49)90121-3.
29. ^ Gonzalez, Alfredo; Budvytyte, Rima; Mosgaard, Lars D; Nissen, Søren; Heimburg, Thomas (10 Sep 2014). "Penetration of Action Potentials During Collision in the Median and Lateral Giant Axons of Invertebrates". Physical Review X. 4. arXiv:1404.3643. Bibcode:2014PhRvX...4c1047G. doi:10.1103/PhysRevX.4.031047.