Finite element modeling of the neuron-electrode interface

Choosing the right electrode representation for modeling real bioelectronic interfaces: a comprehensive guide

Objective.Producing realistic numerical models of neurostimulation electrodes in contact with the electrolyte and tissue, for use in time-domain finite element method simulations while maintaining a reasonable computational burden remains a challenge. We aim to provide a straightforward experimental-theoretical hybrid approach for common electrode materials (Ti, TiN, ITO, Au, Pt, IrOx) that are relevant to the research field of bioelectronics, along with all the information necessary to replicate our approach in arbitrary geometry for real-life experimental applications.Approach.We used electrochemical impedance spectroscopy (EIS) to extract the electrode parameters in the AC regime under different DC biases. The pulsed electrode response was obtained by fast amperometry (FA) to optimize and verify the previously obtained electrode parameters in a COMSOL Multiphysics model. For optimization of the electrode parameters a constant phase element (CPE) needed to be implemented in time-domain.Main results.We find that the parameters obtained by EIS can be used to accurately simulate pulsed response only close to the electrode open circuit potential, while at other potentials we give corrections to the obtained parameters, based on FA measurements. We also find that for many electrodes (Au, TiN, Pt, and IrOx), it is important to implement a distributed CPE rather than an ideal capacitor for estimating the electrode double-layer capacitance. We outline and provide examples for the novel time-domain implementation of the CPE for finite element method simulations in COMSOL Multiphysics.Significance.An overview of electrode parameters for some common electrode materials can be a valuable and useful tool in numerical bioelectronics models. A provided FEM implementation model can be readily adapted to arbitrary electrode geometries and used for various applications. Finally, the presented methodology for parametrization of electrode materials can be used for any materials of interest which were not covered by this work.

Multiscale simulation analysis of passive and active micro/nanoelectrodes for CMOS-based in vitro neural sensing devices

Neuron and neural network studies are remarkably fostered by novel stimulation and recording systems, which often make use of biochips fabricated with advanced electronic technologies and, notably, micro- and nanoscale complementary metal-oxide semiconductor (CMOS). Models of the transduction mechanisms involved in the sensor and recording of the neuron activity are useful to optimize the sensing device architecture and its coupling to the readout circuits, as well as to interpret the measured data. Starting with an overview of recently published integrated active and passive micro/nanoelectrode sensing devices for in vitro studies fabricated with modern (CMOS-based) micro-nano technology, this paper presents a mixed-mode device-circuit numerical-analytical multiscale and multiphysics simulation methodology to describe the neuron-sensor coupling, suitable to derive useful design guidelines. A few representative structures and coupling conditions are analysed in more detail in terms of the most relevant electrical figures of merit including signal-to-noise ratio. This article is part of the theme issue 'Advanced neurotechnologies: translating innovation for health and well-being'.

Lead field theory provides a powerful tool for designing microelectrode array impedance measurements for biological cell detection and observation

Background

Our aim is to introduce a method to enhance the design process of microelectrode array (MEA) based electric bioimpedance measurement systems for improved detection and viability assessment of living cells and tissues. We propose the application of electromagnetic lead field theory and reciprocity for MEA design and measurement result interpretation. Further, we simulated impedance spectroscopy (IS) with two- and four-electrode setups and a biological cell to illustrate the tool in the assessment of the capabilities of given MEA electrode constellations for detecting cells on or in the vicinity of the microelectrodes.

Results

The results show the power of the lead field theory in electromagnetic simulations of cell–microelectrode systems depicting the fundamental differences of two- and four-electrode IS measurement configurations to detect cells. Accordingly, the use in MEA system design is demonstrated by assessing the differences between the two- and four-electrode IS configurations. Further, our results show how cells affect the lead fields in these MEA system, and how we can utilize the differences of the two- and four-electrode setups in cell detection. The COMSOL simulator model is provided freely in public domain as open source.

Conclusions

