Polarization of a Spherical Cell in a Nonuniform Extracellular Electric Field

Restore axonal conductance in a locally demyelinated axon with electromagnetic stimulation

Objective. Axonal demyelination leads to failure of axonal conduction. Current research on demyelination focuses on the promotion of remyelination. Electromagnetic stimulation is widely used to promote neural activity. We hypothesized that electromagnetic stimulation of the demyelinated area, by providing excitation to the nodes of Ranvier, could rescue locally demyelinated axons from conductance failure.Approach. We built a multi-compartment NEURON model of a myelinated axon under electromagnetic stimulation. We simulated the action potential (AP) propagation and observed conductance failure when local demyelination occurred. Conductance failure was due to current leakage and a lack of activation of the nodes in the demyelinated region. To investigate the effects of electromagnetic stimulation on locally demyelinated axons, we positioned a miniature coil next to the affected area to activate nodes in the demyelinated region.Main results. Subthreshold microcoil stimulation caused depolarization of node membranes. This depolarization, in combination with membrane depolarization induced by the invading AP, resulted in sufficient activation of nodes in the demyelinated region and restoration of axonal conductance. Efficacy of restoration was dependent on the amplitude and frequency of the stimuli, and the location of the microcoil relative to the targeted nodes. The restored axonal conductance was due to the enhanced Na+current and reduced K+current in the nodes, rather than a reduction in leakage current in the demyelinated region. Finally, we found that microcoil stimulation had no effect on axonal conductance in healthy, myelinated axons.Significance. Activation of nodes in the demyelinated region using electromagnetic stimulation provides an alternative treatment strategy to restore axonal function under local demyelination conditions. Results provide insights to the development of microcoil technology for the treatment of focal segmental demyelination cases, such as neuropraxia, spinal cord injury, and auditory nerve demyelination.

Quasistatic approximation in neuromodulation

