Vesiculation of healthy and defective red blood cells

Mechanisms underlying the anti-aging activity of bergamot (Citrus bergamia) extract in human red blood cells

Introduction: Aging is a process characterised by a decline in physiological functions. Reactive species play a crucial role in the aging rate. Due to the close relationship between aging and oxidative stress, functional foods rich in phytochemicals are excellent candidates to neutralise age-related changes. Aim: This investigation aims to verify the potential protective role of bergamot (Citrus bergamia, Femminello cultivar) peel and juice extract in a model of aging represented by human red blood cells (RBCs) exposed to D-Galactose (DGal). Methods: Bergamot peel and juice extracts were subjected to RP-HPLC/PDA/MS for determination of their composition in bioactive compounds. Markers of oxidative stress, including ROS production, thiobarbituric acid reactive substances (TBARS) levels -a marker of lipid peroxidation, oxidation of total protein sulfhydryl groups, as well as the expression and anion exchange capability of band 3 and glycated haemoglobin (A1c) production have been investigated in RBCs treated with D-Gal for 24 h, with or without pre-incubation for 15 min with 5 μg/mL peel or juice extract. In addition, the activity of the endogenous antioxidant system, including catalase (CAT) and superoxide dismutase (SOD), as well as the diversion of the RBC metabolism from glycolysis towards the pentose phosphate pathway shunt, as denoted by activation of glucose-6-phosphate dehydrogenase (G6PDH), have been explored. Results: Data shown here suggest that bergamot peel and juice extract i) prevented the D-Gal-induced ROS production, and consequently, oxidative stress injury to biological macromolecules including membrane lipids and proteins; ii) significantly restored D-Gal-induced alterations in the distribution and ion transport kinetics of band 3; iii) blunted A1c production; iv) effectively impeded the over-activation of the endogenous antioxidant enzymes CAT and SOD; and v) significantly prevented the activation of G6PDH. Discussion: These results further contribute to shed light on aging mechanisms in human RBCs and identify bergamot as a functional food rich in natural antioxidants useful for prevention and treatment of oxidative stress-related changes, which may lead to pathological states during aging.

A combined computational and experimental investigation of the filtration function of splenic macrophages in sickle cell disease

Being the largest lymphatic organ in the body, the spleen also constantly controls the quality of red blood cells (RBCs) in circulation through its two major filtration components, namely interendothelial slits (IES) and red pulp macrophages. In contrast to the extensive studies in understanding the filtration function of IES, there are relatively fewer works on investigating how the splenic macrophages retain the aged and diseased RBCs, i.e., RBCs in sickle cell disease (SCD). Herein, we perform a computational study informed by companion experiments to quantify the dynamics of RBCs captured and retained by the macrophages. We first calibrate the parameters in the computational model based on microfluidic experimental measurements for sickle RBCs under normoxia and hypoxia, as those parameters are not available in the literature. Next, we quantify the impact of a set of key factors that are expected to dictate the RBC retention by the macrophages in the spleen, namely, blood flow conditions, RBC aggregation, hematocrit, RBC morphology, and oxygen levels. Our simulation results show that hypoxic conditions could enhance the adhesion between the sickle RBCs and macrophages. This, in turn, increases the retention of RBCs by as much as five-fold, which could be a possible cause of RBC congestion in the spleen of patients with SCD. Our study on the impact of RBC aggregation illustrates a 'clustering effect', where multiple RBCs in one aggregate can make contact and adhere to the macrophages, leading to a higher retention rate than that resulting from RBC-macrophage pair interactions. Our simulations of sickle RBCs flowing past macrophages for a range of blood flow velocities indicate that the increased blood velocity could quickly attenuate the function of the red pulp macrophages on detaining aged or diseased RBCs, thereby providing a possible rationale for the slow blood flow in the open circulation of the spleen. Furthermore, we quantify the impact of RBC morphology on their tendency to be retained by the macrophages. We find that the sickle and granular-shaped RBCs are more likely to be filtered by macrophages in the spleen. This finding is consistent with the observation of low percentages of these two forms of sickle RBCs in the blood smear of SCD patients. Taken together, our experimental and simulation results aid in our quantitative understanding of the function of splenic macrophages in retaining the diseased RBCs and provide an opportunity to combine such knowledge with the current knowledge of the interaction between IES and traversing RBCs to apprehend the complete filtration function of the spleen in SCD.

Evolution of surface area and membrane shear modulus of matured human red blood cells during mechanical fatigue

Mechanical properties of red blood cells (RBCs) change during their senescence which supports numerous physiological or pathological processes in circulatory systems by providing crucial cellular mechanical environments of hemodynamics. However, quantitative studies on the aging and variations of RBC properties are largely lacking. Herein, we investigate morphological changes, softening or stiffening of single RBCs during aging using an in vitro mechanical fatigue model. Using a microfluidic system with microtubes, RBCs are repeatedly subjected to stretch and relaxation as they squeeze into and out of a sudden contraction region. Geometric parameters and mechanical properties of healthy human RBCs are characterized systematically upon each mechanical loading cycle. Our experimental results identify three typical shape transformations of RBCs during mechanical fatigue, which are all strongly associated with the loss of surface area. We constructed mathematical models for the evolution of surface area and membrane shear modulus of single RBCs during mechanical fatigue, and quantitatively developed an ensemble parameter to evaluate the aging status of RBCs. This study provides not only a novel in vitro fatigue model for investigating the mechanical behavior of RBCs, but also an index closely related to the age and inherent physical properties for a quantitative differentiation of individual RBCs.

Membrane Lesions and Reduced Life Span of Red Blood Cells in Preeclampsia as Evidenced by Atomic Force Microscopy

