skip to main content


Title: Modeling and simulation of interstitial fluid flow around an osteocyte in a lacuno-canalicular network
Experiments have shown that external mechanical loading plays an important role in bone development and remodeling. In fact, recent research has provided evidence that osteocytes can sense such loading and respond by releasing biochemical signals (mechanotransduction, MT) that initiate bone degradation or growth. Many aspects on MT remain unclear, especially at the cellular level. Because of the extreme hardness of the bone matrix and complexity of the microenvironment that an osteocyte lives in, in vivo studies are difficult; in contrast, modeling and simulation are viable approaches. Although many computational studies have been carried out, the complex geometry that can involve 60+ irregular canaliculi is often simplified to a select few straight tubes or channels. In addition, the pericellular matrix (PCM) is usually not considered. To better understand the effects of these frequently neglected aspects, we use the lattice Boltzmann equations to model the fluid flow over an osteocyte in a lacuno-canalicular network in two dimensions. We focus on the influences of the number/geometry of the canaliculi and the effects of the PCM on the fluid wall shear stress (WSS) and normal stress (WNS) on an osteocyte surface. We consider 16, 32, and 64 canaliculi using one randomly generated geometry for each of the 16 and 32 canaliculi cases and three geometries for the 64 canaliculi case. We also consider 0%, 5%, 10%, 20%, and 40% pericellular matrix density. Numerical results on the WSS and WNS distributions and on the velocity field are visualized, compared, and analyzed. Our major results are as follows: (1) the fluid flow generates significantly greater force on the surface of the osteocyte if the model includes the pericellular matrix (PCM); (2) in the absence of PCM, the average magnitudes of the stresses on the osteocyte surface are not significantly altered by the number and geometry of the canaliculi despite some quantitative influence of the latter on overall variation and distribution of those stresses; and (3) the dimensionless stress (stress after non-dimensionalization) on the osteocyte surface scales approximately as the reciprocal of the Reynolds number and increasing PCM density in the canaliculi reduces the range of Reynolds number values for which the scaling law holds.  more » « less
Award ID(s):
1951531
NSF-PAR ID:
10327293
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Physics of Fluids
Volume:
34
Issue:
4
ISSN:
1070-6631
Page Range / eLocation ID:
041906
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Transient plasma membrane disruptions (PMD) occur in osteocytes with in vitro and in vivo loading, initiating mechanotransduction. The goal here was to determine whether osteocyte PMD formation or repair is affected by aging. Osteocytes from old (24 months) mice developed fewer PMD (−76% females, −54% males) from fluid shear than young (3 months) mice, and old mice developed fewer osteocyte PMD (−51%) during treadmill running. This was due at least in part to decreased pericellular matrix production, as studies revealed that pericellular matrix is integral to formation of osteocyte PMD, and aged osteocytes produced less pericellular matrix (−55%). Surprisingly, osteocyte PMD repair rate was faster (+25% females, +26% males) in osteocytes from old mice, and calcium wave propagation to adjacent nonwounded osteocytes was blunted, consistent with impaired mechanotransduction downstream of PMD in osteocytes with fast PMD repair in previous studies. Inducing PMD via fluid flow in young osteocytes in the presence of oxidative stress decreased postwounding cell survival and promoted accelerated PMD repair in surviving cells, suggesting selective loss of slower‐repairing osteocytes. Therefore, as oxidative stress increases during aging, slower‐repairing osteocytes may be unable to successfully repair PMD, leading to slower‐repairing osteocyte death in favor of faster‐repairing osteocyte survival. Since PMD are an important initiator of mechanotransduction, age‐related decreases in pericellular matrix and loss of slower‐repairing osteocytes may impair the ability of bone to properly respond to mechanical loading with bone formation. These data suggest that PMD formation and repair mechanisms represent new targets for improving bone mechanosensitivity with aging.

     
    more » « less
  2. ABSTRACT

    Localized apoptosis of osteocytes, the tissue-resident cells within bone, occurs with fatigue microdamage and activates bone resorption. Osteoclasts appear to target and remove dying osteocytes, resorbing damaged bone matrix as well. Osteocyte apoptosis similarly activates bone resorption with estrogen loss and in disuse. Apoptotic osteocytes trigger viable neighbor (ie, bystander) osteocytes to produce RANKL, the cytokine required for osteoclast activation. Signals from apoptotic osteocytes that trigger this bystander RANKL expression remain obscure. Studying signaling among osteocytes has been hampered by lack of in vitro systems that model the limited communication among osteocytes in vivo (ie, via gap junctions on cell processes and/or paracrine signals through thin pericellular fluid spaces around osteocytes). Here, we used a novel multiscale fluidic device (the Macro-micro-nano, or Mμn) that reproduces these key anatomical features. Osteocytes in discrete compartments of the device communicate only via these limited pathways, which allows assessment of their roles in triggering osteocytes RANKL expression. Apoptosis of MLOY-4 osteocytes in the Mμn device caused increased osteocyte RANKL expression in the neighboring compartment, consistent with in vivo findings. This RANKL upregulation in bystander osteocytes was prevented by blocking Pannexin 1 channels as well as its ATP receptor. ATP alone caused comparable RANKL upregulation in bystander osteocytes. Finally, blocking Connexin 43 gap junctions did not abolish osteocyte RANKL upregulation, but did alter the distribution of RANKL expressing bystander osteocytes. These findings point to extracellular ATP, released from apoptotic osteocytes via Panx1 channels, as a major signal for triggering bystander osteocyte RANKL expression and activating bone remodeling. © 2020 American Society for Bone and Mineral Research.

     
    more » « less
  3. Articular cartilage is comprised of two main components, the extracellular matrix (ECM) and the pericellular matrix (PCM). The PCM helps to protect chondrocytes in the cartilage from mechanical loads, but in patients with osteoarthritis, the PCM is weakened, resulting in increased chondrocyte stress. As chondrocytes are responsible for matrix synthesis and maintenance, it is important to understand how mechanical loads affect the cellular responses of chondrocytes. Many studies have examined chondrocyte responses to in vitro mechanical loading by embedding chondrocytes in 3-D hydrogels. However, these experiments are mostly performed in the absence of PCM, which may obscure important responses to mechanotransduction. Here, drop-based microfluidics is used to culture single chondrocytes in alginate microgels for cell-directed PCM synthesis that closely mimics the in vivo microenvironment. Chondrocytes formed PCM over 10 days in these single-cell 3-D microenvironments. Mechanotransduction studies were performed, in which single-cell microgels mimicking the cartilage PCM were embedded in high-stiffness agarose. After physiological dynamic compression in a custom-built bioreactor, microgels exhibited distinct metabolomic profiles from both uncompressed and monolayer controls. These results demonstrate the potential of single cell encapsulation in alginate microgels to advance cartilage tissue engineering and basic chondrocyte mechanobiology. 
    more » « less
  4. Abstract

    Freshwater mussels are dominant ecosystem engineers in many streams throughout North America, yet they remain among the world's most imperiled fauna. Extensive research has quantified the ecological role of mussels in aquatic habitats, but little is known about the interaction between mussels and their surrounding physical and hydrodynamic habitat. Here the physical interactions of mussels with near‐bed flow are investigated in an experimental channel using model mussels. The results show that (1) mussels disrupt the distributions and magnitudes of time‐averaged values of longitudinal flow velocity and Reynolds shear stress depending on mussel density, and (2) at densities of approximately 25 mussels m−2and greater, a hydrodynamic transition occurs where the maximum Reynolds shear stress is displaced from the bed to the height of the mussel canopy, near‐bed longitudinal flow velocity is reduced, and average turbulent shear stresses acting on the mussels are reduced by as much as 64%, thus markedly decreasing the dislodgement potential of the mussels by these stresses. These results provide strong empirical evidence for a positive density‐dependent effect related to flow‐organism interactions and their ecological success, such as enhancing river bed hydrodynamic habitat complexity or decreasing the turbulent shear stresses acting to dislodge mussels from the river bed. This information will improve the understanding of the long‐term persistence of mussel beds and help focus future conservation strategies.

     
    more » « less
  5. Abstract

    Spiral/helical forms of blood flow have been observed in large arteries of the cardiovascular system, but their benefits remain underappreciated. Spiral flow has been postulated to improve near‐wall washout, promoting anti‐atherothrombotic conditions. This research aims to study the washout characteristics of spiral flow, specifically, its ability to increase velocity and wall shear stress (WSS) in atherothrombotic‐prone regions. Using 1.2 cm diameter angled test‐conduits (45°, 90°, 135°) with known recirculation/stasis regions at the bend corners, spiral flow washout potential was evaluated in terms of low velocity and low WSS. Two sub‐studies were conducted: the first utilized a spiral flow‐inducing device to enable qualitative analysis of washout‐potential in both computational fluid dynamic (CFD) simulations and benchtop ultrasound visualization; the second used CFD to study the impact of several induced helical wavelengths on the conduit‐dependent recirculation/stasis zones. Physical models of the angled conduits and spiral flow‐inducer were 3D‐printed to facilitate ultrasound visualization. Compared to straight flow, spiral flow generated by the flow‐inducer significantly cleared the recirculation/stasis zones at the corners of the angled conduits. CFD simulations demonstrated that past a geometry‐dependent threshold, increased helical content improved washout, denoted by decreased regions of low velocity and low WSS. Overall, spiral flow markedly improved washout in difficult to reach areas in the angled conduits. This has several important clinical implications: spiral flow shows great promise in reducing blood‐transport‐related complications and can be used to enhance the performance of future medical devices (eg grafts, mechanical circulatory support devices, hemodialysis access ports).

     
    more » « less