## 1.

## Introduction

Cellular elastic properties are primarily related to their membrane composition and cytoskeletal architecture, which is remodeled by physiological and pathological events during the life cycle. The forces required to observe deformation range from
$\sim 10\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}1000\phantom{\rule{0.3em}{0ex}}\mathrm{pN}$
, depending on the cell type and the area of its application. Accordingly, measured values of cell elasticity vary over orders of magnitude depending on the cell type, health, and technique employed for measurement. For example, malignant fibroblasts have been found with half the elasticity of normal cells as determined by micropipette aspiration,^{1} whereas malignant epithelial cells have a 10-fold lower elastic modulus than normal cells as measured via scanning force microscopy.^{2, 3} Furthermore, metastatic cells taken from patients with suspected lung, breast, and pancreas cancer show 70% reduction in the elastic modulus as measured by atomic force microscopy (AFM).^{4} This study also found that the distribution of cell stiffness for the metastatic cells is much narrower than for healthy cells. Parasitic infection also influences cell mechanics—for example, red blood cells (RBCs) infected by *Plasmodium falciparum* are observed to be significantly stiffer than normal RBCs, as malaria parasites express proteins that stiffen the cell membrane.^{5, 6} Differences in the elastic modulus of normal RBCs and those infected with parasites at different growth stages have shown, in some cases, 10-fold increases in the elastic modulus.^{7, 8} These examples of disease-cell mechanics correlations suggest that cellular mechanical properties may provide a direct route to detecting disease or assessing its progression and motivates the development of a cytometry-based technique for measuring the distribution of elastic parameters.

Traditional bulk methods that subject cell suspensions to shear forces and extract a mean deformability with light scattering^{9, 10} provide average information for large numbers of cells but cannot resolve single-cell properties. In contrast, a range of methods including AFM,^{11, 12} micro-needles,^{13, 14} cell stretching between plates,^{15} micropipette aspiration,^{16} magnetic microbeads,^{17} intracellular methods,^{18, 19} and microfluidic flow-based methods,^{20, 21} can probe the mechanics of small numbers of individual cells. Most recently, optical trapping techniques have been used to manipulate cells in various ways, both with and without attached micro-beads.^{22, 23, 24, 25, 26, 27} In an approach closely related to the method presented here, cells trapped by divergent counterpropagating light sources^{28} have been elongated along the beam axis.^{28, 29} In this technique, however, each cell is sequentially trapped and then stretched, resulting in very low throughput; efforts to increase speed and perform measurements in microfluidic systems have attained
$10\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}100\phantom{\rule{0.3em}{0ex}}\text{cells}\u2215\mathrm{h}$
to date.^{30, 31} However, because both fluid drag on trapped cells and measurement rate are proportional to fluid velocity, extension of this method to flowing systems^{29} to increase throughput would require increasing laser power to optical intensities that result in cell damage.^{32, 33}

To achieve single-cell resolved measurements for large numbers (
${10}^{3}$
or more) of cells at data acquisition rates comparable to current cell cytometers, we have developed a different approach in which an anisotropic optical trap from a single emitter diode laser probes individual cells in a continuous microfluidic flow. In earlier work with linear traps from single-emitter diodes, we demonstrated that these compact light sources can impart significant forces^{34} appropriate for actuating flowing cells. Here, we present both experimental and computational results on simplified new geometries for high-throughput cell-mechanics measurements and first demonstrate that a single anisotropic trap can deform static cells. We then extend this approach to cells in microfluidic flow using long diode bars in a format whose throughput is limited only by the time scale of the cell elastic response. To illustrate the mechanism of force generation and to obtain the cell elastic parameters, we model the cell deformation as a function of applied optical trap power using the immersed boundary method (IBM)^{35} coupled with ray-optics. The IBM has been used extensively to simulate cell–fluid interactions^{35, 36, 37, 38} and cell interactions with ligand-coated surfaces^{37} and can be used to model the deformation in both optical and hydrodynamic force fields. In our studies, the diode laser induces optical forces on the cell surface due to the refractive index mismatch and refraction at the fluid/cell interface. This interaction deforms the cell until its elasticity balances the applied optical forces. In addition, the deformation causes fluid flow that leads to viscous stresses that influence the characteristic times for the induced stretch and its relaxation to equilibrium. This combination of experiment and simulation demonstrates that a single linear diode optical trap provides a convenient geometry for inducing deformation and obtaining the shear modulus of each cell from the observed equilibrium deformation at rates
$>1\phantom{\rule{0.3em}{0ex}}{\mathrm{s}}^{-1}$
.

