The Influence of Electric Field and Confinement on Cell Motility

The ability of cells to sense and respond to endogenous electric fields is important in processes such as wound healing, development, and nerve regeneration. In cell culture, many epithelial and endothelial cell types respond to an electric field of magnitude similar to endogenous electric fields by moving preferentially either parallel or antiparallel to the field vector, a process known as galvanotaxis. Here we report on the influence of dc electric field and confinement on the motility of fibroblast cells using a chip-based platform. From analysis of cell paths we show that the influence of electric field on motility is much more complex than simply imposing a directional bias towards the cathode or anode. The cell velocity, directedness, as well as the parallel and perpendicular components of the segments along the cell path are dependent on the magnitude of the electric field. Forces in the directions perpendicular and parallel to the electric field are in competition with one another in a voltage-dependent manner, which ultimately govern the trajectories of the cells in the presence of an electric field. To further investigate the effects of cell reorientation in the presence of a field, cells are confined within microchannels to physically prohibit the alignment seen in 2D environment. Interestingly, we found that confinement results in an increase in cell velocity both in the absence and presence of an electric field compared to migration in 2D.


Introduction
The asymmetric distribution of ion channels and pumps between the apical and basal surfaces of the endothelial cells surrounding most organs leads to a transendothelial (or transepithelial) potential difference DQ ( = Q apical -Q basal ) of +15 to +60 mV, corresponding to a dc electric field of 0.5-5 V cm 21 [1][2][3]. This is a relatively small field, about six orders of magnitude lower than the threshold field for electroporation of a cell membrane (approx. 2610 6 V cm 21 ) [4]. However, epithelial and endothelial cells are programmed to sense and respond to dc electric fields (dcEFs) in processes such as wound healing, development, and nerve regeneration [1][2][3][5][6][7][8]. Electric fields are also thought to play a role in angiogenesis [3] and metastasis [1].
In cell culture, electric fields influence cell division, polarity, shape, and motility. Many cell types respond to dcEFs, preferentially migrating either to the anode or cathode, a process known as galvanotaxis [2]. In vitro studies of galvanotaxis are usually performed in 2D by analyzing the path traced by individual cells in the presence or absence of an electric field. Most cell types respond to dcEFs of the magnitude of endogenous electric fields, however, the origin of this directionality and the mechanism of galvanotaxis remain unknown.
For many processes of physiological interest cell motion is physically confined. For example, during migration through the extracellular matrix, cells migrate along channels between bundled collagen fibers [9]. Similarly, during wound healing, cells must squeeze between other cells. Various cell migration chambers have been developed for the study of galvanotaxis in 2D [10][11][12][13][14][15], including microfluidic-based platforms [12,14,15]. Our objective is to develop a versatile microfluidic-based platform to study galvanotaxis in 2D and confined geometry.
Here we report on quantitative characterization of the physical and morphological aspects of the motility of NIH 3T3 fibroblasts under an electric field and physical confinement. Although galvanotaxis is usually associated with directional bias towards the cathode or anode we show here that the influence of dcEFs on motility is much more complex. In 2D (no confinement) and in an electric field, cells orient perpendicular to the field vector and migrate preferentially towards the cathode. Unexpectedly we show that the electric field exerts forces on the cells both parallel and perpendicular to the field. These forces are in competition with each other and ultimately govern the trajectories of the cells in the presence of an electric field. At low field, the cells migrate preferentially towards the cathode, however, the perpendicular component of the individual segment vectors is larger than the parallel component. In a larger field, there is a significant increase in average velocity and the parallel component of the individual segment vectors is larger than the perpendicular component. These results suggest that there could be at least two independent signaling pathways that influence cell motility in an electric field. To further probe the effect of perpendicular alignment on directed migration induced by the electric field, 3T3 cells were seeded inside 20 mm channels to physically prevent cell orientation during galvanotaxis. We found that physical confinement results in an increase in cell velocity both in the absence and presence of an electric field compared to migration in 2D. This result could be relevant in understanding galvanotaxis in vivo.

