Haematocrit and Shear Rate Modulate Local Cell-free Layer Thickness and Platelet Margination in Blood Flow Along a Sinusoidal Wall

This study utilizes 3D simulations to demonstrate that hematocrit levels and local shear rates modulate cell-free layer thickness and platelet margination along sinusoidal walls, revealing that low hematocrit promotes platelet accumulation at crests while higher hematocrit leads to uniform distribution, thereby offering mechanistic insights into the coupled evolution of vessel topography and thrombosis.

The Gist

Imagine your bloodstream as a busy, crowded highway. The Red Blood Cells (RBCs) are the massive delivery trucks that take up most of the road, while the Platelets are tiny, nimble motorcycles trying to zip through the traffic to get to the side of the road (the vessel wall) where they might need to stop and patch a pothole (a clot).

This paper is a high-tech simulation of what happens when that highway isn't perfectly flat. Instead, imagine the road has a wavy, sinusoidal pattern—like a series of gentle hills and valleys. This wave pattern represents the bumps and ridges formed by clots that have already started to grow on the vessel wall.

Here is what the researchers discovered, broken down into simple concepts:

1. The "Empty Lane" (The Cell-Free Layer)

In a normal, straight highway, the big trucks (RBCs) tend to stay in the middle of the road because they bounce off the sides. This leaves a thin, empty lane right next to the wall where only the tiny motorcycles (platelets) can drive. Scientists call this the Cell-Free Layer (CFL).

• The Discovery: When the road is wavy, this empty lane changes size depending on the shape of the road.
  • In the Valleys (the dips), the empty lane gets wider. The big trucks avoid the dip, leaving a huge gap.
  • On the Crests (the peaks), the empty lane gets narrower. The trucks get squeezed closer to the wall.

2. The "Traffic Density" (Hematocrit)

The researchers changed how crowded the highway was (this is called Hematocrit).

• Low Traffic (Low Hematocrit): There are fewer trucks. The motorcycles can easily find the empty lane. However, because the empty lane is so wide in the valleys and so narrow on the crests, the motorcycles get stuck on the crests. They pile up on the peaks, making the "bumps" grow taller and sharper.
• High Traffic (High Hematocrit): The highway is packed with trucks. The trucks push each other around so much that they fill up the valleys too. The empty lane becomes very thin everywhere. Now, the motorcycles can't just hide on the peaks; they get pushed into the valleys as well. The result? The clots grow more evenly, creating a smoother, flatter surface rather than sharp spikes.

The Analogy: Think of it like people trying to get to a concert stage.

• If the crowd is thin (Low Traffic), people can easily find the front row, but they only stand on the raised platforms (crests), leaving the floor (valleys) empty.
• If the crowd is packed tight (High Traffic), people are pushed everywhere, filling both the platforms and the floor, creating a uniform wall of fans.

3. The "Wind Speed" (Shear Rate)

The speed of the water flowing past the wall (Shear Rate) is different on the peaks and in the dips.

• On the Crests: The water rushes fast. This is like a strong wind. In real biology, this fast wind triggers a specific "emergency signal" (using a protein called vWF) that tells platelets to grab on.
• In the Valleys: The water slows down significantly. This is like a calm breeze. Here, a different, slower "signal" (using collagen) takes over to help platelets stick.

The Big Picture: Because the road is wavy, you have two different "weather zones" right next to each other. One zone is windy (crests), and one is calm (valleys). This means the body uses different biological tools to build the clot depending on exactly where on the wave the platelet lands.

Why Does This Matter?

The study explains why clots sometimes look like jagged spikes and other times look like smooth mounds.

• If your blood is less crowded, clots will grow tall and spiky on the peaks.
• If your blood is very crowded, clots will grow evenly and smoothly.

The Takeaway for Medicine:
Understanding these "traffic patterns" helps doctors design better drugs. If we know exactly where the platelets are most likely to get stuck (the peaks vs. the valleys) and what kind of "wind" they are facing, we can create targeted medicines that only activate when they hit that specific spot. It's like sending a delivery drone that only drops its package when it detects a specific wind speed, ensuring the medicine works exactly where it's needed without causing side effects elsewhere.

Technical Summary

1. Problem Statement

Platelet aggregation is a critical step in hemostasis and thrombosis. As platelets adhere to damaged endothelium, they form aggregates that alter the local vessel wall topography, creating a feedback loop where surface heterogeneities influence local hemodynamics, which in turn drives further aggregate growth.

• The Gap: While the mechanisms of platelet margination (lateral migration to the wall) are well-studied in straight, smooth channels, it is poorly understood how surface heterogeneities (specifically the wavy geometry of existing aggregates) alter local hemodynamics, Cell-Free Layer (CFL) thickness, and platelet capture sites.
• Objective: To investigate how hematocrit (Ht) and shear rate (via the Capillary number, Ca) modulate RBC dynamics, CFL thickness, and platelet margination near a sinusoidal wall that mimics flow-aligned platelet aggregates.

2. Methodology

The study employs high-fidelity, three-dimensional numerical simulations of particulate blood flow.

