Specific membrane capacitance, cytoplasm conductivity and instantaneous Young’s modulus of single tumour cells

Abstract

As label-free biomarkers, biophysical properties of cells are widely used for cell type classification. However, intrinsic biophysical markers, e.g., specific membrane capacitance (C_{specific membrane}), cytoplasm conductivity (σ_conductivity) and instantaneous Young’s modulus (E_{instantaneous}) measured for hundreds of single cells were not yet reported. In this study, single cells in suspension (adherent cells treated with trypsin) were aspirated through a microfluidic constriction channel at 25 °C, and the entry processes and impedance profiles were recorded and translated to C_{specific membrane}, σ_conductivity and E_{instantaneous}. C_{specific membrane}, σ_conductivity and E_{instantaneous} of five cell types were quantified as 2.10±0.38 μF cm⁻², 0.91±0.15 S m⁻¹ and 5.52±0.95 kPa for H460 cells (n_cell=437); 2.52±0.54 μF cm⁻², 0.83±0.12 S m⁻¹ and 5.54±1.04 kPa for H446 cells (n_cell=410); 2.45±0.57 μF cm⁻², 0.99±0.18 S m⁻¹ and 5.16±1.68 kPa for A549 cells (n_cell=442); 1.86±0.31 μF cm⁻², 1.07±0.18 S m⁻¹ and 3.86±0.81 kPa for 95D cells (n_cell=415); 2.03±0.35 μF cm⁻², 0.99±0.16 S m⁻¹ and 3.49±0.70 kPa for 95C cells (n_cell=290). The database of C_{specific membrane}, σ_conductivity and E_{instantaneous} may serve as a reference for future studies of cellular biophysical properties.

Design Type(s)	cell type comparison design
Measurement Type(s)	Biophysical Chemistry
Technology Type(s)	microfluidics platform
Factor Type(s)	cancer cell line
Sample Characteristic(s)	H460 cell • H446 cell • A549 cell • 95D cell • 95C cell

Machine-accessible metadata file describing the reported data (ISA-Tab format)

Background & Summary

Biophysical properties of single cells include electrical parameters such as specific membrane capacitance (C_{specific membrane}, i.e., capacitance per unit area for the cell membrane modelled as a capacitor) and cytoplasm conductivity (σ_conductivity, i.e., the reciprocal of electrical resistivity or specific electrical resistance for the cytoplasm modelled as a resistor) as well as mechanical parameters such as instantaneous Young’s modulus (E_{instantaneous}, which describes the initial elastic response of a cell, modelled as an incompressible viscoelastic solid, in response to mechanical forces), which are indicators of the status of cytoskeletons and cellular membranes¹. Variations in cellular biophysical properties are closely related to physiological and pathological processes, examples including (1) red blood cells infected by malaria; (2) tumour cells during migration and evasion; (3) leukocytes affected by sepsis; and (4) stem cells undergoing differentiation^2,3.

Conventional techniques for characterizing the electrical properties of cells include patch clamping, electrorotation and dielectrophoresis^4–7. On the other hand, atomic force microscopy, micropipette aspiration and optical tweezers have been used for characterizing the mechanical properties of cells^8–12. Although well established, these approaches suffer from limited throughput (~1 cell per minute). In addition, they are not capable of measuring cellular electrical and mechanical properties simultaneously¹. Therefore, data recording C_{specific membrane}, σ_conductivity and E_{instantaneous} for multiple cell types (hundreds of cells for each cell type) are currently not available in the literature.

Microfluidics is an area focusing on processing small amounts of fluids at the nanoliter to picoliter scale using microfabricated channels with dimensions of tens of micrometers^13–15. The micrometer dimension is comparable to biological cells, and thus microfluidics is well suited for single-cell analysis^16,17. More specifically, advances in microfluidic technology have enabled high-throughput characterization of the biophysical properties of single cells, e.g., through impedance flow cytometry, microfluidic optical stretchers and microfluidic hydrodynamic stretchers^1,8,18. Although powerful, these previously reported microfluidics-based approaches only reported electrical (e.g., impedance flow cytometry) or mechanical (e.g., microfluidic optical stretchers and hydrodynamic stretchers) properties of single cells separately, without demonstrating simultaneous characterization of cellular electrical and mechanical properties.

Microfluidic approaches that enable simultaneous characterization of electrical and mechanical properties of single cells have also been developed, specifically through (1) microcantilever-based electrodes¹⁹, (2) electrodeformation²⁰, (3) micropipette aspiration with impedance spectroscopy²¹, or (4) constriction channels with impedance spectroscopy²². The first three approaches are limited by their low throughput, and thus cannot be used to collect data from hundreds of single cells. In the fourth approach, cells in suspension were aspirated through a constriction channel (with a width and a height that are smaller than the cell diameter) continuously. The deformation and impedance profiles of individual cells as they enter and travel through the constriction channel were recorded as mechanical and electrical signals, respectively, with an estimated throughput of 1 cell per second.

