Dielectrophoresis (DEP) is a technique that has been used to characterize and sort cells based on their distinct dielectric properties [1, 2, 3]. It refers to a force induced on polarisable particles in non-uniform electric fields. This force is proportional to both the field gradient and the real part of the particle’s effective polarisability (per unit volume), described by the Clausius-Mossotti function:
Here,
The effective dielectric properties of a mammalian cell can be modelled by applying the multi-shell method first described by Irimajiri [4]. The nucleoplasm, nuclear membrane, cytoplasm and cytoplasm membrane are represented as concentric spherical shells with distinct radii and dielectric properties (defined by their permittivity and conductivity). By successively merging these shells together an expression for the effective complex permittivity of an entire cell can be produced. Asami
We use the form of figure 1 as a predictive guide to the DEP response of mammalian cells. The frequency regions exhibiting negative and positive DEP force correspond respectively to repulsion from, and attraction to, regions of high electric field concentration that occur close to the edge of electrodes. Between 100 kHz and 1 MHz this force transitions from negative to positive at the first (low frequency) cross-over frequency, fxo1. The value of fxo1 depends upon the cytoplasm membrane capacitance and resistance, the cell diameter and the suspending medium conductivity [2, 6]. Determination of fxo1 for different cells and cell states has been the principle parameter for characterizing and sorting cells by DEP. In figure 1 it can also be seen that, at a frequency between 100 MHz and 1 GHz, a second cross-over, fxo2, is predicted to occur. Above 100 MHz the membrane capacitive reactance should effectively short out the membrane resistance, so that the electric field penetrates into the cell interior [3]. The value of fxo2 for viable cells is thus predicted to depend on the intracellular dielectric properties.
A general expression for determining the DEP cross-over frequencies is derived by equating the real part of equation 1 to zero. This yields equation 2, where εm and σm represent the permittivity and conductivity of the suspending medium; εc and σc the effective permittivity and conductivity of the cell [7]:
The customary practice is to employ the DC approximations derived by Schwan [8] for the low-frequency effective values of σm and εm – namely σm = rGm and εm = rCm, where
From equations 3 and 4 we find that the high frequency DEP cross-over fxo2 is determined by the internal cell dielectric parameters, namely the nucleoplasm permittivity (εnp), nucleoplasm conductivity (σnp), cytoplasm permittivity (εcp), cytoplasm conductivity (σcp) and nucleus volume fraction (v = nucleus volume/cell volume). Based on dielectric studies of mouse lymphocyte suspensions, Asami
Interdigitated electrodes are commonly used in DEP experiments, but understanding their behaviour above 100 MHz requires accurate modelling of the electrical load being driven by the signal source. An overview of the electrical system used to conduct fxo2 measurements is shown in figure 2. The signal generator provides a sine wave between 2 MHz and 500 MHz at up to 2 Vpk-pk. This is fed to a high power amplifier with ~17.5 dB gain over this range when terminated into a 50Ω load. A 50Ω resistor in parallel with the DEP electrodes is used to approximate this termination, connected by a short ¼ metre length of coaxial cable. This prevents reflection issues with longer cables that reduce bandwidth. The interdigitated structure is fabricated from 100nm thick platinum deposited onto a Pyrex substrate, patterned using standard photolithographic techniques and etched by argon plasma milling.
Two interdigitated electrode ‘fingers’ are shown in figure 3, with a signal being applied to the left pad and grounded on the right pad. An example of a current path is displayed in red, flowing from the signal pad and along the top finger track. At some point this current crosses through the suspending medium to the bottom finger, along the remaining track to ground. A continuum of such current paths exist along the overlap length (LO), determined by track resistance, suspending medium conductivity and capacitive coupling between the fingers.
The distributed RCR network of figure 4 has been used to model these parameters. Track resistance (Rt), suspending medium resistance (Rs) and parallel track capacitance (C) components were used for sections of length Δx - LO divided by the number of sections. Sheet resistance from cross bridge resistor measurements was used to estimate Rt. Values for Rs and C per unit length along the fingers were determined by finite element models of a cross section in COMSOL Multiphysics software. The accuracy of this distributed network increases as Δx approaches zero, but beyond five sections the differences between impedance spectra were observed to be marginal. Ten sections were used to provide sufficient detail for voltage distribution plots.
The voltage difference between the adjacent fingers creates the non-uniform electric field
Specifications for two examples of interdigitated electrodes are given in table 1, fabricated from 50 and 100 nm thick platinum.
Model parameters for two electrode design examples.
Design 1 | Design 2 | ||
---|---|---|---|
Number of Fingers | 16 | 16 | |
Finger Overlap (LO) | 20 | 2 | mm |
Finger Width (LW) | 20 | 25 | μm |
Finger Spacing (LS) | 20 | 35 | μm |
Medium Conductivity | 33 | 33 | mS/m |
Sheet Resistance | 3.6 | 1.8 | Ω/□ |
Section Length (Δx) | 2 | 0.2 | mm |
Section Medium Resistance (Rs) | 30 | 340 | kΩ |
Section Track Resistance (Rt) | 28.8 | 14.4 | Ω |
Section Track Capacitance (C) | 1500 | 67 | fF |
Area | 12.8 | 1.9 | mm2 |
The voltage difference profiles shown in figure 5 take the form of catenaries due to the resistive and capacitive loading along the length of the electrode fingers. The bandwidth, defined by a minimum voltage difference of 3.5V, is only 8 MHz and far below the upper range of 320 MHz predicted for fxo2. Electrode design 2 shown in table 1 features an increased finger width and spacing, which reduces the finger overlap, and double the platinum thickness. As shown in figure 6, this flattens the voltage difference profile and increases the DEP operating bandwidth to 400 MHz. However, this is achieved at the cost of significantly reducing the electrode area. Electrode design 2 was used for the DEP studies described here.
Murine myeloma SP2/O cells were grown in RPMI-1640 medium supplemented with 100 units/mL penicillin-streptomycin and 10% foetal calf serum (FCS). The cells were grown in suspension with 5% CO2 in air at 37°C to a density of 1x106 cells/mL. An iso-osmotic cell suspension medium for DEP studies was prepared by adding 3 g/L glucose and 100 g/L sucrose to DI water. Conductivity was adjusted to 33 mS/m with PBS, verified using a conductivity meter (Oakton CON 510) and the pH adjusted to 7.4 with NaOH/HCl. The osmolality of the final medium was 310 mOsm/kg, verified with an osmometer (Advanced Instruments Inc. Model 3300). DEP media of osmolality 250, 390, 435 and 480 mOsm/kg were also prepared by adjusting the sucrose concentration. The SP2/O cells were centrifuged (100g, 5 minutes) and washed twice in 10 mL of the DEP media, before final suspension at a density of 1x107 cells/mL for analysis. The electrodes were mounted in an inverted microscope (Meiji TC5100). By focusing on the electrode plane discrimination between cells undergoing negative DEP (out-of-focus or between the electrodes) and positive DEP (in-focus at electrode edges) was possible.
Figure 7 shows the distribution of values for the high-frequency crossover, fxo2, with a population sample of 418 SP2/O cells suspended in a DEP medium of osmolality 310 mOsm/kg. Under these conditions the vast majority of cells were observed to undergo the transition from negative to positive DEP when sweeping down the frequency range.
The mean value of fxo2 for the cell populations shown in figure 7 is 195 MHz with a standard deviation of 63 MHz. An interesting finding, shown in figure 8, is the reduction of the fxo2 values as a function of time after their suspension in the DEP solution.
A corresponding temporal change in the low-frequency cross-over, fxo1, was not observed, and is not expected unless the cell state is deliberately altered by adding chemical agents to the solution known to induce cell activation or apoptosis for example [2, 9].
The reduction in fxo2 over time was found to be a function of temperature, as shown in figure 9. Cell suspensions were held in a temperature controlled water bath (situated in a refrigerated room for the 10°C run) between measurements. The rate of frequency reduction based on the first two hours clearly increases with temperature, more than doubling in rate (~2.4x) for increases of 10°C between 21°C and 37°C.
At 37°C the rate of fall of fxo2 was so great that accurate determination of it became difficult after two hours. The osmolality of the suspending solution was also found to influence the mean fxo2 value, as shown in figure 10. The higher osmolality solutions tended to result in lower values of fxo2 for the same period of cell suspension in the DEP solution.
The distribution of high frequency cross-over, fxo2, values for SP2/O cells shown in figure 7 was obtained by slowly sweeping down from an applied signal frequency of 400 MHz. At this initial frequency all of the cells were observed to levitate above the electrode plane under the action of negative DEP, but on proceeding below 325 MHz subpopulations of the cells were attracted to the electrode edges by positive DEP. This behaviour confirmed that the predicted high-frequency DEP cross-over shown in figure 1 for mammalian cells can be observed under our experimental conditions, and that the designed electronics and electrodes are capable of exploring this high frequency region of the DEP spectrum. The mean fxo2 value of 195 MHz shown in figure 7 falls well within the 90 MHz to 320 MHz range of values we predicted using the cell dielectric parameters given by Asami
Equations 3 and 4 inform us that the high frequency cross-over, fxo2, is determined by the permittivity and conductivity of both the nucleoplasm and cytoplasm, as well as the nucleus volume fraction. The nucleus volume fraction is known to change with cell cycle [10], and could be one factor responsible for the relatively wide distribution of fxo2 values shown in figure 7. The relationships of the dielectric parameters to specific biological features of a cell are at present not clear, but an insight into the relative sensitivity of fxo2 to them can be obtained.
Table 2 provides a measure of sensitivity in fxo2 to each dielectric parameter based on the modelled values undergoing a +/- 10% variation. A nucleus volume fraction of 0.7 is used here with values approaching ~0.9 being typical of lymphocytes which are associated with our SP2/O cells [11]. Finally, the intracellular conductivity is scaled for a base fxo2 value of 195 MHz to which these variations are normalised.
Sensitivity of fxo2 to dielectric parameters based on Asami
Permittivity | Conductivity | ||
---|---|---|---|
Medium | 31% | 0% | |
Nucleus | 8% | 13% | |
Cytoplasm | 6% | 1% | |
Nucl. Vol. Fraction | 8% |
Fxo2 is unaffected by small changes in medium conductivity although a large increase, beyond ~300 mS/m, results in a collapsed DEP profile with no cross-over occurring. Fxo2 is most sensitive to the medium permittivity, a parameter which is dependent upon temperature but regulated precisely in these experiments [12]. The conductivity, permittivity and volume fraction of the nucleus are the most likely causes for the temporal change displayed in figure 8. A reduction in fxo2 of the magnitude observed requires the effective conductivity of the cell interior to reduce by a factor of nearly 40%. This represents a significant alteration in cell state, and could be directly related to the fact that for our DEP experiments the cells were suspended in low ionic strength solutions, commonly used for the determination of fxo1, for example [2, 6, 9].
Processes that could be responsible for the temporal reduction in fxo2 reported here could include the leakage of ions, such as potassium, down their concentration gradients to the cell exterior, or the osmotic flux of water across the plasma membrane [13]. Such changes would not influence measurements of the lower frequency DEP cross-over at fxo1, unless they lead to significant changes in cell volume, morphology or viability [2, 3, 9]. Determination of the temporal behaviour of fxo2 may prove to be a sensitive indicator of loss of cell viability, and a useful tool to be used when formulating the chemical composition of cell suspending media to be used in DEP studies [14].
Q10 is a widely used quantity in the study of metabolic activity and is the factor by which a 10°C increase in temperature alters a reaction rate [15, 16]. The temperature dependence results shown in figure 9 suggest Q10 is approximately 2.4 for the process in question, between 21°C and 37°C; and within the typical range for biological processes. Based on this data, the Arrhenius plot of figure 11 provides further detail regarding this process.
The linear fit of the data gives an activation energy of ~44 kJ/mol for the underlying process leading to the temporal behaviour of the high frequency DEP cross-over. We note that the temperature dependence of the coefficient of water self-diffusion through lipid multilayers is also described by the Arrhenius law, with an apparent activation energy of ~41 kJ/mol [17], a coincidence which requires further study.
Finally, the trend shown in figure 10, of a decrease in fxo2 with increasing osmolality of the cell suspending medium, suggests an associated decrease in internal cell conductivity and/or permittivity. The mechanism for this may be related to that of the aforementioned temporal reduction in frequency, however the initial rate of reduction appears to be unaffected by the change in osmotic pressure.
An electrical system comprising interdigitated electrodes and signal generator has been designed and tested to observe, for the first time, cells undergoing DEP cross-over in the predicted frequency range above 100 MHz. Electrical modelling of the electrodes based on a distributed RCR network allowed us to determine the limiting factors for high frequency operation. Capacitive loading is the most critical factor, and results in an imposed limitation of the electrode area. The mean fx02 value of 195 MHz observed for murine myeloma SP2/0 cells agrees well with that predicted using the cell interior dielectric properties derived from impedance studies of lymphocytes [5]. This suggests that the internal dielectric properties derived from lymphocytes, and the multi-shell model of a cell, may be applicable for the study of other mammalian cell types using the electrodes and electronics described here.
The theory indicates that the high frequency DEP cross-over at fx02 is sensitive to changes in the dielectric properties of the cell interior, and especially those associated with the nucleus and the relative size of the nucleus to the cell volume. An interesting finding, and one not predicted from the large number of reported studies of the low frequency DEP cross-over at fx01, is the rapid reduction of fx02 with time. At present we do not understand the changes of cell state that are responsible for this behaviour. Experiments to gain a further understanding are currently in progress in our laboratories. One important conclusion to be drawn is that measurements of both the low and high frequency DEP cross-over, at fx01 and fx02, are likely to provide a better indicator of any changes in cell state during DEP experiments than through determination of fx01 alone. This could be of particular importance for DEP studies and manipulations of stem cells [18].