Dynamics of signal propagation and collision in axons

How spatiotemporal dynamics can enhance ecosystem resilience

We study how self-organization in systems showing complex spatiotemporal dynamics can increase ecosystem resilience. We consider a general simple model that includes positive feedback as well as negative feedback mediated by an inhibitor. We apply this model to Posidonia oceanica meadows, where positive and negative feedbacks are well documented, and there is empirical evidence of the role of sulfide accumulation, toxic for the plant, in driving complex spatiotemporal dynamics. We describe a progressive transition from homogeneous meadows to extinction through dynamical regimes that allow the ecosystem to avoid the typical ecological tipping points of homogeneous vegetation covers. A predictable sequence of distinct dynamical regimes is observed as mortality is continuously increased: turbulent regimes, formation of spirals and wave trains, and isolated traveling pulses or expanding rings, the latter being a harbinger of ecosystem collapse, however far beyond the tipping point of the homogeneous cover. The model used in this paper is general, and the results can be applied to other plant-soil spatially extended systems, regardless of the mechanisms behind negative and positive feedbacks.

A Computational Study on the Activation of Neural Transmission in Deep Brain Stimulation

Deep brain stimulation (DBS) is an established treatment for neurodegenerative movement disorders such as Parkinson's disease that mitigates symptoms by overwriting pathological signals from the central nervous system to the motor system. Nearly all computational models of DBS, directly or indirectly, associate clinical improvements with the extent of fiber activation in the vicinity of the stimulating electrode. However, it is not clear how such activation modulates information transmission. Here, we use the exact cable equation for straight or curved axons and show that DBS segregates the signaling pathways into one of the three communicational modes: complete information blockage, uni-, and bi-directional transmission. Furthermore, all these modes respond to the stimulating pulse in an asynchronous but frequency-locked fashion. Asynchrony depends on the geometry of the axon, its placement and orientation, and the stimulation protocol. At the same time, the electrophysiology of the nerve determines frequency-locking. Such a trimodal response challenges the notion of activation as a binary state and studies that correlate it with the DBS outcome. Importantly, our work suggests that a mechanistic understanding of DBS action relies on distinguishing between these three modes of information transmission.

Advances in using ultrasound to regulate the nervous system

Ultrasound is a mechanical vibration with a frequency greater than 20 kHz. Due to its high spatial resolution, good directionality, and convenient operation in neural regulation, it has recently received increasing attention from scientists. However, the mechanism by which ultrasound regulates the nervous system is still unclear. This article mainly explores the possible mechanisms of ultrasound's mechanical effects, cavitation effects, thermal effects, and the rise of sonogenetics. In addition, the essence of action potential and its relationship with ultrasound were also discussed. Traditional theory treats nerve impulses as pure electrical signals, similar to cable theory. However, this theory cannot explain the phenomenon of inductance and cell membrane bulging out during the propagation of action potential. Therefore, the flexoelectric effect of cell membrane and soliton model reveal that action potential may also be a mechanical wave. Finally, we also elaborated the therapeutic effect of ultrasound on nervous system disease such as epilepsy, Parkinson's disease, and Alzheimer's disease.

Thinking about the action potential: the nerve signal as a window to the physical principles guiding neuronal excitability

