Introduction

The rheology of cells determines their response to internal and external mechanical stimuli encountered during normal physiological function. Since the cytoplasm forms the largest part of the cell by volume, its rheological properties are key to understanding cellular shape change in response to mechanical stress. The cytoplasm can be described as a biphasic material that consists of a porous solid phase bathed in an interstitial fluid¹. In such a material, stress relaxation arises from fluid flow through the pores of the solid phase as a result of spatial inhomogeneities in the applied stress field². By studying stress relaxation in micro-indentation experiments combined with perturbations including volumetric deformations together with chemical and genetic treatments, previous work has shown that at sub-second time-scales and micron length-scales, the cytoplasm behaves as a poroelastic material ^1,3. However, the contribution of cytoplasmic poroelastic properties to the mechanical response of cells has not been characterized on the tens of micron length-scale and second-to-minute time-scale relevant to cell physiology. Thus, it is unclear if poroelastic properties must be accounted for in models of cell mechanics.

During mechanotransduction, cells sense external mechanical stresses and translate them into biochemical signals⁴. The intracellular fluid flows and pressure gradients elicited in poroelastic materials in response to local application of stress may represent a stimulus for detection. While deformations are applied at the cellular scale, most mechanosensory processes act at the molecular scale through opening of ion channels and unfolding of proteins. Thus, the spatial extent of the stress and strain fields will determine what proportion of the cell is likely to respond to a stimulus and the temporal evolution of these stress and strain fields will determine the duration over which mechanosensory processes will be stimulated. Thus, understanding how stresses equilibrate in the cytoplasm in response to local application of force is necessary to better understand the physical parameters detected by intracellular mechanotransductory pathways. However, direct characterisation of the spatiotemporal deformation induced by sudden local application of force is currently lacking.

Migrating cells often adopt a polarised morphology with a gradient of non-muscle myosin II (NMII) protein increasing towards the rear^5–7or the front^8,9. Given the poroelastic nature of cytoplasm, theoretical considerations predict that such cortical myosin gradients should result in intracellular pressure gradients driving intracellular fluid flows in the direction opposite to the gradient of NMII^10,11. Intracellular flows have been proposed to participate in protrusion formation and migration but have only been indirectly observed in cells ^12–17. These flows may play a particularly important role during migration in confined environments with low adhesion, where most cell types adopt an amoeboid morphology and extend pressure-driven protrusions at their front^{5–7,18–20}. Although significant evidence points towards a role for myosin-generated pressure gradients in migration⁹, it is unclear if intracellular pressure gradients can be sustained over the minute-long time-scales involved in migration.

Here, we used a combination of cell culture experiments, high-resolution nanoparticle tracking, and finite element (FE) simulations to study how the combination of membrane permeability and cytoplasmic poroelastic properties can lead to steady-state gradients in intracellular pressure. We investigate the dynamics of cellular stress relaxation in response to external and internal mechanical stresses to determine if the poroelastic nature of cytoplasm is relevant to cell physiology on tens of microns length-scales and minute time-scales. We observe a cell-scale mechanical response following application of a local deformation of the cell surface and show that it is due to transient intracellular fluid flows elicited in poroelastic materials. We then reveal experimentally that intracellular pressure gradients lasting several minutes can be sustained in the cytoplasm, signifying that pressure gradients may play a role in migration.

Results and discussion

Whole-cell mechanical equilibration in response to localised deformation necessitates several seconds

Cells are often subjected to localised mechanical forces applied at high strain rates^21–23. In poroelastic materials, application of localised stress pressurises the interstitial fluid. Stress relaxes by flow of water out of the deformed region through the pores of the solid phase. As a consequence, the rate of mechanical equilibration is set by the poroelastic diffusion constant D_р that depends on the elasticity E of the solid phase, the viscosity η of the fluid phase, and the size ξ of the pores through which the interstitial fluid can permeate: D_р ∼Eξ²/ η ²⁴. The time-scale of water flow out of the deformed region is t_р ∼L²/D where L is a length-scale associated with deformation. Previous work has shown that local force application to a cell can result in global changes in the height of the cell surface and revealed the presence of two regimes: one fast due to stress propagation in the cytoskeleton and the other slow, whose origin was unclear²⁵. We hypothesized that the slow response taking place of over tens of microns length-scale and seconds time-scale reflects the poroelastic nature of the cytoplasm.

To examine the role of poroelasticity in the mechanical response of cells, we locally deformed the cell surface with the tip of an AFM cantilever. In our experiments, the deformation had a characteristic length L∼3 µm ( with R∼4 μm the radius of the indenter and δ∼2 μm the indentation depth) applied with a rise time t_r∼100ms (Fig 1A, B). In Hela cells, previous work has reported D_р ∼40µ m² /s ²⁶, leading to an estimate of the characteristic time of fluid efflux t_р ∼225 ms. As this is larger than the rise time, we expect poroelasticity to contribute to mechanical relaxation.

Figure 1:

