Is the surface area of the red cell membrane skeleton locally conserved?

The incompressibility of the lipid bilayer keeps the total surface area of the red cell membrane constant. Local conservation of membrane surface area requires that each surface element of the membrane skeleton keeps its area when its aspect ratio is changed. A change in area would require a flow of lipids past the intrinsic proteins to which the skeleton is anchored. in fast red cell deformations, there is no time for such a flow. Consequently, the bilayer provides for local area conservation. In quasistatic deformations, the extent of local change in surface area is the smaller the larger the isotropic modulus of the skeleton in relation to the shear modulus. Estimates indicate: (a) the velocity of relative flow between lipid and intrinsic proteins is proportional to the gradient in normal tension within the skeleton and inversely proportional to the viscosity of the bilayer; (b) lateral diffusion of lipids is much slower than this flow; (c) membrane tanktreading at frequencies prevailing in vivo as well as the release of a membrane tongue from a micropipette are fast deformations; and (d) the slow phase in micropipette aspiration may be dominated by a local change in skeleton surface.     

Time-dependent elastic extensional RBC deformation by micropipette aspiration: redistribution of the spectrin network?

The time dependence of small elastic extensional RBC deformation by micropipette aspiration has been analyzed. This process shows two-phases which are characterized by time constants of the order of some tenths of seconds and about ten seconds, respectively. The equilibrium tongue length is reached after about 30 s. For the first, fast step we assume that the membrane model of immobilized boundaries holds, i.e., the skeleton is tightly associated with the lipid bilayer and no redistribution of the skeleton with respect to the lipid bilayer is allowed. This lipid-spectrin interaction or anchorage is characterized by some association force density. It has been shown that at a given tongue length the force generated owing to the membrane deformation and acting to redistribute the spectrin, overcomes (in some membrane area) the association force density and results in an additional increase of the sucked membrane length. Equations have been derived to describe this process. From the experimental conditions of an RBC aspiration and the determined tongue length corresponding to the second slow aspiration step, the association force density between the lipid bilayer and the spectrin network may be determined. From literature data and our own results a force density of between 40 and 50 Pa has been estimated. Offprint requests to: D. Lerche     

Stress-free state of the red blood cell membrane and the deformation of its skeleton

The response of a red blood cell (RBC) to deformation depends on its membrane, a composite of a lipid bilayer and a skeleton, which is a closed, twodimensional network of spectrin tetramers as its bonds. The deformation of the skeleton and its lateral redistribution are studied in terms of the RBC resting state for a fixed geometry of the RBC, partially aspirated into a micropipette. The geometry of the RBC skeleton in its initial state is taken to be either two concentric circles, a references biconcave shape or a sphere. It is assumed that in its initial state the skeleton is distributed laterally in a homogeneous manner with its bonds either unstressed, presenting its stress-free state, or prestressed. The lateral distribution was calculated using a variational calculation. It was assumed that the spectrin tetramer bonds exhibit a linear elasticity. The results showed a significant effect of the initial skeleton geometry on its lateral distribution in the deformed state. The proposed model is used to analyze the measurements of skeleton extension ratios by the method of applying two modes of RBC micropipette aspiration.     

Mechanical and functional aspects of membrane skeletons

Membrane skeletons can be characterized as cytoskeletal structures lying parallel to the bilayer part of cellular and organelle membranes. Typical examples are spectrin network and actin-myosin cortex. We approach the problem of elucidating the function of membrane skeletons by theoretically analyzing mechanical models of the cellular behavior. Membranes of different physical and chemical properties are considered. In erythrocytes and some organelles membrane bilayers are smooth and simply underlaid or overlaid by membrane skeletons. It is argued that there the role of a membrane skeleton is, either, to keep the membrane composition laterally homogeneous as it is in the case of the erythrocyte, or, that it is involved in the processes of the lateral separation of integral membrane proteins as it is happening in the case of some intermediate steps of the vesicular membrane trafficking. In the second type of membranes the bilayer part is ruffled and folded, and there the membrane skeletons play a role in the determination of the cortical tension. Here we explore in more detail the mechanical behavior of a cell with such properties of its boundary. The shape transformations are described which occur under the influence (i) of different external forces, i.e., when an originally spherical cell is aspirated into the micropipette or when such a cell is adsorbed on a flat surface, and (ii) of different internal forces on the cell boundary exerted by the cytoskeletal elements.     

