Open Access

Relating membrane potential to impedance spectroscopy


Cite

Introduction

Classification of the dispersions exhibited by impedance spectra of biosystems over a wide frequency range (α, β, γ) was introduced by Schwan [1] to characterize the electric properties of biomaterials. Alpha dispersion is revealed by bioimpedance measurements in the range below 1 Hz and up to 100 kHz. Since mechanisms behind dispersions are often unclear or unknown, the dispersion grouping has been simply based upon a defined frequency range [2].

Even though Gerhard Schwarz seminal paper on low frequency, α dispersion of colloidal particles (related to displacement of counter-ions) was published five decades ago [3], assessment of living cell properties pertaining to α dispersion is still an open issue. Remarkably, such dispersion could basically be explained by various mechanisms formally described by different microscopic models of cell systems, as the ones pioneered by us based either on shape effects (e.g. exhibited by clusters of interconnected cells [4, 5, 6]), or on displacement of counter-ions [7, 8]. Notably, studies relying on the same approaches (and deriving similar results) continue to be published [9, 10]. Moreover, a strong link to the alpha dispersion and counter-ion concepts confirms the classical theories, described for instance in the review article by Foster and Schwan [11].

Theoretical advancement providing a unitary description of both α and β dispersions related to the dielectric behavior of spherical living cells, with emphasis on the relationship between membrane potential and α dispersion was proposed by us in the early nineties [7, 8] However, to our knowledge, experimental non-invasive assessment of membrane potential using impedimetric/dielectric assays has been not reported, as yet.

In this paper we pinpoint the limitations/constraints and suggest solutions enabling measurement of α dispersion on suspensions of living cells, or liposomes consisting of charged lipids and derivation of the related membrane potential.

We briefly introduce the theoretical model, provide the expression for suspension complex permittivity on the parameters of the microscopic model and illustrate the effect on the permittivity and impedance spectra of: I. Parameters pertaining to cell membrane i.e. (i) membrane potential (through the amount of the net charge on the inner side of the membrane), (ii) size of the cells/vesicles, (iii) conductivity of the membrane; II. Parameters of the extra cellular medium: (1) diffusion coefficient and (2) conductivity.

Within the last section of this paper, we emphasize ways and means to check the model as well as a set of parameters pertaining both to cells or vesicles under investigations and to the suspending medium fostering measurable variations of low frequency, α dispersions, exhibited by bioimpedance spectra.

Theory

Aiming to describe the dielectric behavior of a suspension of spherical living cells, we analyze the effect of a uniform ac field of angular frequency ω (ω = 2 π ν) and amplitude E0 on a suspension of shelled spherical particles (of radius R1) with a distribution of (negative) charge, e0 n01, on the inner side of the shell (membrane, of thickness d) [7, 8]; e0 represents the charge of the electron. Diffusion effects are assumed both in the inner (denoted by 1) and external/ suspending medium (denoted by 3), but not in the membrane (denoted by 2). Considering σk as conductivity, and εk as permittivity within the medium k (k=1,2,3) the corresponding complex permittivity and conductivity is given by:

ϵk*=σkiω+εk;σk*=iωϵk*$$\epsilon _{k}^{*}=\frac{{{\sigma }_{k}}}{i\omega }+{{\varepsilon }_{k}};\sigma _{k}^{*}=i\omega \cdot \epsilon _{k}^{*}$$

The external medium is considered as a charge reservoir enabling positive ions (counter-ions) to migrate to the external side of the membrane and balance the negative charge distribution from the inner side (therefore, the total net charge is zero).

By solving the set of equations for the potential (Poisson equation for media 1 and 3, Laplace equation for the shell and the regularity conditions), one obtains the potential in the outer medium (the origin is considered in the centre of a prototype particle) [7]:

