An official website of the United States government
An osteocyte is a bone cell situated inside a hard bone matrix in an interstice (lacuna). It has many dendritic structures called cellular processes that radiate outward from the cell through the bone matrix via cylindrical openings (canaliculi). Osteocytes can sense stress and strain applied by the interstitial fluid flow and respond by releasing biochemical signals that regulate bone remodeling. In vitro experiments have suggested that the stress and strain typically experienced at the macroscale tissue level have to be amplified 10× in order for osteocytes to have a significant response in vivo. This stress and strain amplification mechanism is not yet well understood. Previous studies suggest that the processes are the primary sites for mechanosensation thanks to the tethering elements that attach the process membrane to the canalicular wall. However, there are other potential factors which may also contribute to stress and strain amplification, such as canalicular wall geometry and osteocyte-associated proteins in the interstitial space called pericellular matrix. In this work, we perform computational studies to study how canalicular wall roughness affects stress and strain amplification. Our major finding is that the wall roughness induces significantly greater wall shear stress (WSS) on the process when the wall roughness increases flow resistance; and the roughness has relatively smaller influence on the WSS when the resistance remains the same.
more »
« less
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
- PAR ID:
- 10327293
- 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
-
-
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
-
Nguyen, Thao Vicky; Ethier, C Ross (Ed.)Abstract Cell-laden, scaffold-based tissue engineering methods have been successfully utilized for the treatment of bone fractures and diseases, caused by factors such as trauma, tumors, congenital anomalies, and aging. In such methods, the rate of scaffold biodegradation, transport of nutrients and growth factors, as well as removal of cell metabolic wastes at the site of injury are critical fluid-dynamics factors, affecting cell proliferation and ultimately tissue regeneration. Therefore, there is a critical need to identify the underlying material transport mechanisms and factors associated with cell-seeded, scaffold-based bone tissue engineering. The overarching goal of this study is to contribute to patient-specific, clinical treatment of bone pathology. The overall objective of the work is to establish computational fluid dynamics (CFD) models: (i) to identify the consequential mechanisms behind internal and external material transport through/over porous bone scaffolds designed based on the principles of triply periodic minimal surfaces (TPMS) and (ii) to identify TPMS designs with optimal geometry and flow characteristics for the treatment of bone fractures in clinical practice. In this study, advanced CFD models were established based on ten TPMS scaffold designs for (i) single-unit internal flow analysis, (ii) single-unit external flow analysis, and (iii) cubic, full-scaffold external flow analysis, where the geometry of each design was parametrically created. The influence of several design parameters, such as surface representation iteration, wall thickness, and pore size on geometry accuracy as well as computation time, was investigated in order to obtain computationally efficient and accurate CFD models. The fluid properties (such as density and dynamic viscosity) as well as the boundary conditions (such as no-slip condition, inlet flow velocity, and pressure outlet) of the CFD models were set based on clinical/research values reported in the literature, according to the fundamentals of internal and external Newtonian flow modeling. The main fluid characteristics influential in bone regeneration, including flow velocity, flow pressure, and wall shear stress (WSS), were analyzed to observe material transport internally through and externally over the TPMS scaffold designs. Regarding the single-unit internal flow analysis, it was observed that P.W. Hybrid and Neovius designs had the highest level of not only flow pressure but also WSS. This can be attributed to their relatively flat surfaces when compared to the rest of the TPMS designs. Schwarz primitive (P) appeared to have the lowest level of flow pressure and WSS (desirable for development of bone tissues) due to its relatively open channels allowing for more effortless fluid transport. An analysis of streamline velocity exhibited an increase in velocity togther with a depiction of potential turbulent motion along the curved sections of the TPMS designs. Regarding the single-unit external flow analysis, it was observed that Neovius and Diamond yielded the highest level of flow pressure and WSS, respectively, while Schwarz primitive (P) similarly had a relatively low level of flow pressure and WSS suitable for bone regeneration. Besides, pressure buildup was observed within the inner channels of almost all the TPMS designs due to flow resistance and the intrinsic interaction between the fluid flow and the scaffold walls. Regarding the cubic (full-scaffold) external flow analysis, the Diamond and Schwarz gyroid (G) designs appeared to have a relatively high level of both flow pressure and WSS, while Schwarz primitive (P) similarly yielded a low level of flow pressure and WSS. Overall, the outcomes of this study pave the way for optimal design and fabrication of complex, bone-like tissues with desired material transport properties for cell-laden, scaffold-based treatment of bone fractures.more » « less
-
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
-
Turbulence statistics and blade deformations of flexible emergent canopies impinged by water flows were experimentally investigated across a range of Reynolds numbers Reb=Ubb/ν (where Ub is the bulk incoming flow velocity, b is the blade width, and ν is the water kinematic viscosity) and blade aspect ratios AR=h/b (h is the blade length). Time-resolved particle image velocimetry was used to characterize both the deformation of flexible blades and the surrounding flow fields. Results showed that the blade deformation increased with the growth of both Reb and AR, with higher blade bending causing stronger variations in vertical profiles of streamwise velocities and Reynolds stresses. The drag produced by the presence of flexible canopies was identified as the dominant fluid loading balancing the pressure gradient. This term exhibited distinctive reduction near the water surface region with high blade deformation due to the large local blade inclination angle. Interestingly, in contrast to fully submerged flexible blades where the flow-induced drag increases monotonously with flow speed, a critical Reynolds number Reb,cri was observed, beyond which drag decreased with increasing flow speed until the blade became fully submerged. This phenomenon was explained with theoretical interpretations, which exhibited reasonable agreement with experimental results. Further analysis of unsteady flow dynamics revealed that Reynolds stress within the canopy was dominated by ejection events due to the absence of shear layer at the top of emergent canopy. Additionally, streamwise velocity spectra indicated that flow fluctuations inside the canopy were governed by periodic vortex shedding from blade.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0085299
