Erythrocyte Membrane Model with Explicit Description of the Lipid Bilayer and the Spectrin Network

The membrane of the red blood cell (RBC) consists of spectrin tetramers connected at actin junctional complexes, forming a two-dimensional (2D) sixfold triangular network anchored to the lipid bilayer. Better understanding of the erythrocyte mechanics in hereditary blood disorders such as spherocytosis, elliptocytosis, and especially, sickle cell disease requires the development of a detailed membrane model. In this study, we introduce a mesoscale implicit-solvent coarse-grained molecular dynamics (CGMD) model of the erythrocyte membrane that explicitly describes the phospholipid bilayer and the cytoskeleton, by extending a previously developed two-component RBC membrane model. We show that the proposed model represents RBC membrane with the appropriate bending stiffness and shear modulus. The timescale and self-consistency of the model are established by comparing our results with experimentally measured viscosity and thermal fluctuations of the RBC membrane. Furthermore, we measure the pressure exerted by the cytoskeleton on the lipid bilayer. We find that defects at the anchoring points of the cytoskeleton to the lipid bilayer (as in spherocytes) cause a reduction in the pressure compared with an intact membrane, whereas defects in the dimer-dimer association of a spectrin filament (as in elliptocytes) cause an even larger decrease in the pressure. We conjecture that this finding may explain why the experimentally measured diffusion coefficients of band-3 proteins are higher in elliptocytes than in spherocytes, and higher than in normal RBCs. Finally, we study the effects that possible attractive forces between the spectrin filaments and the lipid bilayer have on the pressure applied on the lipid bilayer by the filaments. We discover that the attractive forces cause an increase in the pressure as they diminish the effect of membrane protein defects. As this finding contradicts with experimental results, we conclude that the attractive forces are moderate and do not impose a complete attachment of the filaments to the lipid bilayer.     

Two-component coarse-grained molecular-dynamics model for the human erythrocyte membrane

We present a two-component coarse-grained molecular-dynamics model for simulating the erythrocyte membrane. The proposed model possesses the key feature of combing the lipid bilayer and the erythrocyte cytoskeleton, thus showing both the fluidic behavior of the lipid bilayer and the elastic properties of the erythrocyte cytoskeleton. In this model, three types of coarse-grained particles are introduced to represent clusters of lipid molecules, actin junctions, and band-3 complexes, respectively. The proposed model facilitates simulations that span large length scales (approximately micrometers) and timescales (approximately milliseconds). By tuning the interaction potential parameters, we were able to control the diffusivity and bending rigidity of the membrane model. We studied the membrane under shearing and found that at a low shear strain rate, the developed shear stress was due mainly to the spectrin network, whereas the viscosity of the lipid bilayer contributed to the resulting shear stress at higher strain rates. In addition, we investigated the effects of a reduced spectrin network connectivity on the shear modulus of the membrane.     

Red Blood Cell Shape and Fluctuations: Cytoskeleton Confinement and ATP Activity

We review recent theoretical work that analyzes experimental measurements of the shape and fluctuations of red blood cells. Particular emphasis is placed on the role of the cytoskeleton and cell elasticity and we contrast the situation of elastic cells with that of fluid-filled vesicles. In red blood cells (RBCs), the cytoskeleton consists of a two-dimensional network of spectrin proteins. Our analysis of the wave vector and frequency dependence of the fluctuation spectrum of RBCs indicates that the spectrin network acts as a confining potential that reduces the fluctuations of the lipid bilayer membrane. However, since the cytoskeleton is only sparsely connected to the bilayer, one cannot regard the composite cytoskeleton membrane as a polymerized object with a shear modulus. The sensitivity of RBC fluctuations and shapes to ATP concentration may reflect the transient defects induced in the cytoskeleton network by ATP.     

Characterization of immunoglobulin binding to isolated human erythrocyte membranes: evidence for selective, temperature-induced binding of naturally occurring autoantibodies to the cytoskeleton