To detect the global response of the cell surface to local force application with high spatial accuracy, we employed defocusing microscopy^3,25. In this technique, collagen-coated fluorescent beads are tethered to the cell surface. The plane of focus is chosen such that the beads are deliberately out of focus and display multiple diffraction rings about their centre (Fig S1A). Bead vertical displacement is monitored with nanometer precision by measuring the temporal evolution of the diameter of the outer diffraction ring in response to the deformation of the cell surface induced by AFM (Fig S1A, B). At steady-state, beads closer than ∼6µm to the AFM tip moved downwards while beads further away moved upwards, as expected for an elastic material subjected to indentation and as previously observed²⁵ (Fig 1C). The amplitude of the steady-state displacement decreased with distance between the bead and the AFM tip (Fig 1C). The temporal evolution of the bead movement displayed a biphasic response at all positions throughout the cell (Fig 1C). The first phase consisted in a fast movement of all the beads independent of their distance from the AFM tip over a time-scale shorter than ∼0.3 s, likely due to elastic effects (Fig 1C inset, Fig S2). In the second phase, the beads relaxed to their final displacement over a characteristic time scale τ_р that increased with increasing distance between AFM tip and the bead (Fig S2, Fig 1C inset, Fig. 1D, Methods). The overall relaxation of the second phase lasted up to ∼5s. Thus, local application of external force leads to a whole-cell mechanical equilibration lasting several seconds, consistent with previous work²⁵.

When stress is applied to a poroelastic material, relaxation takes place through flow of fluid out of the pressurised region at a rate that depends on the hydraulic permeability of the cytoplasm and hence its pore size. If the slow relaxation observed in cells is indeed due to poroelastic effects and intracellular fluid flows occurring at the cellular-scale, changes in the hydraulic pore size ξ should affect the duration of the second phase of bead movement but not the first phase. Pore size can be decreased by increasing the osmolarity of the mediumtreatment ^1,3. Increasing osmolarity by 300mOsm led to a ∼7-fold increase in the characteristic equilibration time τ_р of the slow phase (p=0.0008, Fig 1E) while the duration and amplitude of the first phase of bead movement remained unaffected (Fig S3). Thus, the second-long global changes in cell surface height that occur in response to local application of stress are qualitatively consistent with a poroelastic cytoplasm.

Cytoplasmic flows induced by microinjection follow Darcy’s law

In porous media, pressure imbalances are dissipated by interstitial fluid flows. Though intracellular fluid flows have been inferred to be driven by internal stresses during cell motility^{10–12,14,27} or external stresses during stress relaxation^25,26, a direct evaluation of the relationship between the pressure gradient ∇P_f and the velocity of intracellular fluid flow ν is lacking in cells. In classic porous media, the two are linked by Darcy’s law:

To characterise intracellular flows in response to pressure gradients, we created a fluidic link between a micropipette and an interphase cell using techniques developed for electrophysiology (whole-cell recording) and applied a step pressure increase to the micropipette (Fig 2A, Fig S4A, SI methods). In this configuration, the applied pressure results in fluid injection into the cell, causing the solid phase of the cytoplasm to expand and the height of the cell surface to increase as the fluid progressively permeates through the solid phase of the cytoplasm. In our experiments, we recorded bead displacement using defocusing microscopy for beads at different distances away from the micropipette (Fig 2B-2C). After application of a pressure step, beads moved upwards with an amplitude that decreased with increasing distance from the micropipette for a given time (Fig 2D) and with a time lag that increased linearly with the distance between the bead and the micropipette (see Methods and Fig 2E), indicating an approximately constant velocity of the propagation front. Over longer time-scales (>5s), continued fluid injection led to membrane delamination from the cortex forming several large blebs and, eventually, lysis.

Figure 2:

In our experiments, the duration of the microinjection experiments (∼2s) is short compared to the characteristic time-scale of water flows across the cell membrane (∼20-100s) ²⁸. Therefore, we can neglect fluid flows across the membrane and we expect flows to follow Darcy’s law in response to microinjection. As the velocity of Darcy flows depends on the hydraulic permeability k and the pressure gradient ∇P_f, we plotted the measured velocity v as a function of an estimate of the pressure gradient ΔP/Δx (see SI for estimation of injection pressure). This revealed a linear dependency with a slope k∼1.25 10^-13 m²/(Pa.s) (Fig 2F, r²=0.7028), close to previous estimates^1,29. We then repeated these experiments in metaphase cells in which we could precisely monitor diameter evolution after pressure application (Fig 3A). This revealed a strong dependency of the rate of change in diameter d and the time-lag δt (the onset of diameter increase) on applied pressure (Fig. 3B-D). When we plotted an estimate of the velocity ν ∼ ^d/δt as a function of an estimate of the pressure gradient ∇P_f ∼ΔP/d, this graph again revealed a linear relationship, as expected from Darcy’s law, and yielded an estimate of k∼5 10^-14 m²/(Pa.s).

Figure 3:

Together, these experiments show that pressure gradients lead to intracellular flows whose propagation follows Darcy’s law with a hydraulic permeability consistent with those derived by other experimental approaches ^1,3,29.

Global cellular response to microindendation is compatible with a poroelastic cytoplasm