• Numerical Framework:
  • Solver: Lattice-Boltzmann Method (LBM) with D3Q19 lattice and BGK collision operator to solve Navier-Stokes equations.
  • Coupling: Immersed-Boundary Method (IBM) to couple the fluid (Eulerian) with blood cells (Lagrangian).
  • Cell Modeling:
    • RBCs: Modeled as deformable membranes using a finite-element approach (1280 triangular elements). Elasticity includes in-plane shear/dilation (Skalak's law) and bending resistance, with constraints for constant surface area and volume.
    • Platelets: Modeled as nearly rigid spheres (180 triangular elements, radius 1 µm).
• Geometry:
  • A straight channel with a flat top wall and a sinusoidal bottom wall.
  • The sinusoid mimics streamwise-elongated platelet aggregates observed in microfluidic experiments.
  • Parameters: Amplitude ($A_s$) = 4 µm (1 RBC radius), Wavelength ($\lambda_s$) = 20 µm, Channel height ($L_z$) = 40 µm.
• Simulation Parameters:
  • Hematocrit (Ht): Varied at 0.33, 0.44, and 0.48.
  • Capillary Number (Ca): Varied from 0.1 to 0.6 (representing different shear rates and RBC deformability).
  • Platelet Count: $N_{plt} = N_{RBC}/2$ (higher than physiological to ensure statistical reliability).
  • Duration: Simulations ran for $1.5 \times 10^5$ time steps (~200 advection times) to reach statistical steady state.
• Analysis Metrics:
  • CFL Definition: Region where local RBC volume fraction < 50% of bulk Ht.
  • Platelet Capture: Defined as platelets reaching and remaining within a distance of $1.5 r_{plt}$ (1.5 µm) from the wall.
  • Ordering: RBC nematic order parameter ($Q$) calculated via inertia ellipsoids.

3. Key Contributions

1. Mechanistic Link between Geometry and Margination: Demonstrated that sinusoidal wall geometry creates spatially heterogeneous CFLs and shear rates, directly dictating where platelets accumulate.
2. Role of Ht in Aggregate Morphology: Identified that hematocrit is the primary regulator of platelet distribution along wavy surfaces, shifting accumulation from crests (low Ht) to a more uniform distribution (high Ht).
3. RBC Self-Organization: Revealed that at high Ca and low Ht, RBCs self-organize into ordered, packed structures in the valleys, enhancing hydrodynamic lift and thickening the CFL.
4. Shear-Dependent Adhesion Pathways: Mapped the local shear rate gradients between crests and valleys to specific biological adhesion mechanisms (vWF-mediated vs. collagen-mediated).

4. Key Results

A. Cell-Free Layer (CFL) and RBC Ordering

• Spatial Heterogeneity: The CFL is non-uniform. It is significantly thicker in the valleys and thinner at the crests.
• Dependence on Ht and Ca:
  • CFL thickness at crests and the straight wall decreases markedly with increasing Ht.
  • In valleys, CFL thickness increases with Ca (higher deformability) but decreases with Ht.
• RBC Ordering: At high Ca and low Ht, RBCs in the channel core and valleys form ordered, crystal-like structures. This ordering correlates positively with CFL thickness in the valleys, suggesting that self-organization enhances hydrodynamic lift, pushing RBCs away from the wall and expanding the CFL.

B. Platelet Margination and Capture Sites

• The "Size-Matching" Rule: Platelets preferentially marginate where the local CFL thickness is comparable to the platelet diameter (2 µm).
• Low Ht (0.33): The CFL at the crest is thin (~0.7–1.8 µm), matching the platelet size, while the valley CFL is thick (>2 µm). Consequently, platelets accumulate preferentially at the crests. This suggests a mechanism for the growth of high-amplitude, sharp aggregates.
• High Ht (0.48): Increasing Ht thins the CFL at the crest (making it too small for easy capture) and thins the valley CFL (bringing it closer to the platelet size). This leads to a more uniform platelet distribution along the sinusoid, consistent with experimental observations of smoother aggregate surfaces at high Ht.
• Independence from Ca: Platelet capture probability was found to be largely insensitive to Ca in the tank-treading regime, confirming that Ht is the dominant factor in this geometry.

C. Wall Shear Rate Gradients

• Crest vs. Valley: The sinusoidal geometry creates a pronounced shear rate gradient.
  • Crests: Experience shear rates roughly 2x the reference straight-wall value (reaching >1000 s⁻¹). This regime favors vWF-GPIbα interactions (high-shear adhesion).
  • Valleys: Experience shear rates roughly 0.5x the reference (often <500 s⁻¹). This regime favors collagen/fibrinogen-mediated adhesion (low-shear).
• RBC Influence: The presence of RBCs modulates this gradient. Reduced CFL thickness at the crest lowers the local shear, while RBCs in the valley enhance velocity gradients, increasing local shear. However, the gradient remains substantial.

5. Significance and Implications

• Morphological Evolution: The study provides a mechanistic explanation for how platelet aggregates evolve. Low Ht promotes "spiky" growth at crests due to localized margination, while high Ht promotes "smoother," more uniform growth.
• Therapeutic Targeting: By identifying specific "hotspots" for platelet capture and their associated shear environments, the study offers a rational basis for:
  • Shear-sensitive drug delivery: Designing carriers that target specific shear regimes (e.g., vWF-targeting drugs for crests).
  • Antithrombotic strategies: Tailoring therapies based on local hemodynamic conditions rather than global averages.
• Microfluidic Design: The findings can guide the design of microfluidic assays for studying thrombosis, ensuring that geometric heterogeneities are accounted for to mimic physiological conditions accurately.

In conclusion, this work bridges the gap between micro-scale hemodynamics and macro-scale thrombus morphology, highlighting that hematocrit and local geometry are the primary drivers of where and how platelet aggregates grow, mediated by the interplay between the Cell-Free Layer and shear-dependent adhesion pathways.