Ever since the work of Edgar Adrian, the neuronal action potential has been considered as an electric signal, modeled and interpreted using concepts and theories lent from electronic engineering. Accordingly, the electric action potential, as the prime manifestation of neuronal excitability, serving processing and reliable "long distance" communication of the information contained in the signal, was defined as a non-linear, self-propagating, regenerative, wave of electrical activity that travels along the surface of nerve cells. Thus, in the ground-breaking theory and mathematical model of Hodgkin and Huxley (HH), linking Nernst's treatment of the electrochemistry of semi-permeable membranes to the physical laws of electricity and Kelvin's cable theory, the electrical characteristics of the action potential are presented as the result of the depolarization-induced, voltage- and time-dependent opening and closure of ion channels in the membrane allowing the passive flow of charge, particularly in the form of Na+ and K+ -ions, into and out of the neuronal cytoplasm along the respective electrochemical ion gradient. In the model, which treats the membrane as a capacitor and ion channels as resistors, these changes in ionic conductance across the membrane cause a sudden and transient alteration of the transmembrane potential, i.e., the action potential, which is then carried forward and spreads over long(er) distances by means of both active and passive conduction dependent on local current flow by diffusion of Na+ ion in the neuronal cytoplasm. However, although highly successful in predicting and explaining many of the electric characteristics of the action potential, the HH model, nevertheless cannot accommodate the various non-electrical physical manifestations (mechanical, thermal and optical changes) that accompany action potential propagation, and for which there is ample experimental evidence. As such, the electrical conception of neuronal excitability appears to be incomplete and alternatives, aiming to improve, extend or even replace it, have been sought for. Commonly misunderstood as to their basic premises and the physical principles they are built on, and mistakenly perceived as a threat to the generally acknowledged explanatory power of the "classical" HH framework, these attempts to present a more complete picture of neuronal physiology, have met with fierce opposition from mainstream neuroscience and, as a consequence, currently remain underdeveloped and insufficiently tested. Here we present our perspective that this may be an unfortunate state of affairs as these different biophysics-informed approaches to incorporate also non-electrical signs of the action potential into the modeling and explanation of the nerve signal, in our view, are well suited to foster a new, more complete and better integrated understanding of the (multi)physical nature of neuronal excitability and signal transport and, hence, of neuronal function. In doing so, we will emphasize attempts to derive the different physical manifestations of the action potential from one common, macroscopic thermodynamics-based, framework treating the multiphysics of the nerve signal as the inevitable result of the collective material, i.e., physico-chemical, properties of the lipid bilayer neuronal membrane (in particular, the axolemma) and/or the so-called ectoplasm or membrane skeleton consisting of cytoskeletal protein polymers, in particular, actin fibrils. Potential consequences for our view of action potential physiology and role in neuronal function are identified and discussed.

Bifurcation structure of traveling pulses in type-I excitable media

We study the scenario in which traveling pulses emerge in a prototypical type-I one-dimensional excitable medium, which exhibits two different routes to excitable behavior, mediated by a homoclinic (saddle-loop) and a saddle-node on the invariant cycle bifurcations. We characterize the region in parameter space in which traveling pulses are stable together with the different bifurcations behind either their destruction or loss of stability. In particular, some of the bifurcations delimiting the stability region have been connected, using singular limits, with the two different scenarios that mediated type-I local excitability. Finally, the existence of traveling pulses has been linked to a drift pitchfork instability of localized steady structures.

Traveling pulses in type-I excitable media

We consider a general model exhibiting type-I excitability mediated by a homoclinic and a saddle node on the invariant circle bifurcations. We show how the distinct properties of type-I with respect to type-II excitability confer unique features to traveling pulses in excitable media. They inherit the characteristic divergence of type-I excitable trajectories at threshold exhibiting analogous scalings in the spatial thickness of the pulses. Our results pave the way to identify basic underlying mechanisms behind type-I excitable pulses based solely on the characteristics of the pulse.

Shock and detonation waves at an interface and the collision of action potentials

Action potentials in neurons are known to annihilate each other upon collision, while there are cases where they might penetrate each other. The fate of two waves upon collision is critically dependent on the underlying mechanism of propagation and therefore an understanding of possible outcomes of collision under different conditions is important. Previously, compression waves that travel within the plasma membrane of a neuron have been proposed as a thermodynamic basis for the propagation of action potentials. In this context, it was recently shown that two-dimensional compressive shock waves in the model system of lipid monolayers behave strikingly similar to action potentials in neurons and can even annihilate each other upon head-on collision. However, even a qualitative mechanism remained unclear. To this end, we summarise the fundamentals of shock physics as applied to an interface and recap how it explained the observation of threshold and saturation of shockwaves in the lipid monolayer (all - or - none). We then compare the theory with the soliton model that has the same fundamental premise, i.e. the conservation laws and thermodynamics, and was previously proposed as a model for the nerve pulse propagation. We elaborate on how the two approaches make different predictions with regards to collisions and the detailed structure of the wave-front. As a case study and a new qualitative result, we finally show that previously unexplained annihilation of shock waves in the lipid monolayer is a direct consequence of the nature of state changes, i.e. jump conditions, within these shockwaves.