We define and explain the quasistatic approximation (QSA) as applied to field modeling for electrical and magnetic stimulation. Neuromodulation analysis pipelines include discrete stages, and QSA is applied specifically when calculating the electric and magnetic fields generated in tissues by a given stimulation dose. QSA simplifies the modeling equations to support tractable analysis, enhanced understanding, and computational efficiency. The application of QSA in neuromodulation is based on four underlying assumptions: (A1) no wave propagation or self-induction in tissue, (A2) linear tissue properties, (A3) purely resistive tissue, and (A4) non-dispersive tissue. As a consequence of these assumptions, each tissue is assigned a fixed conductivity, and the simplified equations (e.g. Laplace's equation) are solved for the spatial distribution of the field, which is separated from the field's temporal waveform. Recognizing that electrical tissue properties may be more complex, we explain how QSA can be embedded in parallel or iterative pipelines to model frequency dependence or nonlinearity of conductivity. We survey the history and validity of QSA across specific applications, such as microstimulation, deep brain stimulation, spinal cord stimulation, transcranial electrical stimulation, and transcranial magnetic stimulation. The precise definition and explanation of QSA in neuromodulation are essential for rigor when using QSA models or testing their limits.

State of the Art and Current Challenges on Electroactive Biomaterials and Strategies for Neural Tissue Regeneration

The loss or failure of an organ/tissue stands as one of the healthcare system's most prevalent, devastating, and costly challenges. Strategies for neural tissue repair and regeneration have received significant attention due to their particularly strong impact on patients' well-being. Many research efforts are dedicated not only to control the disease symptoms but also to find solutions to repair the damaged tissues. Neural tissue engineering (TE) plays a key role in addressing this problem and significant efforts are being carried out to develop strategies for neural repair treatment. In the last years, active materials allowing to tune cell-materials interaction are being increasingly used, representing a recent paradigm in TE applications. Among the most important stimuli influencing cell behavior are the electrical and mechanical ones. In this way, materials with the ability to provide this kind of stimuli to the neural cells seem to be appropriate to support neural TE. In this scope, this review summarizes the different biomaterials types used for neural TE, highlighting the relevance of using active biomaterials and electrical stimulation. Furthermore, this review provides not only a compilation of the most relevant studies and results but also strategies for novel and more biomimetic approaches for neural TE.

Towards a more accurate quasi-static approximation of the electric potential for neurostimulation with kilohertz-frequency sources. J Neural Eng 2023;20:10.1088/1741-2552/ad1612. [PMID: 38100821 PMCID: PMC10822676 DOI: 10.1088/1741-2552/ad1612]

Objective.Our goal was to determine the conditions for which a more precise calculation of the electric potential than the quasi-static approximation may be needed in models of electrical neurostimulation, particularly for signals with kilohertz-frequency components.Approach.We conducted a comprehensive quantitative study of the differences in nerve fiber activation and conduction block when using the quasi-static and Helmholtz approximations for the electric potential in a model of electrical neurostimulation.Main results.We first show that the potentials generated by sources of unbalanced pulses exhibit different transients as compared to those of charge-balanced pulses, and this is disregarded by the quasi-static assumption. Secondly, relative errors for current-distance curves were below 3%, while for strength-duration curves these ranged between 1%-17%, but could be improved to less than 3% across the range of pulse duration by providing a corrected quasi-static conductivity. Third, we extended our analysis to trains of pulses and reported a 'congruence area' below 700 Hz, where the fidelity of fiber responses is maximal for supra-threshold stimulation. Further examination of waveforms and polarities revealed similar fidelities in the congruence area, but significant differences were observed beyond this area. However, the spike-train distance revealed differences in activation patterns when comparing the response generated by each model. Finally, in simulations of conduction-block, we found that block thresholds exhibited errors above 20% for repetition rates above 10 kHz. Yet, employing a corrected value of the conductivity improved the agreement between models, with errors no greater than 8%.Significance.Our results emphasize that the quasi-static approximation cannot be naively extended to electrical stimulation with high-frequency components, and notable differences can be observed in activation patterns. As well, we introduce a methodology to obtain more precise model responses using the quasi-static approach, retaining its simplicity, which can be a valuable resource in computational neuroengineering.

Neuron matters: neuromodulation with electromagnetic stimulation must consider neurons as dynamic identities

Neuromodulation with electromagnetic stimulation is widely used for the control of abnormal neural activity, and has been proven to be a valuable alternative to pharmacological tools for the treatment of many neurological diseases. Tremendous efforts have been focused on the design of the stimulation apparatus (i.e., electrodes and magnetic coils) that delivers the electric current to the neural tissue, and the optimization of the stimulation parameters. Less attention has been given to the complicated, dynamic properties of the neurons, and their context-dependent impact on the stimulation effects. This review focuses on the neuronal factors that influence the outcomes of electromagnetic stimulation in neuromodulation. Evidence from multiple levels (tissue, cellular, and single ion channel) are reviewed. Properties of the neural elements and their dynamic changes play a significant role in the outcome of electromagnetic stimulation. This angle of understanding yields a comprehensive perspective of neural activity during electrical neuromodulation, and provides insights in the design and development of novel stimulation technology.

Finding the Location of Axonal Activation by a Miniature Magnetic Coil

Magnetic stimulation for neural activation is widely used in clinical and lab research. In comparison to electric stimulation using an implanted electrode, stimulation with a large magnetic coil is associated with poor spatial specificity and incapability to stimulate deep brain structures. Recent developments in micromagnetic stimulation (μMS) technology mitigates some of these shortcomings. The sub-millimeter coils can be covered with soft, biocompatible material, and chronically implanted. They can provide highly specific neural stimulation in the deep neural structure. Although the μMS technology is expected to provide a precise location of neural stimulation, the exact site of neural activation is difficult to determine. Furthermore, factors that could cause the shifting of the activation site during μMS have not been fully investigated. To estimate the location of axon activation in μMS, we first derived an analytical expression of the activating function, which predicts the location of membrane depolarization in an unmyelinated axon. Then, we developed a multi-compartment, Hodgkin-Huxley (H-H) type of NEURON model of an unmyelinated axon to test the impact of several important coil parameters on the location of axonal activation. The location of axonal activation was dependent on both the parameters of the stimulus and the biophysics properties of the targeted axon during μMS. The activating function analysis predicted that the location of membrane depolarization and activation could shift due to the reversal of the coil current and the change in the coil-axon distance. The NEURON modeling confirmed these predictions. Interestingly, the NEURON simulation further revealed that the intensity of stimulation played a significant role in the activation location. Moderate or strong coil currents activated the axon at different locations, mediated by two distinct ion channel mechanisms. This study reports several experimental factors that could cause a potential shift in the location of neural activation during μMS, which is essential for further development of this novel technology.

Electro-physiology Models of Cells with Spherical Geometry with Non-conducting Center

We study the flow of electrical currents in spherical cells with a non-conducting core, so that current flow is restricted to a thin shell below the cell's membrane. Examples of such cells are fat storing cells (adipocytes). We derive the relation between current and voltage in the passive regime and examine the conditions under which the cell is electro-tonically compact. We compare our results to the well-studied case of electrical current flow in cylinder structures, such as neurons, described by the cable equation. In contrast to the cable, we find that for the sphere geometry (1) the voltage profile across the cell depends critically on the electrode geometry, and (2) the charging and discharging can be much faster than the membrane time constant; however, (3) voltage clamp experiments will incur similar distortion as in the cable case. We discuss the relevance for adipocyte function and experimental electro-physiology.

Modified cable equation incorporating transverse polarization of neuronal membranes for accurate coupling of electric fields

OBJECTIVE

We present a theory and computational methods to incorporate transverse polarization of neuronal membranes into the cable equation to account for the secondary electric field generated by the membrane in response to transverse electric fields. The effect of transverse polarization on nonlinear neuronal activation thresholds is quantified and discussed in the context of previous studies using linear membrane models.

APPROACH

The response of neuronal membranes to applied electric fields is derived under two time scales and a unified solution of transverse polarization is given for spherical and cylindrical cell geometries. The solution is incorporated into the cable equation re-derived using an asymptotic model that separates the longitudinal and transverse dimensions. Two numerical methods are proposed to implement the modified cable equation. Several common neural stimulation scenarios are tested using two nonlinear membrane models to compare thresholds of the conventional and modified cable equations.

MAIN RESULTS

The implementations of the modified cable equation incorporating transverse polarization are validated against previous results in the literature. The test cases show that transverse polarization has limited effect on activation thresholds. The transverse field only affects thresholds of unmyelinated axons for short pulses and in low-gradient field distributions, whereas myelinated axons are mostly unaffected.

SIGNIFICANCE

The modified cable equation captures the membrane's behavior on different time scales and models more accurately the coupling between electric fields and neurons. It addresses the limitations of the conventional cable equation and allows sound theoretical interpretations. The implementation provides simple methods that are compatible with current simulation approaches to study the effect of transverse polarization on nonlinear membranes. The minimal influence by transverse polarization on axonal activation thresholds for the nonlinear membrane models indicates that predictions of stronger effects in linear membrane models with a fixed activation threshold are inaccurate. Thus, the conventional cable equation works well for most neuroengineering applications, and the presented modeling approach is well suited to address the exceptions.

Analysis of upper threshold mechanisms of spherical neurons during extracellular stimulation

Exceeding a certain stimulation strength can prevent the generation of somatic action potentials, as has been demonstrated in vitro with extracellularly stimulated dorsal root ganglion cells as well as retinal ganglion cells. This phenomenon, termed upper threshold, is currently thought to be a consequence of sodium current reversal in strongly depolarized regions. Here we analyze the contribution of membrane kinetics, using spherical model neurons that are stimulated externally with a microelectrode, in more detail. During extracellular pulse application, the electric field depolarizes one part and hyperpolarizes the other part of the cell. Strong transmembrane currents are generated only in the active depolarized region, changing the overall polarization level. The asymmetric membrane voltage distribution caused by the stimulus strongly influences the cell's behavior during and even after the stimulus. Effects on membrane voltage and transmembrane currents during and after the stimulus are shown and discussed in detail. Aside from the sodium current reversal, two more key mechanisms were identified in causing the upper threshold: strong potassium currents and inactivation of sodium channels. The contributions of the mechanisms involved strongly depend on cell properties, stimulus parameters, and other factors such as temperature. The conclusions presented here are based on several retinal ganglion cell models of the Fohlmeister group, a model with original Hodgkin-Huxley membrane, and a pyramidal cell model. NEW & NOTEWORTHY The upper threshold phenomenon in extracellular stimulation is analyzed in detail for spherical cells. Three main mechanisms were identified that prevent the generation of action potentials at high stimulation strengths: 1) strong potassium currents, 2) inactivating sodium ion channels, and 3) sodium current reversal. Ion channel kinetics in retinal ganglion cells, pyramidal cells, and the original Hodgkin-Huxley model were investigated under the influence of an extracellular stimulus.

Shielding effects of myelin sheath on axolemma depolarization under transverse electric field stimulation

Axonal stimulation with electric currents is an effective method for controlling neural activity. An electric field parallel to the axon is widely accepted as the predominant component in the activation of an axon. However, recent studies indicate that the transverse component to the axolemma is also effective in depolarizing the axon. To quantitatively investigate the amount of axolemma polarization induced by a transverse electric field, we computed the transmembrane potential (Vm) for a conductive body that represents an unmyelinated axon (or the bare axon between the myelin sheath in a myelinated axon). We also computed the transmembrane potential of the sheath-covered axonal segment in a myelinated axon. We then systematically analyzed the biophysical factors that affect axonal polarization under transverse electric stimulation for both the bare and sheath-covered axons. Geometrical patterns of polarization of both axon types were dependent on field properties (magnitude and field orientation to the axon). Polarization of both axons was also dependent on their axolemma radii and electrical conductivities. The myelin provided a significant “shielding effect” against the transverse electric fields, preventing excessive axolemma depolarization. Demyelination could allow for prominent axolemma depolarization in the transverse electric field, via a significant increase in myelin conductivity. This shifts the voltage drop of the myelin sheath to the axolemma. Pathological changes at a cellular level should be considered when electric fields are used for the treatment of demyelination diseases. The calculated term for membrane polarization (Vm) could be used to modify the current cable equation that describes axon excitation by an external electric field to account for the activating effects of both parallel and transverse fields surrounding the target axon.

Biophysically realistic neuron models for simulation of cortical stimulation

OBJECTIVE

We implemented computational models of human and rat cortical neurons for simulating the neural response to cortical stimulation with electromagnetic fields.

APPROACH

We adapted model neurons from the library of Blue Brain models to reflect biophysical and geometric properties of both adult rat and human cortical neurons and coupled the model neurons to exogenous electric fields (E-fields). The models included 3D reconstructed axonal and dendritic arbors, experimentally-validated electrophysiological behaviors, and multiple, morphological variants within cell types. Using these models, we characterized the single-cell responses to intracortical microstimulation (ICMS) and uniform E-field with dc as well as pulsed currents.

MAIN RESULTS

The strength-duration and current-distance characteristics of the model neurons to ICMS agreed with published experimental results, as did the subthreshold polarization of cell bodies and axon terminals by uniform dc E-fields. For all forms of stimulation, the lowest threshold elements were terminals of the axon collaterals, and the dependence of threshold and polarization on spatial and temporal stimulation parameters was strongly affected by morphological features of the axonal arbor, including myelination, diameter, and branching.

SIGNIFICANCE

These results provide key insights into the mechanisms of cortical stimulation. The presented models can be used to study various cortical stimulation modalities while incorporating detailed spatial and temporal features of the applied E-field.

Mechanic stress generated by a time-varying electromagnetic field on bone surface

Bone cells sense mechanical load, which is essential for bone growth and remodeling. In a fracture, this mechanism is compromised. Electromagnetic stimulation has been widely used to assist in bone healing, but the underlying mechanisms are largely unknown. A recent hypothesis suggests that electromagnetic stimulation could influence tissue biomechanics; however, a detailed quantitative understanding of EM-induced biomechanical changes in the bone is unavailable. This paper used a muscle/bone model to study the biomechanics of the bone under EM exposure. Due to the dielectric properties of the muscle/bone interface, a time-varying magnetic field can generate both compressing and shear stresses on the bone surface, where many mechanical sensing cells are available for cellular mechanotransduction. I calculated these stresses and found that the shear stress is significantly greater than the compressing stress. Detailed parametric analysis suggests that both the compressing and shear stresses are dependent on the geometrical and electrical properties of the muscle and the bone. These stresses are also functions of the orientation of the coil and the frequency of the magnetic field. It is speculated that the EM field could apply biomechanical influence to fractured bone, through the fine-tuning of the controllable field parameters. Graphical abstract Mechanic stress on bone surface in a time-varying magnetic field.

Possible Effects of Electric Fields on a Pair of Spherical Cells

Electric fields (EF) can induce some physiological or biological effects in neural tissues, which have been explored in many applications such as electroporation. The key to understand the possible underlying mechanisms of such effects tend to be the induced transmembrane potential. Although transmembrane potentials have already been the subject of many theoretical studies, most previous works concerning this topic are mainly focused on the situations of isolated cells. In previous studies, cells are often considered to be three-compartment models with different electroconductivities in different regions (the three compartments are often intracellular regions, membrane, and extracellular regions). In the present paper, we utilize a finite element method (FSM) (with the help of COMSOL^®) to calculate the induced transmembrane potential by the applied EF for a model of two neurons, which may have significant difference on electroporation.

A new theoretical model for transmembrane potential and ion currents induced in a spherical cell under low frequency electromagnetic field

Time-varying electromagnetic fields (EMF) can induce some physiological effects in neuronal tissues, which have been explored in many applications such as transcranial magnetic stimulation. Although transmembrane potentials and induced currents have already been the subjects of many theoretical studies, most previous works about this topic are mainly completed by utilizing Maxwell's equations, often by solving a Laplace equation. In previous studies, cells were often considered to be three-compartment models with different electroconductivities in different regions (three compartments are often intracellular regions, membrane, and extracellular regions). However, models like that did not take dynamic ion channels into consideration. Therefore, one cannot obtain concrete ionic current changes such as potassium current change or sodium current change by these models. The aim of the present work is to present a new and more detailed model for calculating transmembrane potentials and ionic currents induced by time-varying EMF. Equations used in the present paper originate from Nernst-Plank equations, which are ionic current-related equations. The main work is to calculate ionic current changes induced by EMF exposure, and then transmembrane potential changes are calculated with Hodgkin-Huxley model. Bioelectromagnetics. 37:481-492, 2016. © 2016 Wiley Periodicals, Inc.

Neuron matters: electric activation of neuronal tissue is dependent on the interaction between the neuron and the electric field

In laboratory research and clinical practice, externally-applied electric fields have been widely used to control neuronal activity. It is generally accepted that neuronal excitability is controlled by electric current that depolarizes or hyperpolarizes the excitable cell membrane. What determines the amount of polarization? Research on the mechanisms of electric stimulation focus on the optimal control of the field properties (frequency, amplitude, and direction of the electric currents) to improve stimulation outcomes. Emerging evidence from modeling and experimental studies support the existence of interactions between the targeted neurons and the externally-applied electric fields. With cell-field interaction, we suggest a two-way process. When a neuron is positioned inside an electric field, the electric field will induce a change in the resting membrane potential by superimposing an electrically-induced transmembrane potential (ITP). At the same time, the electric field can be perturbed and re-distributed by the cell. This cell-field interaction may play a significant role in the overall effects of stimulation. The redistributed field can cause secondary effects to neighboring cells by altering their geometrical pattern and amount of membrane polarization. Neurons excited by the externally-applied electric field can also affect neighboring cells by ephaptic interaction. Both aspects of the cell-field interaction depend on the biophysical properties of the neuronal tissue, including geometric (i.e., size, shape, orientation to the field) and electric (i.e., conductivity and dielectricity) attributes of the cells. The biophysical basis of the cell-field interaction can be explained by the electromagnetism theory. Further experimental and simulation studies on electric stimulation of neuronal tissue should consider the prospect of a cell-field interaction, and a better understanding of tissue inhomogeneity and anisotropy is needed to fully appreciate the neural basis of cell-field interaction as well as the biological effects of electric stimulation.

Transmembrane potential generated by a magnetically induced transverse electric field in a cylindrical axonal model

During the electrical stimulation of a uniform, long, and straight nerve axon, the electric field oriented parallel to the axon has been widely accepted as the major field component that activates the axon. Recent experimental evidence has shown that the electric field oriented transverse to the axon is also sufficient to activate the axon, by inducing a transmembrane potential within the axon. The transverse field can be generated by a time-varying magnetic field via electromagnetic induction. The aim of this study was to investigate the factors that influence the transmembrane potential induced by a transverse field during magnetic stimulation. Using an unmyelinated axon model, we have provided an analytic expression for the transmembrane potential under spatially uniform, time-varying magnetic stimulation. Polarization of the axon was dependent on the properties of the magnetic field (i.e., orientation to the axon, magnitude, and frequency). Polarization of the axon was also dependent on its own geometrical (i.e., radius of the axon and thickness of the membrane) and electrical properties (i.e., conductivities and dielectric permittivities). Therefore, this article provides evidence that aside from optimal coil design, tissue properties may also play an important role in determining the efficacy of axonal activation under magnetic stimulation. The mathematical basis of this conclusion was discussed. The analytic solution can potentially be used to modify the activation function in current cable equations describing magnetic stimulation.

Transmembrane potential induced on the internal organelle by a time-varying magnetic field: a model study

BACKGROUND

When a cell is exposed to a time-varying magnetic field, this leads to an induced voltage on the cytoplasmic membrane, as well as on the membranes of the internal organelles, such as mitochondria. These potential changes in the organelles could have a significant impact on their functionality. However, a quantitative analysis on the magnetically-induced membrane potential on the internal organelles has not been performed.

METHODS

Using a two-shell model, we provided the first analytical solution for the transmembrane potential in the organelle membrane induced by a time-varying magnetic field. We then analyzed factors that impact on the polarization of the organelle, including the frequency of the magnetic field, the presence of the outer cytoplasmic membrane, and electrical and geometrical parameters of the cytoplasmic membrane and the organelle membrane.

RESULTS

The amount of polarization in the organelle was less than its counterpart in the cytoplasmic membrane. This was largely due to the presence of the cell membrane, which "shielded" the internal organelle from excessive polarization by the field. Organelle polarization was largely dependent on the frequency of the magnetic field, and its polarization was not significant under the low frequency band used for transcranial magnetic stimulation (TMS). Both the properties of the cytoplasmic and the organelle membranes affect the polarization of the internal organelle in a frequency-dependent manner.

CONCLUSIONS

The work provided a theoretical framework and insights into factors affecting mitochondrial function under time-varying magnetic stimulation, and provided evidence that TMS does not affect normal mitochondrial functionality by altering its membrane potential.

Modeling needle stimulation of denervated muscle fibers: voltage-distance relations and fiber polarization effects

A finite-element model of the human thigh was coupled with a 1-D compartment model to simulate the excitation of denervated muscle fibers with a needle electrode. For short electrode-fiber distances, the specific characteristics of the needle geometry determined the areas of lowest threshold values. With increasing distance, these areas shifted toward the needle's center of charge. Comparison of the 1-D model with a 3-D fiber model showed that the assumption of rotational symmetry underlying the 1-D model leads to an overestimation of thresholds. For a 40- micro-diameter fiber stimulated with 50 micross pulses at electrode-fiber distances between 50 microm and 1 mm, the 1-D/3-D threshold ratios were between 1.14 and 1.35 for the muscle fiber model, and between 1.11 and 1.17 for Hodgkin-Huxley membrane properties at 20 degrees C. For both membrane models, the deviation was more pronounced for large fiber diameters and short stimulation pulses. Qualitative results of the 1-D model like voltage-distance relations and predictions of spike initiation sites were correct.

Selective extracellular stimulation of individual neurons in ganglia

Selective control of individual neurons could clarify neural functions and aid disease treatments. To target specific neurons, it may be useful to focus on ganglionic neuron clusters, which are found in the peripheral nervous system in vertebrates. Because neuron cell bodies are found primarily near the surface of invertebrate ganglia, and often found near the surface of vertebrate ganglia, we developed a technique for controlling individual neurons extracellularly using the buccal ganglia of the marine mollusc Aplysia californica as a model system. We experimentally demonstrated that anodic currents can selectively activate an individual neuron and cathodic currents can selectively inhibit an individual neuron using this technique. To define spatial specificity, we studied the minimum currents required for stimulation, and to define temporal specificity, we controlled firing frequencies up to 45 Hz. To understand the mechanisms of spatial and temporal specificity, we created models using the NEURON software package. To broadly predict the spatial specificity of arbitrary neurons in any ganglion sharing similar geometry, we created a steady-state analytical model. A NEURON model based on cat spinal motor neurons showed responses to extracellular stimulation qualitatively similar to those of the Aplysia NEURON model, suggesting that this technique could be widely applicable to vertebrate and human peripheral ganglia having similar geometry.

Transmembrane potential induced in a spherical cell model under low-frequency magnetic stimulation

Time-varying magnetic fields can induce electric fields in the neuronal tissue, a phenomenon that has been recently explored in clinical applications such as peripheral nerve stimulation and transcranial magnetic stimulation. Although the transmembrane potential induced during direct electric stimulation has already been the subject of a number of theoretical studies, an analytical solution for the magnetically induced transmembrane potential change is still unavailable. In addition, although several studies have analyzed the impact of stimulation parameters, including stimulation intensity and frequency, as well as coil design and position, on the amount of tissue polarization, the effects of tissue non-homogeneity on cell polarization have not been fully elucidated. In this study, we have derived an analytical expression for the transmembrane potential induced by a low-frequency magnetic field in a spherical neuronal structure. This model is representative of a spherical cell body or any neuronal structure of a similar shape. The model cell is located in an extracellular medium and possesses a low-conductive membrane and an internal cytoplasm. These three regions represent the basic tissue non-homogeneity of a neuron at a microscopic level. The sensitivity of the induced transmembrane potential to the coil position and to the geometrical and electrical parameters of the model structure was studied in a broad physiologically relevant range. Our results demonstrate that the structure is regionally polarized, with the pattern of polarization depending on the relative positioning between the model cell and the stimulation coil. In addition, both the geometrical and electrical parameters of the structure affect the amount of polarization. These results may be generalized to other neuronal tissues that possess similar non-homogenous properties, but different shapes, such as an axon. Our results support the idea that aside from coil design and position, tissue non-homogeneity could play an important role in determining the effects of magnetic stimulation.

Hybrid finite element method for describing the electrical response of biological cells to applied fields

A novel hybrid finite element method (FEM) for modeling the response of passive and active biological membranes to external stimuli is presented. The method is based on the differential equations that describe the conservation of electric flux and membrane currents. By introducing the electric flux through the cell membrane as an additional variable, the algorithm decouples the linear partial differential equation part from the nonlinear ordinary differential equation part that defines the membrane dynamics of interest. This conveniently results in two subproblems: a linear interface problem and a nonlinear initial value problem. The linear interface problem is solved with a hybrid FEM. The initial value problem is integrated by a standard ordinary differential equation solver such as the Euler and Runge-Kutta methods. During time integration, these two subproblems are solved alternatively. The algorithm can be used to model the interaction of stimuli with multiple cells of almost arbitrary geometries and complex ion-channel gating at the plasma membrane. Numerical experiments are presented demonstrating the uses of the method for modeling field stimulation and action potential propagation.

Combining stimulus direction and waveform for optimization of threshold stimulation of isolated ventricular myocytes

Electric field stimulation is widely used for heart pacing and arrhythmia reversion. In this study, we analysed the influence of waveform and direction of external stimulating electric field on the excitation threshold of isolated ventricular myocytes. The threshold field (E(T)) was lower when the field was applied longitudinally (E(T,L)) rather than transversally (E(T,T)) to the cell major axis. Rheobase was greater for transversal stimulation, but chronaxie and estimated membrane polarization were similar for both directions. The calculated maximal variation in membrane potential at the threshold (DeltaV(T) approximately 15 mV) was insensitive to field direction. As DeltaV(T) values were similar, we assumed that the E(T,T)/E(T,L) ratio might be described solely as the ratio of the major and minor cell semi-axes. Accordingly, the ratio thus estimated was comparable to that determined experimentally. Stimulus waveform significantly affected both E(T) and DeltaV(T), which were greater for monophasic versus biphasic stimuli. Direction and waveform effects were independent. We conclude that (a) direction affects E(T) by its influence on the ability of a given field intensity to cause threshold membrane polarization and (b) threshold-lowering effects of longitudinal stimulation and biphasic waveforms apparently depend on different mechanisms, are additive and thus may be combined to decrease the energy requirement for myocardial stimulation.

Numerical Determination of Transmembrane Voltage Induced on Irregularly Shaped Cells

The paper presents an approach that reduces several difficulties related to the determination of induced transmembrane voltage (ITV) on irregularly shaped cells. We first describe a method for constructing realistic models of irregularly shaped cells based on microscopic imaging. This provides a possibility to determine the ITV on the same cells on which an experiment is carried out, and can be of considerable importance in understanding and interpretation of the data. We also show how the finite-thickness, nonzero-conductivity membrane can be replaced by a boundary condition in which a specific surface conductivity is assigned to the interface between the cell interior (the cytoplasm) and the exterior. We verify the results obtained using this method by a comparison with the analytical solution for an isolated spherical cell and a tilted oblate spheroidal cell, obtaining a very good agreement in both cases. In addition, we compare the ITV computed for a model of two irregularly shaped CHO cells with the ITV measured on the same two cells by means of a potentiometric fluorescent dye, and also with the ITV computed for a simplified model of these two cells.

Field stimulation of cells in suspension: use of a hybrid finite element method

Electric fields are used in a range of applications, including gene transfection, electrochemotherapy of tumors and cardiac defibrillation. Despite the widespread use of electric fields, most of the theoretical and computational studies on discrete cellular tissue have focused on a single cell. In this work, we propose a hybrid finite element method to simulate the effects of external electric fields on clusters of excitable cells. The method can be used to model cells of arbitrary cell geometries and non-linear membrane dynamics. The results show that the response of multiple cell, like a single cell, is a two-stage process consisting of the initial polarization that proceeds with cellular time constant (less than one microsecond) and the actual excitation of the cell membrane that proceeds with the membrane time constant (on the order of milliseconds). The results also show that the stimulation of a given cell depends in part on the arrangement of cells within the field and not simply the location within the field, suggesting that classical approaches that ignores the effect of the cells on the field do not adequately predict the cellular response.