Microfluidic Platform
There are two key features of our platform (Fig. 1). (1) Silver chloride electrodes are integrated into the platform without the need for external salt bridges and solution reservoirs. By integrating the electrodes into the microfluidic platform, our device is sufficiently small to fit in a live cell chamber (Figs. 1D and 1E) so that experiments can be performed under controlled relative humidity, CO 2 partial pressure, and temperature. As described below, by using silver chloride electrodes, the only species involved in the electrochemical reactions at the electrode surface is the consumption or generation of chloride ions. In this way, we can avoid the parasitic reactions associated with noble metal electrodes such as gold or platinum. (2) By controlling the width w of the channels relative to the cell diameter d cell we can study cell migration under an electric field in 2D (w..d cell ) or in confined dimensions (w#d cell ). Such quasi-1D migration is analogous to migration along fibers in 3D gels, such as extracellular matrix [16,17]. Confinement experiments were performed in 20 mm, 50 mm, and 100 mm wide channels. Here we compare cell migration in 2D (w6h, 1000 mm680 mm) and in 20 mm channels (w6h, 20 mm680 mm); experiments in 50 mm and 100 mm channels showed intermediate behavior and hence did not provide any additional insight. The details of platform fabrication are provided in Supporting Information.

Fabrication of Ag/AgCl Electrodes
The electric field was applied between two Ag/AgCl electrodes located at each end of the microfluidic channel ( Fig. 1). Ag/AgCl electrodes have the advantage that the equilibrium at both electrodes is established by the reaction AgCl+e 2 « Ag+Cl 2 and that the current is carried by chloride ion transport. The electrodes were fabricated by electrochemically forming AgCl on two inch lengths of 0.025 inch diameter silver wire (A-M systems, Sequim, WA) using standard procedures [18]. Briefly, the silver wire and a platinum foil cathode were immersed in 1M HCl solution and chloridized for 30 minutes at a current of 5-10 mA cm 22 . The wires were then removed from the HCl solution, rinsed and stored in distilled water (Fig. 1F). The electrodes were coiled to ensure that about two inches was embedded in the agarose gel (see Supporting Material for details). From scanning electron microscope images (Fig. 1F), we determine an average AgCl thickness of 20 mm. A current of about 75 mA is required to maintain a field of about 5.5 V cm 21 within a 1 cm long and 1000 mm wide channel (cross sectional area of 0.08 mm 2 ), and from Faraday's law, we determine that a 5.5 V cm 21 field can be maintained for more than six hours. Experiments were performed in the absence of a field, and with an electric field of 2.2 V cm 21 or 5.5 V cm 21 . The magnitude of the field was confirmed using a four point probe method with two platinum wires inserted at each end of the channel.

Maintenance of Cell Lines
3T3 mouse fibroblast cells (ATCC) were cultured in DMEM media with 10 vol.% calf bovine serum (CBS) and 1 vol.% penicillin-streptomycin (PenStrep). Prior to experiments, cells were washed with phosphate buffer (PBS) and lifted from the surface using trypsin (0.25% trypsin-EDTA in PBS, Sigma). Suspended cells were centrifuged and re-suspended in 500-1000 mL of media before introducing into the platform through the cell injection port using a syringe. The cells were introduced at a density of 300,000 cells mL 21 .

Imaging
Time-lapse phase-contrast images were collected using an inverted microscope with a 10X objective (TiE2000, Nikon, Melville, NY). A live-cell chamber was used to control temperature (37uC), CO 2 (5 vol.%), and relative humidity (95%) (Fig. 1E). The device was located in a custom holder placed on an automatic stage in the live-cell chamber. The holder was customized with electrical feed-throughs to connect the electrodes to an external power supply (Fig. 1E). Time-lapse images were taken at 5 min intervals for three hours using NIS-element software (Advanced Research Edition 3.22.11, Nikon, Melville, NY).

Image Analysis
Time-lapse images were stacked using ImageJ (Version 1.45, NIH, Bethesda, MD) and loaded into Metamorph (Version 7.7, Molecular Devices, Sunnyvale, CA) for analysis of cell paths. Each cell path is made up of individual segments determined by the time between images (Dt = 5 min). Cell paths were only analyzed if the cells did not divide, were not in contact with other cells, and did not leave the field of view. The cell position in each frame was determined from the position of the nucleus. The cell path data was then analyzed using Matlab (MathWorks, Natick, MA). For each cell we determine the persistency, the cell velocity in each segment, the angular orientation of a segment with respect to the field vector, the average directedness, the segment turn angle, and the orientation of the cell with respect to the field vector ( Fig. 2E-G).
(1) The velocity v(t) is determined from the displacement Dz i and elapsed time between successive images (each segment along the cell path). (2) The persistency (or directional persistence) p is a measure of the linearity of a path and is determined from r/gDz i where r is the overall displacement (distance between the start and end points) and gDz i is the length of the path (i.e. the sum of the individual segments along the path). The persistency varies between 0 and 1: p = 0 corresponds to no displacement whereas p = 1 corresponds to a linear path (maximum bias). This term is also known as the directional persistence and in studies of chemotaxis is known as the chemotactic index (CI) [19]. Note that the terminology can be misleading since the persistency does not indicate whether a path is aligned with an external bias, such as electric field or chemoattractant. (3) The directedness d i is a measure of the directional bias of each segment and is given by d i = cos(Q i ), where Q i is the angular orientation of a segment vector ( z z i ) with respect to the field vector and 0u#Q#180u. d av = 0 corresponds to random motion whereas values of 61 correspond to linear motion parallel or antiparallel to the field vector. d av = 1 (Q = 0u) corresponds to the case where cells move in the direction parallel to the field vector, whereas for d av = 2 1 (Q = 180u) cell motion is antiparallel to the field. For all experiments in the presence of an electric field, d av represents the average directedness over one hour during steady state motion after the initial transient response.
(4) The segment turn angle d i is defined as the angle between two successive segment vectors, z z i and z z iz1 (Fig. 2F), and ranges between 0u and 180u: d = 0u corresponds to no change in direction whereas d = 180u corresponds to a complete reversal in direction. By analyzing the distribution of segment turn angles, we can determine the correlation between each step under the influence of a dcEF and under confinement.

Galvanotaxis Alters Cell Morphology and Directs Cell Migration
Many cell types are known to respond to dcEFs with morphological changes and increased motility [20][21][22]. The exposure to dcEFs strongly dictates cell orientation and the direction of migration. In the absence of a dcEF, isolated 3T3 cells display typical fibroblast morphology with an elongated cell body and multiple protrusions ( Fig. 2A and Video S1) [23]. However, on application of an electric field, cells orient perpendicular to the field vector and are characterized by the formation of a broad lamellipodium extending the length of the cell and facing the cathode ( Fig. 2B and Video S2).
In the absence of an electric field, 53% of the cell paths (N = 43) were pointed in the direction of the field vector (towards the cathode) illustrating that there is no significant directional bias (Fig. 2C). However, in the presence of an electric field, the cell paths are strongly biased in the direction of the cathode (Fig. 2D) consistent with previous reports [22,24]. In a 5.5 V cm 21 field, 92% of the cell paths (N = 38) were directed towards the cathode. The average velocity of 3T3 cells in the absence of a field was 0.3860.02 mm min 21 (SE) (Fig. 2H). In a 2.2 V cm 21 field, there was no significant change in cell velocity (0.4060.03 mm min 21 (SE)), but cells preferentially moved toward the cathode and aligned in the direction perpendicular to the field. However, in a 5.5 V cm 21 field, the average velocity increased to 0.6860.03 mm min 21 (SE).
The persistency of the cells in the absence of a field was 0.6060.05 and remained approximately the same in a 2.2 V cm 21 field (0.6460.02 (SE), P = 0.27) (Fig. 2I). Further increasing the field to 5.5 V cm 21 significantly increased the persistency to 0.7160.05 (P,0.01). The average directedness in the absence of a field was 20.0260.06 (SE). The directedness increased to 0.2360.03 (SE) in a 2.2 V cm 21 field and 0.5660.12 (SE) in a 5.5 V cm 21 field (Fig. 2J).
Despite the change in morphology and the increase in directedness of the cell path (Fig. 2J), the average cell velocity under a 2.2 V cm 21 field (0.37 mm min 21 ) is the same as experiments with no field (Fig. 2H). A larger field of 5.5 V cm 21 results in an increase in the average cell velocity (Fig. 2H) and a further increase in the directedness (Fig. 2J). This result reveals a complex interplay between morphology, cell path, and cell speed, and suggests that cell velocity should not be used as the only parameter to probe the influence of electric field on cell motility.
The average cell orientation was strongly influenced by the field. In the absence of a field, the average orientation was 44.8u63.45u (SE), whereas at 2.2 V cm 21 and 5.5 V cm 21 , the average orientations were 71.5u63.1u and 80.3u60.94u, respectively.
We next analyzed the velocity distributions for the individual segments in a cell path. The cell velocities follow an exponential distribution according to f/f 0 = exp(2(v2v 0 )/a v )), where v 0 is the most probable velocity (defined by f 0 ) and a v is the characteristic velocity. Fits were determined for all velocities greater than or equal to the velocity at the peak of the distribution, v 0 (Figs. 3A, D and G). In the absence of an electric field, a v = 0.25 mm min 21 (R 2 = 0.94) (Fig. 3A). Although the persistent random walk model for motility commonly used to describe individual cell motility predicts a Gaussian distribution [19,25,26], exponential velocity distributions have been reported for various cell types [27,28], and it has been suggested that this behavior may be due to the limited production rate of ATP [27]. The distribution remains exponential with a similar exponent in a 2.2 V cm 21 field (R 2 = 0.99) (Fig. 3D), however, in a 5.5 V cm 21 field, a v increased to 0.62 mm min 21 (R 2 = 0.92) (Fig. 3G). This increase is consistent with the observation that the values of v 0 are around 0.20 mm min 21 in each case and that the average velocity remains the same between 0 and a 2.2 V cm 21 field but increases from 0.37 to 0.68 mm min 21 between 0 and a 5.5 V cm 21 field. The coefficient (a v ) is proportional to the average cell speed (Fig.  S2).
To characterize the orientation of individual segments, we analyzed the angular orientation of the segment vectors and the distributions of the parallel (Dx) and perpendicular (Dy) components of the displacements in each segment. In the absence of an electric field, the angular distribution of segment vectors Q i is relatively uniform with an average value of 42u (Fig. 3B). Thus, the corresponding scatter plot of the x-and y-components of each segment is symmetrical about the origin (Fig. S3A). The distributions of the parallel and perpendicular components of each segment overlap and are characterized by a gaussian distribution centered around zero (Fig. 3B inset). In a 2.2 V cm 21 field, a peak around 70u indicates that cells move more in the direction perpendicular to the electric field (Fig. 3E) than parallel to the field. The tendency is further shown in the corresponding scatter plot (Fig. S3B). In the presence of a 2.2 V cm 21 field, the distribution of Dy becomes broader and the frequency of Dx is significantly higher than Dy in the range of 62 mm (Fig. 3E inset), indicating that the average perpendicular component of the displacement is greater than the parallel component. In a 5.5 V cm 21 field, the distribution of Q i becomes strongly biased to small angles with respect to the field vector and can be described by an exponential function according to p/p 0 = exp(2Q/a Q ), where p 0 is the probability at Q = 0u and a Q is the characteristic segment angle. From a fit we obtain a coefficient a Q = 47.6u in a 5.5 V cm 21 field (R 2 = 0.89) (Fig. 3H), indicating that most segment vectors are within 48u of the field vector. A 5.5 V cm 21 electric field also significantly biased the distribution of Dx towards the cathode making it no longer symmetric around zero (Fig. 3H inset).
The distribution of segment turn angles d i can be described by an exponential function according to p/p 0 = exp(-d/a d ), where p 0 is the probability at d = 0 and a d is the characteristic turn angle. In the absence of a field (Fig. 3C), d is weakly exponential with a d = 37.3u (R 2 = 0.79). The exponent decreases progressively with increasing field with a d = 30.3u (R 2 = 0.77) in a 2.2 V cm 21 field and a d = 21.5u (R 2 = 0.79) in a 5.5 V cm 21 field (Figs. 3F and 3I). The segment turn angle and cell velocity distributions remain exponential in the presence of an electric field indicating that the mechanism of motility is the same despite the dramatic change in cell morphology ( Figs. 2A-B). Furthermore, there is no obvious influence of the dominant perpendicular contribution to cell motion seen at 2.2 V cm 21 (Fig. 3E). It has been suggested that the exponential distribution of segment turn angles reported for the migration of eukaryotic cells is due to the random nucleation of an actin-protein complex anchored to the leading edge of a protrusion [29].
As shown above, the angular orientation (Q) of individual segments is random in the absence of a field (Figs. 3B and S3A) but becomes biased in the presence of a field. In a 2.2 V cm 21 field, the distribution of Q becomes biased with a broad peak between 60u and 90u, indicating that the perpendicular component of the bias is larger than the parallel component (Figs. 3E and S3B). Although the cell paths are biased towards the cathode, as seen by the increase in average directedness from 0 to 0.23 (Fig. 2J), the cells are biased to move preferentially perpendicular to the field vector at small fields. Increasing the electric field to 5.5 V cm 21 further biased the trajectories toward the cathode, however, the parallel component of cell motion becomes dominant and the distribution of Q becomes exponential (Fig. 3H). The coefficient a Q = 47.6u implies that most segment vectors are within about 45u of the field vector. This can be seen more clearly in the scatter plot (Fig. S3C). These results suggest a complex response to the electric field. There is a global bias towards the cathode, as well as a mechanism that results in motion perpendicular to the field vector at low fields. At high field, the perpendicular component is dominated by a horizontal component parallel to the field vector, resulting in an exponential distribution of segment angles mostly within a characteristic cone of about 645u around the field vector.
In summary, in response to an electric field, 3T3 cells in 2D are biased towards the cathode. In a small electric field (2.2 V cm 21 ), the average cell velocity and persistency are the same as with no field, however the cells reorient perpendicular to the electric field and are biased towards the cathode. Furthermore, the vertical component of each segment (perpendicular to the field vector) is larger than the horizontal component (parallel to the field). In a larger electric field (5.5 V cm 21 ), the average velocity and directedness increase significantly. In addition, the horizontal component of the segments along a path become larger than the vertical component, and the segment vectors are generally within a cone of about 645u from the field vector. The distributions of cell velocities and segment turn angles are exponential in all cases. These results suggest that the mechanism of motility is due to a random process within the cell. However, the electric field modulates two processes, one that regulates cell motion perpendicular to the field and one that regulates cell motion parallel to the field. Globally, cell motion is not completely random in that a cell has limited ability to make large directional changes due to the orientation of protrusions associated with the spindle like morphology in the absence of a field and the broad cathode-facing lamellipodia in the presence of an electric field.