## 2.

## Materials and Methods

## 2.1.

### Experimental Method

The experimental design is illustrated in Fig. 1
. Rectangular microfluidic channels
$15\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
deep and
$150\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
wide were fabricated in polydimethyl-siloxane (PDMS) using standard soft lithography techniques^{39} and bonded to glass using plasma surface activation. A pair of
$10\times $
(NA 0.25) objectives was used to image the diode laser into the microfluidic channel. A
$200\text{-}\mathrm{mW}$
, single-mode,
$830\text{-}\mathrm{nm}$
,
$1\times 3\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
laser diode was used for measurements on stationary cells, and a
$5\text{-}\mathrm{W}$
, multimode,
$808\text{-}\mathrm{nm}$
,
$1\times 200\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
single-emitter laser was used for measurements in flow. Applied laser power was measured at the system input.

For this initial study, we employ swollen RBCs because they lack internal structure and their elastic properties have been well characterized. In addition and because only the thin, outer phospholipid membrane contributes to the overall cell elasticity, they can be numerically modeled as a simple spherical elastic shell filled with a homogeneous fluid. Experimentally, samples consisted of heparinized bovine RBCs obtained from a slaughterhouse and used within
$24\phantom{\rule{0.3em}{0ex}}\mathrm{h}$
of acquisition. Samples were centrifuged at
$4000\phantom{\rule{0.3em}{0ex}}\mathrm{rpm}$
for
$10\phantom{\rule{0.3em}{0ex}}\mathrm{min}$
, decanted, and diluted in tris buffered saline at
$145\phantom{\rule{0.3em}{0ex}}\mathrm{mOsm}$
to induce swelling. Measurements on spherical cells simplify simulation of the applied optical forces^{40} and eliminate rotation of RBCs in the trap.^{41} RBCs were imaged using Koehler illumination though a
$40\times $
Carl Ziess objective and a
$640\times 480\phantom{\rule{0.3em}{0ex}}\text{pixel}$
CCD camera operating at
$100\phantom{\rule{0.3em}{0ex}}\text{frames}\u2215\mathrm{s}$
. Images were stored and analyzed by first normalizing and subtracting a background frame, smoothing with a
$9\times 9$
Gaussian kernel, and applying a threshold. Contours were then found using OpenCV,^{42} which was also used to determine spatial moments, the centroid, an equivalent ellipse, the ellipse orientation relative to the trapping bar, and the major and minor axes magnitudes corresponding to that ellipse. The changes in the minor
$\left(B\right)$
and major
$\left(A\right)$
axis width during and after stretch relaxation were used to calculate the percentage of minor stretch as
$(B-{B}_{0})\u2215{B}_{0}$
and of major stretch as
$(A-{A}_{0})\u2215{A}_{0}$
. Image analysis of stored video was performed on a computer with a
$2.66\text{-}\mathrm{GHz}$
quad core processor. Image analysis times are equivalent to video length at
$100\phantom{\rule{0.3em}{0ex}}\text{frames}\u2215\mathrm{s}$
and could be done in real time with straightforward program modification.

## 2.2.

### Numerical Method

In the IBM, a finite element model of the cell membrane is used to relate local membrane forces to local membrane deformation. Details of the numerical implementation and validation of the model can be found in previous work.^{36} Here, we modify the IBM to simulate cell deformation due to diode bar optical traps by imposing external optical forces from a ray-optics model of the
$808\text{-}\mathrm{nm}$
light source. We apply these optical forces to find the equilibrium deformation and then remove the forces to simulate the relaxation process. For swollen RBCs, we model the cell as a sphere of initial radius
$a=3\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
with an infinitely thin, hyperelastic, neo-Hookean membrane of negligible bending resistance, characterized solely by the membrane shear modulus
$\mathit{Gh}=3.4\times {10}^{-2}\phantom{\rule{0.3em}{0ex}}\text{dyne}\u2215\mathrm{cm}$
. In this model, the strain energy is given by
${W}_{h}=\mathit{Gh}({\lambda}_{1}^{2}+{\lambda}_{2}^{2}+{\lambda}_{1}^{-2}{\lambda}_{2}^{-2}-3)\u22152$
, where
${\lambda}_{1}$
and
${\lambda}_{2}$
are the planar principal strains. In our simulations, an unstressed spherical cell is first placed in an incompressible Newtonian fluid with the same density
$\rho =1\phantom{\rule{0.3em}{0ex}}\mathrm{g}\u2215{\mathrm{cm}}^{3}$
and viscosity
$\mu =0.8\phantom{\rule{0.3em}{0ex}}\mathrm{cP}$
as the cytoplasmic fluid. The index of refraction is 1.37 for the fluid inside the cell^{43} and 1.33 outside the cell. The fluid domain is a cube with a side that is 8 times the cell radius with periodic boundary conditions. The uniform grid used in our simulations has
${64}^{3}$
nodes, with a grid spacing of
$a\u22158$
, while the finite element cell grid has 20,480 triangular elements. A time step of
${10}^{-5}\phantom{\rule{0.3em}{0ex}}\mathrm{s}$
was used to ensure numerical stability.

