Biophysical Journal: Biophysical Letters
Rheological Behavior of Living Cells Is Timescale-Dependent Dimitrije Stamenovic´,* Noah Rosenblatt,* Martı´n Montoya-Zavala,y Benjamin D. Matthews,y Shaohua Hu,z Be´la Suki,* Ning Wang,§ and Donald E. Ingbery *Department of Biomedical Engineering, Boston University, Boston, Massachusetts; yDepartments of Pathology and Surgery, Children’s Hospital and Harvard Medical School, Boston, Massachusetts; zDepartment of Environmental Health, Harvard School of Public Health, Boston, Massachusetts; and §Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois
ABSTRACT The dynamic mechanical behavior of living cells has been proposed to result from timescale-invariant processes governed by the soft glass rheology theory derived from soft matter physics. But this theory is based on experimental measurements over timescales that are shorter than those most relevant for cell growth and function. Here we report results measured over a wider range of timescales which demonstrate that rheological behaviors of living cells are not timescale-invariant. These findings demonstrate that although soft glass rheology appears to accurately predict certain cell mechanical behaviors, it is not a unified model of cell rheology under biologically relevant conditions and thus, alternative mechanisms need to be considered.
Received for publication 2 July 2007 and in final form 1 August 2007. Address reprint requests and inquiries to D. Stamenovic´, Tel.: 617-353-5902; E-mail:
[email protected].
INTRODUCTION The ability of mammalian cells to deform in response to mechanical forces is essential for behaviors that are critical for life, and thus many seek to identify general principles that govern cell rheology. It has been proposed recently that cell deformability results from timescale-invariant processes governed by soft glass rheology (SGR) theory derived from soft matter physics (1–5). The pivotal idea of SGR is that soft glasses are nonequilibrium systems that flow slowly when stressed such that their material moduli scale with the time and frequency of loading according to a single weak powerlaw (6). This behavior indicates that the rheology of these materials is governed by timescale-invariant processes, regardless of the physical basis for this behavior. SGR has recently gained great interest in biology because similar single power-law behavior has been observed in various cultured cell types with different techniques when analyzed over limited timescales using creep (1022–102 s) (2,3,7) and oscillatory (1022–103 Hz) (1,4,8) measurements. Based on these observations, SGR has been purported to be a fundamental property of living cells (5). The finding that a theory from soft matter physics can describe certain cell behaviors is extremely intriguing. But the timescales explored in those studies are not relevant for most physiological cellular processes (e.g., mitosis, motility, differentiation, apoptosis, and secretion) that take place over minutes to hours within cells that continuously remodel their cytoskeleton. Thus, the biological relevance of SGR remains unclear. To explore whether SGR is applicable at these physiological timescales, we measured the oscillatory response of the same cultured human airway smooth muscle (HASM) cells that were used in the original studies that described the SGR phenomenon (1–3) and with the same oscillatory twistBiophysical Journal: Biophysical Letters
ing magnetometry technique (1), except that we used a wider range of frequencies (1023–103 Hz). The magnitude of the cell dynamic modulus (jG*j) was determined using two protocols: one where jG*j was measured over a high frequency ( f ) range (1021–103 Hz), and the other where the effects of lower frequencies (1023–100 Hz) were studied. Detailed description of experiments is given in Supplementary Material. Interestingly, rather than observing a single power-law response, we found that the jG*j versus f relationships exhibited two distinct power-law regimes—one (100–103 Hz) less frequency-dependent, with exponent a1 ¼ 0.2 6 0.004, and the other (1023–1021 Hz) more so, with exponent a2 ¼ 0.33 6 0.034, separated by a third plateau regime (1021–100 Hz) (Fig. 1); the values of a1 and a2 were significantly different (p , 0.001). Importantly, analysis of the data from 100 to 103 Hz is fully consistent with previous reports claiming that HASM cells exhibit a single power-law (1); however, by including the data from 1023 to 100 Hz, it is obvious that these same cells exhibit multiple regimes that cannot described by a single power-law. Furthermore, when we modulated the cytoskeletal prestress—which is an important determinant of cell rheological behaviors (1–4,9)—by uniform (equibiaxial) stretching the substrate beneath adherent HASM cells (10), we still observed the two power-law regimes (Fig. 1), indicating that this behavior is robust. To rule out the possibility that these results are an artifact of using two different experimental protocols, and to independently confirm the generality of this multi-power-law Editor: Michael Edidin. Ó 2007 by the Biophysical Society doi: 10.1529/biophysj.107.116582 L39
Biophysical Journal: Biophysical Letters
FIGURE 1 The magnitude of the dynamic modulus (jG*j) versus frequency (f ) relationship measured in cultured HASM cells at the baseline (solid symbols) and after cell stretching by an 12% uniform substrate strain (open symbols) using oscillatory magnetic cytometry. Cells displayed two power-law regimes at 100– 103 Hz and at 103–101 Hz, separated by a plateau (101–100 Hz). Each regime was fitted by a function f a (solid lines). Data are means 6SE of median values of jG*j obtained from 180 to 500 magnetic beads distributed in 3–5 different wells for each condition (circles and triangles indicate data from two protocols).
mechanical behavior, we measured the creep response (0.01– 120 s) of a different cell type—bovine capillary endothelial cells—using a different technique (electromagnetic pulling cytometry) that applies sustained tension rather than oscillatory shear (11). Again, we confirmed the existence of two power-law regimes with an intervening plateau (see Supplementary Material), and this is consistent with previous reports (7,11,12). An a priori prediction of the SGR model is that if living cells are glassy systems, then they should deform slowly following a single power-law regime. Our oscillatory and creep data, taken together with past sporadic findings (7,11,12), clearly show that this is not the case. Thus, the underlying processes that govern the rheology of living cells are not timescale-invariant; instead, cells exhibit a transition between power-law regimes that occurs over a well-defined timescale. Moreover, the transition time appears to be cell type-specific: in HASM cells, it occurs between 0.1 and 1 Hz (i.e., 1–10 s) (Fig. 1), whereas it appears at ;10 s and 100 s in bovine capillary endothelial cells (11) (see Supplementary Material) and myocytes (7), respectively. Hence, SGR is not as universal a model as purported (5) when it comes to describing the rheological behavior over the full range of timescales that are biologically relevant for living cells. While SGR has effectively described several dynamic mechanical behaviors of cells (1–5), it has not provided insight into underlying molecular mechanisms, hence its value for cell biology remains limited. In contrast, analysis of cell rheology over a wider range of timescales, as described in this study, could lead to new and more biologically meaningful explanations. For example, theoretical studies of the viscoBiophysical Journal: Biophysical Letters
elasticity of myosin-activated entangled actin networks show that they exhibit multiple power-law behavior during stress relaxation (13) that is consistent with the behaviors observed here. This results from increased longitudinal fluctuations of individual actin filaments due to the mechanical action of molecular myosin motors. However, the longitudinal motion of actin filaments is constrained by multiple crosslinks in living cytoskeletal networks. Moreover, this theory also predicts that increased longitudinal motion of actin filaments will lead to fluidization of the network (13), whereas living cells (1–4) and cross-linked actin gels (14) stiffen when activated by myosin motors. Thus, it seems unlikely that the multiple power-law behavior we observed is due solely to the behavior of activated entangled actin networks in living cells. Alternatively, the emergence of the second power-law regime at low frequencies could be due to cytoskeletal remodeling. However, when we treated HASM cells with jasplakinolide (0.1 mM, 3 min), a drug that induces polymerization of actin filaments, we observed no changes in the low-frequency response (see Supplementary Material). Another possibility is that the high-frequency response is determined primarily by the rheological properties of individual actin filaments (15), whereas the low-frequency response results from these properties combined with the behavior of the higher order (cell-wide) cytoskeletal network; but currently there is no experimental evidence to support this claim. Finally, the high-frequency response might be based on the mechanics of the highly flexible submembranous ankyrinspectrin-actin cytoskeleton, whereas the low-frequency response could be due to the internal microfilament-intermediate filament-microtubule lattice. A novel observation is the existence of a well-defined timescale regime—the plateau (Fig. 1)—where no power-law behavior can be seen. The mechanism responsible for the plateau also remains unclear. This type of plateau behavior can occur when one power-law regime saturates for some reason, such as due to finite system size. Thus, for example, it could be due to a limited length-distribution of actin filaments in the cytoskeleton when analyzed at high frequencies. In contrast, a completely different mechanism may come into play at lower frequencies (e.g., limited length-distributions of intermediate filaments and/or microtubules), and this could generate another power-law regime with a distinct exponent. It is not clear what biological advantage this plateau region might offer from an evolutionary perspective, but considering that energy dissipation is reduced when jG*j is frequencyindependent, it may significantly increase the cell’s functional efficiency. In this context, it is noteworthy that breathing and heartbeat frequencies fall within the range of the plateau region in most species. We also considered the possibility that nonlinearities of the cytoskeleton may contribute to the emergence of the second power law since the nonlinear behavior of HASM cells became apparent below 10 Hz, indicated by the change in the L40
Biophysical Journal: Biophysical Letters
shape of bead displacement versus specific torque loops from elliptic (linear) to sigmoid (nonlinear) (see Supplementary Material). This emerging nonlinear behavior may reflect material nonlinearities of intermediate filaments (16) and actin stress fibers (17), that might become prominent with increasing bead displacements amplitude at low frequencies, or some other mechanisms, and together they could contribute to a new power-law regime. To explore this further, we calculated an index of nonlinearities, the harmonic distortion index (Kd), which characterizes the extent to which the input-output loop deviates from the linear elliptic shape. In a linear system, Kd ¼ 0, and in the presence of nonlinearities Kd . 0 (18). These studies revealed that Kd steadily increased with decreasing f, but this dependence did not become suddenly greater at low frequencies, indicating that the nonlinearities did not correlate with the transition to a new power-law that was observed in the jG*j versus f relationship (Fig. 1) (see Supplementary Material). In summary, regardless of the underlying mechanisms that govern rheological behaviors of cells, our data reveal that these mechanisms are not timescale-invariant, as previously reported (1,2,4). Instead, multiple power-law behavior, often with an intervening plateau, appears to be a property of various types of adherent living cells over the full range of timescales that include physiologically relevant frequencies. These findings demonstrate that SGR is not a unified model of cell rheology and suggest that different micromechanical mechanisms also need to be considered. Whether the SGR model can be integrated into a more general model of cell mechanical behavior remains to be seen. However, an effective model needs to incorporate quantitative descriptions of cell behavior amenable to computational and structural analysis in the context of physiological relevant molecular structures in living cells.
REFERENCES and FOOTNOTES 1. Fabry, B., G. N. Maksym, J. P. Butler, M. Glogauer, D. Navajas, and J. J. Fredberg. 2001. Scaling the microrheology of living cells. Phys. Rev. Lett. 87:148102. 2. Lenormand, G., E. Millet, B. Fabry, J. P. Butler, and J. J. Fredberg. 2004. Linearity and time-scale invariance of the creep function in living cells. J. R. Soc. Interface. 1:91–97. 3. Bursac, P., G. Lenormand, B. Fabry, M. Oliver, D. A. Weitz, V. Viasnoff, J. P. Butler, and J. J. Fredberg. 2005. Cytoskeletal remodeling and slow dynamics in living cells. Nat. Mater. 4:557–561. 4. Deng, L., X. Trepat, J. P. Butler, E. Millet, K. G. Morgan, D. A. Weitz, and J. J. Fredberg. 2006. Fast and slow dynamics of the cytoskeleton. Nat. Mater. 5:636–640. 5. Trepat, X., L. Deng, S. S. An, D. Navajas, D. J. Tschumperlin, W. T. Gerthoffer, J. P. Butler, and J. J. Fredberg. 2007. Universal physical response to stretch in living cells. Nature. 441:592–595. 6. Sollich, P. 1998. Rheological constitutive equation for a model of soft glassy materials. Phys. Rev. E. 58:738–759. 7. Desprat, N., A. Richert, J. Simeon, and A. Asnacios. 2005. Creep function of a single living cell. Biophys. J. 88:2224–2233. 8. Hoffman, B. D., G. Massiera, K. M. Van Citters, and J. C. Crocker. 2006. The consensus mechanics of cultured mammalian cells. Proc. Natl. Acad. Sci. USA. 103:10259–10264. 9. Stamenovic´, D., B. Suki, B. Fabry, N. Wang, and J. J. Fredberg. 2004. Rheology of airway smooth muscle cells is associated with cytoskeletal contractile stress. J. Appl. Physiol. 96:1600–1605. 10. Rosenblatt, N., S. Hu, J. Chen, N. Wang, and D. Stamenovic´. 2004. Distending stress of the cytoskeleton is a key determinant of cell rheological behavior. Biochem. Biophys. Res. Commun. 321:617–622. 11. Overby, D. R., B. D. Matthews, E. Alsberg, and D. E. Ingber. 2005. Novel dynamic rheological behavior of individual focal adhesions measured within single cells using electromagnetic pulling cytometry. Acta Biomater. 1:295–303. 12. Kole, T. P., Y. Tseng, I. Jiang, J. L. Katz, and D. Wirtz. 2005. Intracellular mechanics of migrating fibroblasts. Mol. Biol. Cell. 16:328–338. 13. Liverpool, T. B., A. C. Maggs, and A. Ajdari. 2001. Viscoelasticity of solutions of motile polymers. Phys. Rev. Lett. 86:4171–4174. 14. Mizuno, D., C. Tardin, C. F. Schmidt, and F. C. MacKintosh. 2007. Nonequilibrium mechanics of active cytoskeletal networks. Science. 315:370–373.
SUPPLEMENTARY MATERIAL
15. Rosenblatt, N., A. M. Alencar, A. Majumdar, B. Suki, and D. Stamenovic´. 2006. Dynamics of prestressed semiflexible polymer chains as a model of cell rheology. Phys. Rev. Lett. 97:168101.
To view all of the supplemental files associated with this article, visit www.biophysj.org.
16. Fudge, D. S., K. H. Gardner, V. T. Forsyth, C. Riekel, and J. M. Gosline. 2003. Mechanical properties of hydrated intermediate filaments: insights from hagfish slime threads. Biophys. J. 85:2015–2027.
ACKNOWLEDGMENTS
17. Deguchi, S., T. Ohashi, and M. Sato. 2006. Tensile properties of single stress fibers isolated from cultured vascular smooth muscle cells. J. Biomech. 39:2603–2610.
This work was supported by National Institutes of Health grant Nos. HL33009, GM072744, and CA55833.
18. Suki, B., Z. Hantos, B. Daroczy, G. Alkaysy, and S. Nagy. 1991. Nonlinearity and harmonic distortion of dog lungs measured by lowfrequency forced oscillations. J. Appl. Physiol. 71:69–75.
Biophysical Journal: Biophysical Letters
L41