Next, we sought to quantitively understand how the global cellular response to local indentation might arise from cellular poroelasticity. Previous work has hypothesised that local application of mechanical stress leads to a localised outflow of interstitial fluid from the deformed region that can then propagate through the cell ^1,25. If vertical displacements in indentation experiments were just due to fluid propagation, we would expect that the onset of displacement would occur with a time lag that depends on distance as in the microinjection experiments (Fig 2C-D). Yet, the first phase of displacement shows no time lag and vertical displacements of the cell surface are observed far from the point of indentation (Fig 1C inset). One potential mechanical origin for this may be the submembranous cortex whose mechanics differ markedly from the cytoplasm ³⁰ and which possesses a surface tension γ generated by myosin contractility³¹. This would result in a characteristic length scale over which deformations are transmitted instantaneously in response to local indentation. Using characteristic values of γ∼1 mN/m ³² E_cytoplasm ∼100 Pa ²⁶, and h∼5 μm, we obtain a length scale of ∼7 μm, consistent with the ∼10 μm observed in our experiments (Fig 1D).

The fast vertical displacement is followed by a second slow phase of displacement. With our experimental estimates of k, we can compute the time-scale of relaxation of a vertical displacement arising from a local deformation of a tensed membrane tethered to a poroelastic cytoplasm. For each position outside of the zone of indentation, the time-scale for relaxation is τ_z ∼L²/D_р. D_р ∼kE yielding a numerical estimate of D_р ∼ 10^-11 m²/s and L²∼lδ_z with l∼10 µm the distance over which vertical displacements are observed and a displacement amplitude δ_z∼0.5 µm. With these values, we find τ_z ∼500 ms, qualitatively consistent with the experimentally observed times (Fig 1D).

Thus, the presence of surface tension leads to instantaneous displacement of the cell surface far from the region of indentation (Fig 1C) and these vertical displacements may drive intracellular fluid flows throughout the cell.

Cytoplasmic poroelasticity combined with membrane permeability allows formation of stable intracellular pressure gradients

External stresses applied to the cell surface only give rise to transient intracellular pressure gradients before being dissipated by molecular turnover and intracellular flows. However, spatial inhomogeneities in internal stresses generated by the actomyosin cytoskeleton can be maintained for tens of minutes because they arise as a natural consequence of molecular turnover and contractility³³. In particular, during migration, the presence of gradients in myosin motors increasing from front to rear suggests the existence of a high cortical tension at the cell rear that can be sustained without relaxing^10,11. One implication, given the poroelastic properties of the cytoplasm, is that this surface gradient should result in a sustained pressure gradient continuously driving intracellular fluid flows towards the cell front. Cytoplasmic pressure gradients have been observed in fibroblasts migrating on 2D substrates¹⁴ and in 3D matrices^{9 34}. In addition, intracellular fluid flows have been inferred in cells migrating on 2D substrates^12,14 and suggested to play a motive role in cells migrating through confined environments¹⁷.

To gain a qualitative understanding of how long lasting pressure gradients can be maintained in the cell, we modelled the cell as a pressurised poroelastic material surrounded by a less permeable surface region 250nm in thickness, representing the membrane and the cortex. For this, we calibrated a poroelastic finite element model (Fig 4a-b), adjusting the values of E and D_р such that our model could predict the cell’s response to a localised deformation (Fig 4c) and the dynamics of displacement of the cell surface in response to microinjection (Fig 4d-e). Using this parametrisation, we introduced a sink in a region of the cell periphery (2um µ in diameter) to simulate a region of lower pressure (Fig 4f). In this region, fluid can rapidly leave the cell, which leads to a reduction in cell volume. In turn, this volume reduction will increase intracellular osmolarity and drive a water influx across the membrane in the rest of the cell surface. If the membrane permeability is large enough, efflux through the sink can be compensated by influx through the membrane to maintain a constant cell volume.

Figure 4:

Figure 5:

Following application of a suction pressure to the sink region, the cell reached a steady state presenting a spatial gradient of intracellular pressure with a low pressure close to the sink that increased towards the initial cell pressure far away from this region (Fig 4g-i). To gain insights into the importance of the membrane-cortex layer, we varied its thickness and diffusion constant. While its thickness had little influence ³, its diffusion constant strongly affected the length-scale of the gradient (Fig 4j), with high diffusion constants leading to gradients that reached the initial cell pressure over short length-scales. This is consistent with previous work that showed that localised exposure of cells to high osmolarity medium leads to localised dehydration with a sharp transition in pore size (and hence D_р) between the exposed and non-exposed region ²⁹. This work suggested that membrane permeability was on the order of 100 times lower than cytoplasmic permeability, within the lower range of the parameter values tested here (red curve, Fig 4j). Thus, cytoplasmic poroelasticity combined with passage of water across the membrane can in principle allow maintenance of stable intracellular pressure gradients and pressure compartmentalisation in living cells.

Cells can accommodate intracellular pressure gradients over minute time-scales

To test our predictions experimentally, we introduced a local depressurisation on the cell surface by establishing a fluidic link using a micropipette and bringing its back-pressure to atmospheric pressure. We verified that, in these conditions, a small suction was generated at the tip of the micropipette (Fig S5), signifying that the pressure applied at the pipette tip was lower than the intracellular pressure. Based on the cellular poroelastic properties, the pipette dimensions, and the cell dimensions, steady state is reached for after pressure release, assuming that relaxation is entirely limited by cellular poroelastic properties.

