When subjected to magnetic or electric fields all particles experience torques and forces since they have magnetic and electrical properties related with their containing materials and their shape. Furthermore, due to the particles own electrical charge, magnetization or polarization, they will reveal quite strong mutual interactions. The study of torques and forces exerted on tiny particles (less than approximately 103 micro-meters in diameter) and also the dynamics and mechanics induced by these torques and forces through the action of a magnetic or electric field, is known as particle electromechanics. Lab-
LOC devices apply electrorotation (ER), dielectrophoresis (DEP) and traveling wave dielectrophoresis (twDEP) broadly. All of the mentioned phenomena are created due to the particle induced dipole moment interaction with the electric field. Generalized dielectrophoresis is a suitable title for discussing DEP, ER, and twDEP together in a simultaneous and systematical way 3.
ER is used extensively in manipulation and characterization of particles/cells and is the most accurate method for dielectric characterization among the mentioned phenomena in the literature. It reveals membrane and cells’ interior dielectric properties using the frequencydependent rotational responses. ER is a noninvasive technique, which can be used for finding changes in the physicochemical and morphology properties of microorganisms 4. Applying ER, the influences of different biocide doses on yeast cell suspensions have been studied 5. The specific membrane capacitance value was obtained by Sukhorukov observing rotation of single cells induced by a homogeneous rotating field 6. An essential effort in cell studies is to derive dielectric characteristics of cells from the ER experimental results. In these experiments, the cells inside a medium with permittivity Ɛm and conductivity σm are exposed to an AC field, which rotates uniformly 7. One of the major analyzes in these ER experiments is to determine the frequency where the rate of electrorotation turn into its maximum. This frequency is known as the peak frequency fpk. Since the full expressions for the cell’s effective dielectric properties in terms of fpk are nonlinear equations, they cannot be solved analytically 8. In order to achieve the equivalent cell permittivity ɛp and conductivity σp estimations of low frequency were used by Schwan 9. In the ER experiments with convenient medium conductivity, a couple of peak frequencies may occur. It was suggested to use the high frequency limit for the higher peak frequency 3. In this case, the effective cell permittivity and conductivity can be considered as the cytoplasm permittivity and conductivity, respectively. The estimation of low frequency is useful in the case of measurement inside a less conductive medium and where the peak frequency is adequately low. Furthermore, in the experiments with a higher peak frequency like 100 MHz these estimations are invalid for attaining the dielectric characteristics of cells 10. Thus, this theory is in most situations not suitable for cell dielectric characterization from measured data. Therefore, it is essential to study the cases when the estimations of low frequency would suffer from inefficiency. Discovering the effective substitution for the mentioned theory to investigate the dielectric characteristics of cells using the ER measurements is an essential issue 11.
on-a-chip (LOC) devices which apply electromechanics, have recently become very attractive due to pursuing development of portable, individualized, and implantable devices for point-of-care applications 1. Moreover, LOC devices reveal unique advantages such as, being cost effective, need for very small amount of sample, high level of automation, and high throughput 2.
In this paper, the transformed method in which the comprehensive statements for corresponding permittivity and conductivity of cells are used, is applied 11. Since the current optimization process cannot be used directly for a peak frequency´s nonlinear equation, a distinct method was applied in this study to obtain the cell dielectric characteristics 12. The mentioned transformed theory is credible for the high and the low peak frequencies. On the other hand this method is applicable for all values of conductivity of the medium 11. In order to achieve the cell dielectric properties in this study, genetic algorithm (GA) is applied as a robust means to optimize the equation of modified theory due to its capacity for quick adaptation to the problem and faster convergence without requiring derivative information.
Inhomogeneous particle modeling
Living particles like cells are heterogeneous and a multi shell model could take their heterogeneous structures into account 13. In order to describe non-spherical shells and the plasma membrane dielectric anisotropy, the extended form of this model has been used [14,15].
Maxwell–Wagner polarizations will occur where the various shells join. The temporal dynamics and degrees of these diverse polarizations relate to the dielectric characteristics of the shells constitutive materials.
For a multi shell particle, the Clausius Mossotti function could be achieved by investigating the particle’s effective complex permittivity
in which ɛp and σp represent the particle’s effective permittivity and conductivity respectively.
The terms in Eqs. (3) and (4) can be classified into terms with factor 1/ω2 and terms without this factor. Eqs. (3) and (4) reduce to the forms of Eqs. (5) and (6) keeping each group´s basic term in the numerator and the denominator 11.
Eqs. (5) and (6) are applicable for all angular frequencies, ω. Since the leading term is of two orders greater than the term needed for obtaining the conductance of membrane, the terms of three orders less than the leading term are omitted 11.
In this part, a transformation of the current theory is presented for the cases where the estimation of low frequency limit is not applicable for obtaining the characteristics of cell from the investigations of ER peak frequency. In order to indicate the equivalent permittivity ɛp and conductivity σp , in Eqs. (1), the simplified form of the full expressions Eqs. (3) and (4) are employed in the modified theory 11. Afterwards the nonlinear algebraic equations for fpk is achieved. These equations are implicit equations in terms of cell features, medium permittivity and medium conductivity σm11. The cell characteristics were obtained applying an optimization method to the nonlinear equations, which are mentioned above 12. The equation for the peak rotational angular frequency ωpk and particle/medium dielectric properties is presented below. This equation is obtained setting the derivation of the Clausius Mossotti factor to zero.
where fpk is the peak frequency, ɛm is the effective medium permittivity and σm represents the effective medium conductivity.
Expressing the Eq. (7) explicitly is the essential step in finding the cell characteristics through optimization. Since an explicit solution of fpk from the Eq. (7) cannot be gained, an explicit solution for σm in terms of fpk is used instead. The below equation, which is applied in the modified theory, is obtained by rewriting the Eq. (7) for medium conductivity 11:
As a result, after the shift between the independent and dependent variable’s roles, the optimization is applied. According to Eqs. (3) and (4), (5) and (6), ɛp and σp are not functions of σm, but they both are functions of ω and the cell properties. Therefore, Eq. (8) is an explicit equation for σm.
The number of variables in Eq. (9) is more than that in the recently used theory since the more complete expressions for ɛp and σp shown by Eqs. (5) and (6) are used here. In the present study, d/R and ɛm are assumed to be 1000 and 80 ɛ0, respectively. fpk and R are obtained through the experiments, and the parameters σmem, ɛmem, σint, and ɛint could be obtained via the optimization. This optimization is performed by minimizing the below variance 11,
where applying Eq. (8) would lead to the calculation of σm,theory. The medium conductivity σm,exp is obtained from experiment. The measured data set number, σm and Rfpk,, is shown by N.
In order to optimize Equation (10), GA is applied in this study.
One stochastic method for general search and optimization is genetic algorithm (GA), which is inspired by the biological world. Holland first proposed GA in 1975 18. GA is a population based, iterative, parallel global search algorithm and has been used in many scientific areas such as optimization problems and was found to be very efficient 19. The search is performed by simulating evolution in which more successful organisms in nature appear. Natural selection, which means survival of the fittest is the main idea of GAs. A complicated, long thread of DNA (deoxyribonucleic acid) is called chromosome. Individual´s particular features are determined along these chromosomes. Features are coded by some DNA combination. During reproduction, changes occur in chromosomes. In the crossover process the parents’ chromosomes randomly change thus, the offspring reveals some features of each parent. In addition, some traits are changed during a rarer process called mutation. All the parameters can be altered simultaneously in the GA method; therefore, it can obtain a faster convergence 20.
In this work, a random population was produced to use as the initial population. Then considering the function value of the individual, the fitness was evaluated and best individuals were utilized as parents of the next generation. Genetic transformation such as crossover and mutation applied on these individuals. Two individuals are randomly selected as parents in a crossover process and part of their construction are exchanged to form two new individuals. In the mutation process, one individual is randomly selected from the parent population and its internal representation is changed. In order to improve search performance both of the crossover and mutation rates should be carefully chosen. The produced child population turned to the parent population for the next generation and undergoes the same process until termination conditions, which were determined as stopping criteria, has been satisfied.
The parameters that should be specified in using GA are the number of the initial population and the maximum generation, upper and lower bounds of genes (variables), the chromosomes´ length, the crossover and the mutation probability. The algorithm achieves best results when the values for the first two parameters rise. Nine experimental data sets for medium conductivity and peak rotational frequency obtained by Pethig et al. for insulin secreting cells (primary β-cells and INS-1 insulinoma cells) are applied in optimization in this work 7.
Experimental data sets for medium conductivity and peak rotational frequency.
|Medium conductivity σm (mS/m)||Peak rotational frequency ωpk (MHz)|
For human erythrocytes a radius of equilibrium 2.62 μm and d = 5 nm, d/R=1.91×10-3 is used in this work. In the literature, the following assumptions are used for human erythrocytes and other mammalian cells 21.
The best fitness with the initial population of 60 and the number of generations of 51 was obtained by trial and error. After setting parameters of genetic algorithm in Matlab, the values of fitness vs generation number are shown in Fig. 2.
The mean fitness values are demonstrated in every generation by the points on top of the plot, while the bottom points illustrate the value of the best fitness. The values of the best and mean fitness in 51th generation are equal to 4.7952×10-4 and 4.7956×10-4, respectively. The values of the best fitness’ change decrease at last generations, since the best point is getting closer.
In order to reach optimal cell properties the GA was run many times. Considering the end populations of different runs, the best results for the cell property were demonstrated in Table 2.
the best results for cell properties obtained from genetic optimization
|Membrane Conductivity σmem||7.043 × 10-7 S/m|
|Membrane Permittivity ɛmem||1 × 10-10 F/m|
|Interior Conductivity σint||1.3 S/m|
|Interior Permittivity ɛint||1.507 × 10-9 F/m|
The dielectric properties of cells were obtained applying low frequency approximations. However, transformed theory could be employed in the case the estimations of low frequency are not convenient to obtain the effective permittivity and conductivity of the cell. The transformed theory that is used in this study can lead to logical results for electrorotation experiments performed at the medium with high conductivities. In this theory obtaining dielectric characteristics of cell at every value of medium conductivity using both the high and the low electrorotation peak frequency is possible. Extracting the characteristics of interior and membrane of the cell by measuring electrorotation peak frequency is accomplished in this work. The optimization on medium conductivity´s function in term of peak frequencies resulted in the cell interior and membrane conductivity and permittivity values. GA algorithm is an efficient optimization method, which is used in this study due to its ability of fast convergence and adaptation to the problem.
Kost, G. J. Principles and Practice of Point-of-care Testing. Philadelphia: Lippincott Williams and Wilkins. 2002; Chapter 1: 3-12.
Erickson, D, Li, D. Integrated microuidics devices, Analytica Chimica Acta. 2004; 149: 11-26. https://doi.org/10.1016/j.aca.2003.09.019
Lei, U, Lo, Y.J. Review of the theory of generaliseddielectrophoresis, IET Nanobiotechnology. 2011; 5(3): 86. https://doi.org/10.1049/iet-nbt.2011.0001
Dalton, C, Goater, A. D, Pething, R, Smith, H. V. Viability of Giardia intestinalis Cysts and Viability an Sporulation State of Cyclosporacayetanensis Oocysts Determined by Electrorotation, Applied and Environmental Microbiology, 2001; 586–590
Zhou, XF, Markx, GH, Pethig, R. Effect of biocide concentration on electrorotation spectra of yeast cells, Biochim. Biophy. Acta. 1996; 1281(1): 60-4.
Sukhorukov, V.L, Zimmermann, U. Rotating-Field-Induced Rotation and Measurement of the Membrane Capacitance of Single Mesophyll Cells of A vena sativa, J. Membrane Biol. 1993; 132: 27-40.
Pethig, R, Jakubek, L, Sanger, R.H, Heart, E, Corson, E, Smith, P.J.S. Electrokinetic measurements of membrane capacitance and conductance for pancreatic β-cells, IEE Proc. Nanobiotechnology. 2005; 152(6): 89-193.
Jones, T. B. Electromechanics of Particles, Cambridge University Press, Cambridge. 1995. https://doi.org/10.1017/CBO9780511574498
Arnold, W M, Zimmermann, U. Electro-rotation: Developments of a technique for dielectric measurements on individual cells and particles, J. Electrostatics. 1988; 21: 151-191. https://doi.org/10.1016/0304-3886(88)90027-7
Chung, C, Waterfall, M, Pells, S, Menachery, A, Smith, S, Pethig, R. Dielectrophoretic characterization of mammalian cells above 100 MHz, J. Electr. Bioimp. 2011; 2: 64-71. https://doi.org/10.5617/jeb.196
Lei, U, Sun, P-H, Pethig, R. Refinement of the theory for extracting cell dielectric properties from dielectrophoresis and electrorotation experiments. Biomicrofluidics. 2011; 5 (44109). https://doi.org/10.1063/1.3659282
Huang, Y, Wang, X B, Becker, F F, Gascoyne, P R. Membrane changes associated with the temperature-sensitive P85gag-mos-dependent transformation of rat kidney cells as determined by dielectrophoresis and electrorotation. Biochim. Biophys. Acta, 1996; 1282: 76–84. https://doi.org/10.1016/0005-2736(96)00047-8
Irimajiri, A, Hanai, T, Inouye, A. A dielectric theory of "multistratified shell" model with its application to a lymphoma cell. J Theor Biol. 1979; 21;78(2): 251–269.
Kakutani, T, Shibatani, S, Sugai, M. Electrorotation of non-spherical cells: theory for ellipsoidal cells with an arbitrary number of shells. Bioelectrochem. Bioenerg. 1993; 31: 131. https://doi.org/10.1016/0302-4598(93)80002-C
Sukhorukov, V. L, Meedt, G, Kurschner, M, Zimmermann, U. A single-shell model for biological cells extended to account for the dielectric anisotropy of the plasma membrane. J. Electrost. 2001; 50: 191-204. https://doi.org/10.1016/S0304-3886(00)00037-1
Pethig, R. Dielectrophoresis: An assessment of its potential to aid the research and practice of drug discovery and delivery, Advanced Drug Delivery Reviews. 2013; 65: 1589–1599. https://doi.org/10.1016/j.addr.2013.09.003
Huang, Y, Holzel, R, Pethig, R, Wang, X-B, Differences in the AC electrodynamics of viable and non-viable yeast cells determined through combined dielectrophoresis and electrorotation studies, Phys. Med. Biol. 1992; 37(7): 1499-1517. https://doi.org/10.1088/0031-9155/37/7/003
Holland, J. H. Adaptation in Natural and Artifical Systems, University of Michigan Press, Michigan.1975.
Goldberg, D. E. Genetic algorithm in search, optimization and machine learning. Reading MA Addison Wesley. 1989.
Safak, H, Sahin, M, Gulveren, B, Tomak, M. Efficiency of genetic algorithm and determination of ground state energy of impurity in a spherical quantum dot. Int. J. Mod. Phys. C, 2003; 14: 775-784. https://doi.org/10.1142/S0129183103004917
Gimsa, J, Müller, T, Schnelle, T, and Fuhr, G. Dielectric spectroscopy of single human erythrocytes at physiological ionic strength: dispersion of the cytoplasm. Biophys. J. 1996; 71: 495-506. https://doi.org/10.1016/S0006-3495(96)79251-2