E. Hernández-Balagueraa,
E. López-Doladob and
J. L. Polo*a
aEscuela de Ingeniería Industrial de Toledo, Universidad de Castilla-La Mancha, Avda. Carlos III, s/n, 45071 Toledo, Spain. E-mail: JoseLuis.Polo@uclm.es
bHospital Nacional de Parapléjicos, SESCAM, Finca La Peraleda s/n, 45071 Toledo, Spain
First published on 19th February 2016
We have adapted the current interruption method to obtain electrical equivalent circuits of rat tissues (liver viscera and triceps surae muscle). This paper provides a comprehensive and in-depth explanation of the procedure, and its suitability has been assessed. Experimental data are interpreted using circuit theory and fractional calculus. Histological analysis was also carried out. The technique proposed has a clear physical meaning and we have tried to make this communication accessible to potential readers with multidisciplinary backgrounds. To the best of the authors' knowledge, there is no previous report of adapting the current interruption method to obtain electrical equivalent circuits of tissues using a four-electrode arrangement, circuit theory and fractional calculus.
Fig. 1(A) shows a schematic of single eukaryotic cell. Intracellular (pink color) and extracellular (blue color) spaces are electrolyte solutions containing primarily Na+, K+ and Cl−. Cell membrane (green color), about 7 nm thick, separates intra- and extracellular compartments regulating the ionic concentration differences (cellular homeostasis).2
Let us consider an electric current flowing from left to right through the single cell of Fig. 1(A). In general, a portion of this current will cross the cell membrane and will enter the intracellular fluid. The remaining current will flow through extracellular fluid. Both current paths are implemented in the electrical equivalent circuit (EEC) of Fig. 1(B) using two parallel branches. The upper branch (capacitance CM in series with a resistance RI) models the current path through the cell membrane and the intracellular space. The lower branch (resistance RE) considers the current path through the extracellular compartment. Note that cell membrane acts as an electrical capacitor by involving a charge separation. Intra- and extracellular environments (electrolyte solutions) are represented by the intra- and extracellular resistances, RI and RE, respectively. EEC of Fig. 1(B) was initially proposed by Fricke and Morse.3
Fig. 1(C) shows a portion of tissue which involves a group of similar cells. Thus, the EEC of Fig. 1(D) comprises an array of interconnected EECs (each EEC corresponds to a single cell) which provides a modified Fricke and Morse EEC. The capacitor CM is replaced by a fractional capacitor,4 or more generally known as constant phase element (CPEM) where QM and α are the CPEM parameters (see below).2,5 CPEM takes into account the space distribution of the electrical tissue properties. While the EEC of Fig. 1(B) exhibits one single time-constant, the EEC shown in Fig. 1(D) involves a distribution of time constants.
Fig. 1(E) redraws the EEC of Fig. 1(D). i(t) is the current supplied to the circuit (excitation). Currents and voltages have been labeled in each element. iM(t) and iE(t) are the currents through the CPEM and RE, respectively. vM(t), vI(t) and vE(t) are the voltages across the CPEM, RI, and RE, respectively. Eqn (1) expresses iM(t) as a function of vM(t),
![]() | (1) |
It is important to point out that the CPEM and the resistors RI and RE are LTI (linear time-invariant) circuit elements. Indeed, the parameters QM, α, RI, and RE do not vary with time (time-invariant property) and do not depend on voltages across and/or currents through the corresponding circuit element under consideration (linear property). It follows that the EEC of Fig. 1(E) is a LTI circuit.6
Fig. 1(F) shows the EEC of Fig. 1(E) transformed into the frequency-domain which is usually used in IS. Thus, i(t), iM(t), iE(t), vM(t), vI(t), and vE(t) are replaced by their corresponding phasors Ī, ĪM, ĪE, M,
I, and
E, respectively. The impedances of the circuit elements, ZCPEM(jω), RI, and RE, are also indicated.
Specifically, the impedance of the CPEM (ratio M/ĪM) is written as5
![]() | (2) |
The total impedance between the endpoints of the EEC of Fig. 1(F) is found be
![]() | (3) |
In eqn (3), we find that the coefficient of (jω) – physical frequencies – is the distributed time-constant (τ) of the EEC of Fig. 1(E) when a current is supplied to it; that is,
![]() | (4) |
Observe that if α = 1 (CPEM is an ideal capacitor), eqn (3) yields the impedance of the EEC of Fig. 1(B) and its time constant is written as (RI + RE)CM where (RI + RE) is the resistance seen by CM (a current source is replaced by an open circuit). The inverse of the time constant provides the characteristic frequency ωo (known as frequency pole in first-order (α = 1) LTI circuits)7
![]() | (5) |
It is worth remembering that impedance Z(jω) is a complex number whose real and imaginary parts are the resistance R(ω) and the reactance X(ω), respectively. That is, Z(jω) = R(ω) + jX(ω). Note that at sufficiently low and high frequencies, eqn (3) yields approximately RE and RIRE/(RI + RE), respectively. The same conclusion can be stated from the EEC of Fig. 1(F). Suppose that a sinusoidal current i(t) is connected to the EEC of Fig. 1(E) and a sinusoidal steady-state is reached (refer to Fig. 1(F)). At sufficiently low frequencies, the impedance magnitude of the upper branch |ZCPEM(jω) + RI|, will become large enough to cause |ZCPEM(jω) + RI| ≫ RE, so almost all the current i(t) will flow through the extracellular space (resistor RE), resulting in iE(t) ∼ i(t). At sufficiently high frequencies, RI ≫ |ZCPEM(jω)|, the current will flow through both intra- and extracellular media according to the current-divider rule (circuit consisting of RI and RE).
Interestingly, the expression in eqn (3) is of the same form as that proposed by Cole8
![]() | (6) |
Of specific interest in IS is the Nyquist plot (−X(ω) vs. R(ω) in the complex plane). Nyquist plot of eqn (3) sketches a depressed semicircle (the center of its corresponding circle lies below the real axis).1,2 This depressed semicircle intersects the real axis at two points: RIRE/(RI + RE) and RE at the frequencies ω → ∞ and ω = 0, respectively (see above).
It is easy and intuitive to see that changes in a biological tissue will be reflected in the corresponding EEC. For instance, consider an ischemic tissue. Cell swelling causes a reduction of the extracellular volume which yields an increase of the extracellular resistance RE in the EEC of Fig. 1(D).9 Next, to determine the parameters of the EEC shown in Fig. 1(D), electrical measurements are taken on the tissue under study.
One simple way of obtaining the electrical properties of biological tissues is using the four-electrode or Kelvin method,1,2 as shown in Fig. 2(A). Four electrodes are placed into the tissue in a straight line. An electric current i(t) is injected from electrode 1 (EL1) and collected at electrode (EL4). It determines a potential distribution in the tissue, resulting in a voltage v(t) between the two inner electrodes (EL2 and EL3). This setup (voltage-measuring electrodes placed in line and between the current-injecting electrodes) may usually provide a good signal to noise ratio. Fig. 2(B) shows the EEC of the four-electrode method. EECELi and EECELi/T are the electrical equivalent circuits of the electrode i and (electrode i)/(surrounding tissue) interface, repectively. EECTij is the equivalent circuit of the tissue portion carrying the current i(t) and placed between the two equipotential surfaces (perpendicular to current flow lines) which contain the electrodes i and j. Note that in circuit theory, we find voltages and currents which can also be measured using voltmeters and ammeters. Nevertheless, field theory should be regarded if we want to analyze the potential or current distributions in tissue or in the electrode/tissue interface.10
Simple EECs of ELi and ELi/T interface comprise a resistor and a capacitor (or CPE) in parallel with a resistor, respectively. Electrodes may also involve more complex processes such as anomalous diffusion.11 Typical EECTij is shown in Fig. 1(D). Note that in sinusoidal steady-state analysis, EECELi, EECELi/T, and EECTij are replaced by their corresponding impedances ZELi(jω), ZELi/T(jω), and ZTij(jω), respectively.
Since an ideal voltage-measuring device draws no current (open circuit), it follows from Fig. 2(B) that the measured voltage v(t) between terminal 2 and terminal 3 is equal as that of EECT23. This makes the four-electrode method suitable for characterizing the average electrical properties of the tissue portion which carries the current i(t) and contributes to the voltage v(t) between terminals 2 and 3. Fig. 2(C) shows the EECT23 of the tissue portion which is measured. Note that the measured voltage v(t) is indeed the voltage vE(t) across the extracellular space (resistance RE). We shall henceforth use vE(t) and v(t) interchangeably. Ideally, v(t) does not take into account the effects of electrode polarization (voltages across EECELi/T). Errors in four-electrode configuration have been analyzed elsewhere.12 Note that Fig. 1(E) is repeated here in Fig. 2(C) for the reader's convenience.
As the integrity of the cell membrane is a critical point to maintain the cell's life, the question arises as to what frequency range should be used to characterize its electrical properties. Cellular structure of tissues exhibits a characteristic frequency, ωo/(2π), ranging from several kHz up to MHz. That range is known as β dispersion.13,14 Thus, if a specific tissue has a characteristic frequency ωo/(2π), an ideal situation is to cover an useful frequency range on both sides of ωo to obtain accurate values for the parameters of the EEC. It should be pointed out that the tissues also exhibit other characteristic frequencies (α or γ dispersions) corresponding to different relaxation processes which may be evidenced in an impedance spectrum.13,14
We use entirely the LTI circuit theory.6 Hence, the i(t)–v(t) relationship for the tissue (see Fig. 2(A)) should be described by a linear differential equation with constant coefficients. Although values of current, electric field or voltage have been proposed to ensure linear condition, they result in a large variability depending on the type of tissue and its state.13–15 Additionally, the properties of a biological tissue can vary with time (for instance, during an ischemic event) and, thus, the measurement acquisition time should be sufficiently short to assure that the tissue is “frozen” in time (the electrical properties of the tissue do not change during the measurement acquisition time). This requirement is very critical for some processes such as electroporation.15
An important consideration is the current waveform i(t) to be injected. Characteristic current values should be sufficiently small to ensure the linearity (see above) but sufficiently large to allow an acceptable signal to noise relation without inflicting damage to the tissue (heating, pH changes or accumulation of toxic electrochemical products). The current density (the current per unit area) distribution in the tissue should be adequate to not inflict damage.10
Typically, experimental impedance measurements are obtained in a sequential way: impedance at each frequency is measured using a single sinusoidal-signal which may be time-consuming. The use of excitation signals different from that sinusoidal-signal for analyzing electrical properties of tissue has been barely reported in the scientific literature.15–20 Some advantages such as a drastic reduction of the measurement acquisition time and/or an optimum frequency spectrum which contains sinusoidal components of interest were emphasized.17
EECs of biological tissues have been obtained from the experimental transient-response to a step function.19,20 Basically, a step excitation corresponds to a switch which closes at a specific instant and connects a dc battery to a given circuit.6 For instance, the body composition was obtained from a step-voltage response using the two-electrode method.19 Cole parameters were determined in apple samples using a step-current response, fractional calculus, and also the two-electrode arrangement.20
In this paper, we have obtained EECs of rat tissues using the current interruption method (see below) which has a clear physical meaning. Cell membranes are constant current charged. When the current is cut-off, cell membranes begin to discharge through the surrounding intra/extracellular spaces. The analysis of the discharge voltage curve involves circuit theory and fractional calculus. Importantly, this discharge voltage curve depends only on the tissue's electrical properties.
The current interruption method has been used in electrochemistry to estimate ohmic losses.21 We have adapted this method for obtaining the parameter values of the EEC shown in Fig. 1(D) using a four-electrode arrangement.
![]() | ||
Fig. 3 Current interruption method. (A) Waveform of the current injected into the tissue. (B) Waveform of the resulting extracellular voltage. (C) EEC of Fig. 1(E) at t = 0−. (D) EEC of Fig. 1(E) at t = 0+. (E) EEC of Fig. 1(E) for t > 0. |
Firstly, a dc current IDC (excitation) is injected into the tissue between the two outer electrodes for a sufficiently long time T (T ≫ τ, see below), so that the voltage (response) between the two inner electrodes reaches a steady state (constant value). At time t = 0, the current is abruptly cut off. Fig. 3(C) shows the EEC of Fig. 1(E) at t = 0− (just before the current is interrupted). CPEM is considered to be fully charged (dc steady-state) and acts as an open circuit. Thus all the current IDC flows through the extracellular space (resistor RE), yielding a voltage drop of REIDC. Since no current flows through the cell membrane (upper branch), no voltage appears across the intracellular space (resistor RI). It follows that the voltages across CPEM and RE are equal. Voltages and currents at t = 0− are indicated on the circuit diagram in Fig. 3(C).
Just after interrupting the current (that is, at t = 0+), the energy previously stored in the cell membranes (CPEM) begins to be released and dissipated in the surrounding intra/extracellular spaces. Fig. 3(D) shows the EEC of Fig. 1(E) at t = 0+. It is assumed that the voltage across the CPEM cannot change instantaneously, so that the voltage across it just before (t = 0−) and just after (t = 0+) interrupting the current are equal. At t = 0+, the CPEM voltage drives a current REIDC/(RI + RE) through RE and RI resistors, giving voltages of RE2IDC/(RI + RE) and RIREIDC/(RI + RE), respectively (refer to Fig. 3(D)).
Next, the voltage relaxation (free or natural response) is measured for t > 0. Importantly, the natural response obtained for t > 0 depends only on the tissue's electrical properties. Fig. 3(E) shows the EEC of Fig. 1(E) for t > 0. Let us consider that CPEM is an ideal capacitor (α = 1 and QM = CM). From circuit theory, it can be shown that the expression for v(t) is written as v(t) = [RE2IDC/(RI + RE)]e−t/τ, where the time constant τ is (RI + RE)CM which is as that of eqn (4) with α = 1 and QM = CM. Nevertheless, CPEM involves a fractional derivative – see eqn (1) – and the Mittag–Leffler function is used to include a fractional relaxation.
![]() | (7) |
![]() | (8) |
Note that if α = 1, eqn (8) is the Maclaurin series for the exponential function, resulting in E1(−t/τ) = e−t/τ.
The theoretically predicted waveform of v(t) has been sketched and clearly labeled in Fig. 3(B). Note that RIREIDC/(RI + RE) is the magnitude of the jump discontinuity at t = 0. Fig. 3(B) also shows that the fractional relaxation (eqn (7)) exhibits a faster (slower) decay than the exponential one, for short (large) times.
The values of the parameters of the EEC shown in Fig. 1(D) are determined by performing the following procedure:
(i) RE is found from the previous steady-state (value v(0−)) as
![]() | (9) |
(ii) RI is obtained from the value of v(0+) as
![]() | (10) |
(iii) The transient-voltage response is fitted to eqn (7); α and τ are then determined. Finally, QM is calculated using eqn (4). The fitting procedure was performed using a routine in MATLAB.23
We have assumed that the current is interrupted instantaneously (from t = 0− to t = 0+). Nevertheless a finite (nonzero) fall time (tf) exists in a real-life setup. The tissue will “feel” that the current is abruptly set to zero if tf ≪ τ. It should also be noted that a transient recorder with a sufficiently sampling rate is needed in order to more accurately determine the critical points (jump discontinuity and the initial fast decay of the fractional relaxation) of the experimental waveform of v(t).
Surgical procedures of organ and muscle extraction were performed under intraperitoneal (IP) analgesia with xylazine (10 mg kg−1) and anesthesia with sodium pentobarbital (55 mg kg−1) mixed with atropine (0.05 mg kg−1). Eyes were covered with Lubrithal™ gel to prevent corneal abrasion and dehydration. In all cases, extreme care was taken to minimize muscle tissue damage. When deeply anesthetized, the skin and muscles over the breastbone were cut and the sternal manubrium was grasped to gain access to the thoracoabdominal union. The abdominal cavity was opened by cutting through the abdominal wall in the midline from the tip of the processus xiphoideus to the pecten ossis pubis. The liver was removed and imbibed in saline solution. After that, a right Achilles tenotomy was performed; the triceps surae was dissected, removed and immersed in saline solution just prior to carry out the electrical measurements. Rats were sacrificed immediately after the liver and triceps surae extraction with an IP lethal dose of sodium pentobarbital (110 mg kg−1). Finally, organs were carefully examined looking for any findings of gross damage and properly prepared for subsequent electrical characterization.
We analyzed three samples of livers and other three of triceps surae muscles from three rats. Electrical measurements were taken on each sample using the current interruption method. Just after that, conventional bioimpedance measurements were obtained. This protocol is then repeated twice more in each one of the livers and muscles, to guarantee the reproducibility and repeatability of experimental data. We selected these two particular tissues because of their opposite characteristics related to cellular homogeneity and electrical excitability. So, we chose the liver as an example of high degree of cellular homogeneity but bad electrical excitability and the striated muscle (triceps surae), as an example of low degree of cellular homogeneity but good electrical excitability.
Electrical measurements were performed using four platinum needle electrodes. Each electrode (refer to Fig. 4(A)) was 10 mm long and 0.3 mm in diameter (Grass Technologies, subdermal needle electrode, model F-E2-48). The electrodes were placed into the tissue in a straight line and equally spaced (5 mm). Fig. 4(B) and (C) show the four-electrode arrangement for liver and triceps surae, respectively.
An AutoLab PGSTAT302N Potentiostat/Galvanostat (EcoChemie), combined with the ADC10M module, was used to implement the current interruption method. The instrument was controlled by a computer and driven by NOVA software. A dc current of 100 μA was injected between the two outer electrodes for 1 ms and just after that, the current is abruptly switched off. The voltage between the two inner electrodes was recorded between 1 ms before and 100 or 50 μs after switching off the current. A sampling rate of 10 MHz was used.
A Solartron 1260 frequency response analyzer connected to a Solartron 1294 impedance interface system were used to obtain the impedance measurements. Both instruments were controlled by a computer and driven by SMaRT software. Sinusoidal signals of 100 μA amplitude were injected between the two outer electrodes. The frequency-values were logarithmically spaced (5 steps per decade) between 1 MHz and 1 kHz.
After making electrical measurements, all samples were placed in paraformaldehyde 4% at 4 °C for immersion-fixation and later, three days more in sucrose (30% in phosphate buffered saline, PBS) at 4 °C for cryo-protection. Tissue pieces were mounted on plastic containers, quick-frozen in optimal cutting temperature compound (Tissue Tek, Hatfield, PA) and cut in transverse sections of 10 μm by using a Microm HM550 cryostat. Livers and muscles were initially examined using hematoxylin-van Gieson (Fig. 4(D)) and ATPase (Fig. 4(E)) stainings, respectively. Digital images were collected with a fluorescence Olympus BX51 microscope with 2×, 10×, 20×, and 40× objectives.
Fig. 5(A) and (B) show the resulting voltage for the liver and triceps surae, respectively. The shape of the waveforms is similar to that predicted theoretically (refer to Fig. 3(B)), showing: (i) the previous steady state; (ii) the jump discontinuity. Values of v(0−) and v(0+) have been indicated by arrowheads; and (iii) the voltage relaxation for t > 0. Simulated data are also shown in Fig. 5(A) and (B). There is an excellent agreement between the experimental and simulated data. Average errors of 0.46% and 0.62% were obtained for liver and triceps surae, respectively. The following function was used: , where vexp(ti) and vsim(ti) are the experimental and simulated voltages at time ti, respectively, vmaxexp is the maximum value of vexp(ti), and N is the total number of points.
Table 1 specifies the parameter values of the EEC shown in Fig. 1(D), using the current interruption method. They have been obtained by applying the procedure outlined in Section 2. Additionally, the values of the time constants associated with each of the two tissues have also been listed in Table 1.
Biological tissue | Measurement method | RE (Ω) | RI (Ω) | QM (nF sα−1) | α | τ (μs) |
---|---|---|---|---|---|---|
Liver viscera | Current interruption method | 1160.6 | 106.6 | 361.1 | 0.74 | 30.7 |
Impedance spectroscopy measurements | 1146.1 | 107.8 | 412.6 | 0.73 | 31.5 | |
Triceps surae muscle | Current interruption method | 231.3 | 25.9 | 1432.3 | 0.73 | 19.8 |
Impedance spectroscopy measurements | 238.7 | 25.3 | 1260.3 | 0.72 | 14.8 |
Note that we should check whether end up with a consistent solution. The distributed time constants for the liver and triceps surae are about 30 μs and 20 μs, respectively. Thus we see that a previous steady state is reached because 1 ms ≫ 30 μs and 1 ms ≫ 20 μs. A sampling rate of 10 MHz allowed more accurate values for v(0−) and v(0+) and to measure the time taken for the experimental setup to change the current from the previous steady-state value to zero. This fall time (tf) was of about 1 μs. Thus, the tissue “feels” that the current is abruptly switched off because 1 μs ≪ 30 μs and 1 μs ≪ 20 μs. For clarity's sake, only the samples taken every 1 μs are shown in Fig. 5(A) and (B). The previous steady state is shown during the 20 μs preceding the current interruption.
Typical Nyquist plots for the liver and triceps surae are shown in Fig. 5(C) and (D), respectively. Well defined depressed semicircles can be observed. Experimental impedance data were fitted to the EEC shown in Fig. 1(D) using ZView software. The parameter values obtained are listed in Table 1. Impedance plots also show simulated data which have been obtained using the EEC of Fig. 1(D) and the parameter values given in Table 1. As can be seen there is an excellent agreement between the experimental and simulated results. Average errors of 0.27% and 0.51% were obtained for liver and triceps surae, respectively. We have used the following expression: , where Zexp(jωi) and Zsim(jωi) are the experimental and simulated impedance data at the frequency ωi, respectively, Zmaxexp is the maximum impedance magnitude of experimental data, and N is the total number of points. Note that the calculated average errors are normalized to the maximum value in the data set.24
Table 1 shows that the parameter values obtained using the current interruption method are remarkably close to those found from IS measurements. We note minor discrepancies of less than 3% for RE, RI, and α parameters and of the order of 10% for QM parameter. Determination of QM is more prone to errors because it depends on RE, RI, α, and τ. In turn, α and τ are obtained from the fitting procedure (see Section 2). Even so, a reasonable estimate for QM could be obtained.
The measurement acquisition time (1.1 ms and 1.05 ms for the liver and triceps surae, respectively) using the current interruption method is much shorter than that of the conventional impedance spectroscopy technique (∼35 s). Thus, it may allow monitoring of fast changes in a biological tissue. Nevertheless, as we have seen earlier, specific requirements on current interruption time and sampling rate should be fulfilled.
This journal is © The Royal Society of Chemistry 2016 |