Volume 8, Issue 1 , Pages 93-102, January 2012
The effects of cancer progression on the viscoelasticity of ovarian cell cytoskeleton structures
Article Outline
Abstract
Alterations in the biomechanical properties and cytoskeletal organization of cancer cells in addition to genetic changes have been correlated with their aggressive phenotype. In this study, we investigated changes in the viscoelasticity of mouse ovarian surface epithelial (MOSE) cells, a mouse model for progressive ovarian cancer. We demonstrate that the elasticity of late-stage MOSE cells (0.549 ± 0.281 kPa) were significantly less than that of their early-stage counterparts (1.097 ± 0.632 kPa). Apparent cell viscosity also decreased significantly from early (144.7 ± 102.4 Pa-s) to late stage (50.74 ± 29.72 Pa-s). This indicates that ovarian cells are stiffer and more viscous when they are benign. The increase in cell deformability directly correlates with the progression of a transformed phenotype from a nontumorigenic, benign cell to a tumorigenic, malignant one. The decrease in the level of actin in the cytoskeleton and its organization is directly associated with the changes in cell biomechanical property.
From the Clinical Editor
The authors have investigated changes in the viscoelasticity of mouse ovarian surface epithelial (MOSE) cells and demonstrated that ovarian cells are stiffer and more viscous when they are benign.
Graphical Abstract
The use of Atomic Force Microscopy (top left figure item) for the study of mouse (top right figure item) ovarian cells that have transformed from the benign/healthy state to a malignant/diseased one (middle figure item). Studies have resulted in graphical data along with calculation of key mechanical property parameters in conjunction with cell cytoskeleton component images to link change mechanical behavior with alterations in physiological morphology (bottom left, middle, and right figure items).
Key words: Atomic Force Microscopy, Viscoelasticity, Biomechanics, Cytoskeleton, Ovarian cancer
According to the American Cancer Society studies, cancer is the second most common cause of death in the US, exceeded only by heart disease; it accounts for nearly 1 of every 4 deaths and claims the lives of more than 1,500 Americans every day. The National Cancer Institute estimates overall costs of cancer may reach $158 billion by the year 2020. Deciphering the molecular events that lead to cancer initiation and progression for cancer diagnostic and treatment decisions based on individual genetic fingerprints has been a priority for many researchers in the last decade. However, the multigeneic cause of cancer and the genetic heterogeneity that determines disease severity and response to treatment has been a challenge and affected the success of the investigations.
Technological advances over the past few decades in the fields of nanotechnology, engineering and medicine have encouraged the investigation of the mechanical properties of individual cells to better understand human disease as well as the cell itself. It has been shown that transformed cells are often distinguishable from normal and healthy cells in many aspects, including cell growth, adhesion, morphology and organization of the cytoskeleton structure.1, 2 Therefore, the study of cell biomechanics (elastic and viscous behaviors), which have reportedly been linked to cell structure, is crucial for the development of disease-treating drugs and detection methods.3 As we further characterize the physics of cells, distinct properties could become potentially powerful tools in cancer diagnosis as well as treatment and therapy control.
Due to their “dynamic” nature, biological cells have an intricate microarchitecture that reacts and adapts to changes in both the chemomechanical environment and the disease state.3, 4 To date, various groups have made use of different techniques, including magnetic beads,5 optical tweezers,6 magnetic twisting and pulling cytometry,7, 8 optical stretching rheometry, micropipette aspiration and scanning force microscopy.3 However, the atomic force microscope (AFM) that was invented by Binnig, Quate and Gerber in 1986,9 has especially gained popularity in the task of cell structure characterization because of its unique nanoscale precision and ability to analyze live cultured cells.10 Moreover, because contrast in AFM originates entirely from the interaction force between the cantilever tip and sample, sample staining or fixation are not required.11 The AFM has thus far proved useful in identifying key differences between nontransformed and cancer cells by use of nanoindentation.
Several distinct cell lines in both human and murine species have been assessed in terms of their elastic responses to force stimulation; namely, human bladder,12 prostate,13 breast,14 mesothelial,15 cartilage16 and blood17 as well as murine ovarian18 and fibroblast tissues19 (Table 1). These cell lines were used to investigate the effects of cancer on individual cell structure. All groups have suggested that cancerous cells are “softer” or deform at a higher rate than their healthier, nontransformed counterparts by means of calculating their respective elastic moduli.
Table 1. Summary of previous works on the effects of cancer on cell elasticities via force microscopy
| Cell tissue type | Investigator | Year | Elastic Modulus(kPa) ± SD⁎ | Method | |
|---|---|---|---|---|---|
| Benign | Cancerous | ||||
| Ovarian | |||||
| Zhao et al | 2006 | 1.02 | 0.244 | AFM | |
| Bladder | |||||
| Lekka et al | 1999 | 7.5 ± 3.6 | SFM | ||
| Lekka et al | 1999 | 0.3 ± 0.2 | SFM | ||
| Lekka et al | 1999 | 0.8 ± 0.4 | SFM | ||
| Fibroblast | |||||
| Park et al | 2005 | 1.30 ± 1.0 | AFM | ||
| Park et al | 2005 | 0.6 ± 0.5 | AFM | ||
| Park et al | 2005 | 0.6 ± 0.5 | AFM | ||
| Prostate | |||||
| Faria et al | 2008 | 2.797 ± 0.491 | AFM | ||
| Faria et al | 2008 | 1.401 ± 0.162 | AFM | ||
| Faria et al | 2008 | 2.87 ± 0.052 | AFM | ||
| Breast | |||||
| Li et al | 2008 | 0.41 | AFM | ||
| Li et al | 2009 | 0.425 | AFM | ||
| Li et al | 2008 | 0.73 | AFM | ||
| Li et al | 2009 | 0.85 | AFM | ||
| Nikkhah et al | 2010 | 1.13 ± 0.84 | AFM | ||
| Nikkhah et al | 2010 | 0.51 ± 0.35 | AFM | ||
| Mesothelial | |||||
| S.E. Cross et al | 2008 | 1.97 ± 0.70 | 0.53 ± 0.10 | AFM | |
| Blood | |||||
| Lulevich et al | 2006 | 1000-35000 | AFM | ||
| Rosenbluth et al | 2006 | 0.045 ± 0.035 | AFM | ||
| Rosenbluth et al | 2006 | 0.855 ± 0.670 | AFM | ||
| Rosenbluth et al | 2006 | 0.156 ± 0.087 | AFM | ||
| Articular Cartilage | |||||
| Darling et al | 2007 | 0.343 | AFM | ||
| Darling et al | 2007 | 1.27 ± 0.86 | AFM | ||
| Darling et al | 2007 | 0.787 | AFM | ||
⁎Indentations performed on cell nucleus area. |
Among the various cancer types prevalent in women, ovarian cancer is said to be the fifth most lethal cancer type.20 It is a disease of older women and often diagnosed late, when the disease has already progressed and metastasized, at which point neither surgical nor chemotherapeutic treatment is very effective. Due to this nature, medical sources comment on this cancer type as the silent killer.
Ovarian cancer progression is accompanied by changes in gene expression levels that also affect cellular shape and architecture. Although research in this field is advancing, to date, there are no accounts of early detection for the human ovarian cancer. Moreover, no information exists about the mechanical properties of both malignant and benign human ovarian cells and the time course of changes, due to the lack of available and workable cell lines. Cell viscoelasticity has not sufficiently been determined quantitatively in any ovarian cancer progression model.
Recently, we have developed a primary mouse cell model for progressive ovarian cancer derived from C57BL/6 mice, relying on the spontaneous transformation in mouse ovarian surface epithelial (MOSE) cells in cell culture. During passaging, the cells undergo progressive changes in their growth rate and cell size, lose contact inhibition and acquire the capacity for anchorage-independent growth and tumor formation in vivo. According to their phenotype, they were categorized into early-benign-, intermediate- and late-aggressive stages of the disease.21 Thus, the syngeneic MOSE model represents a valid and novel alternative to human cell lines providing the early-, intermediate- and late-aggressive stages of ovarian cancer. Changes in the expression of cytoskeleton genes and their regulators and an increasingly disorganized cytoskeleton were observed during MOSE cell progression that may be associated with the aggressive phenotype.22 The focus of this fundamental study is to determine whether the differences in the mechanical properties of MOSE cells can be used as a potential marker for cancer detection. We have made use of the AFM method in connection with mathematical models to calculate both the cell elastic and viscoelastic parameters. While performing controlled examinations on MOSE cell lines representing early, intermediate and late stages of ovarian cancer, we have attempted to correlate the changes in elasticity and viscosity to disease state and the subsequent changes in the cell cytoskeleton structure.
Methods
Cell cultivation and sample preparation
The MOSE cell lines representing early (passage no. 15 – 25), intermediate (passage no. 75 – 80), and late stages (passage no.155 – 171) of ovarian cancer were generated as described in our previous studies.21 Before the experiments, cells were cultivated in plastic T-75 (75cm2 area) culture flasks. The environment used to nurture the cells consisted of high glucose Dulbecco's modified eagle's medium (DMEM-HG), which included 40 mL/L (4%) fetal bovine serum (FBS), 3.7g/L of sodium bicarbonate (HCO3), 10 mL/L of insulin-transferin-selenium (ITS), and 10 mL/L (1%) of penicillin-streptomycin solution. For AFM experimentation, cells were grown and harvested in incubators at 37°C in humidified 5% CO2 and later plated with the density of 1x105 cells/slip on 0.15 mm thick, 12 mm2 glass coverslips coated with 0.1 mg/mL collagen type IV (Sigma, St. Louis, Missouri) for 24 – 30 hours. The density of the cells was maintained consistently during all the experimentation to obtain an adequate number of individual cells for the movement range of the AFM cantilever. In addition, 40 μl/3mL of HEPES at 1M concentration was added to the samples to help maintain a physiological pH of 7.2 during the experimentation. The final pH concentration of HEPES in the culture medium was 13.5 mM and was stable for more than 2 hours, which was sufficient for all AFM tests.
Atomic force microscopy
Cell elasticityMeasurements were performed using the force-curve technique on a Multimode V SPM (Veeco Instruments, Santa Barbara, California) with integrated optical microscope. Olympus TR400PSA V-shaped SiNi cantilevers with approximate spring constant values of ∼ 0.02 N/m were employed in all AFM experimentations; exact spring constant values were measured via thermal tuning method. The probes were modified by attaching glass spheres (Duke Scientific, Waltham, Massachusetts) of ∼10μm diameter onto the cantilever free end with two-part epoxy (Miller Stephenson, Sylmar California), which helped the cells remain minimally damaged during contact. In addition, cell-deformation nonlinearity was reduced due to a more homogenous contact between the cells and the probe. The exact diameter of the glass sphere and its attachment location were identified using a HIROX KH-7700 3D Digital Video Microscope. All cell-elasticity measurements were performed in DMEM-HG culture medium by using a MTFML standard contact mode fluid cell at room temperature (∼24°C).
Probes were positioned at the cells' nuclei proximities under optical control, and force curves were acquired at a sample rate of 5 kHz and constant approach velocity of ∼0.5μm/s (Figure 1). At significantly higher speeds (i.e., 10μm/s), the speed-dependent hydrodynamic forces acting on the cantilever will ultimately affect the tip-sample interaction.23 Under the <1μm/s approach velocity setting condition, the movement is slow enough that viscous contributions are small, and force measurements are dominated by elastic behavior.24, 25 A maximum force trigger of 1.5 ± 0.3 nN was implemented for all force curves to maintain a basis for comparison among the cell lines (Figure 2, B).