We verified that, following establishment of a fluidic connection between the cell and the micropipette, no occlusion occurred over time and that actomyosin was not perturbed (Fig S6-7). To confirm fluid efflux from the cell into the micropipette, we labelled cells with a fluorescent dye that becomes cell-impermeant upon cleavage by cellular proteases. We compared fluorescence temporal evolution in pairs of cells, one subjected to depressurisation and a neighbour that was unperturbed (Fig 6). Fluorescence intensity remained constant in the control cells but decreased approximately linearly over the course of 10 minutes in the depressurised cells, indicating the presence of a constant pressure gradient and efflux from the cell. Finally, we asked if the cell volume remained constant during depressurisation by examining the change in radius of prometaphase cells subjected to depressurisation. In these cells, cell radius varied by 10±1% over ten minutes of pressure release with no systematic trend to increase or decrease (N=5 cells, Fig S10). Overall, these experiments indicate that cells have a sufficiently large membrane permeability to allow rapid exchange of fluid across the membrane and that we can experimentally apply a long-lasting localised depressurisation.

Figure 6:

Our simulation results indicate that, for values of D_р measured in the cytoplasm, the low pressure region can be confined to a small zone near the pipette if membrane permeability is sufficiently large (Fig 6c). To experimentally test this, we used cell blebs as pressure gauges. Blebs are quasi-spherical protrusions of the cell membrane that arise due to pressurisation of the cytoplasm by cortical actomyosin contractility^24,35. First, we examined naturally occurring blebs in M2 melanoma cells³⁶. In control conditions, M2 cells bleb profusely with an intracellular pressure P∼400Pa but, when actomyosin contractility is inhibited with a Rho-kinase inhibitor (Y27632), cells no longer bleb ²⁴ and intracellular pressure drops to P∼100Pa (Fig 6A, see SI methods for pressure measurement). Therefore, if D_р is large and pressure is poorly compartmentalised, local depressurisation should lead to fast global decrease in intracellular pressure and blebbing should cease. In contrast, if D_р is low and pressure is strongly compartmentalised, blebbing should continue unperturbed³⁷. In our experiments, following local pressure release (Fig 6B), M2 cells blebbed for several minutes, far longer than necessary for intracellular pressure to reach steady state or for significant exchange of fluid to take place (for t>6min, 19/19 cells still blebbed, Fig 6D, Movie S2). Furthermore, no spatially localised inhibition of blebbing could be noticed close to the pipette (Fig S8). These results suggest that pressure is highly compartmentalised in melanoma blebbing cells and that pressure gradients can be maintained over several minutes. Our pressure measurements indicate that blebs cease to emerge if the intracellular pressure drops below ∼100Pa (Fig 6A). Based on this and assuming a diffusion constant of D_р ∼0.01µm²/s in the membrane-cortex, pressure-release would only induce a sufficiently large depressurisation to stop blebbing in a region ∼1.5µm away from the pipette tip (Fig 6C), consistent with the lack of spatial inhibition of blebbing we observe.

We then verified if HeLa cells could also compartmentalise pressure. In these cells, blebs can be induced within ∼2 minutes by partial depolymerisation of the F-actin cytoskeleton induced by treatment with low doses of latrunculin. In these conditions, blebs emerge because HeLa cells have an intracellular pressure ranging from ∼400Pa in interphase to ∼600Pa in metaphase (Fig 6A) and because, minutes after Latrunculin treatment, cortical actomyosin structures still remain well-defined (Fig S9). When contractility is blocked through inhibition of rho-kinase prior to treatment with latrunculin, blebs no longer emerge³⁸ (Fig S9, S11A). Thus, the growth of latrunculin-induced blebs depends on pressure generated by myosin contractility and they can be used as pressure gauges. When we locally depressurised HeLa cells by establishing a fluidic connection with a micropipette, we found that latrunculin treatment could still induce blebs in both interphase and metaphase cells (Fig 6E, S11B). This effect was independent of the duration over which depressurisation was maintained. When we established a pressure gradient in cells pretreated with Y27632, we did not observe blebs upon latrunculin treatment (Fig 6F, S11B).

As latrunculin triggers blebs indirectly and affects the whole of the actomyosin cytoskeleton, we repeated our experiments using blebs triggered by localised laser ablation of the cortex^35,39. In control conditions, a short pulse of UV laser focused on the cell cortex of a prometaphase cell led to the emergence of a bleb (Fig 6G, S11C). When myosin contractility was inhibited, ablation did not induce blebs (Fig 6G, S11C), consistent with ³⁵. When a pressure gradient was established, laser ablation could still induce blebs after several minutes (Fig 6G, S11C), indicating that pressure remained sufficiently large for bleb growth. Collectively, these experiments show that intracellular pressure is strongly compartmentalised in cells and that stable pressure gradients can be sustained over durations of several minutes.

Discussion