Human plasma contains naturally occurring autoantibodies to the predominant components of the erythrocyte membrane: band 3 and spectrin bands 1 and 2 of the cytoskeleton. The titer of cytoskeletal plasma autoantibodies increases in various hemolytic conditions, suggesting that opsonization of the cytoskeleton may play an important role in the clearance of hemolyzed (not senescent) erythrocytes from the circulation. In this study, we use Alexa Fluor 488 goat anti-human IgG conjugate (Molecular Probes, Eugene, OR, USA), to characterize plasma immunoglobulin binding to erythrocyte membranes from osmotically hemolyzed cells ('ghosts'). The results show that exposure of ghosts to plasma results in 4-fold more immunoglobulin binding to the cytoskeleton than is bound to the proteins contained within the lipid bilayer. Preincubation of the ghosts at 37 degrees C causes 8-fold more immunoglobulin binding to the cytoskeleton compared to bilayer proteins. This temperature-induced change resulted from selective immunoglobulin binding to the cytoskeleton, with no change in immunoglobulin binding to bilayer proteins. However, the rate of increase in cytoskeletal antigenicity at 37 degrees C did correlate with the rate of a conformational change in band 3, a transmembrane protein which serves as a major membrane attachment site for the cytoskeleton. The results of this study suggest that the cytoskeleton is the primary target in the opsonization of hemolyzed erythrocyte membranes by naturally occurring plasma autoantibodies. The conformational changes which occur in ghosts at 37 degrees C are associated with selective exposure of new immunoglobulin binding sites on the cytoskeleton, and with a change in the structure of band 3. We propose a model suggesting that opsonization of the cytoskeleton occurs prior to the decomposition of hemolyzed erythrocytes at 37 degrees C.     

Ultrastructure of the human erythrocyte cytoskeleton and its attachment to the membrane

We attached paraformaldehyde-fixed human erythrocyte ghosts to coated coverslips and sheared them to expose the cytoskeleton. Quick-freeze, deep-etch, rotary-replication, or tannic acid/osmium fixation and plastic embedding revealed the cytoskeleton as a dense network of intersecting straight filaments. Previous negative stain studies on spread skeletons found 5-6 spectrin tetramers intersecting at each actin oligomer, with an estimated 250 such intersections/microns 2 of membrane. In contrast, we found 3-4 filaments at each intersection and approximately 400 intersections/microns 2 of membrane. Immunogold labeling verified that the filaments were spectrin, but their lengths (29-37 nm) were approximately one-third that of extended spectrin dimers. The length and diameter of the filaments were sufficient to accommodate spectrin dimers, but not spectrin tetramers. Our results suggest that, in situ, spectrin dimers may associate as hexamers and octamers, rather than tetramers. We present several explanations that can reconcile our observations on intact cytoskeletons with previous reports on spread material. Extracting sheared ghosts with solutions of low ionic strength removed the cytoskeleton to reveal projections from the cytoplasmic surface of the membrane. These projections contained band 3, as shown by immunogold labeling, and they aggregated to a similar extent as intramembrane particles (IMP) when the cytoskeleton was removed, suggesting a direct relationship between these structures. Quantification indicated a stoichiometry of 2 IMP for each cytoplasmic projection. Cytoplasmic projections presumably contain other proteins besides band 3 since further treatment with high ionic strength solutions extracts peripheral proteins and reduces the diameter of projections by approximately 3 nm.     

Curling and Local Shape Changes of Red Blood Cell Membranes Driven by Cytoskeletal Reorganization

Human red blood cells (RBCs) lack the actin-myosin-microtubule cytoskeleton that is responsible for shape changes in other cells. Nevertheless, they can display highly dynamic local deformations in response to external perturbations, such as those that occur during the process of apical alignment preceding merozoite invasion in malaria. Moreover, after lysis in divalent cation-free media, the isolated membranes of ruptured ghosts show spontaneous inside-out curling motions at the free edges of the lytic hole, leading to inside-out vesiculation. The molecular mechanisms that drive these rapid shape changes are unknown. Here, we propose a molecular model in which the spectrin filaments of the RBC cortical cytoskeleton control the sign and dynamics of membrane curvature depending on whether the ends of the filaments are free or anchored to the bilayer. Computer simulations of the model reveal that curling, as experimentally observed, can be obtained either by an overall excess of weakly-bound filaments throughout the cell, or by the flux of such filaments toward the curling edges. Divalent cations have been shown to arrest the curling process, and Ca²⁺ ions have also been implicated in local membrane deformations during merozoite invasion. These effects can be replicated in our model by attributing the divalent cation effects to increased filament-membrane binding. This process converts the curl-inducing loose filaments into fully bound filaments that arrest curling. The same basic mechanism can be shown to account for Ca²⁺-induced local and dynamic membrane deformations in intact RBCs. The implications of these results in terms of RBC membrane dynamics under physiological, pathological, and experimental conditions is discussed.     