To compare the magnitude of the simulated and measured deformations, we characterize the deformed cell with the Taylor deformation parameter
$\mathit{DF}=(A-B)\u2215(A+B)$
, where
$A$
and
$B$
are the lengths of the major and minor axes of elongated cells. Defined here for the
$x\text{-}y$
plane (Fig. 1),
$\mathit{DF}$
describes the geometrical deformation from perfect spheres
$(\mathit{DF}=0)$
to highly elongated morphologies. To extract the cell shear modulus from experimental data, we simulate the equilibrium deformation
${\mathit{DF}}_{\infty}\left({F}^{*}\right)$
as a function of dimensionless force
${F}^{*}$
for several different applied laser powers, thereby generating a standard curve that is used with measurements of
$\mathit{DF}={\mathit{DF}}_{\infty}-{\mathit{DF}}_{0}$
to extract a corresponding
${F}^{*}$
and
$\mathit{Gh}$
. In this, the dimensionless force is defined as
${F}^{*}={F}_{\mathit{\text{optical}}}\u22153a\mathit{Gh}$
, where
${F}_{\mathit{\text{optical}}}={n}_{1}PQ\u2215c$
, with
$P$
the laser power incident on the cell,
$c$
the speed of light in a vacuum,
${n}_{1}$
the refractive index of the surrounding medium, and
$Q$
a factor describing the amount of momentum transferred, calculated following Ref. 28 for spherical model systems but adapted here for linear optical sources.^{44}
$Q$
is independent of the power applied and depends only on object geometry and reflectance of the medium. In our simulations and under small deformation, we assume constant forces with a constant
$Q$
calculated as 0.011.

## 3.

## Results and Discussion

## 3.1.

### Experimental Results/Static Measurements

We first show in Fig. 2
that a single, anisotropic beam can be used to elongate individual RBCs trapped with a
$1\times 3\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
laser diode under static, no-flow conditions. Laser powers up to
$80\phantom{\rule{0.3em}{0ex}}\mathrm{mW}$
consistently lead to stretching on the order of 5 to 15%, where, at lower applied powers, a clear linear trend is observed. The change in slope and the large error bars at higher applied powers demonstrate a limit to the amount of power an RBC can experience before permanent deformation in a static optical trap if exposed for long times.^{33} If, however, stretching times are kept short, repeated measurements can be made on a single RBC at higher powers, as demonstrated in Fig. 3
and Video 1
. In these measurements, the undeformed RBCs are not spherical; thus, the laser induces both alignment and stretch along the long axis of the laser beam. The major stretch axis of the cell and the laser axis are not perfectly aligned because the
$3\text{-}\mu \mathrm{m}$
“long” axis of the laser is significantly shorter than the RBC diameter
$\left(6\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}8\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}\right)$
. Figure 4
and Video 2
quantify both the measured cell stretch and relaxation as the laser is repeatedly switched on and off. No deviation of measured stretch was seen in
$\sim 50\phantom{\rule{0.3em}{0ex}}\text{cycles}$
. A fit of the measured stretch relaxation to a decaying exponential over multiple measurements indicates that the time scales for the stretching and cell relaxation are on the order of
${\tau}_{r}\sim 0.10\phantom{\rule{0.3em}{0ex}}\mathrm{s}$
, as observed in previous studies.^{45, 46}

## 3.2.

### Experimental Results/Dynamic Measurements