Our experiments revealed that cells could sustain pressure gradients across their cytoplasm over durations of ten minutes, relevant to cell polarisation and migration. In our experiments, we artificially generated an intracellular pressure gradient with a constant high pressure generated by the cell cortex and a low pressure at the tip of a micropipette. Using blebs as reporters of local intracellular pressure, we showed that, despite the presence of an intracellular pressure gradient, pressure over most of the cell periphery remained sufficiently high to continuously generate blebs over several minutes. This result was independent of cell type or cell cycle stage. Computational simulations indicate that the poroelastic properties of the cytoplasm combined with membrane permeability allow the maintenance of stable intracellular pressure gradients.

Global intracellular flows arise from the combination of cortical tension and a poroelastic cytoplasm

Here we show that local application of stress to the cell surface induces intracellular water flows spatially distributed over tens of microns and lasting seconds. When we locally deformed the cell surface with an AFM cantilever, we observed a global change in cell height with two different temporal regimes: one quasi-instantaneous and another that equilibrated more slowly, consistent with previous work²⁵. We attributed the initial fast relaxation to simultaneous motion of the solid and fluid phases and the second slowest relaxation to relative flows between the two phases equilibrating gradients of fluid pressure ²⁶. In line with this picture, when we decreased the cytoplasmic hydraulic pore size ξ, the second phase equilibrated slower taking several seconds. The existence of such spatially distributed intracellular flows may have important consequences for our understanding of mechanotransduction, which detects mechanical changes and converts them into biochemical signals. Mechanosensory processes act at the molecular scale through opening of ion channels or unfolding of proteins. Thus, the spatial extent of the stress and strain fields will determine what proportion of the cell has a stimulus sufficiently large and lasting sufficiently long to trigger these molecular processes. A challenge for future studies will be to link cellular-scale fluid flows and mechanical stresses to molecular-scale activation of mechanosensory processes to determine the strength of mechanical stimulus necessary to elicit a whole-cell mechanosensory response. Interestingly, intracellular flows and pressure gradients might participate in triggering signalling through mechanotransductory pathways located away from the cell membrane. For example, Phospholipase A2 activation at the nuclear membrane of epithelial cells is central to the response of tissues to wounding⁴⁰. Thus, the intracellular fluid flows revealed by our study might stimulate intracellular mechanosensory mechanisms, representing an indirect mechanotransductory mechanism.

Cortical contractility, cytoplasmic poroelastic properties and membrane permeability combine to enable sustained intracellular pressure gradients

Our findings have consequences for our understanding of cell migration in confined environments and low adhesive conditions when cells migrate using blebs^5–7. At steady state, migrating cells form a stable gradient in myosin that increases from front to rear and that powers migration. At the front, protrusion can either consist of a stable bleb when protrusions grow at a rate similar to actin accumulation^6,7 or a succession of blebs when actin accumulation rate is faster^5,18. Disruption of cortical tension at the cell rear by laser ablation leads to cessation of movement as well as localised blebbing in the region of ablation, indicating that myosin accumulation generates a high tension and high pressure at the rear⁵. Our work indicates the cell can sustain long-lasting gradients of pressure and raises the possibility that these may participate in migration. In this picture, the observed gradient in actomyosin distribution might generate an intracellular pressure gradient driving a forward-directed intracellular flow, consistent with some experimental observations^12,13. In support of this, recent work has shown that new cell protrusions emerge in regions of the leading edge where membrane-cortex attachment is the weakest, suggesting a role for pressure as a driving force ⁴¹. Therefore, pressure gradients may play a direct role in the generation of forward protrusion.

Conclusions

In summary, our work shows that steady-state gradients in intracellular pressure can arise through the combination of cytoplasmic poroelastic properties and membrane permeability. Thus, intracellular pressure gradients and intracellular fluid flows may play far more important roles than generally appreciated in cell physiology and poroelastic properties must be considered to gain a quantitative understanding of cellular phenomena such as mechanotransduction, cell shape change, and cell migration.

Materials and methods

Cell Culture

HeLa cells were grown at 37C with 5% CO₂ in DMEM (Life Technologies, UK) supplemented with 10% FBS (Sigma-Aldrich) and 1% Penicillin/Streptomycin. Melanoma (M2) cells were grown in MEM with Earle’s salts (Life Technologies, UK) supplemented with 10% of an 80:20 mix of newborn calf serum: FBS, and 1% Penicillin/Streptomycin. One day before the experiments cells were trypsinized (Trypsin-EDTA, Life Technologies), transferred from tissue culture flasks into glass bottom tissue culture dishes (Willco Wells, The Netherland). Prior to experiments, the medium was replaced with Leibovitz L-15 (Life Technologies, UK) supplemented with 10% FBS.

To visualise the cytoskeleton and nucleic acids, we established previously described stable cell lines: HeLa histone mRFP LifeAct GFP, HeLa LifeAct-Ruby, HeLa MRLC and M2 LifeAct-Ruby using retroviruses and lentiviruses. These were maintained with appropriate selection antibiotics (1mg/mL G418 and/or 250 ng/mL puromycin).

Cells were routinely tested for the presence of mycoplasma using the mycoALERT kit (Lonza). None of the cell lines in this study were found in the database of commonly misidentified cell lines maintained by ICLAC and NCBI Biosample.

Metaphase arrest

