Piezo1 as a force-through-membrane sensor in red blood cells

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: December 14, 2022 (This version)
Accepted: December 1, 2022
Received: August 11, 2022
Preprint posted: August 11, 2022 (Go to version)

1. Of interest
Substrate evaporation drives collective construction in termites

Giulio Facchini, Alann Rathery ... Andrea Perna

Research Article Apr 10, 2024
Further reading

Abstract
Editor's evaluation
Introduction
Results
Discussion
Methods
Data availability
References
Article and author information
Metrics

Abstract

Piezo1 is the stretch activated Ca²⁺ channel in red blood cells that mediates homeostatic volume control. Here, we study the organization of Piezo1 in red blood cells using a combination of super-resolution microscopy techniques and electron microscopy. Piezo1 adopts a non-uniform distribution on the red blood cell surface, with a bias toward the biconcave ‘dimple’. Trajectories of diffusing Piezo1 molecules, which exhibit confined Brownian diffusion on short timescales and hopping on long timescales, also reflect a bias toward the dimple. This bias can be explained by ‘curvature coupling’ between the intrinsic curvature of the Piezo dome and the curvature of the red blood cell membrane. Piezo1 does not form clusters with itself, nor does it colocalize with F-actin, Spectrin, or the Gardos channel. Thus, Piezo1 exhibits the properties of a force-through-membrane sensor of curvature and lateral tension in the red blood cell.

Editor's evaluation

This study reveals the super-resolution localization of native Piezo1 channels in red blood cells. The data provide further evidence to support the force-through membrane mechanism of mechanotransduction. Notably, the authors find that the distribution of proteins on a cell surface can be governed by membrane curvature.

https://doi.org/10.7554/eLife.82621.sa0

Introduction

Red blood cells (RBCs) experience significant mechanical forces in their lifetime as they traverse the mammalian circulatory system, squeezing through capillaries of smaller diameters than themselves (Chien, 1987). Central to the ability of RBCs to deform in response to shear stress is their biconcave disk shape (Lister, 1827; Mohandas and Evans, 1994). RBCs lack organelles and a transcytosolic network and their biconcave shape is instead imparted by the plasma membrane and the membrane skeleton, a 2D quasihexagonal network of actin filament (F-actin) nodes connected by (α₁β₁)₂-spectrin tetramers bound to transmembrane proteins (Fowler, 2013; Gratzer, 1981). In response to shear stress a Ca²⁺ influx is activated in RBCs (Larsen et al., 1981; Dyrda et al., 2010). Ca²⁺ entry into RBCs modulates junctional protein interactions in the membrane skeleton (Nunomura, 2006) in addition to activating the Gardos channel (Gardos, 1958), a K⁺ channel, leading to K⁺ efflux and subsequent cell volume decrease (Danielczok et al., 2017), which together facilitate flexibility of the cell membrane.

Piezo1 was first suggested as the mechanosensitive Ca²⁺ channel in RBCs after identification of gain-of-function mutations in patients with dehydrated xerocytosis (Bae et al., 2013; Albuisson et al., 2013; Andolfo et al., 2013; Zarychanski et al., 2012). Subsequent studies using a hematopoietic PIEZO1 knockout mouse confirmed that the channel is responsible for stretch-induced calcium influx in RBCs (Cahalan et al., 2015). How Piezo1 is organized in the RBC membrane, however, is not known. The unique biconcave shape of the RBC poses interesting questions about how membrane curvature, which locally is known to be important for Piezo1 structure (Guo and MacKinnon, 2017; Lin et al., 2019; Haselwandter and MacKinnon, 2018), might globally impact Piezo1 distribution and function. It has been proposed that due to its intrinsic curvature, Piezo1 might tend to concentrate in the dimple region of RBCs (Svetina et al., 2019), but direct experimental evidence is lacking. More broadly, how Piezo channels are organized in cell membranes is poorly understood. Based on electrophysiological recordings in overexpressing cells (Gottlieb et al., 2012) or fluorescence microscopy of a stable cell line (Ridone et al., 2020), some authors have reported clustering of Piezo1, but channel densities in both cases are likely non-physiological. Extensive patch-clamp studies and simulation of Piezo1 densities in neuro2Acells, which express Piezo1 natively, led others to the conclusion that Piezo1 channels are homogenously expressed across cell membranes and function as independent mechanotransducers (Lewis and Grandl, 2021). Additionally, it is unclear what role the actin cytoskeleton plays in organizing and gating the Piezo1 channel in cells. Piezo1 can be pressure activated in membrane blebs lacking cytoskeleton (Cox et al., 2016) and exhibits spontaneous openings in asymmetric lipid droplet bilayers (Syeda et al., 2016), consistent with a ‘force-from-lipids’ model of activation (Martinac et al., 1990). Others claim that Piezo1 is tethered to the actin cytoskeleton and is gated by a ‘force-from-filament’ mechanism (Wang et al., 2022).

Here, we set out to image and analyze the organization of Piezo1 channels in RBC membranes. The small size of RBCs and their inherent autofluorescence makes them refractory to conventional fluorescence microscopy, so we employed more sensitive super-resolution microscopic methods. Fluorescent Piezo1 spots can be resolved into single Piezo1 channels which, along with electron microscopy studies, confirm that they do not cluster in RBC membranes. Although not clustering, we find a non-uniform distribution of Piezo1 in the RBC membrane with enrichment at the dimple, which can be rationalized through the energetic coupling of the intrinsic curvature of Piezo1 to the curvature of the surrounding membrane. Single particle tracking (SPT) studies of Piezo1 demonstrate its lateral mobility in the plasma membrane. This is consistent with co-immunostaining analysis, suggesting that Piezo1 is not tethered to the membrane skeleton. Mobility of Piezo1 in the RBC membrane permits Piezo1 to explore the curvature landscape of the RBC surface and may allow Piezo1 to respond to local and global changes in membrane curvature as well as membrane tension.

Results

Expression and localization of Piezo1 in red blood cells

To detect Piezo1 in RBCs a Piezo1^HA mouse line was generated (Jackson Laboratories) in which an HA epitope tag was inserted into an extracellular loop of the PIEZO1 gene. The HA-tagged Piezo1 channel functions similar to the WT channel in electrophysiological recordings of HEK cells and could be specifically detected by immunostaining (Figure 1—figure supplement 1). Importantly, only RBCs from the Piezo1^HA mice and not wildtype mice showed fluorescent signal when immunostained with an anti-HA antibody (Figure 1A), confirming clear expression of Piezo1 in RBCs. As shown in 3D reconstructions, single optical sections and volume renderings from reconstructed 3D-structured illumination microscopy (3D-SIM) images (Gustafsson, 2000; Heintzmann and Huser, 2017) of immunostained RBCs, Piezo1 channels appear as spots (Figure 1A, B and Figure 1—figure supplement 2). These spots are distributed on the RBC membrane, here contoured by phalloidin staining of actin, which forms part of the actin–spectrin network bound to the cytosolic face of the RBC membrane (Ballas and Krasnow, 1980; Figure 1—figure supplement 2). 3D reconstructions from z-stacks reveal ~80 Piezo1 puncta per RBC (Figure 1C), suggesting a relatively low copy number of channels per RBC, similar to estimations of Gardos channel number per RBC (Lew et al., 1982; Brugnara et al., 1993) and in contrast to very abundant RBC membrane proteins such as Band3, which is estimated to have a copy number of 1 million per cell (Fairbanks et al., 1971). Supportively, our patch-clamp recordings of RBCs occasionally identified single channel currents that are stretch activated and have a conductance consistent with other reports of Piezo1 (Coste et al., 2010; Del Mármol et al., 2018; Harraz et al., 2022; Figure 1D).

Figure 1 with 2 supplements see all

Download asset Open asset

Detection of Piezo1 in RBCs.

(A) Reconstructed 3D-structured illumination microscopy (3D-SIM) images of red blood cells (RBCs), shown as maximum projections of z-stacks, from control (Ctrl, top row) or *Piezo1^HA* (bottom row) mice labeled using rhodamine phalloidin (magenta) and an antibody directed against the HA tag of Piezo1 (green). Scale bar, 3 μm. (B) 3D volume rendering of a representative labeled Piezo1-HA-KI mouse RBC displayed as XY (left) and XZ (right) axis views, showing the distribution of Piezo1. Scale bar, 2 μm. (C) Quantification of the number of Piezo1 spots per cell from 3D-SIM maximum projections. Mean ± standard deviation (SD) = 80.2 ± 26 spots per cell. Numbers were calculated on individual cells (n = 57) using Imaris spot detection. (D) Representative single channel traces from excised inside-out patch recordings of RBCs in symmetric 150 mM NaMethanesulfonate. The patch was clamped at −80 mV and negative pressure (0 to −30 mm Hg) was applied to the patch pipette. Amplitude histogram (right) enables estimation of single channel conductance. Showing one trace from n = 4 recordings.

Piezo1 does not cluster in RBC membranes

Whether Piezo1 channels cluster in membranes to form mechanosensory domains or are distributed more randomly as independent mechanosensors remains unknown. In 3D-SIM images of labeled intact RBCs, fluorescent Piezo1 spots are relatively uniform in size and do not appear to cluster (Figure 1). However, our SIM data are diffraction limited by an XY resolution of ~120 nm at best and thus each individual spot could feasibly contain a few Piezo1 channels (~23 nm diameter) tightly packed together. To address whether fluorescent spots represent single or multiple Piezo1 channels we therefore employed stimulated emission depletion (STED) microscopy (Hell and Wichmann, 1994), which can achieve resolutions below 50 nm and down to for example 22 nm after deconvolution (Schoonderwoert et al., 2013), depending on the selected dyes. Owing to the absorption of high-power lasers by heme (Schloetel et al., 2019), it was not possible to image intact RBCs by STED. Instead, RBCs adhered to poly-lysine coated coverslips were unroofed by a stream of isotonic buffer (Swihart et al., 2001), washing away hemoglobin and leaving patches of RBC membrane still bound to the actin–spectrin network (Figure 2A). The approximate flat dimensionality of unroofed RBCs also enabled us to use facilitated image acquisition 2D STED mode to achieve maximal XY resolutions of ~40 nm by full width half maximum (FWHM) in single images, closer to the size of a Piezo1 channel. Approximately 26 ± 7 (n = 12 cells) spots were observed per unroofed RBC, which equates to ~0.6 puncta/μm², consistent with the Piezo1 abundance observed in SIM images of intact RBCs. Piezo1 spots in 3D-SIM images of intact RBCs show a range of nearest neighbor distances (Figure 2B) with an average nearest spot distance of 540 ± 37 nm (n = 19 cells). This is similar to the observed distribution of Piezo1 spots in 2D STED images of unroofed RBC membranes, where the nearest spot distance is 544 ± 106 nm (n = 12 cells) and is consistent with fluorescent Piezo1 spots not clustering in the RBC membrane. Further cluster analysis of Piezo1 spots in 2D STED images by 2D spatial analysis (Andrey et al., 2010) shows no statistical difference between the observed spot distribution and a simulation of randomly distributed spots (Figure 2—figure supplement 1).