Electrophysiological-mechanical coupling in the neuronal membrane and its role in ultrasound neuromodulation and general anaesthesia

The current understanding of the role of the cell membrane is in a state of flux. Recent experiments show that conventional models, considering only electrophysiological properties of a passive membrane, are incomplete. The neuronal membrane is an active structure with mechanical properties that modulate electrophysiology. Protein transport, lipid bilayer phase, membrane pressure and stiffness can all influence membrane capacitance and action potential propagation. A mounting body of evidence indicates that neuronal mechanics and electrophysiology are coupled, and together shape the membrane potential in tight coordination with other physical properties. In this review, we summarise recent updates concerning electrophysiological-mechanical coupling in neuronal function. In particular, we aim at making the link with two relevant yet often disconnected fields with strong clinical potential: the use of mechanical vibrations-ultrasound-to alter the electrophysiogical state of neurons, e.g., in neuromodulation, and the theories attempting to explain the action of general anaesthetics. STATEMENT OF SIGNIFICANCE: General anaesthetics revolutionised medical practice; now an apparently unrelated technique, ultrasound neuromodulation-aimed at controlling neuronal activity by means of ultrasound-is poised to achieve a similar level of impact. While both technologies are known to alter the electrophysiology of neurons, the way they achieve it is still largely unknown. In this review, we argue that in order to explain their mechanisms/effects, the neuronal membrane must be considered as a coupled mechano-electrophysiological system that consists of multiple physical processes occurring concurrently and collaboratively, as opposed to sequentially and independently. In this framework the behaviour of the cell membrane is not the result of stereotypical mechanisms in isolation but instead emerges from the integrative behaviour of a complexly coupled multiphysics system.

Computational model of the mechanoelectrophysiological coupling in axons with application to neuromodulation

For more than half a century, the action potential (AP) has been considered a purely electrical phenomenon. However, experimental observations of membrane deformations occurring during APs have revealed that this process also involves mechanical features. This discovery has recently fuelled a controversy on the real nature of APs: whether they are mechanical or electrical. In order to examine some of the modern hypotheses regarding APs, we propose here a coupled mechanoelectrophysiological membrane finite-element model for neuronal axons. The axon is modeled as an axisymmetric thin-wall cylindrical tube. The electrophysiology of the membrane is modeled using the classic Hodgkin-Huxley (H-H) equations for the Nodes of Ranvier or unmyelinated axons and the cable theory for the internodal regions, whereas the axonal mechanics is modeled by means of viscoelasticity theory. Membrane potential changes induce a strain gradient field via reverse flexoelectricity, whereas mechanical pulses result in an electrical self-polarization field following the direct flexoelectric effect, in turn influencing the membrane potential. Moreover, membrane deformation also alters the values of membrane capacitance and resistance in the H-H equation. These three effects serve as the fundamental coupling mechanisms between the APs and mechanical pulses in the model. A series of numerical studies was systematically conducted to investigate the consequences of interaction between the APs and mechanical waves on both myelinated and unmyelinated axons. Simulation results illustrate that the AP is always accompanied by an in-phase propagating membrane displacement of ≈1nm, whereas mechanical pulses with enough magnitude can also trigger APs. The model demonstrates that mechanical vibrations, such as the ones arising from ultrasound stimulations, can either annihilate or enhance axonal electrophysiology depending on their respective directionality and frequency. It also shows that frequency of pulse repetition can also enhance signal propagation independently of the amplitude of the signal. This result not only reconciles the mechanical and electrical natures of the APs but also provides an explanation for the experimentally observed mechanoelectrophysiological phenomena in axons, especially in the context of ultrasound neuromodulation.

Temperature effects on neuronal firing rates and tonic-to-bursting transitions

