Probing nanomechanical responses of cell membranes

Despite the importance in various cellular processes, the nanomechanical responses of the living cell membrane have been elusive due to complexities in the membrane associated with the hidden architecture of multiple molecular components, including the lipid bilayer. Here, combined experimental and theoretical frameworks that can probe and interpret nanomechanical responses of the cell membrane are demonstrated. A magnetic tweezer assay was introduced to apply pico-Newton scale forces to lipids and E-cadherin molecules at the living cell surface. Two unique classes of force-extension curves were identified: one with a deflection transition (Type I) and another with a discontinuous transition (Type II). The repeated observations of these responses, regardless of cell type and targeted cell surface molecule, suggest the Type I and II curves are the primary nanomechanical responses of cell membranes. To reproduce these responses in vitro, a model system using synthetic lipid vesicles was also developed. Together with a finite element model of lipid bilayers, the reproduced responses suggest that the confined fluidity and curvature constraints imposed on the lipid bilayer components of the cell membrane are the main parameters responsible for the generation of these responses. This work provides an insight into how forces on membrane molecules propagate to the lipid bilayer components to generate specific nanomechanical responses. In addition, the consistent results obtained using different methodologies demonstrate that the presented force-probing assays and the theoretical model can serve a combined testbed to investigate nanoscale mechanics of the living cell membrane.