φ(r,θ)=E0rcosθ+E0Hnr2cosθ$$\varphi \left( r,\theta \right)=-{{E}_{0}}\cdot r\cos \theta +\frac{{{E}_{0}}\cdot {{H}_{n}}}{{{r}^{2}}}\cos \theta $$

Hn comprises the entire information of the electrical and geometric parameters of the shelled particles, as revealed by expressions within Eq. 4. Notably, 3×Hn represents the polarizability of the (shelled) particle. Following the Maxwell-Wagner equivalence approach, the complex permittivity of suspension is given by [7, 8]:

εsus*=ε3*(1+p3Hn(R1+d)31pHn(R1+d)3)$$\varepsilon _{sus}^{*}=\varepsilon _{3}^{*}\left( 1+p\cdot \frac{3\cdot \frac{{{H}_{n}}}{{{\left( {{R}_{1}}+d \right)}^{3}}}}{1-p\cdot \frac{{{H}_{n}}}{{{\left( {{R}_{1}}+d \right)}^{3}}}} \right)$$

Where:

Hn=ghσ2+iωε2ε01Q+iωε3ε0Qσ2σ3+iωe0NNR1+d1Qσ3/ghσ2+iωε2ε01+2Q2iωε3ε0Qσ2σ3+iωe0NNR1+d1+2Q+2σ3R1+d3NN=e0n01kBTR12R1+d21+iωR1+d22D3g=1+2δσ1σ1+21δiωσ1ε2ε01+s+σ2σ1iωε1ε0Sδ=R1R1+d3S=σG1R1TanhG1R1iωε1ε0G1R12TanhG1R12G1R1TanhG1R1TanhG1R11S=σ1G1R11iωε1ε0G1R122G1R11Q=σ3iωε3ε01+G3R1+dG32R1+d3+21+G3R1+dGk=σkDkεkε0;k=1,3$$\begin{align}& {{H}_{n}}=\left[ \frac{g}{h}\cdot \left( {{\sigma }_{2}}+i\ \omega \cdot {{\varepsilon }_{2}}{{\varepsilon }_{0}}\left( 1-Q \right)+i\,\omega \cdot {{\varepsilon }_{3}}{{\varepsilon }_{0}}\cdot Q\cdot \frac{{{\sigma }_{2}}}{{{\sigma }_{3}}} \right) \right. \\ & \left. +i\,\omega \cdot {{e}_{0}}\cdot NN\cdot \left( {{R}_{1}}+d \right)\cdot \left( 1-Q \right)-\sigma _{3}^{*} \right]/ \\ & \left[ \frac{g}{h}\cdot \left( {{\sigma }_{2}}+i\,\omega \cdot {{\varepsilon }_{2}}{{\varepsilon }_{0}}\left( 1+2Q \right)-2\,i\,\omega \cdot {{\varepsilon }_{3}}{{\varepsilon }_{0}}\cdot Q\cdot \frac{{{\sigma }_{2}}}{{{\sigma }_{3}}} \right) \right. \\ & \left. +i\,\omega \cdot {{e}_{0}}\cdot NN\cdot \left( {{R}_{1}}+d \right)\cdot \left( 1+2Q \right)+2\sigma _{3}^{*} \right]\cdot {{\left( {{R}_{1}}+d \right)}^{3}} \\ & \ \ \ NN=\frac{{{e}_{0}}\cdot {{n}_{01}}}{{{k}_{B}}T}\frac{\frac{R_{1}^{2}}{{{\left( {{R}_{1}}+d \right)}^{2}}}}{1+i\omega \frac{{{\left( {{R}_{1}}+d \right)}^{2}}}{2\cdot {{D}_{3}}}} \\ & \ \ \ g=\left( 1+2\delta \right){{\sigma }_{1}}\sigma _{1}^{*} \\ & \ \ \ +2\left( 1-\delta \right)\left[ i\omega \cdot {{\sigma }_{1}}{{\varepsilon }_{2}}{{\varepsilon }_{0}}\left( 1+s \right)+{{\sigma }_{2}}\left( {{\sigma }_{1}}-i\omega \cdot {{\varepsilon }_{1}}{{\varepsilon }_{0}}S \right) \right] \\ & \ \ \ \delta ={{\left( \frac{{{R}_{1}}}{{{R}_{1}}+d} \right)}^{3}} \\ & S=\frac{\sigma \cdot \left( {{G}_{1}}{{R}_{1}}-Tanh\ {{G}_{1}}{{R}_{1}} \right)}{i\omega \cdot {{\varepsilon }_{1}}{{\varepsilon }_{0}}\left[ {{\left( {{G}_{1}}{{R}_{1}} \right)}^{2}}\cdot Tanh\ {{G}_{1}}{{R}_{1}}-2\left( {{G}_{1}}{{R}_{1}}-Tanh\,{{G}_{1}}{{R}_{1}} \right) \right]} \\ & Tanh\,{{G}_{1}}{{R}_{1}}\approx 1\Rightarrow S=\frac{{{\sigma }_{1}}\cdot \left( {{G}_{1}}{{R}_{1}}-1 \right)}{i\omega \cdot {{\varepsilon }_{1}}{{\varepsilon }_{0}}\cdot \left[ {{\left( {{G}_{1}}{{R}_{1}} \right)}^{2}}-2\left( {{G}_{1}}{{R}_{1}}-1 \right) \right]} \\ & Q=\frac{{{\sigma }_{3}}}{i\omega \cdot {{\varepsilon }_{3}}{{\varepsilon }_{0}}}\cdot \frac{1+{{G}_{3}}\cdot \left( {{R}_{1}}+d \right)}{G_{3}^{2}\cdot {{\left( {{R}_{1}}+d \right)}^{3}}+2\cdot \left( 1+{{G}_{3}}\cdot \left( {{R}_{1}}+d \right) \right)} \\ & {{G}_{k}}=\sqrt{\frac{\sigma _{k}^{*}}{{{D}_{k}}{{\varepsilon }_{k}}{{\varepsilon }_{0}}}};k=1,3 \\ \end{align}$$