Temperature fluctuations can affect neurological processes at a variety of levels, with the overall output that higher temperatures in general increase neuronal activity. While variations in firing rates can happen with the neuronal system maintaining its homeostatic firing pattern of tonic firing, or bursting, changes in firing rates can also be associated with transitions between the two patterns of firing. Our computer simulations suggest a possible mechanism directly related to the shortening of the duration of the action potential for higher firing rates with temperature increase. Increased temperatures also shorten the period doubling cascade and chaos transition between tonic and burting regimes.

Similarities between action potentials and acoustic pulses in a van der Waals fluid

An action potential is typically described as a purely electrical change that propagates along the membrane of excitable cells. However, recent experiments have demonstrated that non-linear acoustic pulses that propagate along lipid interfaces and traverse the melting transition, share many similar properties with action potentials. Despite the striking experimental similarities, a comprehensive theoretical study of acoustic pulses in lipid systems is still lacking. Here we demonstrate that an idealized description of an interface near phase transition captures many properties of acoustic pulses in lipid monolayers, as well as action potentials in living cells. The possibility that action potentials may better be described as acoustic pulses in soft interfaces near phase transition is illustrated by the following similar properties: correspondence of time and velocity scales, qualitative pulse shape, sigmoidal response to stimulation amplitude (an 'all-or-none' behavior), appearance in multiple observables (particularly, an adiabatic change of temperature), excitation by many types of stimulations, as well as annihilation upon collision. An implication of this work is that crucial functional information of the cell may be overlooked by focusing only on electrical measurements.

The Soliton and the Action Potential - Primary Elements Underlying Sentience

At present the neurological basis of sentience is poorly understood and this problem is exacerbated by only a partial knowledge of how one of the primary elements of sentience, the action potential, actually works. This has consequences for our understanding of how communication within the brain and in artificial brain neural networks (BNNs). Reverse engineering models of brain activity assume processing works like a conventional binary computer and neglects speed of cognition, latencies, error in nerve conduction and the true dynamic structure of neural networks in the brain. Any model of nerve conduction that claims inspiration from nature must include these prerequisite parameters, but current western computer modeling of artificial BNNs assumes that the action potential is binary and binary mathematics has been assumed by force of popular acceptance to mediate computation in the brain. Here we present evidence that the action potential is a temporal compound ternary structure, described as the computational action potential (CAP). The CAP contains the refractory period, an analog third phase capable of phase-ternary computation via colliding action potentials. This would best fit a realistic BNN and provides a plausible mechanism to explain transmission, in preference to Cable Theory. The action potential pulse (APPulse), is made up of the action potential combined with a coupled synchronized soliton pressure pulse in the cell membrane. We describe a model of an ion channel in a membrane where a soliton deforms the channel sufficiently to destroy the electrostatic insulation thereby instigating a mechanical contraction across the membrane by electrostatic forces. Such a contraction has the effect of redistributing the force lengthways thereby increasing the volume of the ion channel in the membrane. Na ions, once attracted to the interior, balance the forces and the channel reforms to its original shape. A refractory period then occurs until the Na ions diffuse from the adjacent interior space. Finally, a computational model of the action potential (the CAP) is proposed with single action potentials significantly including the refractory period as a computational element capable of computation between colliding action potentials.

Dynamics of coupled mode solitons in bursting neural networks

Using an electrically coupled chain of Hindmarsh-Rose neural models, we analytically derived the nonlinearly coupled complex Ginzburg-Landau equations. This is realized by superimposing the lower and upper cutoff modes of wave propagation and by employing the multiple scale expansions in the semidiscrete approximation. We explore the modified Hirota method to analytically obtain the bright-bright pulse soliton solutions of our nonlinearly coupled equations. With these bright solitons as initial conditions of our numerical scheme, and knowing that electrical signals are the basis of information transfer in the nervous system, it is found that prior to collisions at the boundaries of the network, neural information is purely conveyed by bisolitons at lower cutoff mode. After collision, the bisolitons are completely annihilated and neural information is now relayed by the upper cutoff mode via the propagation of plane waves. It is also shown that the linear gain of the system is inextricably linked to the complex physiological mechanisms of ion mobility, since the speeds and spatial profiles of the coupled nerve impulses vary with the gain. A linear stability analysis performed on the coupled system mainly confirms the instability of plane waves in the neural network, with a glaring example of the transition of weak plane waves into a dark soliton and then static kinks. Numerical simulations have confirmed the annihilation phenomenon subsequent to collision in neural systems. They equally showed that the symmetry breaking of the pulse solution of the system leaves in the network static internal modes, sometime referred to as Goldstone modes.