Preeclampsia (PE) presents with maternal de novo hypertension and significant proteinuria and is one of the leading causes of maternal and perinatal morbidity and mortality with unknown etiology. The disease is associated with inflammatory vascular response and severe red blood cell (RBC) morphology changes. This study examined the nanoscopic morphological changes of RBCs from PE women versus normotensive healthy pregnant controls (PCs) and non-pregnant controls (NPCs) applying atomic force microscopy (AFM) imaging. The results revealed that the membrane of fresh PE RBCs differed significantly from healthy ones by the presence of invaginations and protrusions and an increased roughness value (R_rms) (4.7 ± 0.8 nm for PE vs. 3.8 ± 0.5 nm and 2.9 ± 0.4 nm for PCs and NPCs, respectively). PE-cells aging resulted in more pronounced protrusions and concavities, with exponentially increasing R_rms values, in contrast to the controls, where the R_rms parameter decreased linearly with time. The R_rms, evaluated on a 2 × 2 µm² scanned area, for senescent PE cells (13 ± 2.0 nm) was significantly higher (p < 0.01) than that of PCs (1.5 ± 0.2 nm) and NPCs (1.9 ± 0.2 nm). Furthermore, the RBCs from PE patients appeared fragile, and often only ghosts were observed instead of intact cells at 20-30 days of aging. Oxidative-stress simulation on healthy cells led to RBC membrane features similar to those observed for PE cells. The results demonstrate that the most pronounced effects on RBCs in PE patients are related to impaired membrane homogeneity and strongly altered roughness values, as well as to vesiculation and ghost formation in the course of cell aging.

Circulating cell clusters aggravate the hemorheological abnormalities in COVID-19

Microthrombi and circulating cell clusters (CCCs) are common microscopic findings in patients with COVID-19 at different stages in the disease course, implying that they may function as the primary drivers in disease progression. Inspired by a recent flow imaging cytometry study of the blood samples from patients with COVID-19, we perform computational simulations to investigate the dynamics of different types of CCCs, namely white blood cell (WBC) clusters, platelet clusters and red blood cell (RBC) clusters, over a range of shear flows and quantify their impact on the viscosity of the blood. Our simulation results indicate that the increased level of fibrinogen in patients with COVID-19 can promote the formation of RBC clusters at relatively low shear rates, thereby elevating the blood viscosity, a mechanism that also leads to an increase in viscosity in other blood diseases, such as sickle cell disease and type 2 diabetes mellitus. We further discover that the presence of WBC clusters could also aggravate the abnormalities of local blood rheology. In particular, the extent of elevation of the local blood viscosity is enlarged as the size of the WBC clusters grows. On the other hand, the impact of platelet clusters on the local rheology is found to be negligible, which is likely due to the smaller size of the platelets. The difference in the impact of WBC and platelet clusters on local hemorheology provides a compelling explanation for the clinical finding that the number of WBC clusters is significantly correlated with thrombotic events in COVID-19 whereas platelet clusters do not. Overall, our study demonstrates that our computational models based on dissipative particle dynamics can serve as a powerful tool to conduct quantitative investigation of the mechanism causing the pathological alterations of hemorheology and explore their connections to the clinical manifestations in COVID-19.

Morphometric and Nanomechanical Features of Erythrocytes Characteristic of Early Pregnancy Loss

Early pregnancy loss (EPL) is estimated to be between 15 and 20% of all adverse pregnancies. Approximately, half of EPL cases have no identifiable cause. Herein, we apply atomic force microscopy to evaluate the alteration of morphology and nanomechanics of erythrocytes from women with EPL with unknown etiology, as compared to healthy pregnant (PC) and nonpregnant women (NPC). Freshly isolated erythrocytes from women with EPL differ in both the roughness value (4.6 ± 0.3 nm, p < 0.05), and Young’s modulus (2.54 ± 0.6 MPa, p < 0.01) compared to the values for NPC (3.8 ± 0.4 nm and 0.94 ± 0.2 MPa, respectively) and PC (3.3 ± 0.2 nm and 1.12 ± 0.3 MPa, respectively). Moreover, we find a time-dependent trend for the reduction of the cells’ morphometric parameters (cells size and surface roughness) and the membrane elasticity—much faster for EPL than for the two control groups. The accelerated aging of EPL erythrocytes is expressed in faster morphological shape transformation and earlier occurrence of spiculated and spherical-shaped cells, reduced membrane roughness and elasticity with aging evolution. Oxidative stress in vitro contributed to the morphological cells’ changes observed for EPL senescent erythrocytes. The ultrastructural characteristics of cells derived from women with miscarriages show potential as a supplementary mark for a pathological state.

Multiphysics and multiscale modeling of microthrombosis in COVID-19

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. Because of the infectious nature of SARS-CoV-2, patients’ fresh blood samples are limited to access for in vitro experimental investigations. Herein, we employ a novel multiscale and multiphysics computational framework to perform predictive modeling of the pathological thrombus formation in the microvasculature using data from patients with COVID-19. This framework seamlessly integrates the key components in the process of blood clotting, including hemodynamics, transport of coagulation factors and coagulation kinetics, blood cell mechanics and adhesive dynamics, and thus allows us to quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19. Our simulation results show that among the coagulation factors considered, antithrombin and factor V play more prominent roles in promoting thrombosis. Our simulations also suggest that recruitment of WBCs to the endothelial cells exacerbates thrombogenesis and contributes to the blockage of the blood flow. Additionally, we show that the recent identification of flowing blood cell clusters could be a result of detachment of WBCs from thrombogenic sites, which may serve as a nidus for new clot formation. These findings point to potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19. Altogether, our computational framework provides a powerful tool for quantitative understanding of the mechanism of pathological thrombus formation and offers insights into new therapeutic approaches for treating COVID-19 associated thrombosis.