Cells were cultured to reach 70% confluency before being treated with 10 µM s-trityl-L-cysteine (Sigma-Aldrich) overnight to block them in prometaphase. Cells were washed three times using normal medium and then they were released into 20 µM MG132 (Sigma-Aldrich) for 1-2 hours. After this time, many cells were blocked in metaphase. Before the experiments, the medium was replaced with L-15 containing FBS and the same concentration of MG132.

Microscopy and laser ablation

Differential Interference Contrast and epifluorescence imaging was performed on a Nikon TE2000U (Nikon corp., Japan) inverted microscope. Images were captured on an EMCCD camera (Hamamatsu OrcaER, Hamamatsu, Japan) and transferred to a PC running µmanager (Micromanager, CA). Images were acquired using a 100x oil immersion objective lens (NA=1.3, Nikon) with 2x2 binning. The Ruby fluorophore was imaged using 561 nm excitation and collecting emission at 617 nm. GFP was imaged using 488nm excitation and collecting emission at 515nm.

For some experiments, we used an Olympus IX81 inverted microscope equipped with an Olympus FV-1000 scanning laser confocal head. All images were acquired with a 100x oil immersion objective. Imaging of Ruby and mRFP was performed using a 543 nm laser and imaging of GFP was done using a 488 nm laser.

Laser ablation experiments were performed as described in ³⁵ on a scanning confocal microscope (Olympus FV-1000) equipped with two scanning heads. For ablation, the cortex of metaphase HeLa cells was exposed to multiple pulses of a 405nm picosecond pulsed laser (Picoquant). Following induction, blebs grow rapidly before stopping and eventually retracting.

Chemical treatments

Drug treatments were carried out by adding the appropriate concentration to the medium. Drugs were present at all times during imaging and patch clamping experiments. Latrunculin B (250 nM Sigma-Aldrich) was used to induce blebbing in the cells by inducing partial loss of F-actin. Y27632 (50 µM, Sigma-Aldrich) was used to inhibit actomyosin contractility. Sucrose (200 mM, Sigma-Aldrich) was used to increase the medium osmolarity, which resulted in shrinkage of the cells and therefore decreasing the mesh size of cytoplasm. EIPA (50 µM, Sigma-Aldrich), an inhibitor of regulatory volume increase (RVI), was used in whole cell patch clamp experiments to prevent volume increase due to transportation of solutes into cell.

Electrophysiology

The experimental equipment setup consisted of a Digidata 1440A Digitizer and a MultiClamp 700B Amplifier piloted with the pCLAMP 10 Software (all from Molecular Devices, CA). Micropipettes were pulled from thin wall borosilicate capillaires (BF100-78-10, Sutter Instruments, CA) using a Flaming/Brown micropipette puller (Model P-97, Sutter Instruments, CA). Micropipettes had resistance of 6.0 - 6.5 MΩ and a tip diameter of around 2 μm.

For HeLa cells, the pipette was backfilled with a solution was composed of 150 mM K gluconate, 0.005 mM Ca gluconate, 1 mM Mg gluconate, 2 mM K-ATP, 1 mM EGTA, 5mM HEPES, and 5mM Glucose with pH=7.2. For M2 cells, the pipette was backfilled with a solution composed of of 130 mM KCl, 10 mM NaCl, 1 mM MgCl₂, 5mM Na-ATP, 5 mM EGTA, 10 mM HEPES, and 1 mM CaCl₂ with pH=7.2. The bath solution was L15 (Gibco Life Technologies, UK) for both cells types and did not contain FBS because this prevents gigaseal formation.

Data recording and synchronisation

In patch clamp experiments, the pressure transmitter, pinch valves, patch clamp equipment, and the microscope were all connected to the digitizer. This enabled us to control all devices through one platform. The electrophysiology software (pClamp10, Molecular Devices) and the imaging software (Micromanager) were set to communicate to each other. All of the steps to obtain whole cell configuration were performed manually and a macro was created in pClamp 10 to automatically acquire data once whole cell configuration was achieved. At this point by starting the macro in pClamp10, data acquisition, the timing of acquisition of each image, the timing of opening and closing of pinch valves as well as pressure measurement were recorded automatically through one software. This enabled us to synchronize all devices and determine the exact image at which pressure was applied to the micropipette and injection started.

In AFM experiments, for imaging, the camera (Hamamatsu OrcaER, Hamamatsu, Japan) was triggered every 100 ms using one of the output channels of the digitizer. The cantilever displacement and force were also acquired by connecting the corresponding output channels of the AFM to the digitizer using BNC cables. Therefore, similar to patch clamp experiments, we were able to synchronize all equipment by integrating them into one acquisition platform.

Measurement of intracellular pressure

Direct measurements of intracellular pressure were effected using the 900A micropressure system (World Precision Instruments) according to manufacturer’s instructions and as described in ⁹. Briefly, a 0.5-μm micropipette (World Precision Instruments) was filled with a 1-M KCl solution, placed in a microelectrode holder half-cell (World Precision Instruments), and connected to a pressure source regulated by the 900A system. A calibration chamber (World Precision Instruments) was filled with 0.1 M KCl and connected to the 900A system, and the resistance of each microelectrode was set to zero and then secured in a MPC-325 micromanipulator (Sutter Instrument) within an environmental chamber (37°C and 10% CO 2) on an Axiovert 200M microscope (ZEI SS). To measure intracellular pressure, the microelectrode was driven at a 45° angle into the cytoplasm, maintained in place for ≥5 s before being removed. The pressure measurement was calculated as the mean pressure reading during this interval of time.