Figure 2 with 2 supplements see all

Download asset Open asset

Piezo1 does not cluster in red blood cells (RBCs).

(A) Representative 2D stimulated emission depletion (STED) image of an unroofed RBC immunostained with an anti-HA tag antibody for Piezo1 (green) and labeled with phalloidin STAR 580 for F-actin (magenta). Scale bar, 1 μm. (B) 3D reconstruction from SIM imaging of immunostained RBCs, where Piezo1 spots are detected by the spot detector function in Imaris and are colored yellow. The distance from each Piezo1 spot to the nearest Piezo1 spot in this cell was calculated in Imaris and is shown on the right. Scale bar, 1 μm. (C) The same analysis as in (B) was carried out on a 2D STED image of an unroofed RBC, but here Piezo1 spots at the edge of the membrane were masked out from spot detection to prevent edge-based artifacts in estimating nearest neighbor distances. Scale bar, 1 μm. (D) The mean nearest spot distance was calculated for SIM and STED analyses and plotted. For SIM the mean is 540 ± 37 nm (n = 19 cells) and for STED the mean is 544 ± 106 nm (n = 12 cells). (E) (Left) Huygens deconvolution of a Piezo1 puncta, imaged by 2D STED, resolving into a triplet after deconvolution, scale bar is 50 nm. (Bottom right) A top-down view of the Piezo1 channel structure with red asterisks indicating the positions of antibody binding. (Top right) An intensity profile of fluorescent signal from a line-scan of two adjacent bright pixels from the bottom left image, showing a distance of ~25 nm. (F) Negative stain electron microscopy of an unroofed RBC immunostained with an anti-HA tag primary antibody and a secondary antibody conjugated to 18-nm gold. Low magnification image of an unroofed RBC (left) with a region highlighted in red imaged at medium magnification (right). 18-nm gold particles corresponding to labeled Piezo1 channels are highlighted by red circles. As in fluorescence microscopy images, Piezo1 channels do not appear to cluster. (G) Negative stain electron microscopy image at high magnification of a triple-labeled Piezo1 channel.

Optimizing 2D STED image acquisition parameters, including imaging over multiple Z planes and employing adaptive illumination (Heine et al., 2017), as well as image processing through spherical aberration correction and iterative deconvolution (Huygens Professional) (Schoonderwoert et al., 2013) enabled an improvement of the lateral FWHM resolution toward ~25 nm (Figure 2E), approximating the diameter of the Piezo channel. Post-deconvolution, we observed that a significant number of fluorescent spots resolved into triplets of bright spots, consistent with labeling each arm of the trimeric structure of Piezo1 (Saotome et al., 2018; Guo and MacKinnon, 2017; Figure 2E and Figure 2—figure supplement 2). In a single field of view, we can often observe multiple triple-labeled Piezo1 channels (Figure 2—figure supplement 2). Line-scan analysis of the intensity between resolved bright pixels after deconvolution yields a pixel size-limited distance of ~25 nm (Figure 2E). This distance is consistent across analysis of all bright puncta that resolve into triplets, which we would not expect if these puncta represented multiple channels clustering. Further, Piezo1 puncta never resolve beyond triplets, consistent with single Piezo1 channels where all three antibody-binding sites are occupied, and not multiple trimeric Piezo1 channels clustered tightly together. The occurrence of triple-labeled Piezo1 in our STED data suggests we are toward the upper limit of labeling. An important corollary to this conclusion is that, at least if the probability of observing a given number of antibodies bound to each Piezo1 channel is non-cooperative, we should expect only a small fraction of Piezo1 channels to be unlabeled, giving us confidence in our estimation of Piezo1 channel number per RBC (Figure 1C). This settles inconsistencies in the available quantitative proteomics analyses of the RBC proteome (Bryk and Wiśniewski, 2017; Ravenhill et al., 2019) and is in line with a study of highly purified erythrocytes, which estimates approximately one hundred Piezo1 channels per cell (Gautier et al., 2018). We further analyzed Piezo1 distribution in unroofed red blood cells by negative stain electron microscopy using an 18-nm gold-conjugated secondary antibody (Figure 2F, G). Gold particles were distributed over the unroofed RBC membrane without apparent clustering, comparable to the distribution of fluorescent spots in STED images. Gold particles were not observed in RBCs from control mice. By negative stain electron microscopy, we were also able to observe triplets of gold particles (Figure 2G) with an inter-particle distance of ~25 nm, consistent with our 2D STED imaging of labeled trimeric Piezo1 channels. Thus, using multiple imaging approaches, we do not observe Piezo1 clustering in the RBC membrane.

Piezo1 is enriched at the dimple region of RBC membranes

We were curious whether Piezo1 is evenly distributed throughout the RBC membrane or tends to localize to specific regions of the RBC biconcave disk. Because of its intrinsic curvature, it has been hypothesized that Piezo1 concentrates in the curved dimple of the RBC membrane (Svetina et al., 2019). Other force-related proteins such as myosin IIA exhibit a non-uniform membrane density and are enriched at the dimple region (Alimohamadi et al., 2020). Using 2D slices from 3D-SIM images of intact RBCs we segmented the cells into the dimple and rim regions based on F-actin staining of the membrane (Figure 3A) and quantified the number of Piezo1 puncta per μm² surface area of each region using the cell imaging software Imaris (Figure 3B). The dimple accounts for about 17% of the total RBC surface area. We find that the whole RBC and rim have similar Piezo1 puncta densities (0.50 ± 0.16 and 0.46 ± 0.17 μm^–2, respectively) while at the dimple region the density is higher (0.80 ± 0.31 μm^–2). Thus, there is approximately a twofold enrichment of Piezo1 at the dimple region of RBC membranes compared to the rim (Figure 3B). This enrichment of Piezo1 at the dimple region of the RBC membrane is not due to enzymatic deglycosylation of RBCs, which we used to increase antibody binding to HA-tagged Piezo1 (Figure 3—figure supplement 1). To strengthen the hypothesis that it is the intrinsic curvature of Piezo1 that biases it toward the curved dimple region of the RBC we also analyzed the dimple versus rim ratio of two other RBC membrane proteins that are not intrinsically curved: the Gardos channel (Lee and MacKinnon, 2018) and Band3 (Vallese et al., 2022; Xia et al., 2022; Arakawa et al., 2015). We find no statistically significant enrichment of the Gardos channel (KCNN4) or Band3 in the dimple of RBCs (Figure 3C). Further, we find no statistical enrichment of actin at the dimple, indicating that the actin–spectrin meshwork is not more densely packed there (Figure 3C).

Figure 3 with 2 supplements see all

Download asset Open asset

Piezo1 localization is enriched in the curved dimple regions of the red blood cell .

(A) Two dimensional slices in XY and YZ views from SIM images illustrating volume segmentation of RBCs and Piezo1 puncta distribution after labeling with an antibody against Piezo1-HA (green) and rhodamine–phalloidin for F-actin (magenta). Cells were segmented and masked in Imaris as shown. Scale bar, 1 μm. (B) Distribution of Piezo1 puncta in RBCs. The RBC dimple has a higher density of Piezo1 puncta than whole RBCs (Total) (p = <0.0001) and the rim region (p = <0.0001) by the one-way analysis of variance (ANOVA) test with Tukeys post hoc analysis, as indicated on the graph by asterisks where *** indicates p<0.001 and **** indicates p<0.0001. N = 39 cells (C) The ratio of Piezo1 spot density (1.87 ± 0.76, N = 39 cells) in the dimple and rim region of RBCs compared to that of KCNN4 (0.99 ± 0.36, N = 12 cells), Band3 (0.95 ± 0.33, N = 20 cells), and Actin (0.93 ± 0.36, N = 31 cells). The dimple over rim ratio is higher for Piezo1 than KCNN4 (p = 0.0003) and Band3 (p = 0.0003) and Actin (p < 0.0001). (D) Scatterplot of ratio of Piezo1 puncta density in dimple over rim plotted against RBC biconcavity, measured as height ratio a/b (see inset). N = 37 cells. Data are fit by linear regression with equation y = 1.08x + 0.55, R² = 0.34 and p = 0.0001. (E) Probability per unit area for finding a Piezo1 channel at a given RBC membrane location, calculated from a physical model of Piezo1 curvature coupling (see methods) and plotted over the RBC membrane surface. The RBC shape corresponds to Beck’s model of for RBC discocytes (Beck, 1978).

Because we observe a range of relative enrichments of Piezo1 at the dimple we asked whether this correlates with the degree of RBC biconcavity, since we also observe a range of biconcavities via F-actin staining of the RBC membrane (Figure 3—figure supplement 1). The degree of biconcavity was parameterized as the ratio of the maximum to minimum heights of RBCs at a central XZ slice (Figure 3D). Indeed, with increasing RBC biconcavity, we observe a greater enrichment of Piezo1 at the dimple (Figure 3D). If the distribution of Piezo1 in the RBC membrane is determined by the coupling of the channels’ intrinsic curvature to the global curvature of the membrane we would expect an unfavorable distribution of Piezo1 into membrane protrusions where the membrane curvature is opposing to that of Piezo1’s. Treatment of RBCs with NaSalicyate, which preferentially binds into the outer membrane leaflet, generates echinocytes (Li et al., 1999), RBCs that are characterized by ‘thorny’ membrane protrusions. 3D-SIM imaging of labeled echinocytes qualitatively shows a depletion of Piezo1 channels from such membrane protrusions relative to regions of negative membrane curvature, further supporting a mechanism of curvature coupling between the intrinsic curvature of Piezo1 and the surrounding membrane (Figure 3—figure supplement 2).

The relative abundance of Piezo1 in the dimple and positive correlation with the degree of biconcavity indicates that RBC membrane curvature might be the primary determinant of the Piezo1 distribution. The intrinsic curvature of Piezo1 suggests a possible mechanism, based on ‘curvature coupling’, for bias toward the dimple. To examine the plausibility of this mechanism, we developed a simple physical model of Piezo1 curvature coupling in which Piezo1 is represented as a spherical cap with area 450 nm² and approximate radius of curvature 42 nm, the estimated asymptotic value in a planar membrane (Haselwandter et al., 2022a). To describe the RBC shape, we used the model of Beck, 1978 and calculated the mean curvature at each point on the surface (Methods). Then, using the Helfrich functional (Helfrich, 1973), we calculated the minimum shape energy of the membrane surrounding Piezo1 in a vesicle that, if spherical, would have a radius of curvature equal to the reciprocal of the mean curvature at that point on the RBC surface. This procedure associates each point on the RBC surface with a membrane curvature energy which, when applied to the Boltzmann distribution equation, yields a probability distribution for the density of Piezo1 on the membrane surface. Further details of this calculation are outlined in Methods. The surface distribution of Piezo1 calculated according to this method is shown in Figure 3E. The curvature coupling model predicts an approximately two-fold increase in the density of Piezo1 in the dimple relative to the rim, close to the experimentally observed ratio of densities. In the Beck model of RBC shape, the mean curvature at any point on the RBC surface ranges from 0.24 μm⁻¹ (favorable to Piezo1 in the dimple) to about −0.52 μm⁻¹ (unfavorable to Piezo1 in the rim). Because these mean curvatures for the RBC membrane are very small (compare these values to the approximate intrinsic curvature of Piezo1, 1/42 nm ≈24 μm⁻¹), it was not immediately clear to us whether Piezo1 would be sensitive to surface curvature in the RBC. The correspondence between experiment and theory suggests that Piezo1 is indeed sensitive to RBC shape and adopts an equilibrium distribution by a mechanism of curvature coupling. If one thinks about this result, the sensitivity of Piezo1 distribution to the global curvature of the RBC membrane implies that the membrane must be relatively smooth. Indeed, this implication is consistent with electron microscopic images of healthy RBCs (Marikovsky and Danon, 1969; Lewis et al., 1968).