Emerging clinical evidence suggests that thrombosis in the microvasculature of patients with Coronavirus disease 2019 (COVID-19) plays an essential role in dictating the disease progression. We employ a novel multiphysics and multiscale computational framework to investigate the underlying mechanism of the pathological formation of microthrombi and circulating cell clusters in COVID-19. We quantify the contributions of many prothrombotic factors reported in the literature, such as stasis, the derangement in blood coagulation factor levels and activities, inflammatory responses of endothelial cells and leukocytes to the microthrombus formation in COVID-19, through which we identify the potential targets that should be further evaluated, and prioritized in the anti-thrombotic treatment of patients with COVID-19.

Recent Advances in Computational Modeling of Biomechanics and Biorheology of Red Blood Cells in Diabetes

Diabetes mellitus, a metabolic disease characterized by chronically elevated blood glucose levels, affects about 29 million Americans and more than 422 million adults all over the world. Particularly, type 2 diabetes mellitus (T2DM) accounts for 90-95% of the cases of vascular disease and its prevalence is increasing due to the rising obesity rates in modern societies. Although multiple factors associated with diabetes, such as reduced red blood cell (RBC) deformability, enhanced RBC aggregation and adhesion to the endothelium, as well as elevated blood viscosity are thought to contribute to the hemodynamic impairment and vascular occlusion, clinical or experimental studies cannot directly quantify the contributions of these factors to the abnormal hematology in T2DM. Recently, computational modeling has been employed to dissect the impacts of the aberrant biomechanics of diabetic RBCs and their adverse effects on microcirculation. In this review, we summarize the recent advances in the developments and applications of computational models in investigating the abnormal properties of diabetic blood from the cellular level to the vascular level. We expect that this review will motivate and steer the development of new models in this area and shift the attention of the community from conventional laboratory studies to combined experimental and computational investigations, aiming to provide new inspirations for the development of advanced tools to improve our understanding of the pathogenesis and pathology of T2DM.

Computational investigation of blood cell transport in retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). MA leakage or rupture may precipitate local pathology in the surrounding neural retina that impacts visual function. Thrombosis in MAs may affect their turnover time, an indicator associated with visual and anatomic outcomes in the diabetic eyes. In this work, we perform computational modeling of blood flow in microchannels containing various MAs to investigate the pathologies of MAs in DR. The particle-based model employed in this study can explicitly represent red blood cells (RBCs) and platelets as well as their interaction in the blood flow, a process that is very difficult to observe in vivo. Our simulations illustrate that while the main blood flow from the parent vessels can perfuse the entire lumen of MAs with small body-to-neck ratio (BNR), it can only perfuse part of the lumen in MAs with large BNR, particularly at a low hematocrit level, leading to possible hypoxic conditions inside MAs. We also quantify the impacts of the size of MAs, blood flow velocity, hematocrit and RBC stiffness and adhesion on the likelihood of platelets entering MAs as well as their residence time inside, two factors that are thought to be associated with thrombus formation in MAs. Our results show that enlarged MA size, increased blood velocity and hematocrit in the parent vessel of MAs as well as the RBC-RBC adhesion promote the migration of platelets into MAs and also prolong their residence time, thereby increasing the propensity of thrombosis within MAs. Overall, our work suggests that computational simulations using particle-based models can help to understand the microvascular pathology pertaining to MAs in DR and provide insights to stimulate and steer new experimental and computational studies in this area.

How the spleen reshapes and retains young and old red blood cells: A computational investigation

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell’s (RBC) life cycle. Although much progress has been made to elucidate the critical biological processes involved in the maturation of young RBCs (reticulocytes) as well as removal of senescent RBCs in the spleen, the underlying mechanisms driving these processes are still obscure. Herein, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan and investigate the role of the spleen in facilitating the maturation of reticulocytes and in clearing the senescent RBCs. Our simulations reveal that at the beginning of the RBC life cycle, intracellular non-deformable particles in reticulocytes can be biomechanically expelled from the cell upon passage through IES, an insightful explanation of why this peculiar “pitting” process is spleen-specific. Our results also show that immature RBCs shed surface area by releasing vesicles after crossing IES and progressively acquire the biconcave shape of mature RBCs. These findings likely explain why RBCs from splenectomized patients are significantly larger than those from nonsplenectomized subjects. Finally, we show that at the end of their life span, senescent RBCs are not only retained by IES due to reduced deformability but also become susceptible to mechanical lysis under shear stress. This finding supports the recent hypothesis that transformation into a hemolyzed ghost is a prerequisite for phagocytosis of senescent RBCs. Altogether, our computational investigation illustrates critical biological processes in the spleen that cannot be observed in vivo or in vitro and offer insights into the role of the spleen in the RBC physiology.

The spleen, the largest secondary lymphoid organ in humans, not only fulfils a broad range of immune functions, but also plays an important role in red blood cell (RBC) life cycle. In this study, we perform a computational study to simulate the passage of RBCs through interendothelial slits (IES) in the spleen at different stages of their lifespan, a critical biological process that cannot be observed in humans. Our simulation results illustrate a specific role of spleen in shaping young RBCs, which points to a probable missing step in current in vitro RBC culture protocols that fail to generate a majority of typical biconcave RBCs. Our results also reveal that intra-splenic mechanical constraints likely contribute to the final clearance and elimination of aged RBCs. Altogether, we demonstrate that our computational model can provide mechanistic rationales for experimental studies, offer insights into the role of the spleen in the RBC physiology and help the optimization of in vitro RBC culture techniques.

Nanoscale Changes on RBC Membrane Induced by Storage and Ionizing Radiation: A Mini-Review

The storage lesions and the irradiation of blood cellular components for medical procedures in blood banks are events that may induce nanochanges in the membrane of red blood cells (RBCs). Alterations, such as the formation of pores and vesicles, reduce flexibility and compromise the overall erythrocyte integrity. This review discusses the alterations on erythrocytic lipid membrane bilayer through their characterization by confocal scanning microscopy, Raman, scanning electron microscopy, and atomic force microscopy techniques. The interrelated experimental results may address and shed light on the correlation of biomechanical and biochemical transformations induced in the membrane and cytoskeleton of stored and gamma-irradiated RBC. To highlight the main advantages of combining these experimental techniques simultaneously or sequentially, we discuss how those outcomes observed at micro- and nanoscale cell levels are useful as biomarkers of cell aging and storage damage.

Red Blood Cells: Tethering, Vesiculation, and Disease in Micro-Vascular Flow

The red blood cell has become implicated in the progression of a range of diseases; mechanisms by which red cells are involved appear to include the transport of inflammatory species via red cell-derived vesicles. We review this role of RBCs in diseases such as diabetes mellitus, sickle cell anemia, polycythemia vera, central retinal vein occlusion, Gaucher disease, atherosclerosis, and myeloproliferative neoplasms. We propose a possibly unifying, and novel, paradigm for the inducement of RBC vesiculation during vascular flow of red cells adhered to the vascular endothelium as well as to the red pulp of the spleen. Indeed, we review the evidence for this hypothesis that links physiological conditions favoring both vesiculation and enhanced RBC adhesion and demonstrate the veracity of this hypothesis by way of a specific example occurring in splenic flow which we argue has various renderings in a wide range of vascular flows, in particular microvascular flows. We provide a mechanistic basis for membrane loss and the formation of lysed red blood cells in the spleen that may mediate their turnover. Our detailed explanation for this example also makes clear what features of red cell deformability are involved in the vesiculation process and hence require quantification and a new form of quantitative indexing.

Kinetics of self-assembly of inclusions due to lipid membrane thickness interactions

Self-assembly of proteins on lipid membranes underlies many important processes in cell biology, such as, exo- and endo-cytosis, assembly of viruses, etc. An attractive force that can cause self-assembly is mediated by membrane thickness interactions between proteins. The free energy profile associated with this attractive force is a result of the overlap of thickness deformation fields around the proteins which can be calculated from the solution of a boundary value problem. Yet, the time scales over which two inclusions coalesce has not been explored, even though the evolution of particle concentrations on membranes has been modeled using phase-field approaches. In this paper we compute this time scale as a function of the initial distance between two inclusions by viewing their coalescence as a first passage time problem. The mean first passage time is computed using Langevin dynamics and a partial differential equation, and both methods are found to be in excellent agreement. Inclusions of three different shapes are studied and it is found that for two inclusions separated by about hundred nanometers the time to coalescence is hundreds of milliseconds irrespective of shape. An efficient computation of the interaction energy of inclusions is central to our work. We compute it using a finite difference technique and show that our results are in excellent agreement with those from a previously proposed semi-analytical method based on Fourier-Bessel series. The computational strategies described in this paper could potentially lead to efficient methods to explore the kinetics of self-assembly of proteins on lipid membranes.

Tethering, evagination, and vesiculation via cell-cell interactions in microvascular flow

Vesiculation is a ubiquitous process undergone by most cell types and serves a variety of vital cell functions; vesiculation from erythrocytes, in particular, is a well-known example and constitutes a self-protection mechanism against premature clearance, inter alia. Herein, we explore a paradigm that red blood cell derived vesicles may form within the microvascular, in intense shear flow, where cells become adhered to either other cells or the extracellular matrix, by forming tethers or an evagination. Adherence may be enhanced, or caused, by diseased states or chemical anomalies as are discussed herein. The mechanisms for such processes are detailed via numerical simulations that are patterned directly from video-recorded cell microflow within the splenic venous sinus (MacDonald et al. 1987), as included, e.g., as Supplementary Material. The mechanisms uncovered highlight the necessity of accounting for remodeling of the erythrocyte's membrane skeleton and, specifically, for the time scales associated with that process that is an integral part of cell deformation. In this way, the analysis provides pointed, and vital, insights into the notion of what the, often used phrase, cell deformability actually entails in a more holistic manner. The analysis also details what data are required to make further quantitative descriptions possible and suggests experimental pathways for acquiring such.

Valid Presumption of Shiga Toxin-Mediated Damage of Developing Erythrocytes in EHEC-Associated Hemolytic Uremic Syndrome

The global emergence of clinical diseases caused by enterohemorrhagic Escherichia coli (EHEC) is an issue of great concern. EHEC release Shiga toxins (Stxs) as their key virulence factors, and investigations on the cell-damaging mechanisms toward target cells are inevitable for the development of novel mitigation strategies. Stx-mediated hemolytic uremic syndrome (HUS), characterized by the triad of microangiopathic hemolytic anemia, thrombocytopenia, and acute renal injury, is the most severe outcome of an EHEC infection. Hemolytic anemia during HUS is defined as the loss of erythrocytes by mechanical disruption when passing through narrowed microvessels. The formation of thrombi in the microvasculature is considered an indirect effect of Stx-mediated injury mainly of the renal microvascular endothelial cells, resulting in obstructions of vessels. In this review, we summarize and discuss recent data providing evidence that HUS-associated hemolytic anemia may arise not only from intravascular rupture of erythrocytes, but also from the extravascular impairment of erythropoiesis, the development of red blood cells in the bone marrow, via direct Stx-mediated damage of maturing erythrocytes, leading to “non-hemolytic” anemia.

Modelling of Red Blood Cell Morphological and Deformability Changes during In-Vitro Storage

Storage lesion is a critical issue facing transfusion treatments, and it adversely affects the quality and viability of stored red blood cells (RBCs). RBC deformability is a key indicator of cell health. Deformability measurements of each RBC unit are a key challenge in transfusion medicine research and clinical haematology. In this paper, a numerical study, inspired from the previous research for RBC deformability and morphology predictions, is conducted for the first time, to investigate the deformability and morphology characteristics of RBCs undergoing storage lesion. This study investigates the evolution of the cell shape factor, elongation index and membrane spicule details, where applicable, of discocyte, echinocyte I, echinocyte II, echinocyte III and sphero-echinocyte morphologies during 42 days of in-vitro storage at 4 °C in saline-adenine-glucose-mannitol (SAGM). Computer simulations were performed to investigate the influence of storage lesion-induced membrane structural defects on cell deformability and its recoverability during optical tweezers stretching deformations. The predicted morphology and deformability indicate decreasing quality and viability of stored RBCs undergoing storage lesion. The loss of membrane structural integrity due to the storage lesion further degrades the cell deformability and recoverability during mechanical deformations. This numerical approach provides a potential framework to study the RBC deformation characteristics under varying pathophysiological conditions for better diagnostics and treatments.

Cytoskeleton Remodeling Induces Membrane Stiffness and Stability Changes of Maturing Reticulocytes

Reticulocytes, the precursors of erythrocytes, undergo drastic alterations in cell size, shape, and deformability during maturation. Experimental evidence suggests that young reticulocytes are stiffer and less stable than their mature counterparts; however, the underlying mechanism is yet to be fully understood. Here, we develop a coarse-grained molecular-dynamics reticulocyte membrane model to elucidate how the membrane structure of reticulocytes contributes to their particular biomechanical properties and pathogenesis in blood diseases. First, we show that the extended cytoskeleton in the reticulocyte membrane is responsible for its increased shear modulus. Subsequently, we quantify the effect of weakened cytoskeleton on the stiffness and stability of reticulocytes, via which we demonstrate that the extended cytoskeleton along with reduced cytoskeleton connectivity leads to the seeming paradox that reticulocytes are stiffer and less stable than the mature erythrocytes. Our simulation results also suggest that membrane budding and the consequent vesiculation of reticulocytes can occur independently of the endocytosis-exocytosis pathway, and thus, it may serve as an additional means of removing unwanted membrane proteins from reticulocytes. Finally, we find that membrane budding is exacerbated when the cohesion between the lipid bilayer and the cytoskeleton is compromised, which is in accord with the clinical observations that erythrocytes start shedding membrane surface at the reticulocyte stage in hereditary spherocytosis. Taken together, our results quantify the stiffness and stability change of reticulocytes during their maturation and provide, to our knowledge, new insights into the pathogenesis of hereditary spherocytosis and malaria.

A coarse-grained red blood cell membrane model to study stomatocyte-discocyte-echinocyte morphologies

An improved red blood cell (RBC) membrane model is developed based on the bilayer coupling model (BCM) to accurately predict the complete sequence of stomatocyte-discocyte-echinocyte (SDE) transformation of a RBC. The coarse-grained (CG)-RBC membrane model is proposed to predict the minimum energy configuration of the RBC from the competition between lipid-bilayer bending resistance and cytoskeletal shear resistance under given reference constraints. In addition to the conventional membrane surface area, cell volume and bilayer-leaflet-area-difference constraints, a new constraint: total-membrane-curvature is proposed in the model to better predict RBC shapes in agreement with experimental observations. A quantitative evaluation of several cellular measurements including length, thickness and shape factor, is performed for the first time, between CG-RBC model predicted and three-dimensional (3D) confocal microscopy imaging generated RBC shapes at equivalent reference constraints. The validated CG-RBC membrane model is then employed to investigate the effect of reduced cell volume and elastic length scale on SDE transformation, to evaluate the RBC deformability during SDE transformation, and to identify the most probable RBC cytoskeletal reference state. The CG-RBC membrane model can predict the SDE shape behaviour under diverse shape-transforming scenarios, in-vitro RBC storage, microvascular circulation and flow through microfluidic devices.

The fluid membrane determines mechanics of erythrocyte extracellular vesicles and is softened in hereditary spherocytosis

Extracellular vesicles (EVs) are widely studied regarding their role in cell-to-cell communication and disease, as well as for applications as biomarkers or drug delivery vehicles. EVs contain membrane and intraluminal proteins, affecting their structure and thereby likely their functioning. Here, we use atomic force microscopy for mechanical characterization of erythrocyte, or red blood cell (RBC), EVs from healthy individuals and from patients with hereditary spherocytosis (HS) due to ankyrin deficiency. While these EVs are packed with proteins, their response to indentation resembles that of fluid liposomes lacking proteins. The bending modulus of RBC EVs of healthy donors is ~15 k_bT, similar to the RBC membrane. Surprisingly, whereas RBCs become more rigid in HS, patient EVs have a significantly (~40%) lower bending modulus than donor EVs. These results shed light on the mechanism and effects of EV budding and might explain the reported increase in vesiculation of RBCs in HS patients.

Red blood cell disorders are often accompanied by increased release of extracellular vesicles (EVs), but their structural and mechanical properties are not fully understood. Here, the authors show that red blood cell EVs show liposome-like mechanical features and are softened in blood disorder patients.

Erythrocyte Aging, Protection via Vesiculation: An Analysis Methodology via Oscillatory Flow

We demonstrate that erythrocyte deformations, specifically of a type as occur in splenic flow (Zhu et al., 2017), and of the type that promote vesiculation can be caused by simple, yet tailored, oscillatory shear flow. We show that such oscillatory shear flow provides an ideal environment to explore a wide variety of metabolic and biochemical effects that promote erythrocyte vesiculation. Deformation details, typical of splenic flow, such as in-folding and implications for membrane/skeleton interaction are demonstrated and quantitatively analyzed. We introduce a theoretical, essentially analytical, vesiculation model that directly couples to our more complex numerical, multilevel, model that clearly delineates various fundamental elements, i.e., sub-processes, that are involved and mediate the vesiculation process. This analytical model highlights particulary important vesiculation precursors such as areas of membrane/skeleton disruptions that trigger the vesiculation process. We demonstrate, using flow cytometry, that the deformations we experimentally induce on cells, and numerically simulate, do not induce lethal forms of cell damage but do induce vesiculation as theoretically forecasted. This, we demonstrate, provides a direct link to cell membrane/skeletal damage such as is associated with metabolic and aging damage. An additional noteworthy feature of this approach is the avoidance of artificial devices, e.g., micro-fluidic chambers, in which deformations and their time scales are often unrepresentative of physiological processes such as splenic flow.

Plasma Membrane Lipid Domains as Platforms for Vesicle Biogenesis and Shedding?

Extracellular vesicles (EVs) contribute to several pathophysiological processes and appear as emerging targets for disease diagnosis and therapy. However, successful translation from bench to bedside requires deeper understanding of EVs, in particular their diversity, composition, biogenesis and shedding mechanisms. In this review, we focus on plasma membrane-derived microvesicles (MVs), far less appreciated than exosomes. We integrate documented mechanisms involved in MV biogenesis and shedding, focusing on the red blood cell as a model. We then provide a perspective for the relevance of plasma membrane lipid composition and biophysical properties in microvesiculation on red blood cells but also platelets, immune and nervous cells as well as tumor cells. Although only a few data are available in this respect, most of them appear to converge to the idea that modulation of plasma membrane lipid content, transversal asymmetry and lateral heterogeneity in lipid domains may play a significant role in the vesiculation process. We suggest that lipid domains may represent platforms for inclusion/exclusion of membrane lipids and proteins into MVs and that MVs could originate from distinct domains during physiological processes and disease evolution.

Synergistic Integration of Laboratory and Numerical Approaches in Studies of the Biomechanics of Diseased Red Blood Cells

In red blood cell (RBC) disorders, such as sickle cell disease, hereditary spherocytosis, and diabetes, alterations to the size and shape of RBCs due to either mutations of RBC proteins or changes to the extracellular environment, lead to compromised cell deformability, impaired cell stability, and increased propensity to aggregate. Numerous laboratory approaches have been implemented to elucidate the pathogenesis of RBC disorders. Concurrently, computational RBC models have been developed to simulate the dynamics of RBCs under physiological and pathological conditions. In this work, we review recent laboratory and computational studies of disordered RBCs. Distinguished from previous reviews, we emphasize how experimental techniques and computational modeling can be synergically integrated to improve the understanding of the pathophysiology of hematological disorders.

Biophysics of membrane curvature remodeling at molecular and mesoscopic lengthscales

At the micron scale, where cell organelles display an amazing complexity in their shape and organization, the physical properties of a biological membrane can be better-understood using continuum models subject to thermal (stochastic) undulations. Yet, the chief orchestrators of these complex and intriguing shapes are a specialized class of membrane associating often peripheral proteins called curvature remodeling proteins (CRPs) that operate at the molecular level through specific protein-lipid interactions. We review multiscale methodologies to model these systems at the molecular as well as at the mesoscopic and cellular scales, and also present a free energy perspective of membrane remodeling through the organization and assembly of CRPs. We discuss the morphological space of nearly planar to highly curved membranes, methods to include thermal fluctuations, and review studies that model such proteins as curvature fields to describe the emergent curved morphologies. We also discuss several mesoscale models applied to a variety of cellular processes, where the phenomenological parameters (such as curvature field strength) are often mapped to models of real systems based on molecular simulations. Much insight can be gained from the calculation of free energies of membranes states with protein fields, which enable accurate mapping of the state and parameter values at which the membrane undergoes morphological transformations such as vesiculation or tubulation. By tuning the strength, anisotropy, and spatial organization of the curvature-field, one can generate a rich array of membrane morphologies that are highly relevant to shapes of several cellular organelles. We review applications of these models to budding of vesicles commonly seen in cellular signaling and trafficking processes such as clathrin mediated endocytosis, sorting by the ESCRT protein complexes, and cellular exocytosis regulated by the exocyst complex. We discuss future prospects where such models can be combined with other models for cytoskeletal assembly, and discuss their role in understanding the effects of cell membrane tension and the mechanics of the extracellular microenvironment on cellular processes.

Nanoscale membrane architecture of healthy and pathological red blood cells

Red blood cells feature remarkable mechanical properties while navigating through microcirculation vessels and during spleen filtration. An unusual combination of plasma membrane and cytoskeleton physical properties allows red blood cells to undergo extensive deformation. Here we used atomic force microscopy multiparametric imaging to probe how cellular organization influences nanoscale and global mechanical properties of cells in both physiological and pathological conditions. Our data obtained in native conditions confirmed that, compared to healthy cells, cells from patients with hereditary spherocytosis are stiffer. Through vertical segmentation of the cell elasticity, we found that healthy and pathological cells display nanoscale architecture with an increasing stiffness along the direction of the applied force. By decoupling the mechanical response of the plasma membrane from its underlying cytoskeleton, we find that both components show altered properties in pathological conditions. Nanoscale multiparametric imaging also revealed lipid domains that exhibit differential mechanical properties than the bulk membrane in both healthy and pathological conditions. Thanks to correlated AFM-fluorescence imaging, we identified submicrometric sphingomyelin-enriched lipid domains of variable stiffness at the red blood cell surface. Our experiments provide novel insights into the interplay between nanoscale organization of red blood cell plasma membrane and their nanomechanical properties. Overall, this work contributes to a better understanding of the complex relationship between cellular nanoscale organization, cellular nanomechanics and how this 3D organization is altered in pathological conditions.

Change of plasmalogen content of red blood cells in myocardial hypoxia and acidosis

BACKGROUND

The contribution of hypoxic conditions to the chemical composition of membranes is not completely established. Plasmalogens, containing an alkenyl group with aldehydogenic ether linkage, are significant components of membrane lipids and their level can change in oxygen deficiency.

METHODS

Analysis of plasmalogens in red blood cells was performed in patients (n = 17) with coronary heart disease, stable angina (functional class II-III) and coronary atherosclerosis. The control group consisted of 17 healthy volunteers. In addition, isolated blood samples of seven healthy volunteers were analysed before and after 180 min incubation at 37 °C. Fatty acid ethyl esters and diethyl acetals of fatty aldehydes, obtained during sample preparation from red blood cells, were analysed by capillary gas-liquid chromatography. Quantitative assessment of the change of the plasmalogen levels was evaluated as change of the share of fatty aldehyde diethyl acetals in the total sum of fatty aldehyde diethyl acetals and fatty acid ethyl esters.

RESULTS

In comparison with the healthy volunteers, an increase in plasmalogen content of red blood cells and a reduction of the pH of the blood plasma in the group of patients with coronary heart disease were detected. In experimental hypoxia, there was an increase in the plasmalogen content of the red blood cells and a plasma pH decrease in all samples subjected to the incubation.

CONCLUSIONS

The results indicate changes in the physicochemical properties of the cell membrane in hypoxia. One of the most likely reasons of the increase of plasmalogen content in the membranes may be a more significant increase in activity of calcium-dependent phospholipases in comparison with the activity of calcium-independent plasmalogen phospholipases.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

We review recent advances in multiscale modeling of the biomechanical characteristics of red blood cells (RBCs) in hematological diseases, and their relevance to the structure and dynamics of defective RBCs. We highlight examples of successful simulations of blood disorders including malaria and other hereditary disorders, such as sickle-cell anemia, spherocytosis, and elliptocytosis.

Prospects for Human Erythrocyte Skeleton-Bilayer Dissociation during Splenic Flow

Prospects of vesiculation occurring during splenic flow of erythrocytes are addressed via model simulations of RBC flow through the venous slits of the human spleen. Our model is multiscale and contains a thermally activated rate-dependent description of the entropic elasticity of the RBC spectrin cytoskeleton, including domain unfolding/refolding. Our model also includes detail of the skeleton attachment to the fluidlike lipid bilayer membrane, including a specific accounting for the expansion/contraction of the skeleton that may occur via anchor protein diffusive motion, that is, band 3 and glycophorin, through the membrane. This ability allows us to follow the change in anchor density and thereby the strength of the skeleton/membrane attachment. We define a negative pressure between the skeleton/membrane connection that promotes separation; critical levels for this are estimated using published data on the work of adhesion of this connection. By following the maximum range of negative pressure, along with the observed slight decrease in skeletal density, we conclude that there must be biochemical influences that probably include binding of degraded hemoglobin, among other things, that significantly reduce effective attachment density. These findings are consistent with reported trends in vesiculation that are believed to occur in cases of various hereditary anemias and during blood storage. Our findings also suggest pathways for further study of erythrocyte vesiculation that point to the criticality of understanding the biochemical phenomena involved with cytoskeleton/membrane attachment.

OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution

We present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment-simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles-using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation. The algorithm enables the construction of the key spatial searching data structure in our code, i.e., a lattice-free cell list, with a time and space cost linearly proportional to the number of particles in the system. The position and the shape of the cells also adapt automatically to the local density and curvature. The code implements OpenMP parallelization and scales to hundreds of hardware threads. It outperforms a legacy simulator by almost an order of magnitude in time-to-solution and >40 times in problem size, thus providing, to our knowledge, a new platform for probing the biomechanics of red blood cells.

Modeling of the axon membrane skeleton structure and implications for its mechanical properties

Super-resolution microscopy recently revealed that, unlike the soma and dendrites, the axon membrane skeleton is structured as a series of actin rings connected by spectrin filaments that are held under tension. Currently, the structure-function relationship of the axonal structure is unclear. Here, we used atomic force microscopy (AFM) to show that the stiffness of the axon plasma membrane is significantly higher than the stiffnesses of dendrites and somata. To examine whether the structure of the axon plasma membrane determines its overall stiffness, we introduced a coarse-grain molecular dynamics model of the axon membrane skeleton that reproduces the structure identified by super-resolution microscopy. Our proposed computational model accurately simulates the median value of the Young’s modulus of the axon plasma membrane determined by atomic force microscopy. It also predicts that because the spectrin filaments are under entropic tension, the thermal random motion of the voltage-gated sodium channels (Na_v), which are bound to ankyrin particles, a critical axonal protein, is reduced compared to the thermal motion when spectrin filaments are held at equilibrium. Lastly, our model predicts that because spectrin filaments are under tension, any axonal injuries that lacerate spectrin filaments will likely lead to a permanent disruption of the membrane skeleton due to the inability of spectrin filaments to spontaneously form their initial under-tension configuration.

