A comparison of computational models for the extracellular potential of neurons

Finite-element modeling of neuromodulation via controlled delivery of potassium ions using conductive polymer-coated microelectrodes

Objective. The controlled delivery of potassium is an interesting neuromodulation modality, being potassium ions involved in shaping neuron excitability, synaptic transmission, network synchronization, and playing a key role in pathological conditions like epilepsy and spreading depression. Despite many successful examples of pre-clinical devices able to influence the extracellular potassium concentration, computational frameworks capturing the corresponding impact on neuronal activity are still missing.Approach. We present a finite-element model describing a PEDOT:PSS-coated microelectrode (herein, simplyionic actuator) able to release potassium and thus modulate the activity of a cortical neuron in anin-vitro-like setting. The dynamics of ions in the ionic actuator, the neural membrane, and the cellular fluids are solved self-consistently.Main results. We showcase the capability of the model to describe on a physical basis the modulation of the intrinsic excitability of the cell and of the synaptic transmission following the electro-ionic stimulation produced by the actuator. We consider three case studies for the ionic actuator with different levels of selectivity to potassium: ideal selectivity, no selectivity, and selectivity achieved by embedding ionophores in the polymer.Significance. This work is the first step toward a comprehensive computational framework aimed to investigate novel neuromodulation devices targeting specific ionic species, as well as to optimize their design and performance, in terms of the induced modulation of neural activity.

Energy optimisation predicts the capacity of ion buffering in the brain

Neurons store energy in the ionic concentration gradients they build across their cell membrane. The amount of energy stored, and hence the work the ions can do by mixing, can be enhanced by the presence of ion buffers in extra- and intracellular space. Buffers act as sources and sinks of ions, however, and unless the buffering capacities for different ion species obey certain relationships, a complete mixing of the ions may be impeded by the physical conditions of charge neutrality and isotonicity. From these conditions, buffering capacities were calculated that enabled each ion species to mix completely. In all valid buffer distributions, the [Formula: see text] ions were buffered most, with a capacity exceeding that of [Formula: see text] and [Formula: see text] buffering by at least an order of magnitude. The similar magnitude of the (oppositely directed) [Formula: see text] and [Formula: see text] gradients made extracellular space behave as a [Formula: see text]-[Formula: see text] exchanger. Anions such as [Formula: see text] were buffered least. The great capacity of the extra- and intracellular [Formula: see text] buffers caused a large influx of [Formula: see text] ions as is typically observed during energy deprivation. These results explain many characteristics of the physiological buffer distributions but raise the question how the brain controls the capacity of its ion buffers. It is suggested that neurons and glial cells, by their great sensitivity to gradients of charge and osmolarity, respectively, sense deviations from electro-neutral and isotonic mixing, and use these signals to tune the chemical composition, and buffering capacity, of the extra- and intracellular matrices.

Nano-scale solution of the Poisson-Nernst-Planck (PNP) equations in a fraction of two neighboring cells reveals the magnitude of intercellular electrochemical waves

The basic building blocks of the electrophysiology of cardiomyocytes are ion channels integrated in the cell membranes. Close to the ion channels there are very strong electrical and chemical gradients. However, these gradients extend for only a few nano-meters and are therefore commonly ignored in mathematical models. The full complexity of the dynamics is modelled by the Poisson-Nernst-Planck (PNP) equations but these equations must be solved using temporal and spatial scales of nano-seconds and nano-meters. Here we report solutions of the PNP equations in a fraction of two abuttal cells separated by a tiny extracellular space. We show that when only the potassium channels of the two cells are open, a stationary solution is reached with the well-known Debye layer close to the membranes. When the sodium channels of one of the cells are opened, a very strong and brief electrochemical wave emanates from the channels. If the extracellular space is sufficiently small and the number of sodium channels is sufficiently high, the wave extends all the way over to the neighboring cell and may therefore explain cardiac conduction even at very low levels of gap junctional coupling.

A phenomenological computational model of the evoked action potential fitted to human cochlear implant responses

There is a growing interest in biomedical engineering in developing procedures that provide accurate simulations of the neural response to electrical stimulus produced by implants. Moreover, recent research focuses on models that take into account individual patient characteristics.