www.nature.com/scientificreports www.nature.com/scientificreports/ Before making force response measurements for cells incubated with ~65-260 nM of the biotinylated lipids, both the natural fluctuation and induced rotation of the magnetic beads were checked. Beads with minimal fluctuation and wider or non-circular rotations with severe wobbling, which may indicate poor surface targeting, were not sampled during magnetic tweezer experiments (Fig. 1c). Note that sliding of beads on the cell surface was occasionally observed during rotation. Finally, the force-extension curves were recorded by pulling the biotin-conjugated lipids with a constant loading rate and tracking the vertical position of the magnetic beads.
A lot of complexity were observed in the measured responses, including both vertical motion and lateral shifting of the beads even without significant force application. In repeated experiments, however, two types of nanoscale responses with unique mechanistic features were identified. These responses were distinct from the classic tubular structures formed at the cell membrane. Figure 1d shows one representative force-extension curve (referred to as Type I) in the second and fourth loading cycles. This type of force-extension curve starts with a compliant regime and exhibits a biphasic character organized around a deflection or kink in the curve. Figure 1e shows another type of force-extension curve (hereafter referred to as Type II). Without the noticeable deflection, the second, fourth, and fifth cycles of the responses start with an initial stiff barrier and feature a discontinuous transition (i.e. unstable stretching of the membrane) to another stiff regime. In many cases, a small sign for the discontinuous re-transition was also observed during the relaxation. Finally, the microscopic movement of the magnetic beads was observed when the cells were treated with Latrunculin-A. The force-extension response demonstrated the micro-scale extraction of the cell membrane ( Fig. 1f-h). Since Latrunculin-A is known to inhibit actin polymerization 13 , the result suggests an intact cytoskeletal network is important for the generation of Type I and II nanomechanical responses.
Mathematical model of nanomechanical lipid bilayers: comparison between measurements and calculations. Next, whether theories of lipid bilayer mechanics can explain the Type I and II measurements was investigated. To this end, a simple mathematical model of the nanomechanical lipid bilayer that vs. extension curves obtained from successive loading cycles. r ct = 500 nm, r cr = 650 nm (purple); r ct = 500 nm, r cr = 850 nm (gray); r ct = 500 nm, r cr = 1000 nm (dark blue); and r ct = 500 nm, r cr = 1100 nm (blue) were used for calculations in (d). r ct = 160 nm, r cr = 1000 nm (dark blue); r ct = 160 nm, r cr = 1500 nm (blue); r ct = 160 nm, r cr = 1170 nm (purple) were used in e. The calculated membrane shapes (at red marks) were plotted at the top of (d,e). (f) Schematic for the magnetic tweezer experiments with Latrunculin-A treatments. (g) Representative force vs. extension curve with Latrunculin-A. (h) The fraction of magnetic beads that generated membrane tubules (see black arrows in g) with 20 pN in ~15 seconds. (i) Parameterized r ct vs. r cr scatter from multiple bead measurements of three different cell lines (N U2OS = 16, N HUVEC = 10, N IMR90 = 19). Measurements with the discontinous transition (Type II) are shown in red and measurements with the kink (Type I) are shown in blue. The decision boundaries (black line) are: −1.0531 + 0.0351r ct − 0.0105r cr = 0; −4.0420 + 0.038r ct − 0.0091r cr = 0; and −1.3143 + 0.0245r ct − 0.0063r cr = 0 for U2OS, HUVEC, and IMR90, respectively. www.nature.com/scientificreports www.nature.com/scientificreports/ combines a well-known curvature elasticity theory with an expression for membrane surface tension was introduced [14][15][16][17] . The total energy functional in domain Ω Ψ nano membrane was formulated as in Eq. (1).
The first term integrates the energy densities due to the mean (H) and Gaussian curvatures (K) over the membrane area (A). Here k m and k g denote the bending and Gaussian curvature moduli of the membrane, respectively. The second term incorporates the free energy change due to area strain of the membrane. The surface tension (σ) is the function of the area strain of the membrane α = = where A res is the area of the lipid membrane reservoir determined by r cr − r ct , Ω, and r cb (see the next paragraph and Fig. S3 in Supplementary Information online); A res 0 is the area of the lipid reservoir at the resting reference configuration; φ res is the uniform lipid number density; and res 0 φ is res φ at the resting reference configuration. The function integration of σ with respect to α gives the surface strain energy density. Finite element modeling for this theory was performed, whose full description is shown in the Methods (also see an approach in ref. 18 ). The estimates of shape and applied force resulting from the application of a specific displacement of the lipid bilayer can be calculated using this model.
Note the presence of geometric parameters in this model (Fig. 1a): the radius of the lipid reservoir that defines the limit of lateral membrane stretching (r cr ) and the radius of the area where the membrane is tightly associated with the rigid cytoskeleton (r cc ). From r cr and r cc , the radius of the membrane tented (r ct = r cr − r cc ) can be defined. In the rotational axisymmetric configuration, these simple parameters provide a conceptual idea on the physical boundaries responsible for the generation of the nanomechanical responses. Remarkably, with an assumption of point-like pulling, these model-based membrane predictions showed a reasonable fit to the Type I and II observational data. For the Type I data presented in the second and fourth cycles of Fig. 1d, r ct = 500 nm and r cr = 1000 nm were used to fit the calculation to the measurements. Without introducing additional parameter values, a reasonable fit for the Type II data was also achieved by modulating the r ct and generating a sigmoidal-type curve. With a fixed r cr , reducing r ct resulted in the generation of the initial barrier for the force vs. extension calculation. For the fourth and fifth cycle of Fig. 1e, r ct = 160 nm and r cr = 1000 nm were used.
In many cases, using more than one r cr values provided a better interpretation for the single force vs. extension measurement. An additional r cr value (i.e. another limit of the lipid reservoir) was introduced to explain a secondary transition (see blue arrows) followed by the deflection in the first and third cycles of Fig. 1d, and the discontinuous transition in the first, third, and sixth cycles of Fig. 1e (also see other examples in supplementary Fig. S5e). Four r cr values with a shared r ct of 500 nm were used for all four successive cycles of measurements in Fig. 1d. Similarly, three r cr values with a single r ct of 160 nm were used for all six successive cycles of measurements in Fig. 1e.
While measurements were sorted into Type I or Type II responses, a more detailed examination of their force-extension shapes did show some variation. By using the lipid bilayer model, a diagram of Type I and II measurements was generated for Fig. 1i. The data in this diagram were smoothly separated in the r ct -r cr plane, depending on the type of measured force-extension responses. These results suggest that the simple model reasonably explains not only the individual force vs. extension trace but also the whole data set on a higher level. There are also limitations in comparison between the model and the observation. According to the model, for example, beads attached with non-negligible membrane area can also generate Type I and II nanomechanical responses, providing a possibility of an underestimation for the r ct and r cr values (see Supplementary Fig. S3 for the discussion). In addition, same material parameters for the membrane were used for all calculations while they may vary for different region of the real cell surface (see Supplementary Table S1). Nevertheless, predictions made by using theories of lipid bilayers seem to indicate that the bilayer's elastic properties determined by r ct and r cr are important factors for the generation of Type I and II responses. Of note, r ct values in the diagram were similar to the size of membrane compartments determined from imaging molecular hop diffusion 19,20 . Studying whether these share a common biophysical underpinning might be informative.
Pulling lipids at the surface of different cell types. Lipid bilayers are a common factor of the membrane, regardless of the type of the cell. To test observed Type I and II responses in other cell types beyond U2OS cells, the same pulling protocol was applied to living human umbilical vein endothelial cells (HUVEC) and human fetal lung fibroblasts (IMR-90). Because these cell types are distinct from U2OS cells, they may have different underlying cellular substructures. Nevertheless, both Type I and II responses were also observed when pulling the lipids of these cell types. In addition, the model similarly explained these measurements (Supplementary Figs. S4 and S5). While the Type I and II data fell in slightly different regions of the r ct -r cr diagram, the decision boundaries determined for each respective cell line were largely overlapping (Fig. 1i). The result suggests a similar principle of membrane nanomechanics for different cell types, even if their predominant targeting response might be different.
In vitro model of synthetic nanovesicles: repnroduction of Type I and Type II responses. The observed Type I and II nanomechanical responses, likely generated from the lipid bilayer components of living cell membranes, were reproduced in vitro. To this end, synthetic lipid vesicles with a mean diameter of 118 ± 14 (S.D.) nm were made via an extrusion method (Fig. 2d, see Supplementary Information online). These nanoscale vesicles have a limited size of lipid reservoir. The vesicles were immobilized on a glass substrate using biotin-avidin bonding. Magnetic beads were then attached to single vesicles via minimally treated dinitrophenyl-conjugated lipids as demonstrated in Fig. 2a. Before beginning the pulling protocol, rotatory motion of the beads was checked www.nature.com/scientificreports www.nature.com/scientificreports/ as similarly done for the cell membrane. With these nanovesicles, remarkably, the Type I force-extension response was reproduced as shown in Fig. 2e (also see Supplementary Fig. S6a-c).
Analysis performed for the cell membrane measurements suggested increased curvature constraints (i.e., increased r cc or decreased r ct ) can generate the Type II response. To validate this idea, greater levels of rigid body interaction was imposed on the immobilized vesicles by simply increasing the biotin-avidin interactions between the vesicles and the underlying surface. A higher concentration of avidins increased the area of attachment, which led to a flattening of the immobilized vesicles ( Fig. 2b-d, see Supplementary Information online). Finally, when these vesicles of higher avidin condition were pulled, the Type II response was more frequently observed than the Type I (Figs. 2f and 3a, also see Supplementary Fig. S6d,e).
The reproduced vesicle responses were also interpreted with the finite element model (see Methods). The energy functional is expressed in Eq. (2). www.nature.com/scientificreports www.nature.com/scientificreports/ The last term with the Lagrange multiplier (λ) reflects the condition of the fixed volume (V 0 ) spanned by the vesicle. Here, the two geometric parameters introduced for the cell membrane, r cr and r cc , were replaced with r vr (the vesicle radius in a deformation-free configuration) and r vc (the radius of the area of contact with the underlying surface). Figure 2e,f show that calculations with reasonably defined r vr and r vc values explain both Type I and II responses of vesicles (also see Supplementary Fig. S6). The force decomposition analysis shown in Supplementary Fig. S7 provides an interpretation of these biphasic responses 21 . According to the analysis, the initial force-extension curve is dominated by the force due to the membrane curvature development. The deflection (Fig. S7a) and the re-rise of the curve after the turning point ( Fig. S7b) coincide with the point where the force due to membrane stretching becomes significant. The analysis suggests that additional membrane curvature development, to prevent the membrane strain energy from becoming excessive, is responsible for the responses in the larger extension regime.
With a particular interest in the bi-stable system, how the energy landscape of Type II responses interprets the observed discontinuous transitions was investigated in Fig. 3. For a vesicle where Type II force-extension curves were measured from sufficiently repeated loading cycles, the shape was calculated and divided into three partitions (Fig. 3a,b, see Supplementary Information online). Then, the energy of these three partitions (i.e., tip, middle, and bottom partitions) was calculated by taking the vesicle height as a generalized coordinate (Fig. 3c, see Supplementary Fig. S8 and Supplementary Information online). Even with nearly zero tension, the energy landscape of the middle partition showed two minima separated by an energy barrier (Fig. 3c blue arrows). As the tension was increased, the local minimum with a larger extension was shifted downward, becoming a global minimum at around 9 pN (Fig. 3c). Calculations suggest that the transition between these two energy minima appears as the step extension in the Type II force vs. extension curves.
The energy landscape was also evaluated by using a theory widely used in single-molecule force spectroscopy 22 . The membrane calculation suggests that the change of the generalized force to the middle partition is nearly the same with the tip force within the force range examined. Accordingly, transition and re-transition force histograms for the middle partition were generated from magnetic force measurements (Fig. 3a,d,e), where the energy barrier information was estimated ( Fig. 3d-g, see Supplementary Information online). The results of this analysis were consistent with the estimates of the membrane model ( Fig. 3c red marks). Together, the highly controlled measurements and subsequent multidisciplinary analyses performed with synthetic lipid vesicles reveal the biophysical principles that produce the novel mechanical responses of the cell membrane, which include bistability. www.nature.com/scientificreports www.nature.com/scientificreports/ Pulling E-cadherins at the living cell surface. Force is transmitted to the cell surface, in many cases, through membrane receptor proteins. To gain collective insight on how these receptor forces are conveyed to the lipid bilayer, the experimental and theoretical methodologies were also applied to pulling E-cadherin proteins expressed endogenously on the surface of living U2OS cells. Here, the membrane proximal ectodomains of the cadherins were targeted with biotinylated E-cadherin antibodies (Fig. 4a, Supplementary Figs. S2b and S9) 23,24 , so any possible contributions from the extracellular cadherin repeats would be minimized 25 . Notably, force-extension curves with the deflection and the discontinuous transition were also observed in repeated experiments pulling E-cadherins (Fig. 4b,c). Furthermore, when these responses were examined by using the lipid bilayer model, the individual force-extension curves could be fitted with similar r ct and r cr values as for the lipid pulling experiments. Finally, the scatter of these values in the r ct -r cr diagram largely overlaps with the data obtained from the lipid pulling experiments (Fig. 4d). These results suggest the remarkable possibility that lipid bilayer responses can be triggered by the pulling of the E-cadherin. Although it is unclear whether and how the cadherin and any associated molecules contribute to the observed responses, the results demonstrate that mechanical pulling on a membrane-bound protein could produce direct displacements of the lipid bilayer.
To further support the idea of bilayer responses generated via force on a membrane-bound protein, a buffer of cholesterol saturated in 20 mM methyl-β-cyclodextrin (MβCD) was treated while pulling E-cadherin (Fig. 4e). First, the magnetic beads that showed the nanomechanical responses were selected and treated with pure buffer injections. Next, the beads were asked whether similar force-extension responses can be produced even after the buffer flow. In many cases, the force-extension shapes were altered after the flow, likely due to structural changes of the perturbed cell membrane. There were cases, however, without critical modifications in the force vs. extension responses. When these beads were further treated with the cholesterol-MβCD buffer, changes in curves were observed within several cycles of pulling and relaxation of the cadherin. These changes were explained by a decrease of r ct or r cr within the presented model framework (Fig. 4e,f). Cholesterol-MβCD complexes are known to enrich the membrane cholesterol level 26,27 , suggesting that cholesterol is inducing this reduction. Cholesterols are known to modulate rigidity and fluidity of the lipid bilayer 28 . The idea that the changes in the nanomechanical responses are due to changes in the lipid bilayer itself is well-supported by the functional and biochemical properties of membrane cholesterols. www.nature.com/scientificreports www.nature.com/scientificreports/

