The equivalent circuit model for a parallel electrode impedance cytometer is analyzed together with Maxwell’s theory of homogenization of mixtures. It deals with both single shell and double shell cell modeling.
The results show that even small errors in parameter estimation lead to variations in model predictions. Furthermore, the possibility of differentiating different biological states at different analysis frequencies is evident from the modeled signal.
The following report is an extract from a project carried out for the course of Modeling and Simulation of Physiological Systems, held for Medical Engineering at the University of Rome Tor Vergata.
Introduction
The study of cellular properties and the state of life of the cell is fundamental in research. Different analyzes show how cells with different vital states show different dielectric properties that can give rise to different measured signals.
We deal with the measurement using an impedance cytometer and the analysis of a simplified model. This model provides a circuit equivalent and is suitable for modeling the peak signal, with the cell between the two electrodes. The model also makes use of Maxwell’s theory of mixtures to homogenize the dielectric properties in order to estimate a single complex permittivity and therefore the resulting impedance of the system through which to derive the signal.
Despite its simplicity, compared to other more complete approaches such as FEM modeling, this approach allows to obtain very fast and useful results for the distinction between different membrane dielectric properties. The modeling and analysis of how changes in the dielectric properties of the individual constituents are then reflected on the modeled signal are then dealt with.
These properties typically mirror the cellular state of life and it is possible to use this modeling to distinguish signals from different cell families such as apoptodic cells, healthy cells or necrotic cells.
Background
Equivalent circuit modeling challenges Maxwell’s Mixture Theory (MMT) for modeling cellular dielectric properties.
This model allows to describe the signal obtained in an impedance cytometer with parallel electrodes when the cell is between a pair of electrodes. The presence of the cell between the electrodes, together with the perturbation of the electric field, generates a differential current signal which is proportional to the diameter through the Clausis-Mossotti factor.
A simplified version of the signal will be considered in the following report, neglecting additional factors related to the properties of the measurement circuit, that is:
\begin{equation}
S=r^{3} f_{C M}
\end{equation}Circuit equivalence
The circuit in fig. 1 provides for the application of a potential to the upper electrodes and the sampling of the signal from the lower electrodes, as a differential current [1].
By observing the path of the current, it will pass through the electrodes and then through the conductive buffer, possibly even meeting the cell. Each of these materials will contribute its own impedance.
These impedances are in series and therefore the measured signal will be:
\begin{equation}
I_{\mathrm{diff}}=\frac{V_{\mathrm{appl}}}{Z_{\mathrm{med}}+2 Z_{d l}}-\frac{V_{\mathrm{appl}}}{Z_{\mathrm{mix}}+2 Z_{d l}}
\end{equation}Where the impedance ๐mix represents the set of cell and conductive buffer and will be described later through MMT.
Furthermore, in this signal, by rearranging the terms, an impedance difference ฮ๐ = ๐mix – ๐med appears equal to the perturbation induced by the cell. This term is also much lower than the impedance of the medium and this allows us to neglect the quadratic terms and obtain a differential current signal equal to:
\begin{equation}
I_{\mathrm{diff}} \approx \frac{V_{\mathrm{appl}}}{Z_{\mathrm{med}}^{2}} \frac{Z_{\mathrm{mix}}-Z_{\mathrm{med}}}{\left(1+2 \frac{Z_{d l}}{Z_{\text {med }}}\right)^{2}}
\end{equation}The electrode impedance is related to the surface capacity ๐ถ๐๐ and to the size of the channel volume between the pairs. It turns out to be equal to:
\begin{equation}
Z_{d l}=\frac{1}{j \omega C_{d l} w l}
\end{equation}The impedance of the medium is expressed through Ohm’s law, contextualizing it according to the volume of the channel:
\begin{equation}
Z_{\text {med }}=\frac{h}{\sigma^{*} l w k}=\frac{1}{\sigma^{*}} \frac{1}{G}=\frac{1}{j \omega \varepsilon^{*}} \frac{1}{G}
\end{equation}Taking into account also distortion effects of the electric field through a geometric coefficient ๐บ. In general, complex permittivity can be expressed as:
\begin{equation}
\varepsilon^{*}=\varepsilon+\frac{\sigma}{j \omega}
\end{equation}These considerations are similar for the impedance of the mix where, however, the ๐ โ mix must be estimated using Maxwell’s theory of mixtures.
Maxwell’s Mixtures theory
Maxwell’s theory of mixtures allows us to describe the complex permittivity of medium and cell as a function of their properties [2].
In particular, permittivity can be expressed as a function of the Clausius-Mossotti factor as:
\begin{equation}
\varepsilon_{\text {mix }}^{*}=\varepsilon_{\text {med }}^{*} \frac{1+2 \varphi f_{C M}}{1-\varphi f_{C M}}
\end{equation}In turn, this factor is expressed as a function of the properties of the cell:
\begin{equation}
f_{C M}=\frac{\varepsilon_{\text {cell }}^{*}-\varepsilon_{\mathrm{med}}^{*}}{\varepsilon_{\mathrm{cell}}^{*}+2 \varepsilon_{\mathrm{med}}^{*}}
\end{equation}Putting together what is described, we obtain the expression for the differential current in eq. (3), as a function of the Clausius-Mossotti factor:
\begin{equation}
I_{d i f f} \approx-\frac{V_{\text {appl }}}{Z_{\text {med }}} \frac{1}{\left[1+\frac{2 Z_{d l}}{Z_{\text {med }}}\right]^{2}} 3 \varphi f_{C M}
\end{equation}Where ๐ is the volume fraction which takes into account the volume of the cell as:
\begin{equation}
\varphi=\frac{V_{\text {cell }}}{V_{\text {misura }}}=\frac{\frac{4 \pi r^{3}}{3}}{l w h K}
\end{equation}And the measured signal, neglecting the dispersive contributions linked to the electrodes, is proportional to ๐3 (through the volume fraction) and to the Clausius-Mossotti factor.
Shell model
Having seen how to take into account the mixed properties of cell and medium, all that remains is to analyze how to homogenize the properties of the single cell. The cells are homogenized to a single material having ๐โ cell properties and such as to take into account both the properties of the inner part and of the cell membrane.
For cells without a nucleus, such as red blood cells, it is possible to use a single shell model, based therefore on a thin outer membrane and an inner zone whose properties are homogenized as:
\begin{equation}
\varepsilon_{c e l l}^{*} \approx \varepsilon_{i n t}^{*} \frac{\chi}{1+\chi}
\end{equation}Where:
\begin{equation}
\chi=\frac{\varepsilon_{m e m}^{*} / d_{m e m}}{\varepsilon_{i n t}^{*} / r}
\end{equation}For larger and more complex cells it is necessary to consider a more accurate model such as the double shell where four different zones are considered. We consider an internal zone, corresponding to the nucleus, (np) and its membrane ne. To these are added the cyt cytoplasm and the mem cell membrane. In particular, homogenizations are carried out in sequence starting from the central area.
Overall properties
Applying a relationship analogous to eq. (11) we get to estimate the overall properties of the nucleus (nuc). Subsequently these can be joined through the MMT thus describing the cellular extension (int). Finally, homogenization is repeated (eq. (11)) between the internal area and the membrane, thus obtaining a complex permittivity for the entire cell.
| Flag | radius [m] | sgm_mem [S/m] | sgm_int [S/m] | eps_int |
|---|---|---|---|---|
| bead | 3 E-6 | 1 E10 | 2.7 E-3 | 2.5 ฮต0 |
| viable | 9 E-6 | 0 | 0.6 | 60 ฮต0 |
| necrotic | 9 E-6 | 1 E-3 | 0.6 | 60 ฮต0 |
| apoptotic intact | 3 E-6 | 0 | 0.6 | 60 ฮต0 |
| apoptotic permeable | 3 E-6 | 1 E-3 | 0.6 | 60 ฮต0 |
Results
These models are then applied to estimate impedance cytometric signals of different cell families in viable, necrotic and apoptotic state [3]. The parameters used, by De Ninno et al., Are reported in Tab. 1.
In order to represent the physics of the problem, the structure class variables were used within Matlab. These variables allow the addition of some fields and therefore the sequence of fields was used to represent the properties owned by the parent structure. For example, the permittivity of the cell membrane will be the ๐๐๐ field of the ๐๐๐ structure in turn the daughter of ๐๐๐๐.
A graphic scheme is present, for the double shell formulation, in fig. 4, and an example in the appendix. This structure, in addition to being representative of the physical element that possesses this property, allows you to iterate some procedures by taking the individual fields of different structures and to create a modular and easily extensible and updatable code. The code has a main routine main.m and several auxiliary routines.
The data is loaded through the data (flag) routine by passing the flag indicated in Tab. 1 to obtain the relative parameters and the structures ๐๐๐๐ and ๐๐๐. The main routine iterates over the different flags to then obtain the signal in fig. 3. For each iteration, after recalling the material properties, these are passed to the routine ๐๐๐๐๐๐๐๐๐๐๐๐ฒ๐๐๐๐๐๐๐ผ๐๐๐๐.๐ which takes care of modeling the equivalent circuit and returning the signal as in eq. (1) and the homogenized complex permittivity (eq. (7), via ๐๐๐_๐๐๐ก ()). To do this, it requires the calculation of the homogenized properties of the cell and the Clausius-Mossotti factor.
Calculation of the Clausius-Mossotti factor
The Clausius-Mossotti factor is calculated using a corresponding routine. The homogenization of cellular properties is carried out through the MMT.m routine which is divided, via the ๐๐๐๐.๐๐๐๐๐ field, in the single shell or double shell implementation. This routine, in addition to calculating the related complex quantities through the epstile.m routine, calculates the homogenization through the mixture.m subroutine which applies the eq. (11).
From the obtained signal it is possible to obtain the electrical diameter as the cube root of the module and the results are shown in fig. 3a, together with the phase of the signal in fig. 3b. From the analysis of the graphs it is possible to identify different peculiarities depending on the families that you want to distinguish. By carrying out measurements at different frequencies it will be possible to distinguish between different vital states in the analyzed population.
Double-shell model
For the double-shell implementation, simply enter the required fields for cellular properties. The parameters used, known in literature [4], are present in Tab. 2. After inserting the parameters in the properties of the ๐๐๐๐ structure, the MMT() routine takes care of doing the double homogenization and returning the values of the ๐ โ mix.
| cell | .np | .ne | .cyt | .mem |
|---|---|---|---|---|
| .eps ฮต0 | 77 | 40 | 1 E10 | 77 |
| .sgm [S/m] | 1 | 4 E-4 | 1 | 8 E-8 |
These values are graphed through the mix_plot () routine which takes care of normalizing the values and returning the graphs of ๐ and ๐ present in fig. 5.
Influence of variations in individual parameters
In addition, the influence of perturbations on the eight descriptive parameters of the double shell model is analyzed. Variations are made according to Irimajiri, Hanai, and Inouye [4].
The perturbation data is then loaded through the structure ๐๐๐๐๐๐๐๐๐๐๐๐_๐๐๐๐๐๐.๐๐๐๐ and the same procedure is iterated over the eight parameters. The result is present in fig. 6. An example of the iteration procedure is provided in the appendix.
From the graphs it is possible to see how the low frequency behavior is mainly influenced by the variation of the outer layers with predominance for membrane properties. At high frequency, the permittivity is affected little and mainly by varying the properties of the cytoplasm and nucleus, while the conductivity is greatly affected by the variations in the properties of the cytoplasm and nucleus. In the central band the slope of the curves is significantly affected by the variation of ๐cyt. In general, changes in conductivity properties have greater effects on the homogenized complex permittivity.
Conclusions
Such evident variations of the homogenized properties can be used to extract information from cytometric signals in order to analyze the cellular state of life.
Different measured properties provide different signals and it is possible to use this theory to trace the differentiation between different cellular states from the signal.
By analyzing the equivalent circuit model with Maxwell’s theory of mixtures, it is possible to understand at what frequency to perform the measurements depending on the cell populations to be measured.
Code availability
The code is available at the project’s GitHub repository.
References
- [1] Carlos Honrado et al. โSingle-Cell Microfluidic Impedance Cytometry: From Raw Signals to Cell Phenotypes Using Data Analyticsโ. In: Lab on a Chip 21.1 (2021), pp. 22โ54. ISSN: 1473-0197, 1473-0189.
- [2] T. Hanai, N. Koizumi, and A. Irimajiri. โA method for determining the dielectric constant and the conductivity of membrane-bounded particles of biological relevanceโ. en. In: Biophysics of Structure and Mechanism 1.4 (1975), pp. 285โ294. ISSN: 0340-1057, 1432-1017. DOI: 10 . 1007 / BF00537642.
- [3] Adele De Ninno et al. โHigh-throughput label-free characterization of viable, necrotic and apoptotic human lymphoma cells in a coplanar-electrode microfluidic impedance chipโ. en. In: Biosensors and Bioelectronics 150 (Feb. 2020), p. 111887. ISSN: 09565663.
- [4] Akihiko Irimajiri, Tetsuya Hanai, and Akira Inouye. โA dielectric theory of โmulti-stratified shellโ model with its application to a lymphoma cellโ. en. In: Journal of Theoretical Biology 78.2 (May 1979), pp. 251269.
Appendix
Loading data into the structure
Below is a snippet extracted from the data.m procedure within a switch-case on the variable ๐๐๐๐. This section of code allows you to load the structures ๐๐๐๐ and ๐๐๐ with their fields containing the dielectric and geometric properties, together with the double flag which will allow the relative subroutines to carry out the double shell homogenization.
case โdoubleโ
cell.shell=โdoubleโ;
mem.thickness=80*1e-10; %[m]
cell.radius=6.5*1e-6; %[m]
int.radius=cell.radius;
mem.sgm=8e-7*1e-1; %[S/m]
mem.eps= 9*eps0;
int.sgm=0.6;
int.eps=60*eps0;
cell.int=int;
cell.mem=mem;
cell.np.radius=4.67*1e-6;
cell.ne.thickness=400*1e-10;
cell.np.sgm=10*1e-1;
cell.np.eps=77*eps0;
cell.ne.sgm=4e-3*1e-1;
cell.ne.eps=40*eps0;
cell.cyt.sgm=10*1e-1;
cell.cyt.eps=77*eps0;
med.sgm= 1;
med.eps= 77*eps0;
cell.volfrac=0.1;
cell.nuc.volfrac=(cell.np.radius/cell.radius)^3;Perturbation
Example of iteration for the perturbation of the permittivity of the cell membrane. The snippet of the variation of ๐mem is shown. A similar sequence is also used on the other 7 parameters analyzed.
% load pertubation values
pert=perturbation_values.eps.mem;
% string for reference
perturbation=โ{\epsilon}_{mem}โ;
real_value=cell.mem.eps/eps0; % create figure and plot
eps_mem_pert=figure();
cell_pert=cell;
for i=1:length(pert)
% change perturbation properties
cell_pert.mem.eps=pert(i);
[~,~,cell_pert,epstile_mix_pert]=equivalentCircuitModel(med,cell_pert,frequency_span);
mix_plot(eps_mem_pert ,epstile_mix_pert ,frequency_span ,omega ,{colorType1{i},colorType2{i}})
clear epstile_mix_pert
end
% resize figs
eps_mem_pert.Position = [100 100 800 400];
lg=legend(strcat(โ=โ," ",string([pert/eps0,pert/eps0])),โLocationโ,โeastoutsideโ);
lg.Title.String=strcat(perturbation ,โ = โ,num2str(real_value),โ [S/m]โ); title(strcat(โPertubation overโ," ",perturbation))