We present a phenomenological computational model that is customized with the patient’s data provided by the electrically evoked compound action potential (ECAP) for simulating the neural response to electrical stimulus produced by the electrodes of cochlear implants (CIs). The model links the input currents of the electrodes to the simulated ECAP.

Potentials and currents are calculated by solving the quasi-static approximation of the Maxwell equations with the finite element method (FEM). In ECAPs recording, an active electrode generates a current that elicits action potentials in the surrounding auditory nerve fibers (ANFs). The sum of these action potentials is registered by other nearby electrode. Our computational model emulates this phenomenon introducing a set of line current sources replacing the ANFs by a set of virtual neurons (VNs). To fit the ECAP amplitudes we assign a suitable weight to each VN related with the probability of an ANF to be excited. This probability is expressed by a cumulative beta distribution parameterized by two shape parameters that are calculated by means of a differential evolution algorithm (DE). Being the weights function of the current density, any change in the design of the CI affecting the current density produces changes in the weights and, therefore, in the simulated ECAP, which confers to our model a predictive capacity.

The results of the validation with ECAP data from two patients are presented, achieving a satisfactory fit of the experimental data with those provided by the proposed computational model.

The cochlea, found in the inner ear, is the organ where the sound is transformed into an electrical pulse to be transmitted by the neurons to the auditory cortex. Hearing loss can be caused by damage to the hair cells, in which case neuronal excitation is impaired. CIs are devices that replace the normal function of the impaired/damaged Organ of Corti. Computational models allow a better understanding of the mechanisms involved in the electrical stimulation of the auditory nerve. These models can help biomedical engineers to develop new CIs with improved auditory performance. One important aspect of our model is its customization with the patient’s data provided by the recording of the evoked compound action potential (the synchronous firing of a population of electrically stimulated auditory nerve fibers). This phenomenological model allows us to predict the registers of neural stimulation produced when the auditory nerve is stimulated with the CIs. We have validated the proposed model with real data obtained from two patients with CIs.

Effect of extracellular volume on the energy stored in transmembrane concentration gradients

The amount of energy that can be retrieved from a concentration gradient across a membrane separating two compartments depends on the relative size of the compartments. Having a larger low-concentration compartment is in general beneficial. It is shown here analytically that the retrieved energy further increases when the high-concentration compartment shrinks during the mixing process, and a general formula is derived for the energy when the ratio of transported solvent to solute varies. These calculations are then applied to the interstitial compartment of the brain, which is rich in Na^{+} and Cl^{-} ions and poor in K^{+}. The reported shrinkage of this compartment, and swelling of the neurons, during oxygen deprivation is shown to enhance the energy recovered from NaCl entering the neurons. The slight loss of energy on the part of K^{+} can be compensated for by the uptake of K^{+} ions by glial cells. In conclusion, the present study proposes that the reported fluctuations in the size of the interstitial compartment of the brain (expansion during sleep and contraction during oxygen deprivation) optimize the amount of energy that neurons can store in, and retrieve from, their ionic concentration gradients.

Electric discharge of electrocytes: Modelling, analysis and simulation

In this paper, we investigate the electric discharge of electrocytes by extending our previous work on the generation of electric potential. We first give a complete formulation of a single cell unit consisting of an electrocyte and a resistor, based on a Poisson-Nernst-Planck (PNP) system with various membrane currents as interfacial conditions for the electrocyte and a Maxwell's model for the resistor. Our previous work can be treated as a special case with an infinite resistor (or open circuit). Using asymptotic analysis, we simplify our PNP system and reduce it to an ordinary differential equation (ODE) based model. Unlike the case of an infinite resistor, our numerical simulations of the new model reveal several distinct features. A finite current is generated, which leads to non-constant electric potentials in the bulk of intracellular and extracellular regions. Furthermore, the current induces an additional action potential (AP) at the non-innervated membrane, contrary to the case of an open circuit where an AP is generated only at the innervated membrane. The voltage drop inside the electrocyte is caused by an internal resistance due to mobile ions. We show that our single cell model can be used as the basis for a system with stacked electrocytes and the total current during the discharge of an electric eel can be estimated by using our model.

Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries

Mathematical models for excitable cells are commonly based on cable theory, which considers a homogenized domain and spatially constant ionic concentrations. Although such models provide valuable insight, the effect of altered ion concentrations or detailed cell morphology on the electrical potentials cannot be captured. In this paper, we discuss an alternative approach to detailed modeling of electrodiffusion in neural tissue. The mathematical model describes the distribution and evolution of ion concentrations in a geometrically-explicit representation of the intra- and extracellular domains. As a combination of the electroneutral Kirchhoff-Nernst-Planck (KNP) model and the Extracellular-Membrane-Intracellular (EMI) framework, we refer to this model as the KNP-EMI model. Here, we introduce and numerically evaluate a new, finite element-based numerical scheme for the KNP-EMI model, capable of efficiently and flexibly handling geometries of arbitrary dimension and arbitrary polynomial degree. Moreover, we compare the electrical potentials predicted by the KNP-EMI and EMI models. Finally, we study ephaptic coupling induced in an unmyelinated axon bundle and demonstrate how the KNP-EMI framework can give new insights in this setting.

PyPNS: Multiscale Simulation of a Peripheral Nerve in Python

Bioelectronic Medicines that modulate the activity patterns on peripheral nerves have promise as a new way of treating diverse medical conditions from epilepsy to rheumatism. Progress in the field builds upon time consuming and expensive experiments in living organisms. To reduce experimentation load and allow for a faster, more detailed analysis of peripheral nerve stimulation and recording, computational models incorporating experimental insights will be of great help. We present a peripheral nerve simulator that combines biophysical axon models and numerically solved and idealised extracellular space models in one environment. We modelled the extracellular space as a three-dimensional resistive continuum governed by the electro-quasistatic approximation of the Maxwell equations. Potential distributions were precomputed in finite element models for different media (homogeneous, nerve in saline, nerve in cuff) and imported into our simulator. Axons, on the other hand, were modelled more abstractly as one-dimensional chains of compartments. Unmyelinated fibres were based on the Hodgkin-Huxley model; for myelinated fibres, we adapted the model proposed by McIntyre et al. in 2002 to smaller diameters. To obtain realistic axon shapes, an iterative algorithm positioned fibres along the nerve with a variable tortuosity fit to imaged trajectories. We validated our model with data from the stimulated rat vagus nerve. Simulation results predicted that tortuosity alters recorded signal shapes and increases stimulation thresholds. The model we developed can easily be adapted to different nerves, and may be of use for Bioelectronic Medicine research in the future.

A Kirchhoff-Nernst-Planck framework for modeling large scale extracellular electrodiffusion surrounding morphologically detailed neurons

Many pathological conditions, such as seizures, stroke, and spreading depression, are associated with substantial changes in ion concentrations in the extracellular space (ECS) of the brain. An understanding of the mechanisms that govern ECS concentration dynamics may be a prerequisite for understanding such pathologies. To estimate the transport of ions due to electrodiffusive effects, one must keep track of both the ion concentrations and the electric potential simultaneously in the relevant regions of the brain. Although this is currently unfeasible experimentally, it is in principle achievable with computational models based on biophysical principles and constraints. Previous computational models of extracellular ion-concentration dynamics have required extensive computing power, and therefore have been limited to either phenomena on very small spatiotemporal scales (micrometers and milliseconds), or simplified and idealized 1-dimensional (1-D) transport processes on a larger scale. Here, we present the 3-D Kirchhoff-Nernst-Planck (KNP) framework, tailored to explore electrodiffusive effects on large spatiotemporal scales. By assuming electroneutrality, the KNP-framework circumvents charge-relaxation processes on the spatiotemporal scales of nanometers and nanoseconds, and makes it feasible to run simulations on the spatiotemporal scales of millimeters and seconds on a standard desktop computer. In the present work, we use the 3-D KNP framework to simulate the dynamics of ion concentrations and the electrical potential surrounding a morphologically detailed pyramidal cell. In addition to elucidating the single neuron contribution to electrodiffusive effects in the ECS, the simulation demonstrates the efficiency of the 3-D KNP framework. We envision that future applications of the framework to more complex and biologically realistic systems will be useful in exploring pathological conditions associated with large concentration variations in the ECS.

The Role of Potassium Channels in Arabidopsis thaliana Long Distance Electrical Signalling: AKT2 Modulates Tissue Excitability While GORK Shapes Action Potentials