Figure 1.
Optical microscope view screenshot taken during actual AFM experimentation. All indentations were performed on randomly selected cells that had visually “spread” morphologies.

Figure 2.
(A) Explanatory graph showing the full AFM stress relaxation procedure. (B) Force curves obtained via force curve technique on MOSE cells (C) Stress relaxation response curves on MOSE cells. It was visually possible to distinguish between the cell responses because of the precise sensitivity of the AFM.
Cell elastic modulus (Eelastic), was calculated by applying the Hertz model contact theory26, 27 to the AFM force curves. Sneddon's modification to the model was used to describe the behavior of elastic half-space bodies being subject to increasing loads by an infinitely hard surface.28 The model assumes that the surface is continuous, frictionless and incompressible at small deformations,18 which does not hold entirely true for biological cells. However, the Hertz contact model is sufficient to yield a general idea of cell elasticity trends and has been used by different research groups for characterizing cell stiffness.29 AFM indentation was limited to about 500 nm for the force curves to calculate the elastic moduli. This amount generally corresponded to about 15 – 30% of the total cell heights as observed in our scanning electron microscopy results on the MOSE cell line in our previous studies.30 Putting limitations on the indentation depth is an issue that is widely discussed among researchers in this field, stemming from the argument for the validity range of the Hertz model. In general, the findings of researchers have led to one opinion that the Hertz model is valid only for indentation ranges that were equivalent to about < 10% of the cell thickness, and on the other hand, reports have indicated that the same model was reasonably valid for all ranges of indentation from about 4% up to roughly 30%.25,31, 32, 33, 34 In light of these claims, we found it feasible to extend the indentations beyond 10 – 15% of the cell thickness to calculate the biomechanical properties of the MOSE cells.
Because the tip that comes in contact with the cell sample is spherical in shape, the relationship between the force and indentation is described by the following expression:
(1)
(2)
(3)because the Hertz model assumes linearity, this equation can be analyzed as a line in the (F2/3, z-d) plane, from which the Eelastic parameter may be directly extracted from the slope, and the exact contact point (z0, d0) from the intercept value of that expression. The semi-automated mathematical calculation of this point is vital for the accuracy of the Hertz fitting and has been proposed by Guo et al.35 All the force curve analysis which includes the linear curve fitting process as well as the Elastic modulus (Eelastic) calculations were made via MATLAB 2009b software.
Cell viscoelasticity
Viscosity is an important characteristic of a material because, in reality, all materials exhibit some form of time-dependent strain.36 Biological materials, in particular, are known to be able to both store and dissipate applied mechanical energy within their structures through internal frictional interactions, which are highly dependent on the rate of deformation.11 Therefore, to be complete in the structural analysis of biological cells, it is imperative that an analysis or design involving human tissue must incorporate their viscoelastic behavior.
Cell measurements for this parameter were performed by acquiring stress-relaxation responses using the force scripting technique on a Dimension 3100 AFM (Digital Instruments, Plainview, New York) with PicoForce functionality. A closed loop z-position feedback control on the AFM helped to maintain a constant cantilever displacement during experimentation (Figure 2, A). All cell viscosity measurements were performed within a DMEM-HG culture medium fluid meniscus between the AFM fluid probe holder and cell sample at room temperature (24°C).
Indentations were done at the cell nucleus proximity but with a constant approach velocity of 5μm/s, to approximate a step displacement that was appropriate for the stress- relaxation model.34 The data for each stress relaxation curve were sampled at 5 kHz for a total of 60 s upon the cantilever's indentation. A maximum force trigger of 3.0 ± 0.5 nN was implemented to designate the point at which the cantilever indentation halted and the corresponding z-scanner's position was held constant to record the change in deflection of the cantilever as a function of time on the sample (Figure 2, C).
μcell parameter extraction and theoryA viscoelastic behavior of a material can be determined by measuring its stress-relaxation response to a designated step displacement. Specifically, a cell's apparent viscosity (μcell) can be extracted by applying the Hertz contact model for viscoelastic testing of materials.16
Following the work proposed by Darling et al,16 the use of an appended transient response function will give yield to the following expression:
(4)where Er is the relaxed modulus, δ0 is the total deformation of the sample before stress relaxation, and both τσ and τɛ are the relaxation time constants for load and deformation, respectively. Rearranging the terms:
(5)This expression is an analytical equation of the force as a function of time. By fitting Eq. 5 to the force vs. time data, one can directly extract values for Er, τσ, and τɛ. These parameters describe a cell's viscoelastic response as a standard linear solid in the form of a spring-dashpot system in parallel with a spring. Among the several parameters that can be derived from the time constants and relaxed modulus, the apparent viscosity (μ) of the cell sample can be approximated by the following expression:
(6)During stress relaxation, we assumed an approximately constant strain on the cells while the force was recorded over time. This is because the change in cantilever deflection during relaxation was small enough to be considered negligible when compared to the overall indentation up to the maximum force trigger value. The stress relaxation curves analysis, which includes the viscoelastic curve-fitting process, as well as the parameter calculations, were done via MATLAB 2009b software.
Cytoskeletal confocal imaging
Immunofluorescence staining was performed on both the MOSE late and early stages to investigate their cytoskeletal organization. Our group wanted to specifically look into actin filaments in the MOSE cell line at different stages as it has been shown that actin filaments affect cell elasticity and not microtubules primarily.19, 37 For the same reason, the staining of vinculin and tubulin, which are intermediate fibers that also form part of the overall cell structure, was not emphasized in this study. Cells were grown on glass coverslips, fixed with 3% paraformaldehyde followed by permeabilization with 0.5% triton X-100 and quenching with 50 mM glycine. After incubation with the primary antibodies, the cells were incubated with AlexaFluor 488 (Molecular Probes, Eugene, Oregon) or TRITC- (Sigma) conjugated secondary antibodies and mounted with DAPI-containing ProlongGold (Molecular Probes, Eugene, Oregon). Images were taken on a Nikon TE2000 confocal microscope using the NIS Elements software.
Statistical analysis
Every indentable cell was selected at random as long as it had a spread/healthy morphology, irrespective of cancer phenotype, and was subject to a single indentation for either single force curves or stress-relaxation tests. No single cell was indented a second time, as it has been shown that cell morphology will change on contact.38 Shapiro-Wilks tests were conducted to analyze the data for normality distributions. Two sample-independent t-tests were used to see if there were significant differences between the parameter results of each cell line. A 95% confidence interval (P < 0.05) was applied to assess the degrees of difference between the results of all the cell lines.
Results
Several AFM experiments were conducted to best approximate the results with the least amount of variability possible. Cell elasticity measurements on the early, intermediate, and late-stage MOSE cells were performed on samples sizes of 126 – 137 cells (n = 126 – 137). Cell viscoelasticity measurements were performed on sample sizes n = 30 – 33 cells (n = 30 – 33). Both versions of the Hertz contact model theory that were fitted to the acquired respective data gave yield to high correlation coefficients (0.85 ≤ R2 ≤ 0.99) to suggest that theory matched satisfactorily well with experimental data.
The cell elastic property is best described by the elastic modulus, Eelastic. The parameter results were assessed in terms of cell population distributions (Figure 3).