Galvanotaxis is Enhanced Under Confinement
To further study the effects of perpendicular alignment in 2D galvanotaxis on motility, cells are confined within microchannels to prohibit the orientation in the perpendicular direction. In the presence of an electric field and without confinement, cells oriented perpendicular to the field vector with a length of around 100 mm (Fig. 2B), therefore 20 mm wide channels aligned parallel to the field vector prevent the reorientation seen on an unconfined surface. In addition, spatial confinement allows us to impose quasi-1D motion on the cell paths. When confined to 20 mm channels, in the absence of an electric field, 53% of the cells (N = 40) moved in the direction of the field vector although the cell paths were strongly oriented due to the spatial confinement ( Fig. 4B and Video S3). In the absence of a field, the average velocity in the channels was 0.6860.03 mm min 21 (SE), significantly faster than cells with no confinement (P,0.01) (Fig. 4D).
Under confinement, cells were forced to align in the direction parallel to the channels. In most cases, cells migrated with the membrane in contact with the channel bottom and the channel walls. In the presence of a field, instead of a broad and persistent cathode-facing lamellipodia, as seen under no confinement (Fig. 2B), cells tended to contract their trailing edge and migrated with a smaller projected area and with the cell body becoming less flat ( Fig. 4A and Video S4). In the presence of a 5.5 V cm 21 field, 100% of the cells (N = 35) moved in the direction of the field vector. However, as we discuss later, this does not mean that all segments along a path were in the direction of the field vector.
In the presence of a field, cells migrated rapidly in the direction of the cathode (Fig. 4C). For all experiments in 20 mm channels, the average velocity was significantly larger than the corresponding case with no confinement (P,0.01) (Fig. 4D). In a 2.2 V cm 21 field, the average velocity was 0.7060.02 mm min 21 (SE), similar to the velocity in the absence of a field (P = 0.45). For a 5.5 V cm 21 field, the average velocity was 1.0260.07 mm min 21 (SE), significantly larger than for a 2.2 V cm 21 field (P,0.01). The increase in cell velocity in the channels can be due to either the change in cell morphology dictated by the spatial confinement or by the influence of channel walls on motility. An increase in cell speed has also been reported for 3T3 cells on quasi 2 1D lines modified with fibronectin formed by photopatterning and was thought to be due to the suppression of cell spreading and lateral lamellae [30].
The persistency in 20 mm channels was significantly higher than the corresponding case with no confinement (P,0.01). However, the persistency under confinement did not change significantly with electric field. In the absence of a field, the persistency in the channels was 0.7160.03, whereas in a 2.2 V cm 21 field the persistency was 0.7860.03 (P = 0.08). Further increasing the field to 5.5 V cm 21 had no obvious effect on persistency (P = 0.36) (Fig. 4E).
In the absence of an electric field, the average directedness in the channels was 0.0160.1 (SE), very similar to cell motion without confinement (Fig. 4F). The average directedness increased to 0.4860.06 (SE) under a 2.2 V cm 21 field, significantly higher than experiments without confinement (Fig. 4F). Under a 5.5 V cm 21 field, the directedness increased to 0.7160.04 (SE), very similar to experiments without confinement under the same field (P = 0.29), but significantly higher than for a 2.2 V cm 21 field (P,0.05). These results indicate that the influence of confinement is inversely related to the magnitude of the field. In a low field, confinement induces significant directional bias, whereas in a high field the directional bias is dominated by the field. The velocity distributions remain exponential for motility in 20 mm channels both in the absence and presence of an electric field, indicating that the mechanism governing cell speed is the same in all cases, despite the differences in cell morphology. It has been suggested that cell migration under confinement is less dependent on integrin-mediated adhesion but depends largely on microtubule dynamics [31,32] and actin polymerization at the cell membrane [33]. Therefore, although physical confinement can induce cytoskeletal alterations, the rate-limiting step that governs cell speed in unconfined and confined geometry may be the limited production rate of ATP, which gives rise to the exponential distribution of cell speed [27]. The distributions were fit using the same procedure as described previously. The coefficient v 0 was independent of the field, with values of 0.13 mm min 21 (no field), 0.13 mm min 21 (2.2 V cm 21 ) and 0.14 mm min 21 (5.5 V cm 21 ) (Figs. S4A, D, G). The exponents increased with increasing field, with a v = 0.5 mm min 21 (R 2 = 0.88; no field), a v = 0.66 mm min 21 (R 2 = 0.89; 2.2 V cm 21 ), and a v = 1.2 mm min 21 (R 2 = 0.92; 5.5 V cm 21 ) (Figs. S4A, D, G), however, the coefficients are comparable to the values obtained with no confinement (Fig. S4J).
The influence of confinement on cell trajectory was further characterized by examining the angular orientation of the segment vectors and the segment turn angle (Fig. S4). In the absence of a field, the segment orientation angle Q with respect to the field vector was bipolar and symmetrical ( Fig. S4B) since cells have equal probability to travel to the left or the right in the channels. The distribution around zero is exponential with a Q = 5.9u. This is in contrast to cell motion on a 2D surface where the distribution of angles is flat (Fig. 3B). In the presence of a 2.2 V cm 21 field, the distribution becomes asymmetric due to the bias on the cell paths as shown by the higher probability of Q near zero (Fig. S4E). The distribution around zero remains exponential with a Q = 5.1u. In a 5.5 V cm 21 field, the probability of Q near zero further increases and can be characterized by an exponential distribution with a Q = 5.2u (Fig. S4H). However, even though 100% of the cells moved toward the cathode in the presence of a 5.5 V cm 21 field, a small subset of segments along a path were pointed in the direction against the field vector as shown by the small peak near Q = 180u.
The segment turn angle, in the absence of a field, is bipolar but asymmetric with low probability around 180u (Fig. S4C). The distribution around zero is exponential with a d = 7.7u. In the presence of a 2.2 V cm 21 field, the distribution becomes more polarized (Fig. S4F) but the distribution around zero remains exponential with a d = 6.2u. Further increasing the field to 5.5 V cm 21 has no significant effect on the distribution of a d (Fig. S4I) around zero, which remains exponential with a d = 7.0u. The segment turn angle coefficients obtained in 20 mm channels are significantly smaller than the values obtained on 2D surfaces (Fig.  S4K).