Piezo1 is mobile within the red blood cell plasma membrane

The above distribution of Piezo1 according to the curvature coupling model suggests that Piezo1 ought to be laterally mobile so that it can explore the curvature landscape of the RBC membrane. However, the quasihexagonal actin–spectrin meshwork and large junctional protein complexes that make up the RBC membrane layer are known to corral and confine membrane protein diffusion (Sheetz, 1983). Additionally, at least some membrane proteins, including roughly one third of Band3, are essentially immobilized by direct binding to actin–spectrin (Tomishige et al., 1998; Tsuji et al., 1988). We studied Piezo1 mobility in the RBC membrane with SPT using 40-nm gold nanoparticles conjugated to a secondary antibody under differential interference contrast (DIC) microscopy (Wu et al., 2019). The fast time resolution and ability to record over minutes without photobleaching has advantages for studying confined membrane protein diffusion (Tomishige et al., 1998; Daumas et al., 2003). Gold particles could be unambiguously identified on RBC membranes (Figure 4A and Video 1). We imaged at 100 Hz for 2 min. In recordings lasting 2 min we observed a range of Piezo1 diffusion trajectories (Figure 4B–D). Most trajectories exhibited confined diffusion, although sometimes Piezo1 appeared to be freely diffusing. Occasionally, we observed both types of diffusion in a single recording. In all cases, Piezo1 was more mobile than in control cells chemically fixed with paraformaldehyde, when analyzed by a mean squared displacement (MSD) analysis (Figure 4D). For confined trajectories, we find that the confinement diameter is in the range of 200–600 nm. Although limited by the spatial resolution of our recordings, this range is largely consistent with the organization of the actin–spectrin meshwork we observe by 2D STED imaging of unroofed RBCs stained for actin or spectrin (Figure 4—figure supplement 1). In addition to smaller corrals of intact actin–spectrin, we find voids deficient in actin or spectrin staining, approximately 200–500 nm in size (Figure 4—figure supplement 1). Such voids in the actin–spectrin skeleton have also been observed by Stochastic Optical Reconstruction Microscopy (STORM) studies on RBCs (Pan et al., 2018). Plotting the average MSD of the first 5 s of Piezo1 tracks yielded an apparent macroscopic diffusion coefficient of 0.0075 μm² s⁻¹ (Figure 4E), which is similar to that of the mobile population of Band3 (Tomishige et al., 1998) and likely reflects occasional localized periods of confinement within actin–spectrin corrals over time. Microscopic diffusion was also analyzed by plotting the average MSD of the first 50 ms of Piezo1 tracks before the channel will tend to become localized within confinement zones. This analysis yielded an apparent diffusion coefficient of 0.037 μm² s⁻¹, consistent with measurements of freely diffusing Piezo1 in mouse neural stem progenitor cells (Ellefsen et al., 2019). The diffusion coefficient of a free ~50-nm gold nanoparticle in solution has been reported as ~1 μm² s⁻¹ (Giorgi et al., 2019), two orders of magnitude faster than the estimated microscopic diffusion coefficient of gold-labeled Piezo1 in our experiments. It is possible that a 40-nm gold particle may somewhat slow diffusion of membrane-embedded Piezo1, leading to an underestimation of its diffusion coefficient by our calculations. Nevertheless, Piezo1 appears laterally mobile within the RBC membrane.

Figure 4 with 2 supplements see all

Download asset Open asset

Gold-particle tracking of Piezo1 reveals its lateral mobility in the red blood cell membrane.

(A) Example of a RBC imaged under DIC microscopy showing three 40-nm gold particle-labeled Piezo1 channels. Scale bar, 2 μm. (B) Close-up image of the gold particles with trajectories from a 2-min recording at 100 Hz frequency. Scale bar, 1 μm. (C) Representative trajectories corresponding to distinct Piezo1 diffusive behaviors in unfixed conditions, shown in blue (trajectories U1, U2, and U3). The green trajectories correspond to labeled Piezo1 in cells fixed with paraformaldehyde (PFA). (D) MSD against time for the first 60 s of the trajectories U1, U2, and U3 in C corresponding to unfixed conditions (blue) and an average of n = 5 fixed trajectories (green). (E) MSD against time for the first 5 s of an average of 14 trajectories. The data are fit to a straight line with a slope corresponding to a 2D diffusion coefficient of 0.0075 μm² s⁻¹. R² for linear fit to data is 0.99. (F) MSD against time for the first 50 ms of an average of 14 trajectories. The data are fit to a straight line with a slope corresponding to a 2D diffusion coefficient of 0.037 μm² s⁻¹. R² for linear fit to data is 0.98. (G) Image of a RBC with a single gold-labeled Piezo1 channel that was recorded for 1 hr at 1 Hz frequency. Superimposed are the final tracking result (blue) and a circle (red) indicating the approximate segmentation of the RBC dimple. Scale bar, 1 μm.

Video 1

Download asset

posterframe for video — Example video recording of a RBC from a *Piezo1^HA* mouse immunolabeled with three 40-nm gold particles and imaged under DIC.

The recording was made at 100 Hz frequency and lasted 2 min. Scale bar, 1 µm. Video is sped up 2x.

Based on our imaging experiments of fixed RBCs that show an enrichment of Piezo1 in the dimple (Figure 3), we tracked Piezo1’s preference for the dimple versus rim in real time by recording for an hour at 1-Hz frequency (Figure 4G and Video 2). We approximated the dimple as a circle with radius 1.13 μm, comparable to our segmentation based on F-actin staining (Figure 3) and calculated the frequency at which gold-labeled Piezo1 is within 1.13 μm of an origin defined as the center of the RBC (Figure 4G). We find that Piezo1 is observed within the dimple 81% of the time during this 1 hr recording, consistent with a preference of Piezo1 for the dimple over the rim. This observed preference for the RBC dimple was not due to slower diffusion here: analysis of tracks that could easily be segmented as at the dimple or the rim yielded macroscopic diffusion coefficients of 0.0086 and 0.0074 μm² s⁻¹, respectively, which were not statistically different (Figure 4—figure supplement 2). If Piezo1 were to bind to specific sites in the dimple we would expect the diffusion coefficient to be smaller there. As the diffusion coefficient is not smaller, we can understand the observed bias toward the dimple as resulting from curvature coupling between Piezo1’s shape and the membrane’s curvature.

Video 2

Download asset

Piezo1 is not bound to the actin–spectrin meshwork

To further explore the behavior of Piezo1 in RBC membranes we analyzed its spatial proximity to the actin–spectrin meshwork by fluorescence microscopy. Proteins including ankyrin and protein 4.1 bind both to components of the meshwork and to proteins in the membrane, most commonly Band3 proteins, organizing the RBC membrane into functional microdomains composed of channels and transporters tethered to the actin–spectrin skeleton (Lux, 2016). Using co-immunostaining of unroofed RBC membranes, we asked whether Piezo1 is associated with these well-known complexes. Because of the small size of RBCs and the density of actin–spectrin at the membrane, super-resolution microscopy was necessary to address this question. Indeed, others have claimed that Piezo1 does colocalize with spectrin in red blood cells (Dumitru et al., 2021), but this study used diffraction-limited confocal microscopy, whereby every red blood cell membrane protein would be expected to show some degree of colocalization with spectrin. We employed 2D STED to image unroofed RBCs. To analyze spatial proximity of two fluorescent signals in our STED datasets we utilized optimal transport colocalization (OTC) (Tameling et al., 2021) analysis, which is more appropriate for STED microscopy than traditional pixel-based methods. First, we validated our experimental approach by co-immunostaining proteins of known complexes at the RBC membrane. There are three distinct known populations of Band3: one linked to spectrin via ankyrin, another bound to actin via protein 4.1 and other bridging proteins, and finally proteins that remain untethered to the actin–spectrin meshwork (Lux, 2016; Kodippili et al., 2009). A significant fraction of Band3 fluorescent puncta are observed as overlapping or directly adjacent to both F-actin and spectrin puncta by STED microscopy, indicating close spatial proximity (Figure 5A). In contrast, we do not observe a significant spatial proximity of Band3 and Piezo1, suggesting that if Piezo1 is tethered to actin–spectrin it is not as part of the Band3 complex. This is quantified by OTC analysis, where we plot the fraction of pixel intensities that is matched (colocalized) at distances less than the defined threshold distance. Following Nyquist sampling, the pixel size of 12.5 nm served as the lower limit of this threshold and the upper limit was set to 300 nm (Figure 5B). At a threshold of 100 nm, for example, 0.6 or 60% of signals from Band3 complexes are considered proximal enough to specify as colocalized with both actin and spectrin, consistent with the two aforementioned well-known populations of Band3 complexes. By contrast, only 0.15 or 15% of Piezo1 signals are considered proximal enough to specify as colocalized with Band3. We also co-immunostained Piezo1 with F-actin or spectrin to ask whether it is tethered to the RBC membrane skeleton. We do not observe significant spatial proximity between Piezo1 to either F-actin or spectrin by STED microscopy (Figure 5A). At a threshold of 100 nm we observe an OTC score of only 0.1 for Piezo (Figure 5B). Thus, our data show that unlike Band3, Piezo does not tend to exhibit close spatial proximity to actin or spectrin and is unlikely to bind to the membrane skeleton, consistent with its observed mobility in tracking experiments.

Figure 5

Download asset Open asset

Piezo1 does not appear to interact with the membrane skeleton.

(A) Exemplary STED microscopy images assessing interactions of Band3 (upper panels) and Piezo1 (lower panels) with known members of the RBC membrane skeleton as well as with each other and with KCNN4. Close spatial proximity for Band3 with actin and spectrin can be observed in STED images but not for Piezo1. Scale bar, 200 nm. (B) OTCanalysis including 95% confidence error bars of STED microscopy image sections (Band3/F-actin, n = 11; Band3/Spectrin, n = 10; Band3/Piezo1, n = 7; Piezo1/F-actin, n = 10; Piezo1/Spectrin, n = 7; Piezo1/KCNN4, n = 8).

Based on a proposed mechanism of Piezo1-mediated Ca²⁺ influx activating the Gardos channel (Cahalan et al., 2015) (KCNN4), we also analyzed spatial proximity of these two proteins by STED. Piezo1 and KCNN4 do not appear to form complexes in RBCs (Figure 5A, B), but it should be noted, given the small size of the RBC, that colocalization of Piezo1 and KCNN4 is not necessary for a functional coupling of the ion fluxes associated with these two channels (see Discussion).