Figure 3.
MOSE cell line histograms depicting the change in distribution of recorded elastic modulus from one stage of cancer to another.
The AFM indentations on the early-stage MOSE cells resulted in a broad distribution of values for Eelastic and were best described by a log-normal distribution scheme. It was interesting to observe that some cells gave yield to values that were about an order of magnitude larger than some cells in the same experiment, suggesting that this stage of the cell line is still heterogeneous with regard to elasticity because it was possible to detect a very deformable cell in proximity to a more rigid one. Intermediate stage MOSE cells produced a similar but slightly closer distribution, convincing that more cells in this stage correlated to lower elastic moduli. Finally, the late-stage MOSE cells showed a much sharper distribution, giving notice to the large number of cells that corresponded to relatively low (< 1 kPa) Eelastic values. Overall, a noticeable shifting pattern is evident in the distribution of cells from one stage of the cell line to the next. This indicates that an increasing number of low Eelastic cells are present in a cell population when transforming from a benign to an aggressive stage of the same cell line. The peak values that derive from the fitted Gaussian distribution can also be used to notice this change.
Cell viscosity (μcell) was used to represent a cell's viscoelastic response to force stimulation. Indentations and stress-relaxation tests were also analyzed via frequency distributions (Figure 4).

Figure 4.
MOSE cell line histograms depicting the change in distribution of recorded viscosity rates from one stage of cancer to another.
The cell-viscosity rate distribution was similar to that of the elastic modulus results. It is noticeable that higher rates of viscosity are present in the early, benign stage of the MOSE cell line as opposed to an increasingly homogeneous response of cells in the low range of viscosity in the later stages. Again, the peak values that arise from each stage's Gaussian distribution show this trend. The statistical analysis suggested a log-normal distribution scheme that best represented the viscosity trends above.
The elastic modulus and the viscosity rate parameters were also investigated in terms of their average and logarithmic value over each population of cell line (Figure 5). The schematic illustrates the trends from one disease state to another.