Regulation of erythrocyte ghost membrane mechanical stability by chlorpromazine

Chlorpromazine (CPZ), a widely used tranquilizer, is known to induce stomatocytic shape changes in human erythrocytes. However, the effect of CPZ on membrane mechanical properties of erythrocyte membranes has not been documented. In the present study we show that CPZ induces a dose-dependent increase in mechanical stability of erythrocyte ghost membrane. Furthermore, we document that spectrin specifically binds to CPZ intercalated into inside-out vesicles depleted of all peripheral proteins. These findings imply that CPZ-induced mechanical stabilization of the erythrocyte ghost membranes may be mediated by direct binding of spectrin to the bilayer. Membrane active drugs that partition into lipid bilayer can thus induce cytoskeletal protein interactions with the membrane and modulate membrane material properties.     

Bilayer/cytoskeleton interactions in lipid-symmetric erythrocytes assessed by a photoactivable phospholipid analogue

Two mechanisms have been proposed for maintenance of transbilayer phospholipid asymmetry in the erythrocyte plasma membrane, one involving specific interactions between the aminophospholipids of the inner leaflet of the bilayer and the cytoskeleton, particularly spectrin, and the other involving the aminophospholipid translocase. If the former mechanism is correct, then erythrocytes which have lost their asymmetric distribution of phospholipids should display altered bilayer/cytoskeleton interactions. To test this possibility, normal erythrocytes, erythrocytes from patients with chronic myelogenous leukemia or sickle disease, and lipid-symmetric and -asymmetric erythrocyte ghosts were labeled with the radioactive photoactivable analogue of phosphatidylethanolamine, 2-(2-azido-4-nitrobenzoyl)-1-acyl-sn-glycero-3-phospho[14C]ethanolamine ([14C]AzPE), previously shown to label cytoskeletal proteins from the bilayer. The labeling pattern of cytoskeletal proteins in pathologic erythrocytes and lipid-asymmetric erythrocyte ghosts was indistinguishable from normal erythrocytes, indicating that the probe detects no differences in bilayer/cytoskeleton interactions in these cells. In contrast, in lipid-symmetric erythrocyte ghosts, labeling of bands 4.1 and 4.2 and actin, and to a lesser extent ankyrin, by [14C]AzPE was considerably reduced. Significantly, however, labeling of spectrin was unaltered in the lipid-symmetric ghosts, suggesting that its relationship with the bilayer is normal in these lipid-symmetric cells. These results do not support a model in which spectrin is involved in the maintenance of an asymmetric distribution of phospholipids in erythrocytes.     

Intramembrane particle aggregation in erythrocyte ghosts. II. The influence of spectrin aggregation.

Physicochemical properties of mixtures of spectrin and actin extracted from human erythrocyte ghosts have been correlated with ultrastructural changes observed in freeze-fractured erythrocyte membranes. (1) Extracted mixtures of spectrin and actin have a very low solubility (less than 30 mug/ml) near their isoelectric point, pH 4.8. These mixtures are also precipitated by low concentrations of Ca2+, Mg2+, polylysine or basic proteins. (2) All conditions which precipitate extracts of spectrin and actin also induce aggregation of the intramembrane particles in spectrin-depleted erythrocyte ghosts. Precipitation of the residual spectrin molecules into small patches on the cytoplasmic surface of the ghost membrane is thought to be the cause of particle aggregations, implying an association between the spectrin molecules and the intramembrane particles. (3) When fresh ghosts are exposed to conditions which precipitate extracts of spectrin and actin, only limited particle aggregation occurs. Instead, the contraction of the intact spectrin meshwork induced by the precipitation conditions compresses the lipid bilayer of the membrane, causing it to bleb off particle-free, protein-free vesicles. (4) The absence of protein in these lipid vesicles implies that all the proteins of the erythrocyte membrane are immobilized by association with either the spectrin meshwork or the intramembrane particles.     

Effect of Cu2+-ascorbic acid on lipid peroxidation,Mg2+-ATPase activity and spectrin of RBC membrane and reversal by erythropoietin

The effect of erythropoietin (Ep), a glycoprotein hormone, has been studied on lipid peroxidation induced by Cu²⁺ and ascorbate in vitro, Mg²⁺ ATPase activity and spectrin of RBC membrane. Our present investigation reveals that Cu²⁺ and ascorbic acid increases lipid peroxidation of RBC membrane significantly. It has further been observed that under the same experimental condition spectrin, a major cytoskeleton membrane protein, and Mg²⁺-ATPase activity of RBC membrane decrease significantly. However, exogenous administration of Ep completely restores lipid peroxidation and Mg²⁺-ATPase activity and partially recovers spectrin of RBC membrane.     