Here, we present measurements of RBC stretching in rectangular microfluidic channels at flow speeds of
$\sim 50\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}\u2215\mathrm{s}$
with a
$1\times 200\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
emitter operated at
$3.2\phantom{\rule{0.3em}{0ex}}\mathrm{W}$
measured at the imaging plane outside the medium (Fig. 5
). By comparing the cell size before and after exiting the trap, the average stretch was found to be 2 to 3%. Due to conservation of mass, stretch along the minor axis is negative, while stretch along the major axis (parallel to the laser axis) is positive [Fig. 5], and as expected, cells tend to orient
$(\pm 5\phantom{\rule{0.3em}{0ex}}\mathrm{deg})$
along the axis of the trap [Fig. 5]. Simple geometric considerations assuming a constant ellipsoidal volume indicate an average cell deformation along the beam axis of 1.6%. Although this value is reasonable, validation of the constant volume assumption and direct measurement of cell deformation along this axis was not possible with the current experimental setup. One can also recast this data in terms of
$\mathit{DF}$
, as shown in Fig. 6
, where the stretch within the optical trap is seen as a cluster of data above the
${\mathit{DF}}_{\infty}={\mathit{DF}}_{0}$
line. In this single experiment, data for 1275 RBCs were taken in under
$20\phantom{\rule{0.3em}{0ex}}\mathrm{min}$
at a rate of
$\sim 1.2\phantom{\rule{0.3em}{0ex}}\text{cell}\u2215\mathrm{s}$
. This measurement rate is significantly faster than previous single-cell measurement demonstrations.^{29}

## 3.3.

### Computational Results

Next, we present computational results describing the force profiles and deformations generated by the trap laser. The prediction of net optical forces on cells and colloids has been described in the literature;^{47, 48} however, modeling of cell optical stretching forces has only recently been discussed.^{45, 49} In these models, a ray optics approach used to determine the optical forces on large spherical systems is extended to calculation of local stress profiles across the front and back sphere surfaces as a function of refractive index and incident laser beam profile. This approach is valid when the size of the object is much larger than the wavelength of the light,^{50} which is true here as the diameter of the erythrocytes ranges from
$6\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}8\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
. As the input light rays refract, a change in path leads to a change in the momentum carried by the light that is transferred to the interface through conservation. When the interface is an object, its surface absorbs this momentum, and a force proportional to the laser power is created by Newton’s second law. In our calculations, the parallel ray hits the front surface of sphere, and the refracted angle is found using Snell’s law, as detailed in Ref. 28. Multiple reflections are neglected. To integrate the forces imposed by optical stretching into the IBM, the optical forces are first calculated on the nodes of the 3-D membrane finite element grid. As the cell deforms, the applied optical forces are assumed to remain constant but move with the membrane nodes, an approximation valid for small deformations. If the front surface of the cell were significantly deformed, the changing refraction would broaden the stretching profile on the back side of the cell, leading to significant asymmetry in the force distribution. Although simulations that take this into account are certainly feasible, our experiments are performed at lower applied optical powers with smaller induced deformations, where the difference in cell optical forces on the cell is negligible.

Our simulations show that a single diode bar induces stretching forces on both the front and back sides of the cell, with an equilibrium deformation dependent on the applied power. In addition to the steady state deformation induced by the optical stretcher, we observe that a single diode bar produces a net translating force in the $z$ direction, pushing the cell away from the light source. This motion is apparent over the time scale of the applied stretching forces. In an unconfined system, one can estimate the expected $z$ translation away from the optical trap with $\Delta z\left(t\right)\sim {F}_{z}t\u22156\phantom{\rule{0.3em}{0ex}}\mu a$ . Using typical parameters, we estimate ${F}_{z}\sim 5.6\phantom{\rule{0.3em}{0ex}}\mathrm{pN}$ for applied $16\phantom{\rule{0.3em}{0ex}}\mathrm{mW}\u2215\mu \mathrm{m}$ powers, leading to a $z$ translation of $<10\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$ in the $0.1\text{-}\mathrm{s}$ measurement times necessary; however, in our highly confined systems, significant movement in the $z$ direction was not observed. In this case, translation out of the image plane was likely prevented by the presence of the wall, where, with a corresponding ${F}_{z}^{*}<0.02$ , any influence of the physical confinement on deformation is expected to be relatively small.

## 3.4.

### Cell Elasticity Cytometry