Based on custom-developed electromechanical models, these raw biophysical signals obtained from constriction channels with impedance spectroscopy were then translated to size-independent intrinsic biophysical markers including C_{specific membrane}, σ_conductivity and E_{instantaneous}^23–25, enabling the classification of tumour cells with different malignant levels²⁶. However, as proof-of-concept demonstrations, only a limited number of cell types were examined with small populations for each cell type²⁴.

In this study, based on the aforementioned approach, C_{specific membrane}, σ_conductivity and E_{instantaneous} from five types of tumour cells and hundreds of single cells for each cell type were quantified and reported. This study provides a preliminary database for cellular C_{specific membrane}, σ_conductivity and E_{instantaneous}, which may serve as a reference for future studies on characterizing and classifying biological cells based on cellular biophysical properties.

Methods

Working flowchart

The working flowchart for characterizing C_{specific membrane}, σ_conductivity and E_{instantaneous} includes four key steps: device fabrication, cell preparation, device operation and data processing (see Fig. 1). During operation, cells in suspension were aspirated into the microfluidic constriction channels with the deformation and impedance profiles of the cells recorded by a high-speed camera and an impedance analyser, respectively. Raw biophysical data were obtained by processing the images and impedance data captured in experiments, which were then translated to C_{specific membrane}, σ_conductivity and E_{instantaneous}, based on a theoretical electrical model for a cell traveling within the constriction channel and a numerical mechanical model capturing the deformation of a cell when it enters the constriction channel, respectively. The key steps were summarized as follows. The detailed procedures have been reported in a previous publication²⁴.

Figure 1: Working flowchart for simultaneous characterization of C_{specific membrane}, σ_conductivity and E_{instantaneous} of single cells.

Key steps include device fabrication (a), cell preparation (b), device operation (c) and data processing (d) leveraging developed electrical (e) and mechanical (f) models. During device operation, suspended cells were aspirated into the microfluidic constriction channels with the cell deformation and impedance profiles recorded by a high-speed camera and an impedance analyser, respectively. Preliminary biophysical markers including A_{1 kHz}, A_{100 kHz}, L_{instantaneous}, L_transitional and L_elongation were obtained based on processing of image and impedance data, which were then translated to C_{specific membrane}, σ_conductivity and E_{instantaneous}, based on a theoretical electrical model for cells traveling within the constriction channel and a numerical mechanical model capturing cell deformation during the entry process, respectively.

Device fabrication

The microfluidic device consists of a constriction channel (a cross sectional area of 10 μm×10 μm) in polydimethylsiloxane (PDMS) elastomer (Dow Corning Corp., Midland, MI, USA) which was replicated from a SU-8 (MicroChem Corp., Newton, MA, USA) mould master (see Fig. 1a).

Briefly, SU-8 5 was spun coated, prebaked and exposed without development and post exposure bake to form the layer of the constriction channel (10 μm in height). Then SU-8 25 was spin coated on top of the first SU-8 layer, exposed with alignment and developed, forming the cell loading channel with a height of 25 μm. PDMS precursor and curing agents (10:1 in weight) were mixed, poured on channel masters and baked for crosslinking. PDMS channels were then peeled away from the SU-8 masters, punched to form through holes as inlets and outlets, and bonded to glass slides after plasma treatment.

Cell preparation

All cell-culture reagents were purchased from Life Technologies Corporation (Carlsbad, CA, USA). The lung cancer cell lines of H460, H446, A549, 95D and 95C (China Infrastructure of Cell Line Resources, Beijing, China) were cultured at 37 °C in 5% CO₂ in RPMI 1640 medium (11875) supplemented with 10% heat-inactivated fetal bovine serum (10099), 100 units ml⁻¹ penicillin and 100 μg ml⁻¹ streptomycin (15140). Immediately prior to an experiment, cells were trypsinized (25200, 0.25% for 3 min) to form a solution of cell suspension at a concentration of 1×10⁶ cells ml⁻¹ (see Fig. 1b).

Device operation

During operation the device was first filled with culture medium and and a pipette was used to transfer the cell suspension solution to the entrance of the cell loading channel. A pressure calibrator (DPI-610 pressure calibrator, Druck, Billerica, MA, USA) was used to generate a negative pressure (0.5–1.2 kPa), aspirating cells continuously through the constriction channel with silver wires inserted into the inlet and the outlet of the device for impedance profile recording (see Fig. 1c).