Figure 5.
Depiction of change in both elastic modulus and apparent viscosity parameters for the MOSE cell line during cancer progression.
A noticeable downward shift was detected for both the apparent viscosity and elastic modulus parameter values when cells transitioned from a benign state to a more aggressive, tumorigenic one. The average elastic modulus value showed about 50% decrease during the progression of MOSE cells from ∼1.1 kPa in the early MOSE cells to 0.55 kPa in the late MOSE cells. Likewise, the average apparent viscosity value showed a 65% decrease from 147 Pa-s to about 51 Pa-s. For both biomechanical properties, we also noticed a decrease in the standard deviation values, which correlate to the increasing degree of homogeneity of the cells' response in the aggressive stages. The intermediate stage MOSE cell line best served to validate the trends for both cell properties as it seemed to be in between the values for the other stages (0.8 kPa for Eelastic and 103 Pa-s for μcell).
The statistical significances were assessed after transforming the data into the logarithmic scale and we concluded that significant differences were present in terms of average apparent viscosity between the early versus late stage (P < 0.0001), and the intermediate versus late stage (P < 0.0001). Similarly, in terms of average elasticity, significant differences were observed between the early versus intermediate stage (P < 0.0001), the early versus late stage (P < 0.0001), and finally the intermediate versus late stage (P < 0.0001).
A series of additional tests were devoted to the validation of the output of parameter results. A total of 30 indentations (n = 30 cells) were performed for each session. The results showed that 15% – 20% error margins were associated with each cell line's average output, encouraging the hypothesis that there was a gap or difference between cell responses dependent on their stage. For example, there was a clear range between the lowest recorded early-stage cells output (0.886 kPa) and the highest recorded late-stage cells output (0.618 kPa). The intermediate-stage cell line output also situated comfortably between these values.
Noticeable differences were seen for all viscoelastic property parameters for the MOSE cell line (Table 2). The means as well as standard deviations showed a decrease in values, which reflects the observations through the obtained histograms.
Table 2. Elastic and viscoelastic biomechanical properties of the MOSE cell line (mean ± SD) in the progression of malignancy
| Cell property | MOSE early | MOSE-Interm. | MOSE-Late | |
|---|---|---|---|---|
| Elasticity | Eelastic (kPa) | 1.097 ± 0.682 | 0.796 ± 0.441 | 0.549 ± 0.281 |
| Indentation (nm) | ∼500 | ∼500 | ∼500 | |
| Viscoelasticity | Erelaxed (kPa) | 0.371 ± 0.229 | 0.310 ± 0.220 | 0.261 ± 0.090 |
| Einstantaneous(kPa) | 0.459 ± 0.278 | 0.429 ± 0.234 | 0.302 ± 0.096 | |
| Eyoung (kPa) | 0.554 ± 0.349 | 0.472 ± 0.306 | 0.395 ± 0.136 | |
| μ (Pa-s) | 144.7 ± 102.4 | 103.1 ± 53.76 | 50.74 ± 29.72 | |
| τσ (s) | 1.861 ± 0.554 | 1.620 ± 0.449 | 1.299 ± 0.391 | |
| τɛ (s) | 1.437 ± 0.369 | 1.311 ± 0.311 | 1.098 ± 0.280 | |
| Kelvin Spring elements | K1 (Pa) | 371 ± 229 | 310 ± 220 | 261 ± 90 |
| K2 (Pa) | 117.5 ± 85 | 73.36 ± 65 | 49.14 ± 31 | |
| K1+K2 (Pa) | 488.5 ± 313 | 383.4 ± 285 | 310.8 ± 121 |
A set of experiments were devoted to the cytoskeleton imaging of the early- and late-stage MOSE cells (Figure 6).