Discussion

Since its discovery (Coste et al., 2010), Piezo1 has been identified as a primary mechanosensitive molecule in a wide range of force-sensing physiological processes from cell and organ development (Sun et al., 2019; Li et al., 2014; Ranade et al., 2014) to blood pressure regulation (Zeng et al., 2018) and stem cell lineage choice (Pathak et al., 2014). How Piezo1 is organized in different cell types and how this might pertain to the sensing of different types of physiological force is currently not well understood. Here we studied Piezo1 organization in RBCs, both because of the clear importance of Piezo1 in RBC physiology (Bae et al., 2013; Zarychanski et al., 2012; Cahalan et al., 2015) and because of the unique cellular structure of RBCs, particularly their biconcave disk shape.

At the smallest level of Piezo1 organization, that is with itself, we find that Piezo1 does not cluster in the RBC membrane. By 3D-SIM microscopy (Figure 1) we observe diffraction-limited fluorescent spots that could technically represent multiple, closely located, Piezo1 channels. To distinguish these spots as individual channels, we needed to employ 2D STED microscopy (Figure 2). As the power of the depletion laser was increased to improve resolution, the size of each spot decreased correspondingly, but never resolved into multiple spots, even when the achieved resolution was near the size of the Piezo1 channel. Further, after Huygens deconvolution post-processing of our 2D STED data we could resolve triple-antibody-labeled Piezo1 channels, with an inter-subunit diameter consistent with the channel size. Thus, we conclude that each fluorescent spot, whether viewed by 3D-SIM or 2D STED, corresponds to a single Piezo1 channel. The average distance between Piezo1 channels is approximately 550 nm, although occasionally we observe two to three channels within 100-nm distances. This is in good agreement with studies of neuro2A cells, which natively express Piezo1, that also found that Piezo1 channels do not tend to cluster (Lewis and Grandl, 2021). While clustering could affect Piezo1 gating through Piezo’s membrane footprint, the precise coupling of Piezo1 gating to clustering would depend on the specific scenario considered, and clustering does not necessarily facilitate gating (Haselwandter and MacKinnon, 2018). A relatively homogenous membrane distribution of Piezo1 might be well suited to RBCs that are tumbling through solution, experiencing shear forces from all angles, and may contribute to the remarkable plasticity of RBC shapes under stress (Dupire et al., 2012). Of course, our data on RBCs do not rule out the possibility of Piezo1 clustering in other cell types.

Because of the intrinsic curvature of Piezo1 we were curious about its distribution as a function of the curvature of the RBC membrane. Experimentally we were able to ask this question through 3D-SIM microscopy, which yielded sufficiently high-resolution spatial information on Piezo1 channels while enabling easy assignment of the RBC dimple and rim membrane regions (Figure 3). Per μm² area of membrane we find an almost twofold enrichment of Piezo1 in the dimple relative to the rim. Our experimental data can be rationalized by a physical model of curvature coupling between the intrinsic curvature of the Piezo1 dome and the surrounding membrane curvature of the RBC along its biconcave surface (see Methods). Whilst the curved biconcave shape of RBCs is unique, the good agreement between our experimental data and theoretical prediction based on energetic curvature coupling would suggest an inherent relationship between Piezo’s intrinsic curvature and the surrounding membrane that is broadly relevant to other cell types. This is because the curvature coupling theory, which was developed with experiments on lipid bilayer vesicles, depends only on the bending elasticity of membranes and the intrinsic curvature of Piezo1, not on the unique shape of RBCs. It will be interesting to explore whether Piezo1, for example, is enriched in membrane invaginations or other curvature-matched membrane regions in different cell types. Recently, a study analyzing overexpressed Piezo1 in HeLa cells showed depletion of the channel from membrane filopodia, consistent with a mechanism of curvature coupling whereby the channel unfavorably distributes into membrane regions of high opposing curvature (Yang et al., 2022).

Curvature sensing of membrane proteins by curvature matching has been well documented for other intrinsically curved proteins such as Bin/Amphiphysin/Rvs (BAR) (Simunovic et al., 2015; Tsai et al., 2021) proteins. More recently, it has been shown that glycosylation can determine the membrane curvature sensing properties of some proteins (Lu et al., 2022). It is unclear from our current data whether the enrichment of Piezo1 in the RBC dimple is simply a corollary of the curved structure of Piezo1, which determines the channel’s tension sensitivity as described by the membrane dome model (Guo and MacKinnon, 2017; Haselwandter and MacKinnon, 2018), or whether it also has some important consequences for physiological force sensation of RBCs. In this context, it is noteworthy that the force-generating protein myosin IIA also shows an enrichment in the RBC dimple relative to the RBC rim (Alimohamadi et al., 2020).

Consistent with its ability to equilibrate along the RBC surface as a function of membrane curvature, we find that Piezo1 exhibits lateral mobility within the RBC plasma membrane. We used gold nanoparticle tracking under DIC (Figure 4) so we could record for minutes to an hour without being limited by photobleaching and at fast time resolution (100 Hz). This enabled us to record both the confined (<500 nm) diffusion of Piezo1 as well as Piezo1’s ‘hopping’ into more freely diffusing behavior. The spatial scales associated with the confined diffusion regime most likely reflect the compartment sizes of the actin–spectrin meshwork we observe by 2D STED. Such a model of confined diffusion set by actin–spectrin is well documented in RBCs (Tsuji et al., 1988) and other cell types (Fujiwara et al., 2016). Similar to observations by others using STORM microscopy (Pan et al., 2018), we observe ~200–500 nm voids that are deficient in actin or spectrin signal by 2D STED. Similar nanoscale ‘defects’ in the actin–spectrin meshwork have also been observed by atomic force microscopy imaging of RBCs (Nowakowski et al., 2001) and may serve as weak points in the meshwork that facilitate fast RBC remodeling during shear stress. Such defects may also permit faster diffusion of transmembrane proteins such as Piezo1.

Our double labeling experiments suggest that Piezo1 does not form a complex with actin or spectrin, which is consistent with Piezo1’s confined diffusion within the actin–spectrin meshwork compartments, and should be contrasted with the near-immobilization observed for proteins such as CFTR and Band3 (Kodippili et al., 2009, 3; Haggie et al., 2006), which both form complexes with actin. Additionally, our data are not consistent with Piezo1 being tethered to the underlying actin and activated by a ‘force-from-filament’ mechanism as has been traditionally described (Cox et al., 2019; Zhang et al., 2015). Of course, the actin–spectrin meshwork organizes the RBC membrane into compartments whose boundary conditions could impact how changes in local curvature and lateral tension of the membrane gate Piezo1. This could also explain how chemical agents that disrupt the actin cytoskeleton affect Piezo1 currents in other cell types (Gottlieb et al., 2012; Wang et al., 2022), that is through changes in organization of the membrane rather than disruption of a direct tether to the channel itself.

Additionally, we asked whether Piezo1 might be colocalized with the Ca²⁺-activated K⁺ (Gardos) channel. Ca²⁺ flux into the RBC through Piezo1 and subsequent activation of the Gardos channel is thought to be the initial step in mechanically activated RBC volume regulation (Danielczok et al., 2017; Cahalan et al., 2015). We do not, however, observe significant spatial proximity between Piezo1 and Gardos channels. Given the small size of RBCs, approximately 1 μm thick in the dimple and 2–3 μm thick in the rim, and the Piezo1 density of around 0.5 per μm², it seems that colocalization is probably not necessary for rapid opening of the Gardos channel following Piezo activation. For diffusion in three dimensions, the mean dispersion time can be approximated by $τ = \frac{x^{2}}{6 D}$ , $x$ the mean distance and $D$ the diffusion coefficient. The apparent diffusion coefficient for intracellular Ca²⁺ is in the range 13–65 μm²/s (Nakatani et al., 2002), meaning that Ca²⁺ should equilibrate within tens of ms within a RBC. In larger cell types, colocalization of Piezo1 with other channels might be necessary to elicit a functional coupling of ion fluxes. In rat ventricular myocytes, for example, Piezo1 has recently been shown, by diffraction-limited confocal microscopy, to colocalize with transient receptor melastatin (TRPM4) channels (Yu et al., 2022).

The experimentally determined values of the intrinsic curvature and the bending stiffness of Piezo1 in lipid bilayers enable us to estimate changes in Piezo1 curvature at the dimple versus the rim of a RBC membrane (Haselwandter et al., 2022a). We estimate a radius of curvature of 41.4 nm at the dimple and 42.7 nm at the rim and thus the curvature differences across the unperturbed RBC membrane are, alone, insufficient to substantially change the shape of Piezo1, and thus are not expected to open its pore. However, that Piezo1 distributes on the surface of RBCs by sensing the membrane curvature, by corollary, means that it must have the capacity to sense lateral membrane tension by a force-through-membrane mechanism. To see this, consider that Piezo1 has been shown to flex, that is, change its curvature, as a function of local membrane curvature (Haselwandter et al., 2022b). Consider also that membrane tension changes membrane curvature by flattening it. In this manner, curvature sensing and lateral membrane tension sensing are inextricably linked. This linkage lies at the heart of the membrane dome model of Piezo gating (Guo and MacKinnon, 2017; Haselwandter and MacKinnon, 2018; Haselwandter et al., 2022b; Haselwandter et al., 2022a). When RBCs squeeze through small capillaries, the membrane shape as well as the transmembrane pressure change transiently. The resulting changes in membrane curvature and lateral membrane tension can both, according to the membrane dome model of Piezo gating, open Piezo1’s pore to trigger Ca²⁺ entry.

Methods

HA-tag knock in mice

C57BL/6 background mice heterozygous for an HA tag knocked into the coding sequence after amino acid position 893 of Piezo1, generated by CRISPR/Cas9 technology, were purchased from Jackson Laboratories. Heterozygous mice were bred to generate mice homozygous for the HA tag knock in and were genotyped by PCR amplification around the 893 insertion site and gel electrophoresis analysis of the amplified band size.

Preparation of blood

Adult mice (>3 months old) were anesthetized by isoflurane inhalation and retro-orbital bleeds were taken using a heparinized micro-hematrocrit capillary tube (Fisherbrand). Whole blood was washed in Dulbecco’s phosphate-buffered saline (DPBS) three times by centrifugation at 1000 × g to separate RBCs from the buffy coat.

Electrophysiological recordings

Washed RBCs were diluted 1:1000 and plated on uncoated Petri dishes containing bath solution consisting of 140 mM NaMethanesulfonate, 10 HEPES pH 7.4. RBCs were imaged with a Nikon eclipse DIC microscope at ×40 magnification. Pipettes of borosilicate glass (Sutter Instruments; BF150-86-10) were pulled to ~15–20 MΩ resistance with a micropipette puller (Sutter Instruments; P-97) and polished with a microforge (Narishige; MF-83). The pipette was filled with identical bath solution. Recordings were obtained with an Axopatch 200B amplifier (Molecular Devices), filtered at 1 kHz and digitized at 10 kHz (Digidata 1440A; Molecular devices). Gigaseals were obtained by applying a suction pulse (−10 mm Hg) using a high-speed pressure clamp (ALA scientific). Holding at negative voltages often facilitated seal formation. To obtain excised inside-out patches, attached cells were lifted briefly to the air–water interface which about 50% of the time would remove the cell but leave an intact inside-out membrane patch in the pipette.