Spectrin, human erythrocyte shapes, and mechanochemical properties.

Physical studies of human erythrocyte spectrin indicate that isolated spectrin dimers and tetramers in solution are worm-like coils with a persistence length of approximately 20 nm. This finding, the known polyelectrolytic nature of spectrin, and other structural information about spectrin and the membrane skeleton molecular organization have lead us to the hypothesis that the human erythrocyte membrane skeleton constitutes a two-dimensional ionic gel (swollen ionic elastomer). This concept is incorporated in what we refer to as the protein gel-lipid bilayer membrane model. The model accounts quantitatively for red elastic shear modulus and the maximum elastic extension ratio reported for the human erythrocytes membrane. Gel theory further predicts that depending on the environmental conditions, the membrane skeleton modulus of area compression may be small or large relative to the membrane elastic shear modulus. Our analyses show that the ratio between these two parameters affects both the geometry and the stability of the favored cell shapes and that the higher the membrane skeleton compressibility the smaller the values of the gel tension needed to induce cell shape transformations. The main virtue of the protein gel-lipid bilayer membrane model is that it offers a novel theoretical and molecular basis for the various mechanical properties of the membrane skeleton such as the membrane skeleton modulus of area compression and osmotic tension, and the effects of these properties on local membrane skeleton density, cell shape, and shape transformations.     

A power-law rheology-based finite element model for single cell deformation

Physical forces can elicit complex time- and space-dependent deformations in living cells. These deformations at the subcellular level are difficult to measure but can be estimated using computational approaches such as finite element (FE) simulation. Existing FE models predominantly treat cells as spring-dashpot viscoelastic materials, while broad experimental data are now lending support to the power-law rheology (PLR) model. Here, we developed a large deformation FE model that incorporated PLR and experimentally verified this model by performing micropipette aspiration on fibroblasts under various mechanical loadings. With a single set of rheological properties, this model recapitulated the diverse micropipette aspiration data obtained using three protocols and with a range of micropipette sizes. More intriguingly, our analysis revealed that decreased pipette size leads to increased pressure gradient, potentially explaining our previous counterintuitive finding that decreased pipette size leads to increased incidence of cell blebbing and injury. Taken together, our work leads to more accurate rheological interpretation of micropipette aspiration experiments than previous models and suggests pressure gradient as a potential determinant of cell injury.     

Micropipette Aspiration of Substrate-attached Cells to Estimate Cell Stiffness

Growing number of studies show that biomechanical properties of individual cells play major roles in multiple cellular functions, including cell proliferation, differentiation, migration and cell-cell interactions. The two key parameters of cellular biomechanics are cellular deformability or stiffness and the ability of the cells to contract and generate force. Here we describe a quick and simple method to estimate cell stiffness by measuring the degree of membrane deformation in response to negative pressure applied by a glass micropipette to the cell surface, a technique that is called Micropipette Aspiration or Microaspiration.Microaspiration is performed by pulling a glass capillary to create a micropipette with a very small tip (2-50 μm diameter depending on the size of a cell or a tissue sample), which is then connected to a pneumatic pressure transducer and brought to a close vicinity of a cell under a microscope. When the tip of the pipette touches a cell, a step of negative pressure is applied to the pipette by the pneumatic pressure transducer generating well-defined pressure on the cell membrane. In response to pressure, the membrane is aspirated into the pipette and progressive membrane deformation or "membrane projection" into the pipette is measured as a function of time. The basic principle of this experimental approach is that the degree of membrane deformation in response to a defined mechanical force is a function of membrane stiffness. The stiffer the membrane is, the slower the rate of membrane deformation and the shorter the steady-state aspiration length.The technique can be performed on isolated cells, both in suspension and substrate-attached, large organelles, and liposomes.Analysis is performed by comparing maximal membrane deformations achieved under a given pressure for different cell populations or experimental conditions. A "stiffness coefficient" is estimated by plotting the aspirated length of membrane deformation as a function of the applied pressure. Furthermore, the data can be further analyzed to estimate the Young''s modulus of the cells (E), the most common parameter to characterize stiffness of materials. It is important to note that plasma membranes of eukaryotic cells can be viewed as a bi-component system where membrane lipid bilayer is underlied by the sub-membrane cytoskeleton and that it is the cytoskeleton that constitutes the mechanical scaffold of the membrane and dominates the deformability of the cellular envelope. This approach, therefore, allows probing the biomechanical properties of the sub-membrane cytoskeleton.     