Figure 6.
Confocal images of MOSE early- and late-stage cells. A 20x objective was used in air to obtain these results. The images show the early-stage cells as having a much denser concentration of actin stress fibers when compared with their diseased counterparts.
The confocal images confirm there is a progressive dysregulation of the actin cytoskeleton, changing both the thickness and the organization of the actin fibers, which leads to its lesser concentration. These results correlate well with the changes in the general morphology of the early- and late-stage MOSE cells and changes in their cytoskeletal organization in the works of Creekmore et al, where it was found that a 78% reduction in F-actin concentration resulted from the transformation from early stage to late stage.22
Discussion
The acquired AFM experimental data shown above serve best to yield the following claims: Transformed and invasive MOSE cells are more deformable and less viscous in comparison with their benign counterparts. When generalized to large populations, the MOSE cell lines become increasingly more homogeneous in their response to the AFM indentations in the progression of disease. In effect, an increase in cell deformability directly correlates with the progression of a transformed phenotype from a nontumorigenic, benign cell to a tumorigenic, malignant one. Elastic and viscoelastic deformation characteristics are well defined and expressed using the Hertz contact theory model at indentations up to 30% of cell height. This suggests that sample substrate had no effect in the biomechanical property calculations.
Changes in the cytoskeletal organization have been associated with tumorigenesis.39, 40 We have shown previously that gene expression levels are changed during the progression of the MOSE cells with an enrichment of genes in the cytoskeleton functional category.22 Specifically, genes in the actin and actin-binding regulating category exhibited the most drastic changes. This was also reflected in altered protein levels, but, most important, was accompanied by a progressive disorganization of the cytoskeleton. Actin filaments gel the cell periphery, provide the highest resistance to deformation and are the most rigid of the three main components of the cytoskeleton components.3 Because it has been shown that the nucleus region is significantly softer and deformable than the periphery,41 this could directly connect our imaging results and the AFM mechanical measurements obtained as a result of indenting the cell nucleus regions.
Research on the nature of a cell's viscous response generated novel ideas; recent studies on the effect of cytoskeleton disrupting drugs (i.e., colchicine, paclitaxel, and nocodazole) have demonstrated that cells respond with a compensatory mechanism within their structure in the face of exposure and/or stimulation.42, 43, 44 Furthermore, it was revealed that cells that were subject to microtubule dissociating drugs had maintained their structural integrity and showed enhanced rigidity because of increased F-actin polymerization and subsequent increase in F-actin concentration within the cell structure in response to disrupted microtubulin. This phenomenon also resulted in a 30% increase in cell viscosity and cortical tension was increased by about 18%.43 It can be concluded from these findings that increased F-actin concentration is associated with increased cell viscosity. Therefore, the higher concentration of F-actin in the early-stage MOSE cell structure is most likely to be the cause of higher rates of apparent viscosity in comparison with later stages.
Janmey et al conducted a study of the major cytoskeletal fibers' stress resistance and found that actin filaments can sustain much higher stresses and withstand deformations at a much higher rate than other fibers in the cell structure. Furthermore, actin filaments tended to result in larger storage (elastic) moduli with increasing concentration.45 Actin filaments were also shown to have elastic moduli of about ∼2.3 GPa, as found with the method by Wagner et al, suggesting that these fibers were less elastic than rubber (E = 0.001 GPa), but more flexible than bone (E = 20 GPa).46 Because F-actin is thought to provide the highest resistance to deformation until a certain critical strain value3 and has a very rigid composition, these fibers are thought to cause differences in rigidity between the individual cell lines. In effect, a low elastic modulus and apparent viscosity is attributed to a lower concentration of the F-actin in a cancer cell's cytoskeletal structure. Our findings confirm these associations and suggest that the cytoskeleton affects the biomechanical properties of cells. Therefore, changes in these properties can be related to the motility of cancer cells and potentially their invasive potential because it allows cancer cells to deform, squeeze and migrate through size-limiting pores of tissue or vasculature onto other parts of the body.47, 48
It is noteworthy that our findings are a result of the applied model and approach as wwell as the loading rate settings during the AFM experimentations. For example, our studies showed that the biomechanical response changed noticeably as an increased approach velocity, and applied force led to a general increase in all elastic moduli in a proportional manner; these findings were also seen in other groups' works.49, 50 Furthermore, the Hertz contact model is only one of the proposed engineering approaches, and currently no proven model exists to precisely represent cell mechanics without assumptions.
Despite this source of variability, it is beneficial to identify changes in viscoelastic parameters due to cancer progression to aid the understanding of events that cause or support tumor growth and metastasis. To date, several cell lines have been analyzed in terms of their elastic modulus properties. However, their viscoelastic behavior is yet to be fully examined and/or understood. Darling et al tested the viscoelastic behavior in three well-characterized human chondrosarcoma cell lines that each reflected a different degree of invasiveness and malignancy (JJ012, FS090, and 105KC).16 Their finding that viscoelastic behavior changes with increasing malignancy is in general agreement with our own. However, these cell lines do not represent a progressive model like the MOSE model, and established cell lines may greatly adapt their geno- and phenotype to cell culture conditions.
Much more work is needed in the field of viscoelasticity of cells. The development of effective medical devices for the detection and diagnosis of cancer is highly dependent on the degree of realistic modeling of the cell samples. Although the work of previous groups emphasized solely the elastic modulus of cells, research into the individual cell characterization has much promise in the fight of cancer.
The results we obtained indicate that the progressive disorganization and pathophysiological alterations of the cytoskeleton architecture greatly affect their biomechanical properties. This confirms earlier reports and suggests that our results are not specific to our model but could indicate a general trend in the changes of the morphology and mechanics of cells due to cancer progression. The factors that are responsible for the differences in biomechanical properties and invasiveness among the ovarian cell lines are yet to be determined. However, the findings in this study, along with those of other groups, support the idea that highly invasive cancer cells are characterized by unlimited mitotic rates, higher migration rates19,and lower elastic and viscous properties that may be associated with significant differences in cytoskeleton structure. Although further work is needed to enhance our understanding of the complex architecture of living cells and the interaction of molecular and biomechanical events, we believe that our studies will contribute to the development of novel techniques and devices for clinical utility in cancer risk assessment, cancer diagnosis and treatment efficacy.
Acknowledgment
Award No. ECCS-IDR 0925945. The authors would like to thank the staff of the Nanoscale Characterization and Fabrication Laboratory (NCFL) at Virginia Tech for their technical assistance and providing the necessary AFM equipment.
References
- . Interactions of normal, dysplastic, and malignant mammary epithelial cells with fibronectin in vivo and in vitro. J Cancer Res. 1981;41:2115–2125
- . Pleomorphism of human prostatic cancer cells (DU 145) in culture—The role of cytoskeleton. Exp Mol Pathol. 1996;44:235–245
- . Biomechanics and biophysics of cancer cells. Acta Biomaterialia. 2007;3:413–438
- . Molecular biology of the cell. 4th ed.. New York: Garland Science; 2002;
- . A finite element model of cell deformation during magnetic bead twisting. J Appl Physiol. 2002;93:1429–1436
- . Mechanics of the human red blood cell deformed by optical tweezers. J Mech Phys Solids. 2003;51:2259–2280
- . Elastic behavior of cross-linked and bundled actin networks. Science. 2004;304:1301–1305
- . Novel dynamic rheological behavior of individual focal adhesions measured within single cells using electromagnetic pulling cytometry. J Acta Biomaterialia. 2005;1:295–303
- . Atomic force microscope. Phys Rev Letts. 1986;56:930–933
- . Atomic force microscopy in cancer cell research. In: Singh Nalwa H, Webster T editor. Cancer Nanotechnology. 2007;p. 1–17
- . Mechanics, malignancy, and metastatis: the force journey of a tumor cell. J Cancer Metastatis Rev. 2009;28:113–127
- . Elasticity of normal and cancerous human bladder cells studied by scanning force microscopy. Biophysics J. 1999;28:312–316
- Measurement of elastic properties of prostate cancer cells using AFM. Analyst. 2008;133:1498–1500
- . Nanomechanical analysis of cells from cancer patients. Nat Nanotechnol. 2007;2:780–783
- . AFM-based analysis of human metastatic cancer cells. J Nanotechnol. 2008;19:1–8
- . A thin-layer model for viscoelastic, stress-relaxation testing of cells using atomic force microscopy: do cell properties reflect metastatic potential. Biophysics J. 2007;92:1784–1791
- . Force microscopy of nonadherent cells: a comparison of leukemia cell deformability. Biophysics J. 2006;90:2994–3003
- . Elasticity measurement of living cells with an atomic force microscope: data acquisition and processing. Pflugers Arch-Eur J Physiol. 2008;457:551–559
- . The cytoskeletal organization of breast carcinoma and fibroblast cells inside three dimensional (3-D) isotropic silicon microstructures. Biomaterials. 2010;31:1–10
- . Cancer Statistics, 2010. CA Cancer J Clin. 2010;60:1–24
- Sequential molecular and cellular events during neoplastic progression: a mouse syngeneic ovarian cancer model. Neoplasia. 2005;7:944–956
- . Changes in gene expression and cellular organization in a model of progressive ovarian cancer. PLoSOne. 2011;6:1–16
- . Hydrodynamic effects in fast AFM single-molecule force measurements. Eur Biophys J. 2005;34:91–96
- Viscoelastic properties of individual glial cells and neurons in the CNS. Proc Natl Acad Sci USA. 2006;103:17759–17764
- . Endothelial, cardiac muscle and skeletal muscle exhibit different viscous and elastic properties as determined by atomic force microscopy. J Biomechanics. 2001;34:1545–1553
- . Theory of elasticity. 3rd ed.. New York: McGraw-Hill; 1970;
- . On the elastic contact of elastic solids. Journal für die reine und angewandte Mathematik. 1881;92:156–171
- . The relation between load and penetration in the axisymmetric Boussinesq problem for a punch of arbitrary profile. Int J Eng Sci. 1965;3:47–57
- . Atomic force microscopy probing of cell elasticity. J Micron. 2007;38:824–833
- . In: Distinct biomechanical characteristics of mouse ovarian cells as a marker for cancer detection, in 1st International Conference on Applied Bionics and Biomechanics(ICABB), Venice, Italy. 2010;
- . Probing mechanical properties of living cells by atomic force microscopy with blunted pyramidal cantilever tips. J Phys Rev E. 2005;72:1–10
- . Indentation and adhesive probing of a cell membrane with AFM: theoretical model and experiments. Biophysics J. 2005;89:3203–3213
- . Measuring the elastic properties of thin polymer films with the atomic force microscope. Langmuir. 1998;14:3320–3325
- . Viscoelastic properties of human mesenchymally-derived stem cells and primary osteoblasts, chondrocytes, and adipocytes. J Biomech. 2008;41:454–464
- . Packing density and structural heterogeneity of insulin fibrils by AFM nanoindentation. Biomacromolecules. 2006;7:1630–1636
- . Mechanical behavior of materials. Upper Saddle River, New Jersey: Prentice Hall; 1998;
- . Drug-induced changes of cytoskeletal structure and mechanics in fibroblasts - an atomic force microscopy study. Biophysics J. 2000;78:520–535
- . Cell mechanics: mechanical response, cell adhesion, and molecular deformation. Annu Rev Biomed Eng. 2000;2:189–226
- . The microfilament system and malignancy. Semin Cancer Biol. 2008;18:2–11
- . Cytoskeletal network in colon cancer: from genes to clinical application. Int J Biochem Cell Biol. 2004;36:759–765
- . Imaging mechanical properties of living cells by scanning probe microscopy. Curr Nanosci. 2007;3:97–103
- Distinct contributions of microtubule subtypes to cell membrane shape and stability. Nanomedicine. 2007;3:43–52
- . Passive mechanical behavior of human neutrophils: effects of colchicine and paclitaxel. Biophysics J. 1998;74:3282–3291
- . Bio-microrheology: a frontier in microrheology. Biophysics J. 2006;91
- . Viscoelastic properties of vimentin compared with other filamentous biopolymer networks. J Cell Biol. 1991;113:155–160
- . from molecules to motility. 2nd ed.. New York: Garland; 2000;
- . A critical step in metastatis: in vivo analysis of intravasation at the primary tumor. Cancer Res. 2000;60:2504–2511
- . Characterization of solid tumor microvasculature: a three-dimensional analysis using the polymer casting technique. Lab Invest. 1988;58:236–244
- . AFM indentation study of breast cancer cells. Biochem Biophys Res Commun. 2008;374:609–613
- . Rate- and depth-dependent nanomechanical behavior of individual licing chinese hamster ovary cells probed by atomic force microscopy. J Mater Res. 2006;21:1906–1912
Conflict of interest: No conflict of interest was reported by the authors of this paper.
E.M.S. and P.C.R. are supported in part by NIH/NCI CA118846. This study was made possible by the support of the National Science Foundation (NSF) Award No. ECCS-IDR 0925945 as well as the Institute for Critical Technology and Applied Science (ICTAS).
PII: S1549-9634(11)00184-5
doi:10.1016/j.nano.2011.05.012
© 2012 Elsevier Inc. All rights reserved.
Refers to erratum:
- Erratum to "The effects of cancer progression on the viscoelasticity of ovarian cell cytoskeleton structures" [Nanomed Nanotechnol Biol Med. 2012;8:93-102] , 09 March 2012
Volume 8, Issue 1 , Pages 93-102, January 2012

