ABSTRACT ObjectiveTortuous microvessels are characteristic of microvascular remodeling associated with numerous physiological and pathological scenarios. Three‐dimensional (3D) hemodynamics in tortuous microvessels influenced by red blood cells (RBCs), however, are largely unknown, and important questions remain. Is blood viscosity influenced by vessel tortuosity? How do RBC dynamics affect wall shear stress (WSS) patterns and the near‐wall cell‐free layer (CFL) over a range of conditions? The objective of this work was to parameterize hemodynamic characteristics unique to a tortuous microvessel. MethodsRBC‐resolved simulations were performed using an immersed boundary method‐based 3D fluid dynamics solver. A representative tortuous microvessel was selected from a stimulated angiogenic network obtained from imaging of the rat mesentery and digitally reconstructed for the simulations. The representative microvessel was a venule with a diameter of approximately 20 μm. The model assumes a constant diameter along the vessel length and does not consider variations due to endothelial cell shapes or the endothelial surface layer. ResultsMicrovessel tortuosity was observed to increase blood apparent viscosity compared to a straight tube by up to 26%. WSS spatial variations in high curvature regions reached 23.6 dyne/cm2over the vessel cross‐section. The magnitudes of WSS and CFL thickness variations due to tortuosity were strongly influenced by shear rate and negligibly influenced by tube hematocrit levels. ConclusionsNew findings from this work reveal unique tortuosity‐dependent hemodynamic characteristics over a range of conditions. The results provide new thought‐provoking information to better understand the contribution of tortuous vessels in physiological and pathological processes and help improve reduced‐order models.
more »
« less
Angiogenic Microvascular Wall Shear Stress Patterns Revealed Through Three-dimensional Red Blood Cell Resolved Modeling
Abstract The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC “footprints.” Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.
more »
« less
- Award ID(s):
- 2309559
- PAR ID:
- 10465047
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Function
- Volume:
- 4
- Issue:
- 6
- ISSN:
- 2633-8823
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Blood velocity and red blood cell (RBC) distribution profiles in a capillary vessel cross-section in the microcirculation are generally complex and do not follow Poiseuille’s parabolic or uniform pattern. Existing imaging techniques used to map large microvascular networks in vivo do not allow a direct measurement of full 3D velocity and RBC concentration profiles, although such information is needed for accurate evaluation of the physiological variables, such as the wall shear stress (WSS) and near-wall cell-free layer (CFL), that play critical roles in blood flow regulation, disease progression, angiogenesis, and hemostasis. Theoretical network flow models, often used for hemodynamic predictions in experimentally acquired images of the microvascular network, cannot provide the full 3D profiles either. In contrast, such information can be readily obtained from high-fidelity computational models that treat blood as a suspension of deformable RBCs. These models, however, are computationally expensive and not feasible for extension to the microvascular network at large spatial scales up to an organ level. To overcome such limitations, here we present machine learning (ML) models that bypass such expensive computations but provide highly accurate and full 3D profiles of the blood velocity, RBC concentration, WSS, and CFL in every vessel in the microvascular network. The ML models, which are based on artificial neural network and convolution-based U-net models, predict hemodynamic quantities that compare very well against the true data but reduce the prediction time by several orders. This study therefore paves the way for ML to make detailed and accurate hemodynamic predictions in spatially large microvascular networks at an organ-scale.more » « less
-
Abstract ObjectiveFluid shear stress is thought to be a regulator of endothelial cell behavior during angiogenesis. The link, however, requires an understanding of stress values at the capillary level in angiogenic microvascular networks. Critical questions remain. What are the stresses? Do capillaries experience similar stress magnitudes? Can variations explain vessel‐specific behavior? The objective of this study was to estimate segment‐specific shear stresses in angiogenic networks. MethodsImages of angiogenic networks characterized by increased vascular density were obtained from rat mesenteric tissues stimulated by compound 48/80‐induced mast cell degranulation. Vessels were identified by perfusion of a 40 kDa fixable dextran prior to harvesting and immunolabeling for PECAM. Using a network flow‐based segment model with physiologically relevant parameters, stresses were computed per vessel for regions across multiple networks. ResultsStresses ranged from 0.003 to 2328.1 dyne/cm2and varied dramatically at the capillary level. For all regions, the maximum segmental shear stresses were for capillary segments. Stresses along proximal capillaries branching from arteriole inlets were increased compared to stresses along capillaries in more distal regions. ConclusionsThe results highlight the variability of shear stresses along angiogenic capillaries and motivate new discussions on how endothelial cells may respond in vivo to segment‐specific microenvironment during angiogenesis.more » « less
-
How do vessels find optimal radii? Capillaries are known to adapt their radii to maintain the shear stress of blood flow at the vessel wall at a set point, yet models of adaptation purely based on average shear stress have not been able to produce complex loopy networks that resemble real microvascular systems. For narrow vessels where red blood cells travel in a single file, the shear stress on vessel endothelium peaks sharply when a red blood cell passes through. We show that stable shear-stress-based adaptation is possible if vessel shear stress set points are cued to the stress peaks. Model networks that respond to peak stresses alone can quantitatively reproduce the observed zebrafish trunk microcirculation, including its adaptive trajectory when hematocrit changes or parts of the network are amputated. Our work reveals the potential for mechanotransduction alone to generate stable hydraulically tuned microvascular networks.more » « less
-
AbstractSickle cell disease (SCD) is canonically characterized by reduced red blood cell (RBC) deformability, leading to microvascular obstruction and inflammation. Although the biophysical properties of sickle RBCs are known to influence SCD vasculopathy, the contribution of poor RBC deformability to endothelial dysfunction has yet to be fully explored. Leveraging interrelated in vitro and in silico approaches, we introduce a new paradigm of SCD vasculopathy in which poorly deformable sickle RBCs directly cause endothelial dysfunction via mechanotransduction, during which endothelial cells sense and pathophysiologically respond to aberrant physical forces independently of microvascular obstruction, adhesion, or hemolysis. We demonstrate that perfusion of sickle RBCs or pharmacologically-dehydrated healthy RBCs into small venule-sized “endothelialized” microfluidics leads to pathologic physical interactions with endothelial cells that directly induce inflammatory pathways. Using a combination of computational simulations and large venule-sized endothelialized microfluidics, we observed that perfusion of heterogeneous sickle RBC subpopulations with varying deformability, as well as suspensions of dehydrated normal RBCs admixed with normal RBCs, leads to aberrant margination of the less-deformable RBC subpopulations toward the vessel walls, causing localized, increased shear stress. Increased wall stress is dependent on the degree of subpopulation heterogeneity and oxygen tension and leads to inflammatory endothelial gene expression via mechanotransductive pathways. Our multifaceted approach demonstrates that the presence of sickle RBCs with reduced deformability leads directly to pathological physical (ie, direct collisions and/or compressive forces) and shear-mediated interactions with endothelial cells and induces an inflammatory response, thereby elucidating the ubiquity of vascular dysfunction in SCD.more » « less