1. Introduction
Cell screening is a method for observing living cells and evaluating cell physiology to discover cell properties and fate. This method has been widely employed for drug discovery [
1,
2], diagnostic purposes [
3,
4,
5], biological molecules detection [
4], and microbial species detection [
6,
7]. The traditional procedure for screening is based on microtiter plates, in which tests are conducted on 96-well plates. In order to reduce the cost and time of testing and enhance the efficiency, 384-, 1536-, and 9600-well plates were also introduced and employed for monitoring purposes. However, the drawbacks, such as evaporation, sample dispensing, and low amounts of samples, compelled the researchers to develop more efficient platforms [
8]. Lab-on-chip (LOC) microfluidic devices provide the context to monitor cell functions in real-time using a small amount of samples with high levels of controllability. Various microfluidic platforms have been introduced by researchers for cell monitoring, such as continuous-flow-based screening, droplet-based screening, and dielectric spectroscopy (DS) based screening [
8,
9]. Among these methods, DS is a non-invasive, fast, and promising method for cell screening. DS has been widely used for impedance measurements of yeast cells [
10,
11,
12], breast cancer cells [
13], lymphocyte activation profiling [
14], monitoring performance of epithelial tissues [
15], and hepatology research [
16].
DS-based cell screening often utilizes electrical impedance measurements of cell suspensions in the frequency domain. Cells’ dielectric properties, including cell membrane permittivity and cytoplasm conductivity, are extracted using various mathematical models, including the single-shell and double-shell models [
17]. Other models for cells adhered to electrodes and single cells in flow cytometry also exist [
18]. Reliable DS-based cell screening requires the minimization of cell–cell interactions, which were achieved by numerous methods, including the optical, acoustic, and thermophoretic tweezers [
19,
20,
21], manipulation in viscoelastic fluids [
22,
23,
24], and micro-well arrays [
25].
Recently, our group presented a microfluidic DS device that captures individual cells in microwells using positive dielectrophoresis (DEP). This is followed by dielectric measurements and cell release using negative DEP [
25,
26,
27]. The device allows dynamic alterations of the media so that the cell’s response to pH or osmolality variations or cell drug uptake can be measured in real-time. The ability of the device to release the cells after the measurements enable further biological studies [
25]. Using this system, we investigated the responses of PC-3 cells (a highly metastatic prostate cancer cell line) to temporal pH variations [
26] and the anticancer drug, Enzalutamide [
27]. The primary limitation of our previous research was the necessity of using low conductivity buffers (LCB) to avoid the electrode polarization (EP) effects. The EP effects increase with increased buffer conductivity and reduced electrode size, and these overwhelm the impedance data in the β-dispersion regime, where charging of the cell membrane happens [
25]. However, using LCB with conductivities that are ten folds lower than the typical cell growth media (CGM) has major consequences regarding cell shape regulation and pumping of ions out of the cells [
26]. Therefore, an ideal microfluidic DS device should be able to avoid the EP effects, but at the same time, enable accurate measurements in physiological buffers, such as the CGM.
EP arises from charge accumulations on the electrode–electrolyte interface, which is due to the charge carrier prevention to cross the solid–liquid barrier. As a consequence, capacitance on the interface increases, which results in a significant electric potential drop along with impedance rise. EP has a considerable impact on microfluidic devices because of the small electrode surface area. An increase in the interfacial surface area is an effective way to reduce EP, which was achieved by electrodepositing dendritic gold nanostructures on the electrodes’ surfaces [
28,
29]. In this method, oxidation and reduction occur to deposit gold from a solution of gold ions onto the electrode surface. This reduction procedure generates repeated patterns with a fern-leaf-type self-similarity [
30].
In this study, we demonstrate dielectric property measurements of PC-3 cells suspended in CGM as a function of time by using microelectrodes covered with dendritic gold nanostructures that minimize the EP effects. We also demonstrate a methodological approach for extracting the device and cell parameters from the equivalent circuit model. The reported cell membrane capacitance and cytoplasm resistance with 10 s increments show stable measurements for over 15 min, proving that the cells suspended in CGM were not affected by the microfluidic device and the dielectric measurements.
2. Materials and Methods
2.1. Device Fabrication
All chemical and photolithography supplies were purchased from Microchem Corp. (Westborough, MA, USA) and Sigma-Aldrich (St. Louis, MO, USA). To prepare all solutions, 18 MΩ·cm ultrapure water acquired from the Millipore Alpha-Q water system (Bedford, MA, USA) was used. To fabricate electrodes, first, the 2.5 cm × 2.5 cm glass slides were consecutively cleaned using an ultrasonic bath (FB11201, Fisher Scientific, Waltham, MA, USA) with 1 M KOH, acetone, isopropyl alcohol, and deionized (DI) water, for 10 min each at 25 °C and 37 kHz. Then, they were placed inside a 150 °C conventional oven to fully vaporize the water residues. The cleaned glass substrates were spin-coated with the positive photoresist (S1813). The rotation speeds of the spin-coater were 1000 rpm and 4000 rpm for 10 and 30 s, respectively, with 300 rpm/s acceleration/deceleration stages. After the substrates were soft baked for 1 min at 115 °C on a hot plate, and they were exposed to UV light through a transparency mask using a mask aligner (Karl Suss, MJB3) for 11 s. The substrates were submerged into a developer solution (MF-26A) for 17 s, and the UV exposed areas were washed away using DI water. Afterwards, the substrates were sequentially coated with chromium and gold layers for 60 mA-60 s setting using a sputter coater (EMS300TD, Emitech). The substrates were immersed into a PG remover to lift off the gold-coated photoresist regions.
The device in this study contains a pair of electrodes, of which only one has microwells on top of the electrode. To fabricate the microwells, a negative photoresist (SU8) was utilized on the 1 mm × 1 mm square shape gold electrode on the glass slide. An array of 21 × 21 microwells with 30 μm × 30 μm × 30 μm dimensions were formed on the sensing area. The dendritic nanostructures were fabricated on the bottom side of the microwells using the electrochemical deposition method in the optimal condition, as shown in our previous study [
29], at −0.7 volt for 60 min. A three-electrode potentiostat/galvanostat system was utilized for electrochemical deposition (EZstatPro, Nuvant, Crown Point, IN, USA). Nanostructures were electrochemically built on the planar gold electrodes in 0.5 mg/mL Sodium Tetrachloroaurate (III) (AuCl
4Na
2H
2O) solution (Sigma-Aldrich, Budapest, Hungary) using the deposition setup in the potentiostatic mode [
18].
The inlet and outlet fluidic holes were drilled on the glass slides using a diamond drill bit by a mini electric drill (Dremel 4000). A 70 μm thick double-sided tape was cut using a craft cutter (Silver Bullet) and used to make a microchannel. The microchannel was attached to the top electrode, and then the electrode pair were aligned using the mask aligner. The two electrodes were clamped to enhance adhesion between the slides, and copper tapes were connected to the electrodes for electrical connection. A schematic view of the microfluidic device is shown in
Figure 1a, while the images of the microwells and SEM pictures showing the dendritic structures on the electrodes are shown in
Figure 1b.
2.2. Experimental Setup
The experimental setup is shown in
Figure 1c. The PC-3 cells with 10
5 cells/mL concentration suspended in RPMI-1640 growth medium (Sigma-Aldrich) were injected into the channel from the inlet, and the cells were allowed to settle down in the wells for 5 min. This step was repeated three times, then the number of filled wells was counted under a microscope. To check the consistency of the gravity-based cell capture procedure, this method was applied to ten different devices. The mean value of the percentage of filled wells was 49.89% with ±2.60% standard deviation. The mean value corresponds to 220 filled wells out of 441 microwells. A syringe pump is used to supply growth media at a flow speed of 3 µL/min to wash the rest of the cells from the device. The cell suspension was collected from the outlet in a drain. The microfluidic chip was connected to the impedance analyzer (HP Agilent 4194A) to record the impedance spectrum. The impedance was measured in the 10 kHz–40 MHz frequency range at 20 mV AC voltage. Each impedance measurement sweep took six seconds for 314 discrete, logarithmically spaced frequencies, and the sweep was repeated every ten seconds. Then, the data from the impedance measurements were transferred to a PC and processed using a self-developed MATLAB code.
2.3. Cell Preparation
The PC-3 cancer cells were obtained from the Urology Lab at UT Southwestern, Dallas, Texas. The approximate diameter of the PC-3 cells is 22.0 ± 4 μm. The cells were cultured in RPMI 1640 supplemented with 5% fetal bovine serum (FBS), penicillin (100 IU/mL), and streptomycin (100 μg/mL). The cells were grown in an incubator (Thermo Scientific) at 37 °C with a 5% CO2 atmosphere. After 72 h, about 90% of the flask area was covered with cells, the growth medium was extracted, and the cells were washed with 1 × PBS, and 3 mL TrypLE was added to detach the cells from the surface. After incubating the cells with TrypLE for five minutes, 6 mL of complete growth medium was added, and the suspension was transferred to 15 mL centrifuge tubes. The cell suspension was centrifuged at 2000 rpm for four minutes, the supernatant was extracted, and a fresh medium was added to obtain 105 cells/mL suspension for the experiments.
3. DS Method
DS is a method used to characterize the dielectric properties of biological cells by applying small AC voltage over a desired frequency range. The impedance can be obtained from the measured response current of an applied voltage. An equivalent circuit model with fundamental electrical components is used to analyze the impedance measurements, and the dielectric properties of cells are extracted from the impedance spectrum.
Figure 1d shows the elements used for the equivalent circuit model for empty and filled microwells. We assume all EP effects to be identical for filled and empty wells, and these are characterized with the constant phase element model (CPE) with the impedance of Z
CPE. EP is experienced at the electrolyte–electrode interface. Accumulation of charge at the interface causes large potential drop and low measurement sensitivity. This effect becomes dominant in the low-frequency range at high voltages and in highly conductive buffers. This unfavorable behavior of EP is modeled using the CPE model, defined as follows:
where
and
are the exponent and CPE coefficient, respectively. The angular frequency is denoted by
, and
is the imaginary unit, where
. The parameter
varies between 0 to 1, where 1 represents purely capacitive and 0 corresponds to purely resistive behavior. The value of
increases and
decreases with increasing the surface area [
29].
Each filled channel has channel resistance, R
sol, while each empty and filled well has its own corresponding resistance, R
w,e and R
w,f, respectively. The microwells made out of SU8 are assumed to have capacitive behavior denoted with C
SU8, while the biological cells in each filled well are modeled using the cell membrane capacitance, C
mem, and cytoplasm resistance, R
cyt. The stray capacitance in the entire system is shown with C
f. It is important to note that each filled and empty well is assumed to be connected in parallel to each other. Details of this equivalent circuit model can be found in [
26]. The total impedance can be written as:
where
and
represent the impedance of empty and filled microwells, respectively, and n
f and n
e show the number of filled and empty microwells. Considering the equivalent circuit model, the impedance of empty microwells is written as:
The filled microwell impedance is given as:
To calculate
as an SU8 capacitance, the following equation can be used:
where
is the relative permittivity of SU8, which is 3.2, and
and t are the SU8 surface area and SU8 thickness, respectively.
Data Fitting Procedure
In order to find each electrical component from the circuit model and experimental data, real and imaginary parts of the impedance measurement were fitted to the real and imaginary parts of the equivalent circuit model using a nonlinear least-square method [
26]. Overall, the following parameters are fit to the experimental data:
,
,
,
, and
,
, and
.
The following procedure was used to find the unknown parameters, and an example of extracted parameters for each step is shown in
Table 1. First, the impedance of the growth media (GM) was measured. Then, the nonlinear least square method was applied to the GM data to find the unknown electrical components:
,
,
,
, and
(
Table 1, Step 1).
and
represent solution resistance and resistance of the empty well, which should remain constant for both GM and CGM measurements. Thus, the boundaries of these two parameters were selected to be in the range of the GM in all data fits. After finding these values, PC-3 cells were injected into the wells, which resulted in filling approximately 220 wells, and the impedance of the system was measured. The curved fit was obtained, and all unknown parameters were calculated. It was observed that the new
, from CGM were close to the values obtained for GM curve fitting, as expected before. However, due to the presence of cells in the wells, other parameters had different values compared with GM values (
Table 1, Step 2). This curve fit procedure for CGM was performed for all time steps, which proved that
,
, and
remained approximately constant through time (
Table 1, Step 3). Therefore, we can conclude that impedance data changes through time were mostly related to the dielectric properties of cells (
and
).
4. Results and Discussion
DS measurements are usually performed in kHz to MHz frequency ranges, and the EP effect results in the reduction in measurement sensitivity in sub-MHz frequencies, especially in high conductive buffers, such as growth media. In this study, in order to decrease this undesirable effect that can overshadow the biological impedance response, the electrode–electrolyte interfacial area is increased by electrodepositing dendritic gold nanostructures on electrode surfaces on the bottom of the microwells, using three electrodes electrochemical deposition.
Figure 1b shows an SEM image of the microwells array with a zoomed-in view of a single microwell and with details of the gold nanostructured electrode (GNE) morphology inside a well. In order to prove the benefits of the GNE electrode for the measurement of the biological cells, firstly, the impedance spectra of the growth media using both gold planar electrode (GPE) and GNE are measured, and the comparison is shown in
Figure 2a,b. According to this figure, the measured impedance range for GPE with growth media is 420 (Ω) in 10 kHz and 100 (Ω) in 40 MHz. However, with electrode surface area modification, the impedance in 10 kHz frequency dropped to 250 (Ω), while it stayed the same as the GPE at 40 MHz. This clearly shows that GNE decreases the EP effect in sub-MHz frequency ranges and enhances the usable frequency bandwidth for biological measurements.
As mentioned before, in order to extract cell properties, such as membrane capacitance and cytoplasmic resistance, the experimental results are fitted into an equivalent circuit model using the data fitting methodology presented in the previous section. Lines in
Figure 2a represent the fitted data on the experimental results for the device with growth media. However, the data fit error for GPE is higher than GNE (the R square for the corresponding fitted curve of GPE and GNE are approximately 0.8 and 0.99, respectively). The reason for data fit deviation from experimental results is the overshadowing of the impedance spectra due to the EP effect. Thus, the extracted parameters for GPE are not reliable. However, the results from the mathematical model fitted well on the experimental results of GNE due to the reduction in EP, which results in the extraction of acceptable parameters. In addition to the devices containing only growth media, the mathematical model is also applied to the acquired experimental data from GNE and GPE devices filled with PC-3 cells.
Figure 2c demonstrates an accurate fit between model and experiment in device with GNE, which exhibits EP reduction. The extracted parameters from all graphs in
Figure 2 are summarized in
Table 2.
One of the important parameters in the equivalent circuit model is the number of filled vs empty microwells. The number of filled wells in the GNE devices could not be directly measured due to the opaqueness induced by the dendritic electrodes. As a result, the gravitational fill studies, as discussed in the Experimental Setup section, were performed using identical sized devices without the dendritic electrodes. Systematic experiments have shown that typically 220 micro-wells out of 441 were filled on average with a standard deviation of approximately ±11.5 wells. In order to verify the importance of the number of filled vs empty wells on the circuit model, we fit the circuit parameters for 220/441 (~50%), 233/441 (~53%), and 208/441 (~47%) fill ratios. The extracted model parameters are shown in
Table 3. The effects of the well fill ratio on the cell membrane capacitance and cytoplasm resistance are 3% and 1%, respectively, while the filled well resistance varies less than 1%. These results prove the robustness of the equivalent circuit model.
One of the challenges in this study was the selection of channel width. Therefore, the effects of the channel width on the equivalent circuit model are also considered.
Figure 3 shows the impedance spectra obtained from different sized channels with their corresponding curve fits based on the circuit model. The length of the channels is 15 mm for all devices, while the widths of channels are 2.1, 1.7, and 1.4 mm for device 1, device 2, and device 3, respectively. The extracted parameters from the fitted impedance curves are also summarized in
Table 2. Cell properties such as membrane capacitance and cytoplasm resistance are merely independent of the channel width. Additionally, the solution resistance and resistance of empty wells from the three different channels are almost equal to each other, as expected from a consistent model. These results confirm that different channel sizes do not impose an effect on the extracted properties of the cells. This also shows the consistency of the electrical circuit model and the sequential data fitting procedure explained in the previous section. In addition to the cell properties,
Table 2 also includes K and α values, which represent the CPE coefficients of interfacial impedance. According to this table, the electrode with a dendritic structure has a lower α (lower capacitance effect) compared with GNE, which shows DS measurement sensitivity enhancement in sub-MHz frequencies. However, GPE has higher α (higher capacitance effect) and lower K indicating the higher EP effect. It should also be noted that the values of K and α for the GPE electrode lacks accuracy, as shown in
Figure 2. Moreover, by increasing the width of the channel, α decreased and K increased due to the reduction in the total impedance of the system.
A significant characteristic of a promising microfluidic device for cell screening is its potential of maintaining the cells in a stable condition during the measurements. To verify this capability of the device, the membrane capacitance and cytoplasm resistance of PC-3 cells are measured every 10s for 1000 s while constantly feeding the system with high conductive growth media (RPMI-1640) at 37 °C.
Figure 4 shows the stability of the extracted cell parameters as a function of time with C
mem = 100.30 ± 0.87 pF and R
cyt = 100.0020 ± 0.0001 kΩ. The results show that constant cell properties can be maintained in the device for over 16 min. We extended this period to 2.5 h and have not observed variations in the cell properties (data not shown for brevity).
The employed frequency range to extract cell properties in this study was 10 kHz−40 MHz. The fitted data in the previous figures and tables to find were calculated using 314 data points distributed logarithmically. In order to reduce the measurement time, the number of data points for property extraction was reduced to 88, 46, and 27 in the mathematical model. For these three data point numbers, the frequency ranges of 1MHz–10 MHz, 3MHz–10 MHz, and then 5 MHz–10 MHz were considered, respectively. The values of Rsol, Rw,e, K, α, Cf were obtained from the 10 kHz–10 MHz range of frequency using the 314 data points, and these parameters were assumed to be constants. Thereafter, values were found for the three mentioned frequency ranges. The results proved that using 314, 88, 46, and 27 data points to find the, parameters led to the same results. As a consequence, a smaller range of frequency with a lower number of data points can be measured during the experiment to find the cell properties, which means lower measurement time and faster data analysis.