Fast responses to an external threat depend on the rapid transmission of signals through a plant. Action potentials (APs) are proposed as such signals. Plant APs share similarities with their animal counterparts; they are proposed to depend on the activity of voltage-gated ion channels. Nonetheless, despite their demonstrated role in (a)biotic stress responses, the identities of the associated voltage-gated channels and transporters remain undefined in higher plants. By demonstrating the role of two potassium-selective channels in Arabidopsis thaliana in AP generation and shaping, we show that the plant AP does depend on similar Kv-like transport systems to those of the animal signal. We demonstrate that the outward-rectifying potassium-selective channel GORK limits the AP amplitude and duration, while the weakly-rectifying channel AKT2 affects membrane excitability. By computational modelling of plant APs, we reveal that the GORK activity not only determines the length of an AP but also the steepness of its rise and the maximal amplitude. Thus, outward-rectifying potassium channels contribute to both the repolarisation phase and the initial depolarisation phase of the signal. Additionally, from modelling considerations we provide indications that plant APs might be accompanied by potassium waves, which prime the excitability of the green cable.

Electrodiffusion phenomena in neuroscience: a neglected companion

The emerging technological revolution in genetically encoded molecular sensors and super-resolution imaging provides neuroscientists with a pass to the real-time nano-world. On this small scale, however, classical principles of electrophysiology do not always apply. This is in large part because the nanoscopic heterogeneities in ionic concentrations and the local electric fields associated with individual ions and their movement can no longer be ignored. Here, we review basic principles of molecular electrodiffusion in the cellular environment of organized brain tissue. We argue that accurate interpretation of physiological observations on the nanoscale requires a better understanding of the underlying electrodiffusion phenomena.

An Evaluation of the Accuracy of Classical Models for Computing the Membrane Potential and Extracellular Potential for Neurons

Two mathematical models are part of the foundation of Computational neurophysiology; (a) the Cable equation is used to compute the membrane potential of neurons, and, (b) volume-conductor theory describes the extracellular potential around neurons. In the standard procedure for computing extracellular potentials, the transmembrane currents are computed by means of (a) and the extracellular potentials are computed using an explicit sum over analytical point-current source solutions as prescribed by volume conductor theory. Both models are extremely useful as they allow huge simplifications of the computational efforts involved in computing extracellular potentials. However, there are more accurate, though computationally very expensive, models available where the potentials inside and outside the neurons are computed simultaneously in a self-consistent scheme. In the present work we explore the accuracy of the classical models (a) and (b) by comparing them to these more accurate schemes. The main assumption of (a) is that the ephaptic current can be ignored in the derivation of the Cable equation. We find, however, for our examples with stylized neurons, that the ephaptic current is comparable in magnitude to other currents involved in the computations, suggesting that it may be significant-at least in parts of the simulation. The magnitude of the error introduced in the membrane potential is several millivolts, and this error also translates into errors in the predicted extracellular potentials. While the error becomes negligible if we assume the extracellular conductivity to be very large, this assumption is, unfortunately, not easy to justify a priori for all situations of interest.

From Maxwell's equations to the theory of current-source density analysis

Despite the widespread use of current‐source density (CSD) analysis of extracellular potential recordings in the brain, the physical mechanisms responsible for the generation of the signal are still debated. While the extracellular potential is thought to be exclusively generated by the transmembrane currents, recent studies suggest that extracellular diffusive, advective and displacement currents—traditionally neglected—may also contribute considerably toward extracellular potential recordings. Here, we first justify the application of the electro‐quasistatic approximation of Maxwell's equations to describe the electromagnetic field of physiological origin. Subsequently, we perform spatial averaging of currents in neural tissue to arrive at the notion of the CSD and derive an equation relating it to the extracellular potential. We show that, in general, the extracellular potential is determined by the CSD of membrane currents as well as the gradients of the putative extracellular diffusion current. The diffusion current can contribute significantly to the extracellular potential at frequencies less than a few Hertz; in which case it must be subtracted to obtain correct CSD estimates. We also show that the advective and displacement currents in the extracellular space are negligible for physiological frequencies while, within cellular membrane, displacement current contributes toward the CSD as a capacitive current. Taken together, these findings elucidate the relationship between electric currents and the extracellular potential in brain tissue and form the necessary foundation for the analysis of extracellular recordings.