Discussion
Numerous cellular events begin with force applied across the cell membrane. In addition to asking for how applied forces are transduced into a certain cellular phenotype 2,5,29 , an important subject in recent studies of mechanobiological processes is to investigate how the mechanical responses propagate within the force-bearing elements 30 . Detailed mechanical responses had been identified for various systems at the single-molecular level in vitro 25,[31][32][33][34][35] . However, the nanomechanical responses of cell membranes that can serve an important role in the propagation of mechanical energy across the cell surface had been largely elusive.
Here, by targeting cell surface molecule lipids and E-cadherins, two primary nanomechanical responses of living cell membranes were identified, one with a deflection (Type I) and the other with a discontinuous transition (Type II) of the force-extension curve. These responses were universal regardless of tested cell type and targeted molecule in the cell membrane. In vitro vesicle and theoretical analyses performed for the observed responses suggested that, while force receptors can vary on the cell surface, the shape of the mechanical responses can be determined by the lipid bilayer components in which the molecule is inserted (Fig. 4g). The modulation of the fluidity and flexibility of the lipid bilayer are important determinants of the properties of these responses. Questions remain to be answered for how cells take advantage or disadvantage of these unique mechanical characteristics of the lipid bilayer for specific physiological functions. Overall, the work provides standardized living, synthetic, and theoretical frameworks to study nanoscale mechanics of the cell membrane. Furthermore, the work demonstrates how mechanical cues can propagate from the receptor proteins to the lipid bilayer components to contribute to nanomechanical responses of the cell membrane, and thus provides a coupled force transmission paradigm for various processes at the cell surface.