Molecular dynamics simulations showing 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) membrane mechanoporation damage under different strain paths

Continuum finite element material models used for traumatic brain injury lack local injury parameters necessitating nanoscale mechanical injury mechanisms be incorporated. One such mechanism is membrane mechanoporation, which can occur during physical insults and can be devastating to cells, depending on the level of disruption. The current study investigates the strain state dependence of phospholipid bilayer mechanoporation and failure. Using molecular dynamics, a simplified membrane, consisting of 72 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) phospholipids, was subjected to equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial tensile deformations at a von Mises strain rate of 5.45 × 10⁸ s^?1, resulting in velocities in the range of 1 to 4.6 m·s^?1. A water bridge forming through both phospholipid bilayer leaflets was used to determine structural failure. The stress magnitude, failure strain, headgroup clustering, and damage responses were found to be strain state-dependent. The strain state order of detrimentality in descending order was equibiaxial, 2:1 non-equibiaxial, 4:1 non-equibiaxial, strip biaxial, and uniaxial. The phospholipid bilayer failed at von Mises strains of .46, .47, .53, .77, and 1.67 during these respective strain path simulations. Additionally, a Membrane Failure Limit Diagram (MFLD) was created using the pore nucleation, growth, and failure strains to demonstrate safe and unsafe membrane deformation regions. This MFLD allowed representative equations to be derived to predict membrane failure from in-plane strains. These results provide the basis to implement a more accurate mechano-physiological internal state variable continuum model that captures lower length scale damage and will aid in developing higher fidelity injury models.     

Simulations of the erythrocyte cytoskeleton at large deformation. II. Micropipette aspiration.

Coarse-grained molecular models of the erythrocyte membrane's spectrin cytoskeleton are presented in Monte Carlo simulations of whole cells in micropipette aspiration. The nonlinear chain elasticity and sterics revealed in more microscopic cytoskeleton models (developed in a companion paper; Boey et al., 1998. Biophys. J. 75:1573-1583) are faithfully represented here by two- and three-body effective potentials. The number of degrees of freedom of the system are thereby reduced to a range that is computationally tractable. Three effective models for the triangulated cytoskeleton are developed: two models in which the cytoskeleton is stress-free and does or does not have internal attractive interactions, and a third model in which the cytoskeleton is prestressed in situ. These are employed in direct, finite-temperature simulations of erythrocyte deformation in a micropipette. All three models show reasonable agreement with aspiration measurements made on flaccid human erythrocytes, but the prestressed model alone yields optimal agreement with fluorescence imaging experiments. Ensemble-averaging of nonaxisymmetrical, deformed structures exhibiting anisotropic strain are thus shown to provide an answer to the basic question of how a triangulated mesh such as that of the red cell cytoskeleton deforms in experiment.     

Elastic properties of the red blood cell membrane that determine echinocyte deformability

The natural biconcave shape of red blood cells (RBC) may be altered by injury or environmental conditions into a spiculated form (echinocyte). An analysis is presented of the effect of such a transformation on the resistance of RBC to entry into capillary sized cylindrical tubes. The analysis accounts for the elasticity of the membrane skeleton in dilation and shear, and the local and nonlocal resistance of the bilayer to bending, the latter corresponding to different area strains in the two leaflets of the bilayer. The shape transformation is assumed to be driven by the equilibrium area difference (A₀, the difference between the equilibrium areas of the bilayer leaflets), which also affects the energy of deformation. The cell shape is approximated by a parametric model. Shape parameters, skeleton shear deformation, and the skeleton density of deformed membrane relative to the skeleton density of undeformed membrane are obtained by minimization of the corresponding thermodynamic potential. Experimentally, A₀ is modified and the corresponding discocyte–echinocyte shape transition obtained by high-pressure aspiration into a narrow pipette, and the deformability of the resulting echinocyte is examined by whole cell aspiration into a larger pipette. We conclude that the deformability of the echinocyte can be accounted for by the mechanical behavior of the normal RBC membrane, where the equilibrium area difference A₀ is modified.     