Lead field theory can be successfully applied in MEA design for the IS based assessment of biological cells providing the necessary visualization and insight for MEA design. The proposed method is expected to enhance the design and usability of automated cell and tissue manipulation systems required for bioreactors, which are intended for the automated production of cell and tissue grafts for medical purposes. MEA systems are also intended for toxicology to assess the effects of chemicals on living cells. Our results demonstrate that lead field concept is expected to enhance also the development of such methods and devices.

A novel bio-mimicking, planar nano-edge microelectrode enables enhanced long-term neural recording

Our inability to accurately monitor individual neurons and their synaptic activity precludes fundamental understanding of brain function under normal and various pathological conditions. However, recent breakthroughs in micro- and nano-scale fabrication processes have advanced the development of neuro-electronic hybrid technology. Among such devices are three-dimensional and planar electrodes, offering the advantages of either high fidelity or longer-term recordings respectively. Here, we present the next generation of planar microelectrode arrays with “nano-edges” that enable long-term (≥1 month) and high fidelity recordings at a resolution 15 times higher than traditional planar electrodes. This novel technology enables better understanding of brain function and offers a tremendous opportunity towards the development of future bionic hybrids and drug discovery devices.

Effect of planar microelectrode geometry on neuron stimulation: finite element modeling and experimental validation of the efficient electrode shape

BACKGROUND

Microelectrode arrays have been used successfully for neuronal stimulation both in vivo and in vitro. However, in most instances currents required to activate the neurons have been in un-physiological ranges resulting in neuronal damage and cell death. There is a need to develop electrodes which require less stimulation current for neuronal activation with physiologically relevant efficacy and frequencies.

NEW METHOD

The objective of the present study was to examine and compare the stimulation efficiency of different electrode geometries at the resolution of a single neuron. We hypothesized that increasing the electrode perimeter will increase the maximum current density at the edges and enhance stimulation efficiency. To test this postulate, the neuronal stimulation efficacy of common circular electrodes (smallest perimeter) was compared with star (medium perimeter), and spiral (largest perimeter with internal boundaries) electrodes. We explored and compared using both a finite element model and in vitro stimulation of neurons isolated from Lymnaea central ganglia.

RESULTS

Interestingly, both the computational model and the live neuronal stimulation experiments demonstrated that the common circular microelectrode requires less stimulus to activate a cell compared to the other two electrode shapes with the same surface area. Our data further revealed that circular electrodes exhibit the largest sealing resistance, stimulus transfer, and average current density among the three types of electrodes tested.

COMPARISON WITH EXISTING METHODS

Average current density and not the maximum current density at the edges plays an important role in determining the electrode stimulation efficiency.

CONCLUSION

Circular shaped electrodes are more efficient in inducing a change in neuronal membrane potential.

A CMOS IC-based multisite measuring system for stimulation and recording in neural preparations in vitro

In this report, we describe the system integration of a complementary metal oxide semiconductor (CMOS) integrated circuit (IC) chip, capable of both stimulation and recording of neurons or neural tissues, to investigate electrical signal propagation within cellular networks in vitro. The overall system consisted of three major subunits: a 5.0 × 5.0 mm CMOS IC chip, a reconfigurable logic device (field-programmable gate array, FPGA), and a PC. To test the system, microelectrode arrays (MEAs) were used to extracellularly measure the activity of cultured rat cortical neurons and mouse cortical slices. The MEA had 64 bidirectional (stimulation and recording) electrodes. In addition, the CMOS IC chip was equipped with dedicated analog filters, amplification stages, and a stimulation buffer. Signals from the electrodes were sampled at 15.6 kHz with 16-bit resolution. The measured input-referred circuitry noise was 10.1 μ V root mean square (10 Hz to 100 kHz), which allowed reliable detection of neural signals ranging from several millivolts down to approximately 33 μ Vpp. Experiments were performed involving the stimulation of neurons with several spatiotemporal patterns and the recording of the triggered activity. An advantage over current MEAs, as demonstrated by our experiments, includes the ability to stimulate (voltage stimulation, 5-bit resolution) spatiotemporal patterns in arbitrary subsets of electrodes. Furthermore, the fast stimulation reset mechanism allowed us to record neuronal signals from a stimulating electrode around 3 ms after stimulation. We demonstrate that the system can be directly applied to, for example, auditory neural prostheses in conjunction with an acoustic sensor and a sound processing system.