Microinjection and pressure release setup

A pressure sensor (IMPRESS sensors and systems, IMP-LR-C0238-7A4-BAV-00-000) was used to measure the magnitude and temporal evolution of applied pressure. A glass Erlenmeyer was used as a pressure reservoir and was connected to a plastic tube which could be opened by a computer-controlled pinch valve. The tube was then connected to the pressure sensor and the micropipette holder. Two digital manometers were used to monitor the pressure in the reservoir and just before the pipette holder. The length of all tubing was approximately 1 meter and the time for pressure to propagate through the tubes, plus the response time of switches was around 30 ms (data not shown), which is less than our frame interval (66.7 or 100 ms). The pressure transmitter, pinch valves, patch clamp equipment, and the microscope were all connected to the digitizer. Before the start of the experiments, pressure in the reservoir was set. Once the whole-cell patch clamping configuration was achieved, a pulse of pressure was applied by opening the valve, resulting in injection of fluid into the cell. In pressure release experiments, after forming a whole-cell configuration, opening the valve resulted in connecting the pipette to the open atmosphere with pressure of zero (P_atm = 0).

Statistical analysis

All statistical analysis was carried using Microsoft Excel and MATLAB (Mathworks Inc, Cambridge, UK). Custom written codes in MATLAB were used to calculate desired parameters such as: area, speed, time delays, displacements, p value, etc. In all graphs, error bars indicate standard deviation. In box plots, the whiskers represent range of data.

Image processing

Fiji was used for producing kymographs and preparation of images for the figures.

Functionalisation of fluorescent beads

Yellow-Green carboxylate-modified fluorescent nanobeads with a diameter of 500nm (FluoSpheres, Molecular Probes, Invitrogen) were coated with collagen-I following the manufacturer’s protocol. To attach nanobeads to the cell membrane prior to experiments, HeLa cells were incubated for ∼30 min with a dilute solution (1:100 dilution) containing the fluorescent collagen coated nanobeads. Unattached beads were then washed out prior to experimentation.

Defocusing microscopy

Collagen coated fluorescent beads were added to cell culture dishes 30 minutes before the start of the experiments. Some of the beads attached to the cells. The cells were washed in L15 three times to remove unattached beads. Imaging solution was Leibovitz L-15 (Life Technologies, UK) supplemented with 10% FBS. Defocusing microscopy was implemented as described in ²⁵. A cell with 2 or 3 attached beads was selected. The motion of integrin-bound fluorescent beads was tracked in three-dimensions using defocusing fluorescence microscopy. By focusing a few microns above the bead plane, each bead appeared as a set of concentric rings. The distance between the beads and image plane is directly related to the radius of the outer ring and is used to determine relative z displacements. Time-lapse images were acquired every 67 or 100 ms and stored as a stack for each bead.

A code was written in MATLAB to automatically track the motion of the beads. Briefly, for a selected bead, a line intensity profile along the diameter of the concentric circles was obtained for each image in the stack. The radius of the outer ring was obtained from the Gaussian fit to the first and last peak in each image. For calibration, changes in radius were then converted into changes Z by moving a bead by a known distance using a piezoelectric stage ³ and acquiring images.

Atomic Force Microscopy and Data Analysis

Indentations of cells by AFM were performed using a JPK NanoWizard-1 AFM (JPK, Berlin, Germany) mounted on an inverted microscope (IX-81, Olympus, Berlin, Germany). The day prior to experimentation, cells were plated onto 35mm glass bottom Petri dishes. Experiments were performed at room temperature and cells were maintained in Leibovitz L15 medium (Life Technologies) supplemented with 10% FBS (Sigma-Aldrich) and MG132 (10µM). Before each experiment, the spring constant of the cantilever was calibrated using the thermal noise method implemented in the AFM software (JPK SPM). The sensitivity of the cantilever was measured from the slope of force-distance curves acquired on glass. For apparent stiffness measurements, we used soft cantilevers with V-shaped tips (BioLever OBL-10, Bruker; nominal spring constant of 0.006 N m^-1).

For each measurement, the cantilever was first aligned above the cell of interest using the optical microscope. Then, it was lowered towards the cell with an approach speed of 10 µm/s until reaching a force setpoint of 5nN and then kept the cantilever at a constant height.

Finite element modelling

We conducted finite element simulations to model the mechanical response at the cell scale, specifically focusing on the local deformations of the living cell surface and pressure gradients driving intracellular cytosolic flows. These simulations aimed to capture the deformation behaviour of a poroelastic cell when subjected to changes in effective pressure, mimicking scenarios such as fluid injection or pressure release in one region of the cell surface. Finite element models were developed using ABAQUS (version 2018). We used nonlinear geometry and unstructured mesh in our FE simulations, and also took into consideration a neo-Hookean isotropic porohyperelastic model ³. This was due to the large mechanical deformation range observed during our AFM, microinjection, and pressure release experiments. The best mesh and domain sizes were determined via mesh convergence studies, and a tolerance for the maximum pore pressure change per increment was calculated for SOILS analysis in our simulations.