Immunofluorescence of intact RBCs using 3D-SIM

Washed RBCs were deglycosylated by incubation with 10% PNGase F (New England Biolabs) at 37°C for 2 hr, which we found increased antibody binding. RBCs were then diluted 1:50 in DPBS with 4% PFA and incubated overnight at room temperature to fix. Fixed cells were washed three times in DPBS by centrifugation at 1000 × g for 5 min and then permeabilized in DPBS + 0.3% Triton X-100 for 15 min. Permeabilized cells were then blocked in DPBS + 4% bovine serum albumin (BSA) + 1% normal goat serum (blocking buffer) overnight at 4°C. Permeabilized and blocked RBCs were then incubated with a rabbit anti-HA primary antibody (3724; Cell Signaling Technology) diluted in blocking buffer for 1 hr at room temperature, washed three times in blocking buffer and then incubated in 1:500 dilution of Alexa-488-conjugated goat anti-rabbit secondary (ab150077; Abcam) mixed with ×1 rhodamine–phalloidin (Invitrogen) for 1 hr, followed by washing three times in blocking buffer. For 3D-SIM imaging of KCNN4 and Band3 in RBCs, 1:100 dilution of rabbit anti-Band3 primary (Proteintech) and 1:500 dilution of rabbit anti-KCN44 (76647, Invitrogen) were used. 1:500 dilution of Alexa-488-conjugated goat anti-rabbit (ab150077; Abcam) was used as the secondary along with 1× rhodamine–phalloidin (Invitrogen).

3D-SIM images were acquired using a DeltaVision OMX V4/Blaze system (Cytiva) fitted with an Olympus ×100/1.40 NA UPLSAPO oil objective and Photometrics Evolve EMCCD cameras. 488 and 568 nm laser lines were used for excitation and the corresponding emission filters sets were 528/48 and 609/37 nm, respectively. Image stacks were acquired with an optical section spacing of 125 nm. SI reconstruction and Image Registration were performed with softWoRx v 6.1 software using Optical Transfer Functions (OTFs) generated from Point Spread Functions (PSFs) acquired from 100 nm green and red FluoSpheres and alignment parameters refined using 100 nm TetraSpeck beads (Invitrogen).

Imaris version 9.9 was used for analyzing 3D-SIM images. Detection of Piezo1 fluorescent spots was performed with Imaris spot detection with an estimated XY diameter of 120 nm, based on measurements of fluorescent spots in 2D slices and the approximate XY resolution limit of 3D-SIM experiments. Spots were mostly uniform in size with a maximum XY diameter of 180 nm. The intensity threshold for spot detection was set to between 3000 and 5000 based on excluding detection of weak background fluorescent signal on the coverslip surface.

Echinocyte formation

Washed RBCs that were degclyosylated were incubated in DPBS + 10 mM NaSalicylate for 30 min at 37°C to generate echinocytes before fixation and immunostaining as described above.

Immunofluorescence of unroofed RBCs using STED microscopy

Washed RBCs that were deglycosylated as above were plated onto poly-D-lysine coated high-performance coverslips and allowed to adhere at room temperature for 15 min. Plated RBCs were unroofed by applying a stream of 10 ml DPBS through a 20 gauge needle at a ~20^o angle (Swihart et al., 2001). Unroofed RBCs were fixed in DPBS with 4% PFA for 15 min and then washed three times in DPBS before blocking in blocking buffer for 30 min. This was followed by immunostaining with 1:500 dilution of rabbit anti-HA primary (3724; Cell Signaling Technology) for 1 hr at room temperature, washing three times in blocking buffer and then incubating in 1:500 dilution of anti-rabbit STAR RED secondary (Abberior) mixed with 1× phalloidin STAR 580 (Abberior) for 1 hr. Immunostained unroofed cells were then, again, washed three times in blocking buffer and mounted in uncured Prolong Diamond Antifade mountant (Invitrogen) before imaging. For co-immunostaining experiments, additional antibodies used were as follows: 1:100 dilution of rabbit anti-Band3 primary (Proteintech), 1:100 dilution of mouse anti-alpha 1 Spectrin (ab11751), and 1:500 dilution of rabbit anti-KCNN4 (76647, Invitrogen). For co-immunostaining experiments with rabbit primary antibodies, rat anti-HA (Roche) primary was used to label Piezo1 and in all cases the species-appropriate STAR RED and STAR 580 (Abberior) secondaries were used.

STED microscopy was performed using a Facility Line STED microscope (Abberior Instruments) equipped with Olympus IX83 stand, Olympus UPLXAPO ×100/1.45 NA oil objective, pulsed excitation lasers with time gating (405, 488, 561, and 640 nm) and a 775-nm pulsed STED depletion laser, 4 Avalance Photo Diode detectors, adaptive illumination packages (DyMIN/ RESCue), and a deformable mirror for correction of spherical aberrations. Abberior Imspector software version 16.3.14287-w2129 with Lightbox interface was used for image acquisition. Fluorophore excitation was facilitated at 580 nm (Phalloidin STAR 580) and 640 nm (STAR RED), whereas depletion was achieved using a wavelength of 775 nm. A pixel size of 12.5 nm was used in 2D STED mode. Excitation laser power, depletion laser power, line averaging/accumulation, and pixel dwell time were optimized for balancing signal-to-noise ratio (SNR), and resolution, while minimizing photobleaching. Adaptive illumination techniques RESCue (Staudt et al., 2011) and DyMIN (Heine et al., 2017) were employed to reduce photobleaching and increase SNR.

Acquired STED images were saved as OBF files, which were deconvolved using theoretical PSF modeled from microscopy parameters in Huygens Professional version 22.04 (Scientific Volume Imaging). Imported images were cropped as needed, microscopic parameters were edited for refractive index mismatch correction, the excitation fill factor was set to 1.38 per software developer recommendations, and the STED saturation factor was tested within a range of 40–80, as per vendor recommendations, for optimal resolution. The deconvolution wizard was used to estimate background manually in the raw images and the acuity option was set for optimal image sharpness. Deconvolved images were exported as 16-bit TIFF files.

Cluster analysis on 2D STED images was performed using the spatial statistics 2D/3D ImageJ plugin (Andrey et al., 2010).

Image analysis in Imaris

Imaris version 9.9 was used for analyzing fluorescent images. Detection of Piezo1 fluorescent spots was performed with Imaris spot detection with an estimated XY diameter of 120 nm, based on measurement of fluorescent spots in 2D slices and the approximate XY resolution limit of 3D-SIM experiments.

Negative stain electron microscopy of unroofed RBCs

Washed and deglycosylated RBCs were diluted 1:100 and added to carbon coated 400 mesh copper grids (CF400-CU, Electron Microscopy Sciences) that had been glow discharged and pre-coated in poly-D-lysine solution. After 30 min incubation, RBCs were unroofed in DPBS as described above. Unroofed RBCs were blocked for 30 min at room temperature in blocking buffer and then incubated with 1:200 rabbit anti-HA primary antibody (3724; Cell Signaling Technology) for 1 hr at room temperature. Unroofed cells were then washed three times in blocking buffer and incubated with 1:20 diluted 18-nm gold-conjugated goat anti-rabbit secondary antibody (Jackson Laboratories) for 1 hr at room temperature. After rinsing three times in blocking buffer, grids were then washed three times in 150 NaCl, 20 HEPES 7.5 to remove phosphate before staining in 1% uranyl acetate for 1 min. Images were collected on a Tecnai G2 Spirit BioTWIN Transmission Electron Microscope.

Probability distribution calculation of Piezo1 along the RBC membrane

The Beck RBC model was used to describe the surface of a biconcave disk with diameter 8.0 μm, dimple thickness 1.0 μm, maximum rim thickness 2.5 μm, and area 138 μm² (Beck, 1978). The mean curvature at each point on this surface was calculated as $H (x) = \frac{1}{2} (c_{1} + c_{2})$ , where $c_{1}$ and $c_{2}$ are the principal curvatures, which are functions of the position $x$ on the surface (Carmo, 2016). The following calculation was carried out to assign a (relative) membrane energy, and associated probability density, to Piezo1 at each position $x$ on the RBC surface. The aim of this calculation is not to provide a detailed description of RBC membrane shape but, rather, to test whether Piezo1 localization can couple to curvatures on the RBC surface that are 1–2 orders of magnitude smaller than the Piezo dome curvature. In particular, Piezo1, modeled as a spherical cap (dome) with area 450 nm² and radius of curvature 42 nm (Haselwandter et al., 2022a), was placed into a vesicle whose radius of curvature, if it were spherical, equals $1 / H (x)$ . We note that the above conditions specify the area of the vesicle, which equals the free membrane area plus 450 nm². For positions on the RBC surface at which the orientation (sign) of the RBC curvature matches the orientation of the Piezo1 intrinsic curvature (concave side facing outside the cell) Piezo1 was placed into the vesicle with an inside-out orientation, that is, matching the orientation of the vesicle curvature. Otherwise, Piezo1 was placed into the vesicle with an outside-out orientation. The energy of the vesicle free membrane (outside the perimeter of the Piezo dome) was calculated through minimization of the Helfrich membrane bending energy

G_{M} = \frac{K_{b}}{2} \int d A {(c_{1} + c_{2})}^{2},

where the lipid bilayer bending modulus $K_{b} \approx 20 k_{B} T$ and the integral is carried out over the entire vesiclefree membrane, subject to the boundary conditions at the Piezo dome-free membrane interface mandated by the shape of the Piezo dome (Helfrich, 1973; Haselwandter et al., 2022b; Haselwandter et al., 2022a). The resulting energy as a function of RBC curvature (and, thus, as a function of the position $x$ on the RBC surface) was applied to the Boltzmann distribution equation to calculate the probability per unit area of observing Piezo1 across the RBC surface. This probability density is graphed on the surface of the Beck RBC model in Figure 3. Three major assumptions underlie our calculation of the curvature coupling energy. First, Piezo1 is treated as a spherical cap. Second, for each position on the RBC surface the two principal curvatures, $c_{1}$ and $c_{2}$ , are used to define a radius of curvature equal to the reciprocal of the mean curvature. Both the first and second assumption simplify the problem by imposing radial symmetry on the system of the Piezo dome and its surrounding membrane. The third assumption is that Piezo1 is sensing the curvature associated with the overall shape of the RBC, rather than the curvature of a locally specified surface (small membrane compartment) with nearby fixed boundaries.

We also carried out a joint energy minimization of the free membrane as described above while simultaneously permitting Piezo1 to flex, that is, to change its radius of curvature, as part of the energy minimization. This calculation was made possible through our recently determined bending modulus of the Piezo dome itself (Haselwandter et al., 2022a). This calculation leads to the conclusion that Piezo1 in an unperturbed RBC should change its mean dome radius of curvature by less than approximately 2 nm as a function of position on the RBC surface.