Transient Response in Directedness and Orientation
In addition to the steady state behavior, we also studied the transient response to an electric field. The key parameters in the transient response are the directedness and the overall cell orientation. In the absence of confinement, on application of a 5.5 V cm 21 electric field, cell re-orientation occurs over about 2 hours (Fig. 5A). The time-lapse images show that cells first send out protrusions in the desired locations to align perpendicular to the field vector before developing stable cathode-facing lamellipodia. The time dependence of the average directedness of the cells follows an exponential increase (Fig. 5B) with a time constant t of 43.5 minutes (Fig. 5D). Similarly, the average cell orientation angle (h avg ) increases exponentially (Fig. 5C) with a time constant t of 40 minutes (Fig. 5D). The progressively smaller error bars indicate that cell orientation gradually converges to a steady state perpendicular to the field vector. Cell orientation also increased exponentially in the presence of a 2.2 V cm 21 field (Fig. S5); however, cells aligned substantially slower as suggested by a greater time constant (t = 95 minutes).

Conclusions
Through detailed analysis of cell paths we show that the influence of electric field on motility is much more complex than simply imposing a directional bias towards the cathode or anode. External electric fields with magnitude similar to endogenous electric field modulate several processes within 3T3 cells, involving overall cell reorientation, and motion both parallel and perpendicular to the field. Detailed analysis of cell paths reveal forces both parallel and perpendicular to the electric field that drives migration in a field dependent manner. Confinement in 2D channel prevents cell orientation perpendicular to the field and results in increased cell velocity. These results provide new insight into the biophysical response of cells to electric fields that can guide further research into the signaling pathways that regulate galvanotaxis.  Figure S3 Scatter plots of the x-and y-components of each segment in no field (A), a 2.2 V cm 21 field (B), and a 5.5 V cm 21 field (C). In the absence of a field, the scatter plot is symmetrical about the origin (A). In a 2.2 V cm 21 field, Dy has a greater effect on cell motion and therefore, the x-and ycomponents of each segment are scattered along the y-axis with a higher frequency between 60u to 90u and 260u to 290u (B). In a 5.5 V cm 21 field, horizontal bias toward the cathode becomes dominant and the majority of the x-and y-components are within 645u (C). (TIF) Figure S4 Distributions of average cell velocity, segment orientation (h), and segment turn angle (d) in 20 mm channels with no field (A, B, C), in a 2.2 V cm 21 field (D, E, F), and in a 5.5 V cm 21 field (G, H, I). The distribution of velocity is exponential in the absence of an electric field (R 2 = 0.88) (A), in a 2.2 V cm 21 (R 2 = 0.89) (D), and in a 5.5 V cm 21 field (R 2 = 0.92) (G). (B) In the absence of a field, the distribution of segment angles is bipolar and symmetrical with peaks at h = 0u and 180u. (C) In the absence of a field, the segment turn angle remains bipolar but with a large exponential distribution around 0u and a small distribution around 180u. (D) Distribution of velocity remains exponential in a 2.2 V cm 21 field. (E) In the presence of a field, the approximately symmetrical distribution of d becomes strongly biased towards small angles. (F) There is relatively little change in the distribution of segment turn angles in the presence of a field. (G) Further increasing the electric filed to 5.5 V cm 21 , significantly increases the strength of the exponential, a v increases from 0.5 to 1.