Computational analysis of the tether-pulling experiment to probe plasma membrane-cytoskeleton interaction in cells

Primary cilium: a paradigm for integrating mathematical modeling with experiments and numerical simulations in mechanobiology

Primary cilia are non-motile, solitary (one per cell) microtubule-based organelles that emerge from the mother centriole after cells have exited the mitotic cycle. Identified as a mechanosensing organelle that responds to both mechanical and chemical stimuli, the primary cilium provides a fertile ground for integrative investigations of mathematical modeling, numerical simulations, and experiments. Recent experimental findings revealed considerable complexity to the underlying mechanosensory mechanisms that transmit extracellular stimuli to intracellular signaling many of which include primary cilia. In this invited review, we provide a brief survey of experimental findings on primary cilia and how these results lead to various mathematical models of the mechanics of the primary cilium bent under an external forcing such as a fluid flow or a trap. Mathematical modeling of the primary cilium as a fluid-structure interaction problem highlights the importance of basal anchorage and the anisotropic moduli of the microtubules. As theoretical modeling and numerical simulations progress, along with improved state-of-the-art experiments on primary cilia, we hope that details of ciliary regulated mechano-chemical signaling dynamics in cellular physiology will be understood in the near future.

A Bundle of Mechanisms: Inner-Ear Hair-Cell Mechanotransduction

In the inner ear, the deflection of hair bundles, the sensory organelles of hair cells, activates mechanically-gated channels (MGCs). Hair bundles monitor orientation of the head, its angular and linear acceleration, and detect sound. Force applied to MGCs is shaped by intrinsic hair-bundle properties, by the mechanical load on the bundle, and by the filter imparted by the environment of the hair bundle. Channel gating and adaptation, the ability of the bundle to reset its operating point, contribute to hair-bundle mechanics. Recent data from mammalian hair cells challenge longstanding hypotheses regarding adaptation mechanisms and hair-bundle coherence. Variations between hair bundles from different organs in hair-bundle mechanics, mechanical load, channel gating, and adaptation may allow a hair bundle to selectively respond to specific sensory stimuli.

Excess area dependent scaling behavior of nano-sized membrane tethers

Thermal fluctuations in cell membranes manifest as an excess area ([Formula: see text]) which governs a multitude of physical process at the sub-micron scale. We present a theoretical framework, based on an in silico tether pulling method, which may be used to reliably estimate [Formula: see text] in live cells. We perform our simulations in two different thermodynamic ensembles: (i) the constant projected area and (ii) the constant frame tension ensembles and show the equivalence of our results in the two. The tether forces estimated from our simulations compare well with our experimental measurements for tethers extracted from ruptured GUVs and HeLa cells. We demonstrate the significance and validity of our method by showing that all our calculations performed in the initial tether formation regime (i.e. when the length of the tether is comparable to its radius) along with experiments of tether extraction in 15 different cell types collapse onto two unified scaling relationships mapping tether force, tether radius, bending stiffness κ, and membrane tension σ. We show that [Formula: see text] is an important determinant of the radius of the extracted tether, which is equal to the characteristic length [Formula: see text] for [Formula: see text], and is equal to [Formula: see text] for [Formula: see text]. We also find that the estimated excess area follows a linear scaling behavior that only depends on the true value of [Formula: see text] for the membrane, based on which we propose a self-consistent technique to estimate the range of excess membrane areas in a cell.

Model for probing membrane-cortex adhesion by micropipette aspiration and fluctuation spectroscopy

We propose a model for membrane-cortex adhesion that couples membrane deformations, hydrodynamics, and kinetics of membrane-cortex ligands. In its simplest form, the model gives explicit predictions for the critical pressure for membrane detachment and for the value of adhesion energy. We show that these quantities exhibit a significant dependence on the active acto-myosin stresses. The model provides a simple framework to access quantitative information on cortical activity by means of micropipette experiments. We also extend the model to incorporate fluctuations and show that detailed information on the stability of membrane-cortex coupling can be obtained by a combination of micropipette aspiration and fluctuation spectroscopy measurements.

The local forces acting on the mechanotransduction channel in hair cell stereocilia

In hair cells, mechanotransduction channels are located in the membrane of stereocilia tips, where the base of the tip link is attached. The tip-link force determines the system of other forces in the immediate channel environment, which change the channel open probability. This system of forces includes components that are out of plane and in plane relative to the membrane; the magnitude and direction of these components depend on the channel environment and arrangement. Using a computational model, we obtained the major forces involved as functions of the force applied via the tip link at the center of the membrane. We simulated factors related to channels and the membrane, including finite-sized channels located centrally or acentrally, stiffness of the hypothesized channel-cytoskeleton tether, and bending modulus of the membrane. Membrane forces are perpendicular to the directions of the principal curvatures of the deformed membrane. Our approach allows for a fine vectorial picture of the local forces gating the channel; membrane forces change with the membrane curvature and are themselves sufficient to affect the open probability of the channel.

Dynamics of the primary cilium in shear flow