Methods
Finite element modeling for lipid membranes. Weak form of the problem. Finite element methods for the membrane theory are provided here in detail. Since Eqs. (1) and (2) share the identical formalisms without the third term of Eq. (2), descriptions are adapted for Eq. (2). As shown in the results section, the total energy functional in (2) is expressed with three terms: (1) the Canham-Helfrich curvature energy term, (2) a simple area strain energy term, and (3) a term for the fixed volume V 0 with a Lagrange multiplier λ. Note that λ = 0 can be defined for Eq. (1).
The differential forms of the surface are described using the Monge gauge with respect to the parametric coordinate of the model. The mean curvature H and Gaussian curvature K for the rotational axisymmetric membrane geometry can be expressed as follows 36,37 = . and π = r dV 2 (−h s )ds are the vesicle axisymmetric area and volumetric elements, respectively. The constants k m and k g in Eq. (2) are the bending and Gaussian curvature modulus of the membrane, respectively. Here k g = 0.5 k m is used. An expression for the membrane surface tension σ was determined from Eq. (4) 16,17  where A res is the area of the vesicle membrane determined by r vc , Ω, and r cb (see Fig. S3 in Supplementary Information online); A res 0 is the resting reference area of the vesicle in the spherical configuration; φ res is the uniform lipid number density; and res 0 φ is res φ at the resting reference configuration. Therefore, the function integration of surface tension with respect to the area strain α from initial zero strain α 0 to the strain under consideration α c estimates surface strain energy density as denoted in the second term of Eq. (2). See Supplementary Table S1 for the membrane parameters used in this study.
In the variational context, finite element solutions can be found by taking the first variation of Eq. (2) equal to zero-a condition necessary to minimize free energy. To this end, the variational equation is defined in (5) as follows: . Here ψ H is integrand associated with the curvature energy and ψ λ is integrand for the term of volumetric constraints when Eq. (2) is expanded with respect to parametric domain s. The expression T α is defined as derived in the Appendix A in Supplementary Information online. Here, the Gaussian curvature energy contribution can be omitted because surface integration of the Gaussian curvature is invariant for the axisymmetric deformation. Therefore, their variations are identically zero, according to the Gauss-Bonnet theorem 38 . The functions δh and δs denote variations of the membrane shape functions h and r, respectively.
Since Eq. (5) is expressed by the maximum second order of the membrane shape functions, the collection of trial solutions S and admissible variations V with boundary conditions up to its first derivative are defined as follows: s s s s here, h 0 and r 0 are the prescribed displacements at s = 0, as well as h s,0 and r s,0 are the prescribed slope at s = 0. Similarly, h 1 , r 1 , and h s,1 , r s,1 are the prescribed displacements and slopes at s = 1. Finally, the weak form of the presented boundary value problem can be stated, "Find (h, r) ∈ S and λ ∈  satisfying Eq. (5) for ∀(δh, δr) ∈ V and ∀δλ ∈ ". Galerkin form. The Galerkin form of the problem can be stated by approximating function spaces S (6) and V (7) in the weak statement to the finite dimensional space S h and V h . Here, the parameterized functions for membrane shape (h h , r h ) and their variation (δh h , δr h ) belong to S h and V h , respectively, i.e., (h h , r h ) ∈ S h ⊂ S and (δh h , δr h ) ∈ V h ⊂ V. Therefore, the Galerkin approximated solution of Eq. (5) can be found by solving (8) as follows: B-spline-based approximation. With the given Galerkin form of the problem, the structure of the functions δh h , δr h , u h , and v h in space V h , as well as the structure of the given functions w h and x h that satisfy the boundary conditions need to be defined. Since δh h , δr h , u h , and v h belong to H 2 functions, C 1 conforming elements are required. The model parameterizes axisymmetric membrane shape with parametric B-spline functions. Although quadratic functions are used here, the numerical framework introduced in this work can be extended for any type of spline function family.
Instead of parameterizing δh h and δr h (and u h and v h ) independently to provide full two-dimensional degrees of freedom (DOFs) for the finite element nodes, the motion of each node was constrained into the one-dimensional normal direction of a given reference curve. As previously discussed, such treatment was effective for avoiding the so-called zero-energy mode in numerical methods for Helfrich-type membrane models 18 . Since the model solves the shape of the membrane by applying infinitesimal displacement steps during nonlinear calculations, here, the reference curve was simply defined from the membrane shape calculated in the previous displacement step.
Based on the approach summarized above, the parameterized variation of the shape δh h and δr h can be expressed as follows by setting that a spline basis function N(s) i is associated with the i th node in the discretized parametric coordinate s. To model the parameterized shape function h h , w h is defined as in Eq. (11) using the same spline basis function N. In Eqs. (11) and (12), N −1 , N 0 , N n+1 , and N n+2 are B-spline basis functions for the boundary region, and dh −1 , dh 0 , dh n+1 , and dh n+2 are the given fixed values that define the boundary conditions. d 1  The first and second derivatives with respect to s can be derived from the given structure of δh h , h h , δr h , and r h . Substituting those into Eq. (8) as well as the arbitrariness of c i in Eqs. (9) and (10) and δλ results in coupled n + 1 nonlinear simultaneous equations, as denoted by the vector notation in (15).
a n 1 Here, the residual vector G has a total of n + 1 rows with its a th row component G a (for 1 ≤ a ≤ n) being defined in Eq. (16). There, each equation contains five nodal unknowns in using the quadratic B-spline basis function.
Here, A, B, C, D and E represent parameterized A h , B h , C h , D h , and E h in Eq. (8), respectively. The n + 1 th element of the residual vector that is associated with a Lagrange multiplier can be defined as follows: Here, similarly, F represents the parameterized F h in Eq. (8).
Linearization: jacobian matrix for the newton-raphson method. Given n + 1 nonlinear equations in (15), a tangential operator (i.e., n + 1 by n + 1 Jacobian matrix) can be derived to use Newton's method to iteratively obtain solutions of the nonlinear equation system (i.e., d 1 , d 2 , …, d n and λ). For this purpose, the Jacobian matrix J of the residual vector G is defined as follows: Here, the elements in the a th -row and the b th -column of the matrix (for 1 ≦ a ≦ n, 1 ≦ b ≦ n) can be found from Eq. Finally, by substituting an initial guess for the solutions d 0 (i.e. d 1 , d 2 , …, d n and λ) to the Jacobian matrix J and the residual vector G, the solution vector d for j + 1 th Newton's iteration can be calculated from Here fixed values for the d were substituted into the terms associated with the area strain T α to avoid the potent numerical oscillation and thus the divergence of the iterative process. The strain www.nature.com/scientificreports www.nature.com/scientificreports/ value for k th displacement step was predicted from k − 1 th step and corrected for k + 1 th step. Such treatments might be supported by assuming the infinitesimal lateral oscillation of the membrane area from the outside of the parametric domain for different displacement steps. To account some large deformation of the cell membrane, the infinitesimal shifting of the r cr was also allowed to satisfy a condition that the tented membrane area of the k th stretching step is greater than or equal with that of the k − 1 th step and by assuming the curvature deformation limit. Gaussian quadrature and standard mapping techniques were used to compute element arrays for the Jacobian matrix J and the residual vector G. With given boundary values, this iterative process is continued until the Euclidean norm of the difference of two subsequent solution vectors converges to a certain tolerance. See Supplementary Information and Algorithm S1 online for simulation details.