Collision of two action potentials in a single excitable cell

BACKGROUND

It is a common incident in nature, that two waves or pulses run into each other head-on. The outcome of such an event is of special interest, because it allows conclusions about the underlying physical nature of the pulses. The present experimental study dealt with the head-on meeting of two action potentials (AP) in a single excitable plant cell (Chara braunii internode).

METHODS

The membrane potential was monitored with multiple sensors along a single excitable cell. In control experiments, an AP was excited electrically at either end of the cell cylinder. Subsequently, stimuli were applied simultaneously at both ends of the cell in order to generate two APs that met each other head-on.

RESULTS

When two action potentials propagated into each other, the pulses did not penetrate but annihilated (N=26 experiments in n=10 cells).

CONCLUSIONS

APs in excitable plant cells did not penetrate upon meeting head-on. In the classical electrical model, this behavior is specifically attributed to relaxation of ion channel proteins. From an acoustic point of view, annihilation can be viewed as a result of nonlinear material properties (e.g. a phase change).

GENERAL SIGNIFICANCE

The present results suggest that APs in excitable animal and plant cells belong to a similar class of nonlinear phenomena. Intriguingly, other excitation waves in biology (intracellular waves, cortical spreading depression, etc.) also annihilate upon collision and are thus expected to follow the same underlying principles as the observed action potentials.

Branching morphology determines signal propagation dynamics in neurons

Computational modeling of signal propagation in neurons is critical to our understanding of basic principles underlying brain organization and activity. Exploring these models is used to address basic neuroscience questions as well as to gain insights for clinical applications. The seminal Hodgkin Huxley model is a common theoretical framework to study brain activity. It was mainly used to investigate the electrochemical and physical properties of neurons. The influence of neuronal structure on activity patterns was explored, however, the rich dynamics observed in neurons with different morphologies is not yet fully understood. Here, we study signal propagation in fundamental building blocks of neuronal branching trees, unbranched and branched axons. We show how these simple axonal elements can code information on spike trains, and how asymmetric responses can emerge in axonal branching points. This asymmetric phenomenon has been observed experimentally but until now lacked theoretical characterization. Together, our results suggest that axonal morphological parameters are instrumental in activity modulation and information coding. The insights gained from this work lay the ground for better understanding the interplay between function and form in real-world complex systems. It may also supply theoretical basis for the development of novel therapeutic approaches to damaged nervous systems.

Solitonic conduction of electrotonic signals in neuronal branchlets with polarized microstructure

A model of solitonic conduction in neuronal branchlets with microstructure is presented. The application of cable theory to neurons with microstructure results in a nonlinear cable equation that is solved using a direct method to obtain analytical approximations of traveling wave solutions. It is shown that a linear superposition of two oppositely directed traveling waves demonstrate solitonic interaction: colliding waves can penetrate through each other, and continue fully intact as the exact pulses that entered the collision. These findings indicate that microstructure when polarized can sustain solitary waves that propagate at a constant velocity without attenuation or distortion in the absence of synaptic transmission. Solitonic conduction in a neuronal branchlet arising from polarizability of its microstructure is a novel signaling mode of electrotonic signals in thin processes (<0.5 μm diameter).

Bifurcation transitions in gap-junction-coupled neurons

Here we investigate transitions occurring in the dynamical states of pairs of distinct neurons electrically coupled, with one neuron tonic and the other bursting. Depending on the dynamics of the individual neurons, and for strong enough coupling, they synchronize either in a tonic or a bursting regime, or initially tonic transitioning to bursting via a period doubling cascade. Certain intrinsic properties of the individual neurons such as minimum firing rates are carried over into the dynamics of the coupled neurons affecting their ultimate synchronous state.