Red blood cell membrane fluctuations and shape controlled by ATP-induced cytoskeletal defects

We show theoretically how adenosine 5'-triphosphate (ATP)-induced dynamic dissociations of spectrin filaments (from each other and from the membrane) in the cytoskeleton network of red blood cells (RBC) can explain in a unified manner both the measured fluctuation amplitude as well as the observed shape transformations as a function of intracellular ATP concentration. Static defects can be induced by external stresses such as those present when RBCs pass through small capillaries. We suggest that the partially freed actin at these defect sites may explain the activation of the CFTR membrane-bound protein and the subsequent release of ATP by RBCs subjected to deformations. Our theoretical predictions can be tested by experiments that measure the correlation between variations in the binding of actin to spectrin, the activity of CFTR, and the amount of ATP released.     

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     

Barrier-free paths of directed protein motion in the erythrocyte plasma membrane.

A model is presented for the steric interaction between a plasma membrane protein and the membrane cytoskeleton in the human erythrocyte. The cytoskeleton is treated as a network of polymer chains attached to a flat bilayer, and the membrane protein is a hemisphere of effective radius R(e) with center on the bilayer edge. The simulation is used to investigate the barrier-free path L for linear guided motion of a protein in the bilayer plane. It is shown that the barrier-free paths of small proteins can be used to extract the effective in-plane diameter of cytoskeletal components. For example, the in-plane diameter of an ankyrin attachment site is found to be approximately 12 nm in the simulation, or twice the computational spectrin diameter. The barrier-free paths of large proteins (R(e) > 23 nm) vanish when the proteins are corralled by the cytoskeleton. For intermediate size proteins, L decreases approximately as L is directly proportional to S-1.4 where S is proportional to the sum of the protein and cytoskeleton chain radii.     

Triton shells of intact erythrocytes

About 40% of human erythrocyte membrane protein is resistant to solubilization in 0.5% Triton X-114. These components comprise a structure called a Triton shell roughly similar in size and shape to the original erythrocyte and thus constitute a cytoskeleton. With increasing concentrations of Triton the lipid content of the Triton shell decreases dramatically, whereas the majority of the protein components remain constant. Exceptions to this rule include proteins contained in band 3, the presumed anion channel, and in band 4 which decrease with increasing Triton concentration. The Triton-insoluble complex includes spectrin (bands 1 and 2), actin (band 5), and bands 3′ and 7. Component 3′ has an apparent molecular weight of 88,000 daltons as does 3; but unlike 3, it is insensitive to protease treatment of the intact cell, has a low extinction coefficient at 280 nm, and is solubilized from the shells in alkaline water solutions. Component 7 also has a low extinction coefficient at 280 nm. Spectrin alone is solubilized from the Triton shells in isotonic media. The solubilized spectrin contains no bound Triton and coelectrophoreses with spectrin eluted in hypotonic solutions from ghosts. Electron micrographs of fixed Triton shells stained with uranyl acetate show the presence of numerous filaments which appear beaded and are 80–120 Å in diameter. The filaments cannot be composed mainly of actin, but enough spectrin is present to form the filaments. Triton shells may provide an excellent source of material useful in the investigation of the erythrocyte cytoskeleton.     

CASK and protein 4.1 support F-actin nucleation on neurexins

Rearrangements of the actin cytoskeleton are involved in a variety of cellular processes from locomotion of cells to morphological alterations of the cell surface. One important question is how local interactions of cells with the extracellular space are translated into alterations of their membrane organization. To address this problem, we studied CASK, a member of the membrane-associated guanylate kinase homologues family of adaptor proteins. CASK has been shown to bind the erythrocyte isoform of protein 4.1, a class of proteins that promote formation of actin/spectrin microfilaments. In neurons, CASK also interacts via its PDZ domain with the cytosolic C termini of neurexins, neuron-specific cell-surface proteins. We now show that CASK binds a brain-enriched isoform of protein 4.1, and nucleates local assembly of actin/spectrin filaments. These interactions can be reconstituted on the cytosolic tail of neurexins. Furthermore, CASK can be recovered with actin filaments prepared from rat brain extracts, and neurexins are recruited together with CASK and protein 4.1 into these actin filaments. Thus, analogous to the PDZ-domain protein p55 and glycophorin C at the erythrocyte membrane, a similar complex comprising CASK and neurexins exists in neurons. Our data suggest that intercellular junctions formed by neurexins, such as junctions initiated by beta-neurexins with neuroligins, are at least partially coupled to the actin cytoskeleton via an interaction with CASK and protein 4.1.     