Single particle tracking and diffusion analysis

Washed RBCs were adhered to poly-lysine coated 35 mm glass-bottom dishes (MatTeK) and blocked in blocking buffer before incubation with 1:500 dilution of rabbit anti-HA primary (3725; Cell Signaling Technology) for 1 hr at room temperature. Cells were then washed in blocking buffer and incubated with 1:20 diluted 40-nm gold-conjugated goat anti-rabbit secondary antibody (Jackson Laboratories) for 1 hr at room temperature. Cells were washed in blocking buffer and then DPBS before imaging. The movement of colloidal gold particles was observed at room temperature using DIC microscopy on an eclipse Ti2 inverted microscope (Nikon) using a ×40/0.65 objective lens and coupled to an Orca Fusion CMOS camera (Hamamatsu). Video sequences were recorded on Nikon NIS-Elements AR5.4.1 software.

Detection and tracking were carried out in ImageJ imaging analysis software using the TrackMate plugin (Tinevez et al., 2017; Ershov et al., 2021). First, trainable WEKA segmentation (Arganda-Carreras et al., 2017) was employed to generate a model of gold nanoparticles that could be used for spot detection. Detection was then carried out in TrackMate using the Linear Assignment Problem (LAP) tracker. Tracks were exported as XYT coordinates and further analyzed in MATLAB. The MATLAB class msdanalyzer (Tarantino et al., 2014) was used for MSD analysis of trajectories and estimation of diffusion coefficients.

Spatial proximity analysis by OTC

Co-immunostained STED images were analyzed using an open-source R code package (Tameling et al., 2021; Tameling and Naas, 2021) based on OTC. Sections of 64 × 64 pixels (800 × 800 nm) were picked only if they contained fluorescent signal and were analyzed in parallel to get 95% confidence bands.

Data availability

Data generated and analyzed are included in main figures and figure supplements of manuscript. Source files for raw data used for tracking analysis in Figure 4 and Figure 4—figure supplement 1 have been deposited to Dryad: https://doi.org/10.5061/dryad.c2fqz61c7. Source files for raw data used for optimal transport colocalization analysis in Figure 5 have been deposited to Dryad: https://doi.org/10.5061/dryad.1ns1rn8x6.

The following data sets were generated

1. Mackinnon R
(2022) Dryad Digital Repository
Data from: Piezo1 as a force-through-membrane sensor in red blood cells.
https://doi.org/10.5061/dryad.c2fqz61c7
1. Mackinnon R
(2022) Dryad Digital Repository
Data from: Piezo1 as a force-through-membrane sensor in red blood cells.
https://doi.org/10.5061/dryad.1ns1rn8x6

References

(2013) Dehydrated hereditary stomatocytosis linked to gain-of-function mutations in mechanically activated Piezo1 ion channels
Nature Communications 4:1884.
https://doi.org/10.1038/ncomms2899
- PubMed
- Google Scholar
(2020) Non-Uniform distribution of myosin-mediated forces governs red blood cell membrane curvature through tension modulation
PLOS Computational Biology 16:e1007890.
https://doi.org/10.1371/journal.pcbi.1007890
- PubMed
- Google Scholar
1. Andolfo I
2. Alper SL
3. De Franceschi L
4. Auriemma C
5. Russo R
6. De Falco L
7. Vallefuoco F
8. Esposito MR
9. Vandorpe DH
10. Shmukler BE
11. Narayan R
12. Montanaro D
13. D’Armiento M
14. Vetro A
15. Limongelli I
16. Zuffardi O
17. Glader BE
18. Schrier SL
19. Brugnara C
20. Stewart GW
21. Delaunay J
22. Iolascon A
(2013) Multiple clinical forms of dehydrated hereditary stomatocytosis arise from mutations in Piezo1
Blood 121:3925–3935.
https://doi.org/10.1182/blood-2013-02-482489
- PubMed
- Google Scholar
1. Andrey P
2. Kiêu K
3. Kress C
4. Lehmann G
5. Tirichine L
6. Liu Z
7. Biot E
8. Adenot PG
9. Hue-Beauvais C
10. Houba-Hérin N
11. Duranthon V
12. Devinoy E
13. Beaujean N
14. Gaudin V
15. Maurin Y
16. Debey P
(2010) Statistical analysis of 3D images detects regular spatial distributions of centromeres and chromocenters in animal and plant nuclei
PLOS Computational Biology 6:e1000853.
https://doi.org/10.1371/journal.pcbi.1000853
- PubMed
- Google Scholar
1. Arakawa T
2. Kobayashi-Yurugi T
3. Alguel Y
4. Iwanari H
5. Hatae H
6. Iwata M
7. Abe Y
8. Hino T
9. Ikeda-Suno C
10. Kuma H
11. Kang D
12. Murata T
13. Hamakubo T
14. Cameron AD
15. Kobayashi T
16. Hamasaki N
17. Iwata S
(2015) Crystal structure of the anion exchanger domain of human erythrocyte band 3
Science 350:680–684.
https://doi.org/10.1126/science.aaa4335
- PubMed
- Google Scholar
(2017) Trainable weka segmentation: a machine learning tool for microscopy pixel classification
Bioinformatics 33:2424–2426.
https://doi.org/10.1093/bioinformatics/btx180
- PubMed
- Google Scholar
(2013) Xerocytosis is caused by mutations that alter the kinetics of the mechanosensitive channel piezo1
PNAS 110:E1162–E1168.
https://doi.org/10.1073/pnas.1219777110
- PubMed
- Google Scholar
1. Ballas SK
2. Krasnow SH
(1980)
Structure of erythrocyte membrane and its transport functions

Annals of Clinical and Laboratory Science 10:209–219.
- PubMed
- Google Scholar
1. Beck JS
(1978) Relations between membrane monolayers in some red cell shape transformations
Journal of Theoretical Biology 75:487–501.
https://doi.org/10.1016/0022-5193(78)90358-2
- PubMed
- Google Scholar
(1993) Ca (2+) -activated K+ transport in erythrocytes. Comparison of binding and transport inhibition by scorpion toxins
The Journal of Biological Chemistry 268:8760–8768.
https://doi.org/10.1016/S0021-9258(18)52939-2
- PubMed
- Google Scholar
1. Bryk AH
2. Wiśniewski JR
(2017) Quantitative analysis of human red blood cell proteome
Journal of Proteome Research 16:2752–2761.
https://doi.org/10.1021/acs.jproteome.7b00025
- PubMed
- Google Scholar
1. Cahalan SM
2. Lukacs V
3. Ranade SS
4. Chien S
5. Bandell M
6. Patapoutian A
(2015) Piezo1 links mechanical forces to red blood cell volume
eLife 4:e07370.
https://doi.org/10.7554/eLife.07370
- PubMed
- Google Scholar
Book
1. Carmo M
(2016)
Differential Geometry of Curves & Surfaces