Modeling extracellular electrical stimulation: II. Computational validation and numerical results

The validity of approximate equations describing the membrane potential under extracellular electrical stimulation (Meffin et al 2012 J. Neural Eng. 9 065005) is investigated through finite element analysis in this paper. To this end, the finite element method is used to simulate a cylindrical neurite under extracellular stimulation. Laplace's equations with appropriate boundary conditions are solved numerically in three dimensions and the results are compared to the approximate analytic solutions. Simulation results are in agreement with the approximate analytic expressions for longitudinal and transverse modes of stimulation. The range of validity of the equations describing the membrane potential for different values of stimulation and neurite parameters are presented as well. The results indicate that the analytic approach can be used to model extracellular electrical stimulation for realistic physiological parameters with a high level of accuracy.

Nanoscale Nonlinear dynamic characterization of the neuron-electrode junction

Extracellular recordings from neurons using microelectrode and field effect transistor arrays suffer from many problems including low signal to noise ratio, signal attenuation due to counter-ion diffusion from the bulk extracellular medium and a modification of the shape of the cell-generated potentials due to the presence of a highly dispersive dielectric medium in the cell-electrode cleft. Attempts to date to study the neuron-electrode interface have focused on point or area contact linear-equivalent-circuit models. We present here the results obtained from a 'data-true' nonlinear dynamic characterization of the neuron-electrode junction using Volterra-Wiener modeling. For the characterization, NG108-15 cells were cultured on microelectrode arrays and stimulated with broadband Gaussian white noise under voltage clamp mode. A Volterra-Wiener model was then estimated using the input signal and the extracellular signal recorded on the microelectrode. The existence of the second order Wiener kernel confirmed that the recorded extracellular signal had a nonlinear component. The verification of the estimated model was carried out by employing the intracellular action potential as an input to the Volterra-Wiener model and comparing the predicted extracellular response with the corresponding extracellular signal recorded on the microelectrode. We believe that a 'data-true' Volterra-Wiener model of the neuron-electrode junction shall not only facilitate a direct insight into the physicochemical processes taking place at the interface during signal transduction but will also allow one to evolve strategies for engineering the neuron-electrode interface using surface chemical modification of the microelectrodes.

Patch-clamping of primary cardiac cells with micro-openings in polyimide films

Patch-clamping is a powerful method for investigating the function and regulation of ionic channels. Currently, great efforts are being made to automate this method. As a step towards this goal, the feasibility of patch-clamping primary cells with a microscopic opening in a planar substrate was tested. Using standard microfabrication and ion beam technology, small-diameter openings (2 and 4 microm) were formed in polyimide films (thickness 6.5 microm). Single cells (sheep Purkinje heart cells, Chinese hamster ovary cells) in a suspension were positioned on top of the opening and sucked towards the opening to improve adhesion of the cell to the planar substrate, hence increasing the seal resistance. Voltage/current measurements yielded a median seal resistance of 1.3 Mohms with 4 microm openings (n=24) and 26.0 Mohms with 2 microm openings (n = 75), respectively. With 2 microm openings, successful loose-patch recordings of TTX-sensitive inward currents and action potentials in sheep Purkinje heart cells (n = 18) were made. In rare cases, gigaseals (n = 4) were also measured, and a whole-cell configuration (n = 1) could be established. It was concluded that the simple planar patch approach is suitable for automated loose-patch recordings from cells in suspension but will hardly be suitable for high-throughput whole-cell patch-clamping with high-resistance seals.