Dataset Information

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells.

ABSTRACT: We compare the predictive capability of two mathematical models for red blood cells (RBCs) focusing on blood flow in capillaries and arterioles. Both RBC models as well as their corresponding blood flows are based on the dissipative particle dynamics (DPD) method, a coarse-grained molecular dynamics approach. The first model employs a multiscale description of the RBC (MS-RBC), with its membrane represented by hundreds or even thousands of DPD-particles connected by springs into a triangular network in combination with out-of-plane elastic bending resistance. Extra dissipation within the network accounts for membrane viscosity, while the characteristic biconcave RBC shape is achieved by imposition of constraints for constant membrane area and constant cell volume. The second model is based on a low-dimensional description (LD-RBC) constructed as a closed torus-like ring of only 10 large DPD colloidal particles. They are connected into a ring by worm-like chain (WLC) springs combined with bending resistance. The LD-RBC model can be fitted to represent the entire range of nonlinear elastic deformations as measured by optical-tweezers for healthy and for infected RBCs in malaria. MS-RBCs suspensions model the dynamics and rheology of blood flow accurately for any vessel size but this approach is computationally expensive for vessel diameters above 100?m. Surprisingly, the much more economical suspensions of LD-RBCs also capture the blood flow dynamics and rheology accurately except for small-size vessels comparable to RBC diameter. In particular, the LD-RBC suspensions are shown to properly capture the experimental data for the apparent viscosity of blood and its cell-free layer (CFL) in tube flow. Taken together, these findings suggest a hierarchical approach in modeling blood flow in the arterial tree, whereby the MS-RBC model should be employed for capillaries and arterioles below 100?m, the LD-RBC model for arterioles, and the continuum description for arteries.

SUBMITTER: Pan W

PROVIDER: S-EPMC3149761 | biostudies-literature | 2011 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells.

Pan Wenxiao W Fedosov Dmitry A DA Caswell Bruce B Karniadakis George Em GE

Microvascular research 20110527 2

We compare the predictive capability of two mathematical models for red blood cells (RBCs) focusing on blood flow in capillaries and arterioles. Both RBC models as well as their corresponding blood flows are based on the dissipative particle dynamics (DPD) method, a coarse-grained molecular dynamics approach. The first model employs a multiscale description of the RBC (MS-RBC), with its membrane represented by hundreds or even thousands of DPD-particles connected by springs into a triangular net ...[more]

PMID: 21640731

Similar Datasets

Project description: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.

Project description:In this work, we compared the dynamics of motion in a linear shear flow of individual red blood cells (RBCs) from healthy and pathological donors (Sickle Cell Disease (SCD) or Sickle Cell-β-thalassemia) and of low and high densities, in a suspending medium of higher viscosity. In these conditions, at lower shear rates, biconcave discocyte-shaped RBCs present an unsteady flip-flopping motion, where the cell axis of symmetry rotates in the shear plane, rocking to and fro between an orbital angle ±ϕ observed when the cell is on its edge. We show that the evolution of ϕ depends solely on RBC density for healthy RBCs, with denser RBCs displaying lower ϕ values than the lighter ones. Typically, at a shear stress of 0.08 Pa, ϕ has values of 82 and 72° for RBCs with average densities of 1.097 and 1.115, respectively. Surprisingly, we show that SCD RBCs display the same ϕ-evolution as healthy RBCs of same density, showing that the flip-flopping behavior is unaffected by the SCD pathology. When the shear stress is increased further (above 0.1 Pa), healthy RBCs start going through a transition to a fluid-like motion, called tank-treading, where the RBC has a quasi-constant orientation relatively to the flow and the membrane rotates around the center of mass of the cell. This transition occurs at higher shear stresses (above 0.2 Pa) for denser cells. This shift toward higher stresses is even more remarkable in the case of SCD RBCs, showing that the transition to the tank-treading regime is highly dependent on the SCD pathology. Indeed, at a shear stress of 0.2 Pa, for RBCs with a density of 1.097, 100% of healthy RBCs have transited to the tank-treading regime vs. less than 50% SCD RBCs. We correlate the observed differences in dynamics to the alterations of RBC mechanical properties with regard to density and SCD pathology reported in the literature. Our results suggest that it might be possible to develop simple non-invasive assays for diagnosis purpose based on the RBC motion in shear flow and relying on this millifluidic approach.

Dataset Information

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells.

Publications

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets