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.

Red blood cell motion and deformation in a curved microvessel

The flow of cells through curved vessels is often encountered in various biomedical and bioengineering applications, such as red blood cells (RBCs) passing through the curved arteries in circulation, and cells sorting through a shear-induced migration in a curved channels. Most of past numerical studies focused on the cell deformation in small straight microvessels, or on the flow pattern in large curved vessels without considering the cell deformation. However, there have been few attempts to study the cell deformation and the associated flow pattern in a curved microvessel. In this work, a particle-based method, smoothed dissipative particle dynamics (SDPD), is used to simulate the motion and deformation of a RBC in a curved microvessel of diameter comparable to the RBC diameter. The emphasis is on the effects of the curvature, the type and the size of the curved microvessel on the RBC deformation and the flow pattern. The simulation results show that a small curved shape of the microvessel has negligible effect on the RBC behavior and the flow pattern which are similar to those in a straight microvessel. When the microvessel is high in curvature, the secondary flow comes into being with a pair of Dean vortices, and the velocity profile of the primary flow is skewed toward the inner wall of the microvessel. The RBC also loses the axisymmetric deformation, and it is stretched first and then shrinks when passing through the curved part of the microvessel with the large curvature. It is also found that a pair of Dean vortices arise only under the condition of De>1 (De is the Dean number, a ratio of centrifugal to viscous competition). The Dean vortices are more easily observed in the larger or more curved microvessels. Finally, it is observed that the velocity profile of primary flow is skewed toward the inner wall of curved microvessel, i.e., the fluid close to the inner wall flows faster than that close to the outer wall. This is contrary to the common sense in large curved vessels. This velocity skewness was found to depend on the curvature of the microvessel, as well as the viscous and inertial forces.

Control performance of paper-based blood analysis devices through paper structure design

In this work, we investigated the influence of paper structure on the performance of paper-based analytical devices that are used for blood analysis. The question that we aimed to answer is how the fiber type (i.e., softwood and hardwood fibers) influences the fiber network structure of the paper, which affects the transport of red blood cells (RBCs) in paper. In the experimental design, we isolated the influence of fiber types on the paper structure from all other possible influencing factors by removing the fines from the pulps and not using any additives. Mercury porosimetry was employed to characterize the pore structures of the paper sheets. The results show that papers with a low basis weight that are made with short hardwood fibers have a higher porosity (i.e., void fraction) and simpler pore structures compared with papers made with long softwood fibers. RBC transport in paper carried by saline solution was investigated in two modes: lateral chromatographic elution and vertical flow-through. The results showed that the complexity of the paper's internal pore structure has a dominant influence on the transport of RBCs in paper. Hardwood fiber sheets with a low basis weight have a simple internal pore structure and allow for the easy transport of RBCs. Blood-typing sensors built with low basis weight hardwood fibers deliver high-clarity assays. Softwood fiber papers are found to have a more complex pore structure, which makes RBC transport more difficult, leading to blood-typing results of low clarity. This study provides the principle of paper sheet design for paper-based blood analysis sensors.

Micro-scale blood particulate dynamics using a non-uniform rational B-spline-based isogeometric analysis

The current research presents a novel method in which blood particulates - biconcave red blood cells (RBCs) and spherical cells are modeled using isogeometric analysis, specifically Non-Uniform Rational B-Splines (NURBS) in 3-D. The use of NURBS ensures that even with a coarse representation, the geometry of the blood particulates maintains an accurate description when subjected to large deformations. The fundamental advantage of this method is the coupling of the geometrical description and the stress analysis of the cell membrane into a single, unified framework. Details on the modeling approach, implementation of boundary conditions and the membrane mechanics analysis using isogeometric modeling are presented, along with validation cases for spherical and biconcave cells. Using NURBS - based isogeometric analysis, the behavior of individual cells in fluid flow is presented and analyzed in different flow regimes using as few as 176 elements for a spherical cell and 220 elements for a biconcave RBC. This work provides a framework for modeling a large number of 3-D deformable biological cells, each with its own geometric description and membrane properties. To the best knowledge of the authors, this is the first application of the NURBS - based isogeometric analysis to model and simulate blood particulates in flow in 3D.

Equivalent mechanical properties of biological membranes from lattice homogenization

The goal of this manuscript is to set up a novel methodology for the calculation of the effective mechanical properties of biological membranes viewed as repetitive networks of elastic filaments, based on the discrete asymptotic homogenization method. We will show that for some lattice configurations, flexional effects due to internal structure mechanisms at the unit cell scale lead to additional flexional effects at the continuum scale, accounted for by an internal length associated to a micropolar behavior. Thereby, a systematic methodology is established, allowing the prediction of the overall mechanical properties of biological membranes for a given network topology, as closed form expressions of the geometrical and mechanical micro-parameters. The peptidoglycan and the erythrocyte have been analyzed using this methodology, and their effective moduli are calculated and recorded versus the geometrical and mechanical lattice parameters. A classification of lattices with respect to the choice of the equivalent continuum model is proposed: The Cauchy continuum and a micropolar continuum are adopted as two possible effective medium, for a given beam model. The relative ratio of the characteristic length of the micropolar continuum to the unit cell size determines the relevant choice of the equivalent medium. In most cases, the Cauchy continuum is sufficient to model membranes in most of their configurations. The peptidoglycan network may exhibit a re-entrant hexagonal lattice, for which micropolar effects become important. This is attested by the characteristic length becoming larger than the beam length for such configurations. The homogenized moduli give accurate results for both membranes, as revealed by comparison with experimental measurements or simulation results from the literature at the network scale. A first insight into the nonlinear mechanical behavior of the hexagonal and triangular networks is lastly investigated using a perturbative method.