Dover Publications, INC.
- Google Scholar
1. Chien S
(1987) Red cell deformability and its relevance to blood flow
Annual Review of Physiology 49:177–192.
https://doi.org/10.1146/annurev.ph.49.030187.001141
- PubMed
- Google Scholar
1. Coste B
2. Mathur J
3. Schmidt M
4. Earley TJ
5. Ranade S
6. Petrus MJ
7. Dubin AE
8. Patapoutian A
(2010) Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels
Science 330:55–60.
https://doi.org/10.1126/science.1193270
- PubMed
- Google Scholar
1. Cox CD
2. Bae C
3. Ziegler L
4. Hartley S
5. Nikolova-Krstevski V
6. Rohde PR
7. Ng C-A
8. Sachs F
9. Gottlieb PA
10. Martinac B
(2016) Removal of the mechanoprotective influence of the cytoskeleton reveals Piezo1 is gated by bilayer tension
Nature Communications 7:10366.
https://doi.org/10.1038/ncomms10366
- PubMed
- Google Scholar
1. Cox CD
2. Bavi N
3. Martinac B
(2019) Biophysical principles of ion-channel-mediated mechanosensory transduction
Cell Reports 29:1–12.
https://doi.org/10.1016/j.celrep.2019.08.075
- PubMed
- Google Scholar
(2017) Red blood cell passage of small capillaries is associated with transient Ca2+-mediated adaptations
Frontiers in Physiology 8:979.
https://doi.org/10.3389/fphys.2017.00979
- PubMed
- Google Scholar
1. Daumas F
2. Destainville N
3. Millot C
4. Lopez A
5. Dean D
6. Salomé L
(2003) Confined diffusion without fences of a G-protein-coupled receptor as revealed by single particle tracking
Biophysical Journal 84:356–366.
https://doi.org/10.1016/S0006-3495(03)74856-5
- PubMed
- Google Scholar
(2018) Piezo1 forms a slowly-inactivating mechanosensory channel in mouse embryonic stem cells
eLife 7:e33149.
https://doi.org/10.7554/eLife.33149
- PubMed
- Google Scholar
1. Dumitru AC
2. Stommen A
3. Koehler M
4. Cloos A-S
5. Yang J
6. Leclercqz A
7. Tyteca D
8. Alsteens D
(2021) Probing Piezo1 localization upon activation using high-resolution atomic force and confocal microscopy
Nano Letters 21:4950–4958.
https://doi.org/10.1021/acs.nanolett.1c00599
- PubMed
- Google Scholar
(2012) Full dynamics of a red blood cell in shear flow
PNAS 109:20808–20813.
https://doi.org/10.1073/pnas.1210236109
- PubMed
- Google Scholar
1. Dyrda A
2. Cytlak U
3. Ciuraszkiewicz A
4. Lipinska A
5. Cueff A
6. Bouyer G
7. Egée S
8. Bennekou P
9. Lew VL
10. Thomas SLY
(2010) Local membrane deformations activate Ca2+-dependent K+ and anionic currents in intact human red blood cells
PLOS ONE 5:e9447.
https://doi.org/10.1371/journal.pone.0009447
- PubMed
- Google Scholar
1. Ellefsen KL
2. Holt JR
3. Chang AC
4. Nourse JL
5. Arulmoli J
6. Mekhdjian AH
7. Abuwarda H
8. Tombola F
9. Flanagan LA
10. Dunn AR
11. Parker I
12. Pathak MM
(2019) Myosin-ii mediated traction forces evoke localized piezo1-dependent ca2+ flickers
Communications Biology 2:298.
https://doi.org/10.1038/s42003-019-0514-3
- PubMed
- Google Scholar
Preprint
1. Ershov D
2. Phan MS
3. Pylvänäinen JW
4. Rigaud SU
5. Le Blanc L
6. Charles-Orszag A
7. Conway JRW
8. Laine RF
9. Roy NH
10. Bonazzi D
11. Duménil G
12. Jacquemet G
13. Tinevez JY
(2021) Bringing TrackMate into the Era of Machine-Learning and Deep-Learning
bioRxiv.
https://doi.org/10.1101/2021.09.03.458852
- Google Scholar
(1971) Electrophoretic analysis of the major polypeptides of the human erythrocyte membrane
Biochemistry 10:2606–2617.
https://doi.org/10.1021/bi00789a030
- PubMed
- Google Scholar
1. Fowler VM
(2013) The human erythrocyte plasma membrane
Current Topics in Membranes 72:39–88.
https://doi.org/10.1016/B978-0-12-417027-8.00002-7
- Google Scholar
1. Fujiwara TK
2. Iwasawa K
3. Kalay Z
4. Tsunoyama TA
5. Watanabe Y
6. Umemura YM
7. Murakoshi H
8. Suzuki KGN
9. Nemoto YL
10. Morone N
11. Kusumi A
(2016) Confined diffusion of transmembrane proteins and lipids induced by the same actin meshwork lining the plasma membrane
Molecular Biology of the Cell 27:1101–1119.
https://doi.org/10.1091/mbc.E15-04-0186
- PubMed
- Google Scholar
1. Gardos G
(1958) The function of calcium in the potassium permeability of human erythrocytes
Biochimica et Biophysica Acta 30:653–654.
https://doi.org/10.1016/0006-3002(58)90124-0
- PubMed
- Google Scholar
1. Gautier E-F
2. Leduc M
3. Cochet S
4. Bailly K
5. Lacombe C
6. Mohandas N
7. Guillonneau F
8. El Nemer W
9. Mayeux P
(2018) Absolute proteome quantification of highly purified populations of circulating reticulocytes and mature erythrocytes
Blood Advances 2:2646–2657.
https://doi.org/10.1182/bloodadvances.2018023515
- PubMed
- Google Scholar
1. Giorgi F
2. Coglitore D
3. Curran JM
4. Gilliland D
5. Macko P
6. Whelan M
7. Worth A
8. Patterson EA
(2019) The influence of inter-particle forces on diffusion at the nanoscale
Scientific Reports 9:12689.
https://doi.org/10.1038/s41598-019-48754-5
- PubMed
- Google Scholar
(2012) Gating the mechanical channel Piezo1: a comparison between whole-cell and patch recording
Channels 6:282–289.
https://doi.org/10.4161/chan.21064
- PubMed
- Google Scholar
1. Gratzer WB
(1981) The red cell membrane and its cytoskeleton
The Biochemical Journal 198:1–8.
https://doi.org/10.1042/bj1980001
- PubMed
- Google Scholar
1. Guo YR
2. MacKinnon R
(2017) Structure-Based membrane dome mechanism for Piezo mechanosensitivity
eLife 6:e33660.
https://doi.org/10.7554/eLife.33660
- PubMed
- Google Scholar
1. Gustafsson MGL
(2000) Surpassing the lateral resolution limit by a factor of two using structured illumination microscopy
Journal of Microscopy 198:82–87.
https://doi.org/10.1046/j.1365-2818.2000.00710.x
- PubMed
- Google Scholar
1. Haggie PM
2. Kim JK
3. Lukacs GL
4. Verkman AS
(2006) Tracking of quantum dot-labeled CFTR shows near immobilization by C-terminal PDZ interactions
Molecular Biology of the Cell 17:4937–4945.
https://doi.org/10.1091/mbc.e06-08-0670
- PubMed
- Google Scholar
(2022) Piezo1 is a mechanosensor channel in central nervous system capillaries
Circulation Research 130:1531–1546.
https://doi.org/10.1161/CIRCRESAHA.122.320827
- PubMed
- Google Scholar
1. Haselwandter CA
2. MacKinnon R
(2018) Piezo’s membrane footprint and its contribution to mechanosensitivity
eLife 7:e41968.
https://doi.org/10.7554/eLife.41968
- PubMed
- Google Scholar
(2022a) Elastic properties and shape of the Piezo dome underlying its mechanosensory function
PNAS 119:e2208034119.
https://doi.org/10.1073/pnas.2208034119
- PubMed
- Google Scholar
(2022b) Quantitative prediction and measurement of piezo’s membrane footprint
PNAS 119:e2208027119.
https://doi.org/10.1073/pnas.2208027119
- PubMed
- Google Scholar
1. Heine J
2. Reuss M
3. Harke B
4. D’Este E
5. Sahl SJ
6. Hell SW
(2017) Adaptive-illumination STED nanoscopy
PNAS 114:9797–9802.
https://doi.org/10.1073/pnas.1708304114
- PubMed
- Google Scholar
1. Heintzmann R
2. Huser T
(2017) Super-Resolution structured illumination microscopy
Chemical Reviews 117:13890–13908.
https://doi.org/10.1021/acs.chemrev.7b00218
- PubMed
- Google Scholar
1. Helfrich W
(1973) Elastic properties of lipid bilayers: theory and possible experiments
Zeitschrift für Naturforschung C 28:693–703.
https://doi.org/10.1515/znc-1973-11-1209
- Google Scholar
1. Hell SW
2. Wichmann J
(1994) Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy
Opt Lett 19:780–782.
https://doi.org/10.1364/ol.19.000780
- PubMed
- Google Scholar
1. Kodippili GC
2. Spector J
3. Sullivan C
4. Kuypers FA
5. Labotka R
6. Gallagher PG
7. Ritchie K
8. Low PS
(2009) Imaging of the diffusion of single band 3 molecules on normal and mutant erythrocytes
Blood 113:6237–6245.
https://doi.org/10.1182/blood-2009-02-205450
- PubMed
- Google Scholar
(1981) Physiological shear stresses enhance the Ca2+ permeability of human erythrocytes
Nature 294:667–668.
https://doi.org/10.1038/294667a0
- PubMed
- Google Scholar
1. Lee C-H
2. MacKinnon R
(2018) Activation mechanism of a human SK-calmodulin channel complex elucidated by cryo-EM structures
Science 360:508–513.
https://doi.org/10.1126/science.aas9466
- PubMed
- Google Scholar
(1982) Properties of the Ca2+-activated K+ channel in one-step inside-out vesicles from human red cell membranes
Nature 296:742–744.
https://doi.org/10.1038/296742a0
- PubMed
- Google Scholar
(1968) Demonstration of an internal structure within the red blood cell by ion etching and scanning electronmicroscopy
Nature 220:614–616.
https://doi.org/10.1038/220614a0
- PubMed
- Google Scholar
1. Lewis AH
2. Grandl J
(2021) Piezo1 ion channels inherently function as independent mechanotransducers
eLife 10:e70988.
https://doi.org/10.7554/eLife.70988
- PubMed
- Google Scholar
1. Li A
2. Seipelt H
3. Müller C
4. Shi Y
5. Artmann GM
(1999) Effects of salicylic acid derivatives on red blood cell membranes
Pharmacology & Toxicology 85:206–211.
https://doi.org/10.1111/j.1600-0773.1999.tb02010.x
- PubMed
- Google Scholar
1. Li J
2. Hou B
3. Tumova S
4. Muraki K
5. Bruns A
6. Ludlow MJ
7. Sedo A
8. Hyman AJ
9. McKeown L
10. Young RS
11. Yuldasheva NY
12. Majeed Y
13. Wilson LA
14. Rode B
15. Bailey MA
16. Kim HR
17. Fu Z
18. Carter DA
19. Bilton J
20. Imrie H
21. Ajuh P
22. Dear TN
23. Cubbon RM
24. Kearney MT
25. Prasad RK
26. Evans PC
27. Ainscough JF
28. Beech DJ
(2014) Piezo1 integration of vascular architecture with physiological force
Nature 515:279–282.
https://doi.org/10.1038/nature13701
- PubMed
- Google Scholar
1. Lin Y-C
2. Guo YR
3. Miyagi A
4. Levring J
5. MacKinnon R
6. Scheuring S
(2019) Force-Induced conformational changes in Piezo1
Nature 573:230–234.
https://doi.org/10.1038/s41586-019-1499-2
- PubMed
- Google Scholar
1. Lister JJ
(1827) Xxvi. Notice of some miscroscopic observations of the blood and animal tissues
The Philosophical Magazine 2:130–138.
https://doi.org/10.1080/14786442708674422
- Google Scholar
1. Lu C-H
2. Pedram K
3. Tsai C-T
4. Jones T
5. Li X
6. Nakamoto ML
7. Bertozzi CR
8. Cui B
(2022) Membrane curvature regulates the spatial distribution of bulky glycoproteins
Nature Communications 13:3093.
https://doi.org/10.1038/s41467-022-30610-2
- PubMed
- Google Scholar
1. Lux SE
(2016) Anatomy of the red cell membrane skeleton: unanswered questions
Blood 127:187–199.
https://doi.org/10.1182/blood-2014-12-512772
- PubMed
- Google Scholar
1. Marikovsky Y
2. Danon D
(1969) Electron microscope analysis of young and old red blood cells stained with colloidal iron for surface charge evaluation
The Journal of Cell Biology 43:1–7.
https://doi.org/10.1083/jcb.43.1.1
- PubMed
- Google Scholar
(1990) Mechanosensitive ion channels of E. coli activated by amphipaths
Nature 348:261–263.
https://doi.org/10.1038/348261a0
- PubMed
- Google Scholar
1. Mohandas N
2. Evans E
(1994) Mechanical properties of the red cell membrane in relation to molecular structure and genetic defects
Annual Review of Biophysics and Biomolecular Structure 23:787–818.
https://doi.org/10.1146/annurev.bb.23.060194.004035
- PubMed
- Google Scholar
(2002) Calcium diffusion coefficient in rod photoreceptor outer segments
Biophysical Journal 82:728–739.
https://doi.org/10.1016/S0006-3495(02)75435-0
- PubMed
- Google Scholar
(2001) Imaging erythrocytes under physiological conditions by atomic force microscopy
Biochimica et Biophysica Acta (BBA) - Biomembranes 1514:170–176.
https://doi.org/10.1016/S0005-2736(01)00365-0
- Google Scholar
1. Nunomura W
(2006) Regulation of protein 4.1R interactions with membrane proteins by ca2+ and calmodulin
Frontiers in Bioscience 11:1522.
https://doi.org/10.2741/1901
- Google Scholar
1. Pan L
2. Yan R
3. Li W
4. Xu K
(2018) Super-resolution microscopy reveals the native ultrastructure of the erythrocyte cytoskeleton
Cell Reports 22:1151–1158.
https://doi.org/10.1016/j.celrep.2017.12.107
- Google Scholar
1. Pathak M
2. Nourse JL
3. Tran T
4. Hwe J
5. Arulmoli J
6. Le DTT
7. Bernardis E
8. Flanagan LA
9. Tombola F
(2014) Stretch-activated ion channel Piezo1 directs lineage choice in human neural stem cells
PNAS 111:16148–16153.
https://doi.org/10.1073/pnas.1409802111
- Google Scholar
1. Ranade SS
2. Qiu Z
3. Woo SH
4. Hur SS
5. Murthy SE
6. Cahalan SM
7. Xu J
(2014) Piezo1, a mechanically activated ion channel, is required for vascular development in mice
PNAS 111:10347–10352.
https://doi.org/10.1073/pnas.1409233111
- Google Scholar
1. Ravenhill BJ
2. Kanjee U
3. Ahouidi A
4. Nobre L
5. Williamson J
6. Goldberg JM
7. Antrobus R
8. Dieye T
9. Duraisingh MT
10. Weekes MP
(2019) Quantitative comparative analysis of human erythrocyte surface proteins between individuals from two genetically distinct populations
Communications Biology 2:350.
https://doi.org/10.1038/s42003-019-0596-y
- PubMed
- Google Scholar
1. Ridone P
2. Pandzic E
3. Vassalli M
4. Cox CD
5. Macmillan A
6. Gottlieb PA
7. Martinac B
(2020) Disruption of membrane cholesterol organization impairs the activity of Piezo1 channel clusters
The Journal of General Physiology 152:e201912515.
https://doi.org/10.1085/jgp.201912515
- PubMed
- Google Scholar
(2018) Structure of the mechanically activated ion channel Piezo1
Nature 554:481–486.
https://doi.org/10.1038/nature25453
- PubMed
- Google Scholar
(2019) Guided STED nanoscopy enables super-resolution imaging of blood stage malaria parasites
Scientific Reports 9:4674.
https://doi.org/10.1038/s41598-019-40718-z
- PubMed
- Google Scholar
(2013) Huygens STED deconvolution increases signal-to-noise and image resolution towards 22 nm
Microscopy Today 21:38–44.
https://doi.org/10.1017/S1551929513001089
- Google Scholar
1. Sheetz MP
(1983)
Membrane skeletal dynamics: role in modulation of red cell deformability, mobility of transmembrane proteins, and shape