AFM microindentation simulations

We ran simulations on two-dimensional rectangular sections to computationally investigate the influence of poroelasticity in the temporal mechanical response of cells to indentation (Fig 4a). To minimize edge effects, the cell was represented as a cylindrical disk (20μm radius, 20μm thickness) indented by 2μm with an infinitely rigid indenter representing the AFM cantilever tip. Frictionless and impermeable contact between the indenter and the cell was assumed, and a no-slip condition was imposed on the bottom surface of the cylindre. Pore pressure was set to zero except at the indenter contact surface to simulate fluid drainage. The simulation domain was discretized using quadratic quadrilateral CAX8P elements. Mesh sensitivity checks were performed to ensure independence of results on element size. The simulation consisted of a ramp step followed by a hold step. The ramp step entails indenting the material surface instantaneously (∼0.01s) with specific indentation displacement (∼2um). The phase in which the force reaction is released is known as the hold step. The main parameter of interest was displacement of thr cell surface in response to localized indentation. The force-relaxation curve and final displacement of the cell surface as a function of distance from the AFM tip were analyzed to extract the elastic modulus and poroelastic diffusion coefficient by fitting to the experimental data (For a thorough approach and curve fitting techniques, see ⁴²).

Fluid injection and pressure release

We ran FE simulations of a poroelastic cell that responds to the application of effective pressure changes at its top boundary where the micropipette contact its surface. The initial cell shape was idealised based on confocal profiles (Fig.1b). We modelled the cell as an axisymmetric elliptical cap with a diameter of 40 µm and a thickness of 4.5 µm attached to a substrate. The pipette was in contact with the top surface in the centre of the cell sufficiently long to enable equilibration of effective pressure with the same value as the endogenous cell pressure, P_eff = P_app - p_in. The difference between internal (P_in) and applied (P_app) pressure induces pressure gradients (P_eff) across the pipette-cell boundary and causes fluid to flow across this boundary. For fluid injection P_eff > 0→ P_app > p_in, and for fluid releasing P_eff < 0→ P_app < p_in were considered. Fig. 4b shows that the pressure increase within the pipette at t=0s led to fluid injection into the cell causing a rise in surface height and pressure changes until it reaches steady state, P_eff = P_app - P_in = 0. Our simulation had a ramp time (t=2s) over which applied pressure increased the internal pressure. The poroelastic material domain was discretised using the quadratic quadrilateral CPE8P element and the sensitivity of the FE simulations to domain size and mesh element numbers were checked. When the injection is applied, the cell surface reacts with a time lag (δ_t) that increases with the the distance from the micropipette. δ_t is defined as the time for which surface movement reaches 10% of the steady state displacement. In the FE simulations, the elastic modulus—which was extracted from AFM microindentation tests—was taken into account to fit the experimental fluid injection curves and derive P_eff and D. To show the impact of the membrane-cortex layer on intracellular pressure gradients and pressure compartmentalisation in living cells, we modelled cells as a double layer poroelastic material ²⁹. The cytoplasmic material is parameterized by its elasticity E₂, poroelastic diffusion constant D₂, and pressure P_in and it is surrounded by a less permeable thin layer representing the cortex parameterised by E₁, D₁, and P_in. A no-slip condition was imposed between the two layers.

Parameterisation of the Porohyperelastic model

As in previous studies, we used a fixed Poisson ratio of υ=0.3 ³ and extracted elastic modulus and diffusion constant by fitting experimental steady-state vertical displacement of beads tethered to the cell surface in response to localized indentation with FE simulation (Fig 4c). Considering the surface displacement curve, we extracted the following elastic modulus and diffusion constant, E=1.8kPa and D=28 μm²s^-1, consistent with our experimental results and previous work. To determine effective pressure (P_eff), elastic modulus (E), and poroelastic diffusion constant (D) for the fluid microinjection of the HeLa cells, we used optimisation processes to fit the vertical displacement of nanobeads positioned on the cell surface after 2s microinjection with FE simulations (Fig. 4d). The first step in the optimisation process was to conduct the simulations to match the experimental curves and get a first approximations of D and E taking into account P_eff =500Pa based on experimental measurements. Next, we adjusted E and D allowed to minimise error compared to the experimental curves and this yielded D= 13 µm².s^-1 and E=1.2 kPa.

Acknowledgements

The authors are grateful to members of the GC lab for feedback and suggestions on the manuscript. The authors wish to acknowledge Prof Guillaume Salbreux and Dr Pragya Srivastava for discussions. MM and GC were supported by grant WT092825 from the Wellcome Trust. GC was supported by a University Research Fellowship from the Royal Society.

Additional information

Author Contributions

MM and GC conceived the experiments. MM performed and analysed all experiments, except the intracellular pressure measurements that were performed by RP and LK. MHE and EM developed the model and performed analysis of the simulations. MM and MHE prepared the figures. GC, MM, MHE, and EM wrote the manuscript.

Additional files

Supplementary figures