Estimating the sensitivity of mechanosensitive ion channels to membrane strain and tension

Bone adapts to its environment by a process in which osteoblasts and osteocytes sense applied mechanical strain. One possible pathway for the detection of strain involves mechanosensitive channels and we sought to determine their sensitivity to membrane strain and tension. We used a combination of experimental and computational modeling techniques to gain new insights into cell mechanics and the regulation of mechanosensitive channels. Using patch-clamp electrophysiology combined with video microscopy, we recorded simultaneously the evolution of membrane extensions into the micropipette, applied pressure, and membrane currents. Nonselective mechanosensitive cation channels with a conductance of 15 pS were observed. Bleb aspiration into the micropipette was simulated using finite element models incorporating the cytoplasm, the actin cortex, the plasma membrane, cellular stiffening in response to strain, and adhesion between the membrane and the micropipette. Using this model, we examine the relative importance of the different cellular components in resisting suction into the pipette and estimate the membrane strains and tensions needed to open mechanosensitive channels. Radial membrane strains of 800% and tensions of 5 10(-4) N.m(-1) were needed to open 50% of mechanosensitive channels. We discuss the relevance of these results in the understanding of cellular reactions to mechanical strain and bone physiology.     

Kinematics of red cell aspiration by fluorescence-imaged microdeformation.

Maps of fluorescing red cell membrane components on a pipette-aspirated projection are quantitated in an effort to elucidate and unify the heterogeneous kinematics of deformation. Transient gradients of diffusing fluorescent lipid first demonstrate the fluidity of an otherwise uniform-density bilayer and corroborate a "universal" calibration scale for relative surface density. A steep but smooth and stable gradient in the densities of the skeleton components spectrin, actin, and protein 4.1 is used to estimate large elastic strains along the aspirated skeleton. The deformation fields are argued to be an unhindered response to loading in the surface normal direction. Density maps intermediate to those of the compressible skeleton and fluid bilayer are exhibited by particular transmembrane proteins (e.g., Band 3) and yield estimates for the skeleton-connected fractions. Such connected proteins appear to occupy a significant proportion of the undeformed membrane surface and can lead to steric exclusion of unconnected integral membrane proteins from regions of network condensation. Consistent with membrane repatterning kinematics in reversible deformation, final vesiculation of the projection tip produces a cell fragment concentrated in freely diffusing proteins but depleted of skeleton.     

The effect of rimantadine on the structure of model and biological membranes

The action of the antiviral drug rimantadine on the structure of bilayer lipid membranes (BLM) and RBC membranes was investigated. Structural changes in BLM were recorded by ionophore conductivity changes and by changes in the third harmonic of capacity current signal due to lateral compression of BLM in an electric field. It was shown that the adsorption of rimantadine on BLM results in an increase in ionophore mobility in bilayer membranes of dioleolyllecithin (DOL) and common lipids of bovine brain (CL) and in a decrease in those of azolectin (A). Relative changes in the third harmonic signal also depend on the membrane composition and have different signs. The results may be explained by the rimantadine action on the lipid bilayer structure: "rigidification" of A-membranes and "fluidization" of BLM from DOL and CL. Structural reorganization of RBC membranes as investigated by the ability of the cells to enter a micropipette (inner diameter greater than or equal to 3 microns) thereby undergoing deformation. It was shown that rimantadine influences RBC deformability due to drug induced inhomogenous mechanical membrane properties. Also, rimantadine accelerated the process of artificially induced aggregation of erythrocytes. The relation of the effects on artificial and biological membranes, and the structural changes in the lipid phase of membrane are discussed.     

The human erythrocyte membrane skeleton may be an ionic gel

In the first paper in this series (Stokke et al. Eur Biophys J 1986, 13:203-218) we developed the general theory of the mechanochemical properties and the elastic free energy of the protein gel--lipid bilayer membrane model. Here we report on an extensive numerical analysis of the human erythrocyte shapes and shape transformations predicted by this new cell membrane model. We have calculated the total elastic free energy of deformation of four different cell shape classes: disc-shaped cells, cup-shaped cells, crenated cells, and cells with membrane invaginations. We find that which of these shape classes is favoured depends strongly on the spectrin gel osmotic tension, IIGu, and the surface tensions, IIEu and IIPu, of the extracellular and protoplasmic halves of the membrane lipid bilayer, respectively. For constant ratio IIEu/IIPu greater than O large negative or positive values of IIGu favour respectively the crenated and invaginated cell shape classes. For small absolute values of IIGu, IIEu, and IIPu, biconcave or cup-shaped cells are the stable ones. Our numerical analysis shows that the higher the membrane skeleton compressibility is, the smaller are the values of IIGu needed to induce cell shape transformation. We find that the stable and metastable shapes of discocytes and stomatocytes generally depend both on the shape of the stressfree membrane skeleton and the membrane skeleton compressibility.     