Seminars in Hematology 20:175–188.
- PubMed
- Google Scholar
(2015) When physics takes over: BAR proteins and membrane curvature
Trends in Cell Biology 25:780–792.
https://doi.org/10.1016/j.tcb.2015.09.005
- PubMed
- Google Scholar
1. Staudt T
2. Engler A
3. Rittweger E
4. Harke B
5. Engelhardt J
6. Hell SW
(2011) Far-Field optical nanoscopy with reduced number of state transition cycles
Optics Express 19:5644–5657.
https://doi.org/10.1364/OE.19.005644
- PubMed
- Google Scholar
1. Sun W
2. Chi S
3. Li Y
4. Ling S
5. Tan Y
6. Xu Y
7. Jiang F
8. Li J
9. Liu C
10. Zhong G
11. Cao D
12. Jin X
13. Zhao D
14. Gao X
15. Liu Z
16. Xiao B
17. Li Y
(2019) The mechanosensitive Piezo1 channel is required for bone formation
eLife 8:e47454.
https://doi.org/10.7554/eLife.47454
- PubMed
- Google Scholar
(2019) A model of Piezo1-based regulation of red blood cell volume
Biophysical Journal 116:151–164.
https://doi.org/10.1016/j.bpj.2018.11.3130
- PubMed
- Google Scholar
(2001) Atomic force microscopy of the erythrocyte membrane skeleton
Journal of Microscopy 204:212–225.
https://doi.org/10.1046/j.1365-2818.2001.00960.x
- PubMed
- Google Scholar
1. Syeda R
2. Florendo MN
3. Cox CD
4. Kefauver JM
5. Santos JS
6. Martinac B
7. Patapoutian A
(2016) Piezo1 channels are inherently mechanosensitive
Cell Reports 17:1739–1746.
https://doi.org/10.1016/j.celrep.2016.10.033
- PubMed
- Google Scholar
Software
1. Tameling C
2. Naas J
(2021) Ctameling/OTC: optimal transport colocalization, version v1.0
Zenodo.

https://doi.org/10.5281/zenodo.4553632
1. Tameling C
2. Stoldt S
3. Stephan T
4. Naas J
5. Jakobs S
6. Munk A
(2021) Colocalization for super-resolution microscopy via optimal transport
Nature Computational Science 1:199–211.
https://doi.org/10.1038/s43588-021-00050-x
- PubMed
- Google Scholar
1. Tarantino N
2. Tinevez J-Y
3. Crowell EF
4. Boisson B
5. Henriques R
6. Mhlanga M
7. Agou F
8. Israël A
9. Laplantine E
(2014) Tnf and IL-1 exhibit distinct ubiquitin requirements for inducing NEMO-IKK supramolecular structures
The Journal of Cell Biology 204:231–245.
https://doi.org/10.1083/jcb.201307172
- PubMed
- Google Scholar
(2017) TrackMate: an open and extensible platform for single-particle tracking
Methods 115:80–90.
https://doi.org/10.1016/j.ymeth.2016.09.016
- PubMed
- Google Scholar
(1998) Regulation mechanism of the lateral diffusion of band 3 in erythrocyte membranes by the membrane skeleton
The Journal of Cell Biology 142:989–1000.
https://doi.org/10.1083/jcb.142.4.989
- PubMed
- Google Scholar
1. Tsai F-C
2. Simunovic M
3. Sorre B
4. Bertin A
5. Manzi J
6. Callan-Jones A
7. Bassereau P
(2021) Comparing physical mechanisms for membrane curvature-driven sorting of BAR-domain proteins
Soft Matter 17:4254–4265.
https://doi.org/10.1039/d0sm01573c
- PubMed
- Google Scholar
1. Tsuji A
2. Kawasaki K
3. Ohnishi S
4. Merkle H
5. Kusumi A
(1988) Regulation of band 3 mobilities in erythrocyte ghost membranes by protein association and cytoskeletal meshwork
Biochemistry 27:7447–7452.
https://doi.org/10.1021/bi00419a041
- PubMed
- Google Scholar
1. Vallese F
2. Kim K
3. Yen LY
4. Johnston JD
5. Noble AJ
6. Calì T
7. Clarke OB
(2022) Architecture of the human erythrocyte ankyrin-1 complex
Nature Structural & Molecular Biology 29:706–718.
https://doi.org/10.1038/s41594-022-00792-w
- PubMed
- Google Scholar
1. Wang J
2. Jiang J
3. Yang X
4. Zhou G
5. Wang L
6. Xiao B
(2022) Tethering piezo channels to the actin cytoskeleton for mechanogating via the cadherin-β-catenin mechanotransduction complex
Cell Reports 38:110342.
https://doi.org/10.1016/j.celrep.2022.110342
- PubMed
- Google Scholar
1. Wu Y
2. Ali MRK
3. Chen K
4. Fang N
5. El-Sayed MA
(2019) Gold nanoparticles in biological optical imaging
Nano Today 24:120–140.
https://doi.org/10.1016/j.nantod.2018.12.006
- Google Scholar
1. Xia X
2. Liu S
3. Zhou ZH
(2022) Structure, dynamics and assembly of the ankyrin complex on human red blood cell membrane
Nature Structural & Molecular Biology 29:698–705.
https://doi.org/10.1038/s41594-022-00779-7
- PubMed
- Google Scholar
Preprint
1. Yang S
2. Miao X
3. Arnold S
4. Li B
5. Ly AT
6. Wang H
7. Wang M
8. Guo X
9. Pathak MM
10. Zhao W
11. Cox CD
12. Shi Z
(2022) Membrane Curvature Governs the Distribution of Piezo1 in Live Cells
bioRxiv.
https://doi.org/10.1101/2022.06.22.497259
- Google Scholar
1. Yu Z-Y
2. Gong H
3. Kesteven S
4. Guo Y
5. Wu J
6. Li JV
7. Cheng D
8. Zhou Z
9. Iismaa SE
10. Kaidonis X
11. Graham RM
12. Cox CD
13. Feneley MP
14. Martinac B
(2022) Piezo1 is the cardiac mechanosensor that initiates the cardiomyocyte hypertrophic response to pressure overload in adult mice
Nature Cardiovascular Research 1:577–591.
https://doi.org/10.1038/s44161-022-00082-0
- Google Scholar
(2012) Mutations in the mechanotransduction protein Piezo1 are associated with hereditary xerocytosis
Blood 120:1908–1915.
https://doi.org/10.1182/blood-2012-04-422253
- PubMed
- Google Scholar
1. Zeng W-Z
2. Marshall KL
3. Min S
4. Daou I
5. Chapleau MW
6. Abboud FM
7. Liberles SD
8. Patapoutian A
(2018) PIEZOs mediate neuronal sensing of blood pressure and the baroreceptor reflex
Science 362:464–467.
https://doi.org/10.1126/science.aau6324
- PubMed
- Google Scholar
1. Zhang W
2. Cheng LE
3. Kittelmann M
4. Li J
5. Petkovic M
6. Cheng T
7. Jin P
8. Guo Z
9. Göpfert MC
10. Jan LY
11. Jan YN
(2015) Ankyrin repeats convey force to gate the NOMPC mechanotransduction channel
Cell 162:1391–1403.
https://doi.org/10.1016/j.cell.2015.08.024
- PubMed
- Google Scholar

Article and author information

Author details

George Vaisey

Laboratory of Molecular Neurobiology and Biophysics, Howard Hughes Medical Institute, The Rockefeller University, New York, United States

Contribution
Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-8359-1314
Priyam Banerjee

Bio-Imaging Resource Center, The Rockefeller University, New York, United States

Contribution
Resources, Methodology, Writing – review and editing, Design of imaging experiments

Competing interests
No competing interests declared
Alison J North

Bio-Imaging Resource Center, The Rockefeller University, New York, United States

Contribution
Resources, Methodology, Writing – review and editing, Design of imaging experiments

Competing interests
No competing interests declared
Christoph A Haselwandter

Department of Physics and Astronomy and Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, United States

Contribution
Formal analysis, Writing – review and editing, Developed calculations of probability distribution along RBC membrane

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5012-5640
Roderick MacKinnon

Laboratory of Molecular Neurobiology and Biophysics, Howard Hughes Medical Institute, The Rockefeller University, New York, United States

Contribution
Conceptualization, Formal analysis, Supervision, Funding acquisition, Writing – original draft, Writing – review and editing, Developed calculations of probability distribution along RBC membrane

For correspondence
mackinn@rockefeller.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7605-4679

Funding

National Center for Research Resources (S10RR031855)

Priyam Banerjee
Alison J North

National Science Foundation (DMR-2051681)

Christoph A Haselwandter

National Science Foundation (DMR-1554716)

Christoph A Haselwandter

National Institutes of Health (GM043949)

Roderick MacKinnon

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

Super-resolution microscopy work was performed in the Bio-Imaging Resource Center at Rockefeller University, RRID:SCR_017791. The OMX 3D-SIM system was funded by Award Number S10RR031855 from the National Center for Research Resources. This work was also supported at USC by NSF Grant No. DMR-2051681 and by NSF Grant No. DMR-1554716 (to CAH) and at Rockefeller University by NIH grant GM43949 (to RM). RM is an investigator of the Howard Hughes Medical Institute.

Ethics

All animal procedures were reviewed and approved by the Institutional Animal Care and Use Committee at The Rockefeller University, protocol # 22056.

Version history

Preprint posted: August 11, 2022 (view preprint)
Received: August 11, 2022
Accepted: December 1, 2022
Version of Record published: December 14, 2022 (version 1)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.