An inverted microscope (IX71, Olympus Inc., Tokyo, Japan) in connection with a high-speed camera (M320S, Phantom Inc., Bublin, OH, USA) was used to capture the process of cellular entry and travelling in the constriction channel at 200 frames per second. Impedance data at both 1 and 100 kHz was measured by a lock-in amplifier (7270, Signal Recovery, Oak Ridge, TN, USA) with a sampling rate of 25 points per second. All the characterization experiments were conducted within 30 min trypsinization of the cells at the room temperature (25 °C).

Data processing

At the stage when a cell enters the constriction channel, two preliminary parameters were quantified based on image processing: L_{instantaneous} as the aspiration length when a cell instantaneously jumps into the channel and L_transitional as the aspiration length at the point where the creep deformation of the cell ends (see Fig. 1d)²⁵. During the stage when the cell is travelling within the constriction channel, the impedance profiles and images of the elongated cell were analysed which produced three preliminary parameters: A_{1 kHz} and A_{100 kHz}, i.e., the ratios between the impedance amplitudes with and without a travelling cell at 1 and 100 kHz, respectively, as well as L_elongation, i.e., aspiration length of the cell when travelling within the constriction channel (see Fig. 1d)²³.

A theoretical model²³ was previously developed to model the electrical response of a cell as it is traveling within the constriction channel, thereby enabling the conversion of the measured A_{1 kHz}, A_{100 kHz} and L_elongation to C_{specific membrane} and σ_conductivity as size-independent intrinsic electrical markers for single cells. Briefly, the electrical response of a cell was represented by R_cytoplasm and C_membrane in series where R_leak indicated sealing properties during the cellular travelling process in the constriction channel. R_leak was derived from impedance data at 1 kHz (A_{1 kHz}) and C_membrane and R_cytoplasm were derived from impedance at 100 kHz (A_{100 kHz}), respectively, which, were further translated to C_{specific membrane} and σ_cytoplasm based on geometrical information of the constriction channel (see Fig. 1e). Note that in this study, the lumped electrical model was used for data interpretation where the membrane portion was represented by C_membrane and the cytoplasm portion was represented by R_cytoplasm, which, to an extent, neglects the potential effects of ion channels in the cell membrane and the membranes of cytosolic organelles in cytoplasm. Future studies may develop more accurate cellular electrical models.

On the other hand, a numerical mechanical model²⁵ was developed to model the cell deformation as it enters the constriction channel, where the channel walls were modelled as rigid surfaces and the cell was modelled as an incompressible viscoelastic solid. The key mechanical parameter to be extracted was E_{instantaneous}. Based on numerical simulations, the relationships between L_{instantaneous}, L_transitional and E_{instantaneous} were obtained as follows:

$$\begin{array}{l}{L}_{ins\tan \tan eous}/D_{channel}=\left(44.27f_c^2-37.24f_c+13.70\right)\times P_{aspiration}/E_{ins\tan \tan eous}+\left(-5.31f_c^2+2.84f_c-0.59\right)\\ L_{transitional}/D_{channel}=\left(-60.40f_c^2+40.05f_c-8.68\right)\times P_{aspiration}/E_{ins\tan \tan eous}+\left(1.99f_c^2+0.03f_c+1.60\right)\end{array}$$

where D_channel represents channel geometrical information, ƒ_c represents the friction on cell-wall interfaces and P_aspiration represents pressure applied to aspirate cells into the constriction channel (see Fig. 1f). The two unknown parameters, i.e., E_{instantaneous} and ƒ_c, were solved from these two equations since all the other parameters can be measured from experiments.

Data Records

Raw impedance and image data of five cell types are available at Dryad repository (Data Citation 1, Raw Impedance and Image Data of Single Tumour Cells). The file is an Excel composed of five sheets named after five cell lines of H460, H446, A549, 95D and 95C, respectively. In each sheet of a specific cell type, the first row is the cell number and within each cell number, there are seven columns which are time (collected from impedance analyser), amplitude at 1 kHz, phase at 1 kHz, amplitude at 100 kHz, phase at 100 kHz, time (collected from high-speed camera) and aspiration length.

C_{specific membrane}, σ_conductivity and E_{instantaneous} of five cell types are available at Dryad repository (Data Citation 1, Biophysical Properties of Single Tumour Cells). The file is an Excel composed of five sheets named after five cell lines of H460, H446, A549, 95D and 95C, respectively. In each sheet of a specific cell type, there are nine columns which are cell number (column A), A_{1 kHz} (column B), A_{100 kHz} (column C), L_{instantaneous} (column D), L_transitional (column E), L_elongation (column F), C_{specific membrane} (column G), σ_conductivity (column H) and E_{instantaneous} (column I). In this file, each row represents the biophysical parameters of a specific cell, including A_{1 kHz}, A_{100 kHz}, L_{instantaneous}, L_transitional, L_elongation, C_{specific membrane}, σ_conductivity and E_{instantaneous}.