The human erythrocyte membrane skeleton may be an ionic gel. III. Micropipette aspiration of unswollen erythrocytes

We have carried out a theoretical analysis of micropipette aspiration of unswollen erythrocytes using the protein-gel-lipid-bilayer membrane model and taking into account that the modulus of area compression of the membrane skeleton may depend on the environmental conditions. Our analysis shows that the aspiration pressure needed to obtain a certain membrane projection length is strongly dependent on the ratio between the membrane skeleton modulus of area compression and the elastic shear modulus. Our analysis therefore predicts that micropipette aspiration of unswollen erythrocytes may be a sensitive method for detection of changes in this ratio. The analysis thus also shows that micropipette aspiration of unswollen erythrocytes can not be used to determine the membrane shear modulus unless something is known about the membrane skeleton modulus of area compression.     

Remodeling the shape of the skeleton in the intact red cell.

The role of the membrane skeleton in determining the shape of the human red cell was probed by weakening it in situ with urea, a membrane-permeable perturbant of spectrin. Urea by itself did not alter the biconcave disk shape of the red cell; however, above threshold conditions (1.5 M, 37 degrees C, 10 min), it caused an 18% reduction in the membrane elastic shear modulus. It also potentiated the spiculation of cells by lysophosphatidylcholine. These findings suggest that the contour of the resting cell is not normally dependent on the elasticity of or tension in the membrane skeleton. Rather, the elasticity of the skeleton stabilizes membranes against deformation. Urea treatment also caused the projections induced both by micropipette aspiration and by lysophosphatidylcholine to become irreversible. Furthermore, urea converted the axisymmetric conical spicules induced by lysophosphatidylcholine into irregular, curved and knobby spicules; i.e., echinocytosis became acanthocytosis. Unlike controls, the ghosts and membrane skeletons obtained from urea-generated acanthocytes were imprinted with spicules. These data suggest that perturbing interprotein associations with urea in situ allowed the skeleton to evolve plastically to accommodate the contours imposed upon it by the overlying membrane.     

Membrane bending is critical for the stability of voltage sensor segments in the membrane

The interaction between membrane proteins and the surrounding membrane is becoming increasingly appreciated for its role in regulating protein function, protein localization, and membrane morphology. In particular, recent studies have suggested that membrane deformation is needed to stably accommodate proteins harboring charged amino acids in their transmembrane (TM) region, as it is energetically prohibitive to bury charge in the hydrophobic core of the bilayer. Unfortunately, current computational methods are poorly equipped for describing such deformations, as atomistic simulations are often too short to observe large-scale membrane reorganization and most continuum approaches assume a flat membrane. Previously, we developed a method that overcomes these shortcomings by using elasticity theory to characterize equilibrium membrane distortions in the presence of a TM protein, while using traditional continuum electrostatic and nonpolar energy models to determine the energy of the protein in the membrane. Here, we linked the elastostatics, electrostatics, and nonpolar numeric solvers to permit the calculation of energies for nontrivial membrane deformations. We then coupled this procedure to a robust search algorithm that identifies optimal membrane shapes for a TM protein of arbitrary chemical composition. This advance now permits us to explore a host of biological phenomena that were beyond the scope of our original method. We show that the energy required to embed charged residues in the membrane can be highly nonadditive, and our model provides a simple mechanical explanation for this nonadditivity. Our results also predict that isolated voltage sensor segments do not insert into rigid membranes, but membrane bending dramatically stabilizes these proteins in the bilayer despite their high charge content. Additionally, we use the model to explore hydrophobic mismatch with regard to nonpolar peptides and mechanosensitive channels. Our method is in quantitative agreement with molecular dynamics simulations at a tiny fraction of the computational cost.