Native ultrastructure of the red cell cytoskeleton by cryo-electron tomography

Erythrocytes possess a spectrin-based cytoskeleton that provides elasticity and mechanical stability necessary to survive the shear forces within the microvasculature. The architecture of this membrane skeleton and the nature of its intermolecular contacts determine the mechanical properties of the skeleton and confer the characteristic biconcave shape of red cells. We have used cryo-electron tomography to evaluate the three-dimensional topology in intact, unexpanded membrane skeletons from mouse erythrocytes frozen in physiological buffer. The tomograms reveal a complex network of spectrin filaments converging at actin-based nodes and a gradual decrease in both the density and the thickness of the network from the center to the edge of the cell. The average contour length of spectrin filaments connecting junctional complexes is 46 ± 15 nm, indicating that the spectrin heterotetramer in the native membrane skeleton is a fraction of its fully extended length (∼190 nm). Higher-order oligomers of spectrin were prevalent, with hexamers and octamers seen between virtually every junctional complex in the network. Based on comparisons with expanded skeletons, we propose that the oligomeric state of spectrin is in a dynamic equilibrium that facilitates remodeling of the network as the cell changes shape in response to shear stress.     

Membrane tethers formed from blood cells with available area and determination of their adhesion energy

Fundamental to all mammalian cells is the adherence of the lipid bilayer membrane to the underlying membrane associated cytoskeleton. To investigate this adhesion, we physically detach the lipid membrane from the cell by mechanically forming membrane tethers. For the most part these have been tethers formed from either neutrophils or red cells. Here we do a simple thermodynamic analysis of the tether formation process using the entire cell, including tether, as the control volume. For a neutrophil, we show that the total adhesion energy per unit area between lipid membrane and cytoskeleton depends on the square of the tether force. For a flaccid red cell, we show that the total adhesion energy minus the tension in the spectrin cytoskeleton depends also on the square of the tether force. Finally, we discuss briefly the viscous flow of membrane. Using published data we calculate and compare values for the various adhesion energies and viscosities.     

Control of band 3 lateral and rotational mobility by band 4.2 in intact erythrocytes: release of band 3 oligomers from low-affinity binding sites.

Band 4.2 is a human erythrocyte membrane protein of incompletely characterized structure and function. Erythrocytes deficient in band 4.2 protein were used to examine the functional role of band 4.2 in intact erythrocyte membranes. Both the lateral and the rotational mobilities of band 3 were increased in band 4.2-deficient erythrocytes compared to control cells. In contrast, the lateral mobility of neither glycophorins nor a fluorescent phospholipid analog was altered in band 4.2-deficient cells. Compared to controls, band 4.2-deficient erythrocytes manifested a decreased ratio of band 3 to spectrin, and band 4.2-deficient membrane skeletons had decreased extractability of band 3 under low-salt conditions. Normal band 4.2 was found to bind to spectrin in solution and to promote the binding of spectrin to ankyrin-stripped inside-out vesicles. We conclude that band 4.2 provides low-affinity binding sites for both band 3 oligomers and spectrin dimers on the human erythrocyte membrane. Band 4.2 may serve as an accessory linking protein between the membrane skeleton and the overlying lipid bilayer.     

Regulation of Protein Mobility via Thermal Membrane Undulations

The in-plane diffusivelike motion of membrane bound proteins on the surface of cells is considered. We suggest, on the basis of theoretical arguments and simulation, that thermally excited undulations of the lipid bilayer may serve as a mechanism for proteins to hop between adjacent regions on the cell surface separated by barriers composed of internal cellular structure (e.g., the cytoskeleton). We specifically investigate the mobility of band 3 dimer on the surface of red blood cells where the spectrin cytoskeletal meshwork defines a series of “corrals” on the cell surface known to hinder protein motion. Previous models of this system have postulated that the cytoskeleton must deform to allow passage of membrane bound proteins out of these corral regions and have ignored fluctuations of the bilayer. Our model provides a complementary mechanism and we posit that the mobility of real proteins in real cells is likely the result of several mechanisms acting in parallel.