experimental Methods
The mechanical pulling experiments were carried out using a previously developed magnetic tweezer apparatus 9,10 . After introducing polystyrene reference beads to the bottom surface of the channel slides purchased from ibidi, cells were seeded and cultured on the slides until they formed an epithelial monolayer. Then, streptavidin-coated magnetic beads (Dynabeads ® M-280 streptavidin) were introduced to the upper surface of the cellular layer after treatment and washing of ~65-260 nM biotinylated lipids (Fig. 1), and ~0.5 μg/ml of biotinylated E-Cadherin monoclonal antibody (Fig. 4). The force vs. extension responses were measured with a constant loading rate in 2.1-2.6 pN/s. For the cholesterol flow experiments, a sufficient amount of cholesterol powders was mixed with 20 mM methyl-β-cyclodextrin. After overnight incubation at 37 °C with vortexing, undissolved cholesterol remains were filtered from the buffer. This cholesterol solution was applied with a flow velocity of 3.3 ul/sec during the buffer exchange process.
For the vesicle experiments, a channel of ~20 μl volume where the bottom coverslip was coated with polyethylene glycol (PEG) polymer chains was constructed. Two types of PEG-coated coverslips were prepared with two different ratios used in mixing the PEG polymer chains and the biotin-conjugated PEG polymer chains. Vesicles were made by following a typical extrusion protocol 39,40 using polycarbonate membrane filters with 100 nm pores.
The anti-dinitrophenyl antibody was conjugated to the magnetic beads (Dynabeads ® M-270 Carboxylic Acid) to target dinitrophenyl-conjugated lipids in the vesicle. Atomic force microscopic images for the vesicle were acquired with NanoWizard Ultra Speed (JPK Instruments) in Quantitative Imaging (QI) mode. Triangular Si 3 N 4 cantilevers with the spring constant of 0.02 N/m were used. Both vesicle and live-cell experiments were performed at room temperature (22-25 °C). See Supplementary Information online for experimental method details.