Super-resolution microscopy has suggested that the actin cytoskeleton structure differ between various neuronal subcompartments. To determine the possible implication of the differing actin cytoskeleton structure, we determined the stiffness of the plasma membrane of neuronal subcompartments using atomic force microscopy (AFM). We found that axons are almost ~6 fold stiffer than the soma and ~2 fold stiffer than dendrites. By using a particle-based model for the surface membrane skeleton of the axon that comprises actin rings connected with spring filaments to represent the axonal structure, we show that regions neighboring actin rings are stiffer than areas between these rings. In these in between sub-regions, the spectrin filaments determine stiffness. Our modeling also shows that because the spectrin filaments are under tension, the thermal jitter of the actin-associated ankyrin particles, connected to the middle area of spectrin filaments, is minimal. As a result, we propose that the sodium channels bound to ankyrin particles will maintain an ordered distribution along the axon. We also predict that laceration of the spectrin filaments due to injury will cause a permanent damage to the axon since spontaneous repair of the spectrin network is not possible as spectrin filaments are under entropic tension.

MD/DPD Multiscale Framework for Predicting Morphology and Stresses of Red Blood Cells in Health and Disease

Healthy red blood cells (RBCs) have remarkable deformability, squeezing through narrow capillaries as small as 3 microns in diameter without any damage. However, in many hematological disorders the spectrin network and lipid bilayer of diseased RBCs may be significantly altered, leading to impaired functionality including loss of deformability. We employ a two-component whole-cell multiscale model to quantify the biomechanical characteristics of the healthy and diseased RBCs, including Plasmodium falciparum-infected RBCs (Pf-RBCs) and defective RBCs in hereditary disorders, such as spherocytosis and elliptocytosis. In particular, we develop a two-step multiscale framework based on coarse-grained molecular dynamics (CGMD) and dissipative particle dynamics (DPD) to predict the static and dynamic responses of RBCs subject to tensile forcing, using experimental information only on the structural defects in the lipid bilayer, cytoskeleton, and their interaction. We first employ CGMD on a small RBC patch to compute the shear modulus, bending stiffness, and network parameters, which are subsequently used as input to a whole-cell DPD model to predict the RBC shape and corresponding stress field. For Pf-RBCs at trophozoite and schizont stages, the presence of cytoadherent knobs elevates the shear response in the lipid bilayer and stiffens the RBC membrane. For RBCs in spherocytosis and elliptocytosis, the bilayer-cytoskeleton interaction is weakened, resulting in substantial increase of the tensile stress in the lipid bilayer. Furthermore, we investigate the transient behavior of stretching deformation and shape relaxation of the normal and defective RBCs. Different from the normal RBCs possessing high elasticity, our simulations reveal that the defective RBCs respond irreversibly, i.e., they lose their ability to recover the normal biconcave shape in successive loading cycles of stretching and relaxation. Our findings provide fundamental insights into the microstructure and biomechanics of RBCs, and demonstrate that the two-step multiscale framework presented here can be used effectively for in silico studies of hematological disorders based on first principles and patient-specific experimental input at the protein level.

Membrane Compartmentalization Reducing the Mobility of Lipids and Proteins within a Model Plasma Membrane

The cytoskeleton underlying cell membranes may influence the dynamic organization of proteins and lipids within the bilayer by immobilizing certain transmembrane (TM) proteins and forming corrals within the membrane. Here, we present coarse-grained resolution simulations of a biologically realistic membrane model of asymmetrically organized lipids and TM proteins. We determine the effects of a model of cytoskeletal immobilization of selected membrane proteins using long time scale coarse-grained molecular dynamics simulations. By introducing compartments with varying degrees of restraints within the membrane models, we are able to reveal how compartmentalization caused by cytoskeletal immobilization leads to reduced and anomalous diffusional mobility of both proteins and lipids. This in turn results in a reduced rate of protein dimerization within the membrane and of hopping of membrane proteins between compartments. These simulations provide a molecular realization of hierarchical models often invoked to explain single-molecule imaging studies of membrane proteins.

Modeling of band-3 protein diffusion in the normal and defective red blood cell membrane

We employ a two-component red blood cell (RBC) membrane model to simulate lateral diffusion of band-3 proteins in the normal RBC and in the RBC with defective membrane proteins. The defects reduce the connectivity between the lipid bilayer and the membrane skeleton (vertical connectivity), or the connectivity of the membrane skeleton itself (horizontal connectivity), and are associated with the blood disorders of hereditary spherocytosis (HS) and hereditary elliptocytosis (HE) respectively. Initially, we demonstrate that the cytoskeleton limits band-3 lateral mobility by measuring the band-3 macroscopic diffusion coefficients in the normal RBC membrane and in a lipid bilayer without the cytoskeleton. Then, we study band-3 diffusion in the defective RBC membrane and quantify the relation between band-3 diffusion coefficients and percentage of protein defects in HE RBCs. In addition, we illustrate that at low spectrin network connectivity (horizontal connectivity) band-3 subdiffusion can be approximated as anomalous diffusion, while at high horizontal connectivity band-3 diffusion is characterized as confined diffusion. Our simulations show that the band-3 anomalous diffusion exponent depends on the percentage of protein defects in the membrane cytoskeleton. We also confirm that the introduction of attraction between the lipid bilayer and the spectrin network reduces band-3 diffusion, but we show that this reduction is lower than predicted by the percolation theory. Furthermore, we predict that the attractive force between the spectrin filament and the lipid bilayer is at least 20 times smaller than the binding forces at band-3 and glycophorin C, the two major membrane binding sites. Finally, we explore diffusion of band-3 particles in the RBC membrane with defects related to vertical connectivity. We demonstrate that in this case band-3 diffusion can be approximated as confined diffusion for all attraction levels between the spectrin network and the lipid bilayer. By comparing the diffusion coefficients measured in horizontal vs. vertical defects, we conclude that band-3 mobility is primarily controlled by the horizontal connectivity.