p denotes volume concentration of the suspended particles, T the absolute temperature, and kB Boltzmann constant.

Impedance spectra are derived based on εsus*$\varepsilon _{sus}^{*}$using:

Zsus*=GFiωεsus*$$Z_{sus}^{*}=\frac{GF}{i\omega \cdot \varepsilon _{sus}^{*}}$$

where GF represents the Geometric Factor of the measurement chamber.

Essentially, a suspension of shelled particles with a distribution of charges on the inner side, exhibits two types of dispersions:

α, described by counter-ions displacement. When no field is applied (i.e., at equilibrium) the dipole moment of the system (shelled particle and counter-ions) is zero. When applying the AC field, the counter-ions move along the external side of the membrane, their distribution changes and therefore the dipole moment of the system becomes non-negligible. By increasing the angular frequency of the AC field, the capacity of counter-ions to follow the field is gradually exceeded and therefore dispersion occurs usually at frequencies lower (hence the notation α ) than the ones due to interfacial polarization (denoted by β). Notably, the dispersion described by this mechanism originates from the net charge asymmetry across the shell, directly related to the potential difference i.e., the membrane potential [7].

β, described by the Maxwell-Wagner interfacial polarization.

Effect of cell membrane parameters on the permittivity and impedance spectra of a suspension of living cells

The surface charge distribution n01 is related both to the membrane potential Δφ0=dR1ε2ε0(d+R1)e0n01$\Delta {{\varphi }_{0}}=-\frac{d\cdot {{R}_{1}}}{{{\varepsilon }_{2}}{{\varepsilon }_{0}}\left( d+{{R}_{1}} \right)}\cdot {{e}_{0}}{{n}_{01}}$and to the suspension complex permittivity εsus*$$\varepsilon _{sus}^{*}$$via Hn, (Eq. 3) through the term NN (Eq. 4).