To prepare a general relation between measured deformation and applied laser power, a set of simulations was performed at varying applied forces [Fig. 7
] with a characteristic time for deformation of
${t}_{0}=\mu {a}^{2}\u2215{F}_{\mathit{\text{optical}}}$
and
${t}_{c}=\mu a\u22153\mathit{Gh}$
for stretch relaxation. These results demonstrate that cells rapidly stretch but reach a maximum deformation at longer times. Upon removal of the applied forces, the cell relaxes back to initial shape, where as expected, the equilibrium deformation increases with the applied force, as does the response time needed to reach equilibrium. We have tested this approach by simulating the deformation of model cells of different size and elastic modulus and found that the resulting deformations fall on the same master curve. Figure 7 shows the dimensionless force
${F}^{*}$
as a function of
$\Delta {\mathit{DF}}_{\infty}$
, a plot (curve fit) that allows for determination of the cell shear modulus
$\mathit{Gh}$
from the experimental laser diode measurements. In this, cells are stretched until equilibrium
$\Delta {\mathit{DF}}_{\infty}$
at a given power
$P$
and the corresponding
${F}^{*}$
found from the measured deformation, as illustrated in Fig. 7.
$\mathit{Gh}$
is then obtained from
$\mathit{Gh}={n}_{1}PQ\u22153ac{F}^{*}$
. Figure 8
demonstrates this approach on the data presented in Figs. 5 and 6 under high-throughput conditions, where the resulting scatter in
$\mathit{Gh}$
is apparent in the population of stretched cells and a median value of
$\mathit{Gh}=3.14\times {10}^{-3}\phantom{\rule{0.3em}{0ex}}\text{dyne}\u2215\mathrm{cm}$
is measured. This value falls within the range of values reported for the shear modulus of human RBCs of
$2.5\times {10}^{-3}\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}1.3\times {10}^{-2}\phantom{\rule{0.3em}{0ex}}\text{dyne}\u2215\mathrm{cm}$
^{22, 26, 28} measured under varying osmotic conditions. Note that for the individual measurements in which
$\Delta \mathit{DF}$
is negative, our approach results in a negative cell stiffness. Despite the inaccuracy inherent in each individual measurement using bright-field imaging of the cells, we obtain results that cluster around a physically reasonable median value. Higher contrast methods that resolve the cell sizes more precisely would eliminate these points. Furthermore, the reasonable agreement between our median value and published values of
$\mathit{Gh}$
implies that even relatively simple instrumentation with low-precision size measurements could be used to create practical devices.

The theoretical maximum measurement rate for stretching in flow is set by the response time of the cells. In particular, the cells must remain stretched for times within the trap on the order of the relaxation time constant ${\tau}_{r}$ , for which we use the value of ${\tau}_{r}=0.10\phantom{\rule{0.3em}{0ex}}\mathrm{s}$ determined in the static measurements of Fig. 4. Defining a characteristic cell dimension $\sigma $ that includes average spacing between cells allows prediction of maximum theoretical device throughput $\theta =L\u2215\sigma \tau $ . With a diode bar of length $L=200\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$ and using a characteristic spacing of $20\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$ , measurements of $100\phantom{\rule{0.3em}{0ex}}\text{cells}\u2215\mathrm{s}$ are achievable given the current system and experimental parameters. Significantly higher throughput will be achievable with stiffer cells, which have faster relaxation times, and by parallelizing the design, which is a desirable and achievable goal with the microfluidic platform used in this study.

## 4.

## Conclusions

In conclusion, we demonstrate that a simple asymmetric optical trap created by an inexpensive diode laser can be used to optically stretch blood cells in both static and dynamic, flowing environments. This approach is both practical and scalable, making it a useful technology for rapid measurements of cell deformability. Of particular advantage is that the method is optically based and is therefore inherently compatible with current optical cell-characterization technologies. In addition, we simulate the transient cell deformation induced by a single linear diode bar optical stretcher and compare the deformation and relaxation by varying laser power on cells with neo-Hookean membrane properties. Our simulations show that the forces imposed by a single diode bar optical stretcher can be used to both deform and translate cells and that a dimensionless single master curve relating laser power to cell deformation enables one to determine the cell elasticity of significant numbers of individual cells.

## Acknowledgments

The authors gratefully acknowledge support from the National Science Foundation under Grant No. DBI-045468, the National Institutes of Health under Grant No. RO1-AI063366, and from the Butcher Foundation. RJ is a staff member in the Quantum Physics Division of NIST.

## References

*Plasmodium falciparum*,” Proc. Natl. Acad. Sci. U.S.A., 104 (22), 9213 –9217 (2007). https://doi.org/10.1073/pnas.0703433104 0027-8424 Google Scholar

*Plasmodium falciparum*–infected erythrocytes,” Proc. Natl. Acad. Sci. U.S.A., 100 14618 –14622 (2003). https://doi.org/10.1073/pnas.2433968100 0027-8424 Google Scholar

*Escherichia coli*in optical traps,” Biophys. J., 77 2856 –2863 (1999). https://doi.org/10.1016/S0006-3495(99)77117-1 0006-3495 Google Scholar