In this work, the equilibrium shape and dynamics of a primary cilium under flow are investigated by using both theoretical modeling and experiment. The cilium is modeled as an elastic beam that may undergo large deflection due to the hydrodynamic load. Equilibrium results show that the anchoring effects of the basal body on the cilium axoneme behave as a nonlinear rotational spring. Details of the rotational spring are elucidated by coupling the elastic beam with an elastic shell. We further study the dynamics of cilium under shear flow with the cilium base angle determined from the nonlinear rotational spring, and obtain good agreement in cilium bending and relaxing dynamics when comparing between modeling and experimental results. These results potentially shed light on the physics underlying the mechanosensitive ion channel transport through the ciliary membrane.

Stereocilia membrane deformation: implications for the gating spring and mechanotransduction channel

In hair cells, although mechanotransduction channels have been localized to tips of shorter stereocilia of the mechanically sensitive hair bundle, little is known about how force is transmitted to the channel. Here, we use a biophysical model of the membrane-channel complex to analyze the nature of the gating spring compliance and channel arrangement. We use a triangulated surface model and Monte Carlo simulation to compute the deformation of the membrane under the action of tip link force. We show that depending on the gating spring stiffness, the compliant component of the gating spring arises from either the membrane alone or a combination of the membrane and a tether that connects the channel to the actin cytoskeleton. If a bundle is characterized by relatively soft gating springs, such as those of the bullfrog sacculus, the need for membrane reinforcement by channel tethering then depends on membrane parameters. With stiffer gating springs, such as those from rat outer hair cells, the channel must be tethered for all biophysically realistic parameters of the membrane. We compute the membrane forces (resultants), which depend on membrane tension, bending modulus, and curvature, and show that they can determine the fate of the channel.

Internal forces, tension and energy density in tethered cellular membranes

We analyze tethered cellular membranes by considering the membrane resultants, tension and densities of two modes of energy, bending and adhesion. These characteristics are determined based on a computational (finite-difference) analysis of membrane shape. We analyze the relative contribution and distribution of the membrane characteristics in four typical zones of the membrane surface. Using an axisymmetric model, we found that the meridional and circumferential components of the resultant are different near the tether body and they converge to the value of membrane tension farther from the tether. At the beginning of the area of membrane detachment from the cytoskeleton, the density of bending energy is on the same order of magnitude as membrane tension (resultant). Away from the tether, the bending energy density quickly decreases and becomes of the same order as that of the adhesion energy in the membrane-cytoskeleton attachment area. In that area, both modes of energy are significantly smaller than the membrane tension. We also consider the effect of the membrane bending modulus on the distribution of the membrane characteristics. An increase in the bending modulus results in changing the length and position on the membrane surface of zone 1 characterized by significant evolution of the resultant components. It also results in shortening zone 2 that covers the rest of the area of membrane detachment. The obtained results can help in a better interpretation of the measurements of membrane mechanical properties as well as in analyses of proteins and channels in curved membranes.

Constitutively active ezrin increases membrane tension, slows migration, and impedes endothelial transmigration of lymphocytes in vivo in mice

ERM (ezrin, radixin moesin) proteins in lymphocytes link cortical actin to plasma membrane, which is regulated in part by ERM protein phosphorylation. To assess whether phosphorylation of ERM proteins regulates lymphocyte migration and membrane tension, we generated transgenic mice whose T-lymphocytes express low levels of ezrin phosphomimetic protein (T567E). In these mice, T-cell number in lymph nodes was reduced by 27%. Lymphocyte migration rate in vitro and in vivo in lymph nodes decreased by 18% to 47%. Lymphocyte membrane tension increased by 71%. Investigations of other possible underlying mechanisms revealed impaired chemokine-induced shape change/lamellipod extension and increased integrin-mediated adhesion. Notably, lymphocyte homing to lymph nodes was decreased by 30%. Unlike most described homing defects, there was not impaired rolling or sticking to lymph node vascular endothelium but rather decreased migration across that endothelium. Moreover, decreased numbers of transgenic T cells in efferent lymph suggested defective egress. These studies confirm the critical role of ERM dephosphorylation in regulating lymphocyte migration and transmigration. Of particular note, they identify phospho-ERM as the first described regulator of lymphocyte membrane tension, whose increase probably contributes to the multiple defects observed in the ezrin T567E transgenic mice.

Cell cytoskeleton and tether extraction

We perform a detailed investigation of the force × deformation curve in tether extraction from 3T3 cells by optical tweezers. Contrary to conventional wisdom about tethers extracted from cells, we find that actin filaments are present within them, so that a revised theory of tether pulling from cells is called for. We also measure steady and maximum tether force values significantly higher than previously published ones for 3T3 cells. Possible explanations for these differences are investigated. Further experimental support of the theory of force barriers for membrane tube extension is obtained. The potential of studies on tether pulling force × deformation for retrieving information on membrane-cytoskeleton interaction is emphasized.

Models of dynamic extraction of lipid tethers from cell membranes

When a ligand that is bound to an integral membrane receptor is pulled, the membrane and the underlying cytoskeleton can deform before either the membrane delaminates from the cytoskeleton or the ligand detaches from the receptor. If the membrane delaminates from the cytoskeleton, it may be further extruded and form a membrane tether. We develop a phenomenological model for this process by assuming that deformations obey Hooke's law up to a critical force at which the cell membrane locally detaches from the cytoskeleton and a membrane tether forms. We compute the probability of tether formation and show that tethers can be extruded only within an intermediate range of force loading rates and pulling velocities. The mean tether length that arises at the moment of ligand detachment is computed as are the force loading rates and pulling velocities that yield the longest tethers.