We consider a suspension of shelled particles with R1 ~ 2 μm (p = 0.06, σ1 ~ 0.2 S/m, σ3 ~ 0.377 S/m, ε2 ~ 12, σ2 = 10-6 S/m, d ~ 10-2μm, D1 ~ 2×10-10 m2/s, D3 ~ 2×10-9 m2/s) with membrane potential -150 mV, -75 mV and 0 V respectively (fig 1 A, B and C).

Fig. 1A

Permittivity spectra: Δϕ0 = -150 mV (solid), Δϕ0 = -75 mV (dashed), Δϕ0 ~ 0 V (dot-dashed line)

Fig. 1B

Spectra of Impedance magnitude relative to the value at 1 kHz: Δϕ0=-150 mV (solid), Δϕ0 = -75 mV (dashed), Δϕ0 ~ 0 V (dot-dashed line)

Fig. 1C

Phase Spectra: Δϕ0=-150 mV (solid), Δϕ0=-75mV(dashed), Δϕ0 ~0V(dot-dashed line)

As revealed by fig. 1A, the permittivity spectra present a clear dependency of α dispersion on the membrane potential. However, for cells with radius ~2 μm, the spectra of impedance magnitude relative to the value at 1 kHz (impedance level prior to β dispersion) reveal (fig. 1B) very small decrements related to α dispersion, ΔZr ≤ 5×10-3 % raising tough experimental constraints. The same challenge is related to phase variations in the α dispersion (fig 1C) where changes Δθ ≤ 2×10-3 degrees are emphasized.

When considering suspensions of larger cells (R1 ~ 0.5 mm), impedance variation significantly enhances, but α dispersion occurs at lower frequencies (<10 mHz) which considerably increases measurement time and experimental constraints (due to electrode polarization, probe stability etc.).

Figure 2 (A, B and C) reveals the dependency of α and β dispersions on membrane conductivity for all parameters as previously considered but R1 ~ 0.5 mm and membrane conductivity σ2 with 10-6 S/m and 10-5 S/m assigned values.

Fig.2.A

Permittivity spectra for suspensions of large cells (R1 ~ 0.5 mm), Δϕ0 = -150 mV; σ2 = 10-6 S/m (solid); σ2 = 10-5 S/m (dashed line)

Fig. 2.B

Spectra of Impedance magnitude relative to the value at 100 Hz for suspensions of large cells (R1 ~ 0.5 mm), Δϕ0 = -150 mV; σ2 = 10-6 S/m (solid line); σ2 = 10-5 S/m (dashed)

Fig. 2.c

Phase Spectra for suspensions of large cells (R1 ~ 0.5 mm), Δϕ0 = -150 mV; σ2 = 10-6 S/m (solid line); σ2 = 10-5 S/m (dashed)

Figure 2A presents the permittivity spectra, whereas Figures 2 B and 2 C illustrate the spectra corresponding to impedance magnitude and phase. Notably, impedance variation related to α dispersion relative to its value at 100 Hz (impedance level prior to β dispersion), increases in this case to ΔZr ~ 0.5 %.

Effect of the outer medium parameters on the permittivity and impedance spectra of a suspension of charged vesicles (liposomes)

The dependency of α dispersion on viscosity (diffusion) and conductivity of the outer medium is described by expressions in Eq. 4.

We consider a suspension of charged vesicles with R1 ~ 50 nm (p = 0.25, charge density: e0 n01 ~ 1.6×10-2 C/m2, σ1 = 0.2 S/m, ε2 = 4, σ2 = 10-6 S/m, d = 4 nm, D1 = 2×10-10 m2/s).

When assuming: D3 ~ 5×10-9 m2/s, σ3 ~ 0.377 S/m and σ3 ~ 0.1 S/m we find that α and β dispersions superpose: due to the relatively high value of the diffusion coefficient in the suspending medium, the α dispersion is shifted over the β dispersion.