Technical Validation

Figure 2 shows a sequence of microscopic images (a–d), showing the entry process into the constriction channel for a single cell. In addition, the impedance data for the same cell was shown in Fig. 2e and the aspiration length as a function of time was obtained by image processing and shown in Fig. 2f. The detailed data for individual cells are included in (Data Citation 1, Raw Impedance and Image Data of Single Tumour Cells).

Figure 2: Raw imaging and impedance data for a typical cell to enter and travel in the constriction channel (length: 200 μm, width: 10 μm and height: 10 μm) with preliminary biophysical parameters quantified.

(a–d) Microscopic pictures of a cell’s entry and travelling process in the constriction channel with raw impedance data of the same cell shown in e and processed aspiration length versus time shown in f. Based on data processing, five preliminary biophysical parameters including A_{1 kHz}, A_{100 kHz}, L_{instantaneous}, L_transitional, and L_elongation were obtained.

Based on the processing of impedance profiles and image processing, five preliminary biophysical markers of A_{1 kHz}, A_{100 kHz}, L_{instantaneous}, L_transitional, and L_elongation were obtained. When the cell entered the constriction channel, an instantaneous jump into the constriction channel was initially observed (see Fig. 2a), enabling the quantification of L_{instantaneous} in Fig. 2e. Then, a gradual increase in aspiration length (see Fig. 2b,e) was observed due to viscoelastic creep. The creep deformation ended when L_transitional was reached (see Fig. 2c,e). Figure 2d shows an image of the cell travelling within the constriction channel where L_elongation was quantified. Furthermore, impedance ratios at A_{1 kHz} and A_{100 kHz} were also quantified during the cellular travelling process in the constriction channel (see Fig. 2e). The averages and standard deviations of five preliminary biophysical markers (A_{1 kHz}, A_{100 kHz}, L_{instantaneous}, L_transitional and L_elongation) for five types of tumour cells were summarized in Table 1. The detailed data for individual cells are included in (Data Citation 1, Biophysical Properties of Single Tumour Cells).

Table 1 A summary of biophysical parameters of five types of tumour cells.

Figure 3a–e show the distributions of C_{specific membrane}, σ_conductivity and E_{instantaneous} for five types of tumour cells, respectively. For each type of tumour cells, the scatter plots of (i) C_{specific membrane} versus D_cell, (ii) σ_conductivity versus D_cell, (iii) E_{instantaneous} versus D_cell and (iv) C_{specific membrane} versus σ_conductivity versus E_{instantaneous} were included. Note that D_cell represents the diameter of the cell under measurement, which was calculated from L_elongation based on the assumption of volume conservation when the cell deforms. As shown in Table 1, C_{specific membrane}, σ_conductivity and E_{instantaneous} of five cell types were found to be 2.10±0.38 μF cm⁻², 0.91±0.15 S m⁻¹ and 5.52±0.95 kPa for H460 cells (n_cell=437); 2.52±0.54 μF cm⁻², 0.83±0.12 S m⁻¹ and 5.54±1.04 kPa for H446 cells (n_cell=410); 2.45±0.57 μF cm⁻², 0.99±0.18 S m⁻¹ and 5.16±1.68 kPa for A549 cells (n_cell=442); 1.86±0.31 μF cm⁻², 1.07±0.18 S m⁻¹ and 3.86±0.81 kPa for 95D cells (n_cell=415); 2.03±0.35 μF cm⁻², 0.99±0.16 S m⁻¹ and 3.49±0.70 kPa for 95C cells (n_cell=290). The detailed data of C_{specific membrane}, σ_conductivity and E_{instantaneous} for individual cells are included in (Data Citation 1, Biophysical Properties of Single Tumour Cells).

Figure 3: Scatter plots of C_{specific membrane}, σ_conductivity and E_{instantaneous} of five tumour cells.

The number of measured cells for each type was 437 for H460 (a), 410 for H446 (b), 442 for A549 (c), 415 for 95D (d) and 290 for 95C cells (e). In addition, for each cell type, scatter plots of (i) C_{specific membrane} versus D_cell, (ii) σ_conductivity versus D_cell, (iii) E_{instantaneous} versus D_cell and (iv) C_{specific membrane} versus σ_conductivity versus E_{instantaneous} were included.