Figures 3 A, B and C show permittivity, impedance magnitude and phase spectra respectively, emphasizing that separation of the frequency ranges of α and β dispersions for suspensions of charged liposomes is achieved by lowering the diffusion coefficient (increasing the viscosity) of the suspending medium to D3 = 5×10-10 m2/s. Also lowering the conductivity of the external medium (to σ3 ~ 0.1 S/m) one notices an increase of α dispersion magnitude as revealed by both permittivity and impedance spectra.

Fig. 3.A

Permittivity spectra for suspensions of charged liposomes: σ3 = 0.377 S/m & D3 = 5×10-9 m2/s (solid); σ3 = 0.377 S/m & D3 = 5×10-10 m2/s (dashed); σ3 = 0.1 S/m & D3 = 5×10-10 m2/s (dot-dashed)

Fig. 3.B

Spectra of Impedance magnitude relative to the value at 100 Hz for suspensions of charged liposomes σ3 = 0.377 S/m & D3 = 5 ×10-9 m2/s (solid); σ3 = 0.377 S/m & D3 = 5×10-10 m2/s (dashed); σ3 = 0.1 S/m & D3 = 5×10-10 m2/s (dot-dashed)

Fig. 3.C

Phase Spectra for suspension of charged liposomes σ3 = 0.377 S/m & D3 = 5×10-9 m2/s (solid); σ3 = 0.377 S/m & D3 = 5×10-10 m2/s (dashed); σ3 ~ 0.1 S/m & D3 = 5×10-10 m2/s (dot-dashed)

Discussion

This study emphasises the dependency of α dispersion (related to displacement of counter-ions) on membrane potential and on several parameters pertaining to the membrane (shell) and to the suspending media, highlighting a set of experimental constraints for measuring impedance (magnitude and phase) spectra exhibiting both α and β dispersions on suspensions of living cells or charged liposomes. Aiming for increased generality to avoid spectra dependency on experimental set-up (i.e. geometric factor of the measurement chamber) the impedance magnitude spectra where plotted relative to the value at a frequency corresponding either to an almost constant value (plateau) between the α and β dispersions, or to the level of impedance magnitude prior to the α dispersion when a plateau region preceding the β dispersion cannot be defined.

Whereas the dependency on membrane potential (modeled based on a distribution of net charge on the inner side of the membrane) is clearly shown by permittivity spectra, impedance spectra reveal very tiny variations pertaining to α dispersion emphasizing the challenges to perform such measurements. We show that larger size of the cells/vesicles is accompanied by an increased variation of the impedance; however, the dispersion shifts towards lower frequencies (< 1 mHz) posing significant experimental challenges (due to electrode polarization and specimen stability).

A slight increase of membrane conductivity is accompanied by an increase of the α dispersion as revealed by both the permittivity and impedance (magnitude and phase) spectra. Therefore, we stress on possible biosensing application of bioimpedance assays related to assessment of pore formation effect of various stimuli.

Analysis of the effects of diffusion and conductivity of the suspending medium on impedimetric (dielectric) behavior of suspensions of charged liposomes reveals that α and β dispersions superpose when assuming the parameters of a normal (physiologic like) medium. However, we show that an increase of medium viscosity (decrease of diffusion coefficient) provides a means to shift α dispersion towards lower frequencies enabling both dispersion to be distinctively appraised. Notably, the decrease of diffusion coefficient is normally accompanied by a decrease of the conductivity (providing an additional increase of impedance variation in respect to α dispersion). Nevertheless, when coping with biological cells, isotonic conditions have to be assured, therefore the limitations for lowering the conductivity of the suspending media beyond physiological values are evident.

Advances of the experimental set-up and appropriate choice of both the type of cells (or liposomes) and of parameters of the suspending medium could lead to detailed interrogation and improvement of the theoretical model, as well as open new applications in biomedicine and related fields.