Lactation alters fluid flow and solute transport in maternal skeleton: A multiscale modeling study on the effects of microstructural changes and loading frequency

Human cortical bone intrinsic permeability distribution based on 3D canalicular morphology

Bone permeability is a key parameter that drives osteocyte-based mechanobiological modelling and remodelling. While previous experimental and numerical studies have estimated bone permeability based on the morphology of the lacuno-canalicular network, these studies often relied on simplified geometries. In the current study, bone permeability was characterized using more realistic canalicular geometry for the morphological data. Bone samples harvested from 27 human femoral bones were investigated using synchrotron radiation-based nano-computed tomography with a voxel size of 100 nm. After segmenting the canaliculi and lacunae, each canaliculus was investigated individually by applying a distance map and watershed algorithms. Bone permeability based on canalicular morphology was then assessed using the Kozeny relation, which defines the permeability of a porous medium with capillary-like pores. An averaged intrinsic permeability value of 8.8 10-18 m2 was obtained. It should be noted that this study considered an empty canalicular network, however in vivo, both cellular and peri-cellular matrices decrease space for interstitial fluid flow and thus permeability. Furthermore, a voxel size of 100 nm does not allow for the detection of smaller canaliculi, which may also modify average permeability. With the current data set and the analytic process applied, the results showed a heterogeneous permeability distribution within bone tissue, both when comparing osteonal and interstitial tissues and within an individual osteon. A difference was observed between male and female samples, and permeability appeared to significantly decrease with age. Finally, a significant correlation was found between permeability and canalicular length density, defined as canalicular length per unit bone volume. This study proposes a new form of the Kozeny law to express bone canalicular permeability as a proportional relationship with the canalicular length density. Importantly, this parameter can be directly quantified through confocal fluorescence microcopy, which is more convenient than synchrotron radiation-based nano-computed tomography. In conclusion, the current study confirms that confocal microscopy can be serve as a reliable tool to estimate bone permeability. However, the permeability values calculated here are solely based on canalicular morphology and do not consider cellular and peri-cellular intra-canalicular features.

Multiscale interstitial fluid computation modeling of cortical bone to characterize the hydromechanical stimulation of lacunar-canalicular network

Bone tissue is a biological composite material with a complex hierarchical structure that could continuously adjust its internal structure to adapt to the alterations in the external load environment. The fluid flow within bone is the main route of osteocyte metabolism, and the pore pressure as well as the fluid shear stress generated by it are important mechanical stimuli perceived by osteocytes. Owing to the irregular multiscale structure of bone tissue, the fluid stimulation that lacunar-canalicular network (LCN) in different regions of the tissue underwent remained unclear. In this study, we constructed a multiscale conduction model of fluid flow stimulus signals in bone tissue based on the poroelasticity theory. We analyzed the fluid flow behaviors at the macro-scale (whole bone tissue), macro-meso scale (periosteum, interstitial bone, osteon and endosteum), and micro-scale (lacunar-osteocyte-canalicular) levels. We explored how fluid stimulation at the tissue level correlated with that at the cellular level in cortical bone and characterized the distributions of the pore pressure, fluid velocity and fluid shear stress that the osteocytes experienced across the entire tissue structure. The results showed that the initial conditions of intramedullary pressure had a significant impact on the pore pressure of Haversian systems, but had a relatively small influence on the fluid velocity. The osteocyte which were located at different positions in the bone tissue received very distinct fluid stimuli. Osteocytes in the vicinity of the Haversian Canals experienced higher fluid shear stress stimulation. When the permeability of the LCN was within the range from 10-21 m2 to 10-18 m2, the distribution of pressure, fluid velocity and fluid shear stress within the osteon near the periosteum and endosteum was significantly different from that in other parts of the bone. However, when the permeability was less than 10-22 m2, such a difference did not exist. Particularly, the flow velocity at the lacunae was markedly higher than that in the canaliculi. Meanwhile, the pore pressure and fluid shear stress were conspicuously lower than those in the canaliculi. In this study, we considered the interconnections of different biofunctional units at different scales of bone tissue, construct a more complete multiscale model of bone tissue, and propose that osteocytes at different locations receive different fluid stimuli, which provides a reference for a deeper understanding of bone mechanotransduction.

Numerical study of interstitial fluid flow behavior in osteons under dynamic loading

BACKGROUND

The porous structure in bone tissue is essential for maintaining the physiological functions and overall health of intraosseous cells. The lacunar-canalicular net (LCN), a microscopic porous structure within osteons, facilitates the transport of nutrients and signaling molecules through interstitial fluid flow. However, the transient behavior of fluid flow within these micro-pores under dynamic loading conditions remains insufficiently studied.

METHODS

The study constructs a fluid-solid coupling model including the Haversian canal, canaliculi, lacunae, and interstitial fluid, to examine interstitial fluid flow behavior within the LCN under dynamic loading with varying frequencies and amplitudes. The relationship between changes of LCN pore volume and fluid velocity, and pressure is researched.

RESULTS

The results demonstrate that increasing strain amplitude leads to significant changes of LCN pore volume within osteons. In a complete loading cycle, with the increase of compressive strain, the pore volume in the osteon gradually shrinks, and the pressure gradient in the LCN increases, which promotes the increase of interstitial fluid velocity. When the compressive strain reaches the peak value, the flow velocity also reaches the maximum. In the subsequent unloading process, the pore volume began to recover, the pressure gradient gradually decreased, the flow rate decreased accordingly, and finally returned to the steady state level. At a loading amplitude of 1000 µε, the pore volume within LCN decreases by 1.1‰. At load amplitudes of 1500 µε, 2000 µε, and 2500 µε, the pore volume decreases by 1.6‰, 2.2‰ and 2.7‰ respectively, and the average flow velocity at the center of the superficial lacuna is 1.36 times, 1.77 times, and 2.14 times that at 1000 µε, respectively. Additionally, at a loading amplitude of 1000 µε under three different loading frequencies, the average flow velocities at the center of the superficial bone lacuna are 0.60 μm/s, 1.04 μm/s, and 1.54 μm/s, respectively. This indicates that high-frequency and high-amplitude dynamic loading can promote more vigorous fluid flow and pressure fluctuations with changes in LCN pore volume.

CONCLUSIONS

Dynamic mechanical loading can significantly enhance the interstitial fluid flow in LCN by the changes of LCN pore volume. and dynamic loading promoted fluid flow in shallow lacunae significantly higher than that in deep lacunae. The relationship between changes of LCN pore volume and interstitial fluid flow behavior has implications for drug delivery and bone tissue engineering research.

Osteocyte Dendrites: How Do They Grow, Mature, and Degenerate in Mineralized Bone?

Osteocytes, the most abundant bone cells, form an extensive cellular network via interconnecting dendrites. Like neurons in the brain, the long-lived osteocytes perceive mechanical and biological inputs and signal to other effector cells, leading to the homeostasis and turnover of bone tissues. Despite the appreciation of osteocytes' vital roles in bone biology, the initiation, growth, maintenance, and eventual degradation of osteocyte dendrites are poorly understood due to their full encasement by mineralized matrix. With the advancement of imaging modalities and genetic models, the architectural organization and molecular composition of the osteocyte dendrites, as well as their morphological changes with aging and diseases, have begun to be revealed. However, several long-standing mysteries remain unsolved, including (1) how the dendrites are initiated and elongated when a surface osteoblast becomes embedded as an osteocyte; (2) how the dendrites maintain a relatively stable morphology during their decades-long life span; (3) what biological processes control the dendrite morphology, connectivity, and stability; and (4) if these processes are influenced by age, sex, hormones, and mechanical loading. Our review of long, thin actin filament (F-actin)-containing processes extending from other cells leads to a working model that serves as a starting point to investigate the formation and maintenance of osteocyte dendrites and their degradation with aging and diseases.

Impact of gravity on fluid flow and solute transport in the bone lacunar-canalicular system: a multiscale numerical simulation study

Different gravity fields have important effects on the structural morphology of bone. The fluid flow caused by loadings in the bone lacunar-canalicular system (LCS), converts mechanical signals into biological signals and regulates bone reconstruction by affecting effector cells, which ensures the efficient transport of signaling molecules, nutrients, and waste products. In this study, the fluid flow and mass transfer effects of bone lacunar-canalicular system at multi-scale were firstly investigated, and a three-dimensional axisymmetric fluid-solid coupled finite element model of the LCS within three continuous osteocytes was established. The changes in fluid pressure field, flow velocity field, and fluid shear force variation on the surface of osteocytes within the LCS were studied comparatively under different gravitational fields (0 G, 1 G, 5 G), frequencies (1 Hz, 1.5 Hz, 2 Hz) and forms of cyclic compressive loading. The results showed that different frequencies represented different exercise intensities, suggesting that high-intensity exercise may accelerate the fluid flow rate within the LCS and enhance osteocytes activity. Hypergravity enhanced the transport of solute molecules, nutrients, and signaling molecules within the LCS. Conversely, the mass transfer in the LCS may be inhibited under microgravity, which may cause bone loss and eventually lead to the onset of osteoporosis. This investigation provides theoretical guidance for rehabilitative training against osteoporosis.

Aging decreases osteocyte peri-lacunar-canalicular system turnover in female C57BL/6JN mice

Osteocytes engage in bone resorption and mineralization surrounding their expansive lacunar-canalicular system (LCS) through peri-LCS turnover. However, fundamental questions persist about where, when, and how often osteocytes engage in peri-LCS turnover and how these processes change with aging. Furthermore, whether peri-LCS turnover is associated with natural variation in cortical tissue strain remains unexplored. To address these questions, we utilized confocal scanning microscopy, immunohistochemistry, and scanning electron microscopy to characterize osteocyte peri-LCS turnover in the cortical (mid-diaphysis) and cancellous (metaphysis) regions of femurs from young adult (5 mo) and early-old-age (22 mo) female C57BL/6JN mice. LCS bone mineralization was measured by the presence of perilacunar fluorochrome labels. LCS bone resorption was measured by immunohistochemical marker of bone resorption. The dynamics of peri-LCS turnover were estimated from serial fluorochrome labeling, where each mouse was administered two labels between 2 and 16 days before euthanasia. Osteocyte participation in mineralizing their surroundings is highly abundant in both cortical and cancellous bone of young adult mice but significantly decreases with aging. LCS bone resorption also decreases with aging. Aging has a greater impact on peri-LCS turnover dynamics in cancellous bone than in cortical bone. Lacunae with recent peri-LCS turnover are larger in both age groups. While peri-LCS turnover is associated with variation in tissue strain between cortical quadrants and intracortical location for 22 mo mice, these associations were not seen for 5 mo mice. The impact of aging on decreasing peri-LCS turnover may have significant implications for bone quality and mechanosensation.

Numerical simulation on mass transfer in the bone lacunar-canalicular system under different gravity fields

The bone lacunar-canalicular system (LCS) is a unique complex 3D microscopic tubular network structure within the osteon that contains interstitial fluid flow to ensure the efficient transport of signaling molecules, nutrients, and wastes to guarantee the normal physiological activities of bone tissue. The mass transfer laws in the LCS under microgravity and hypergravity are still unclear. In this paper, a multi-scale 3D osteon model was established to mimic the cortical osteon, and a finite element method was used to numerically analyze the mass transfer in the LCS under hypergravity, normal gravity and microgravity and combined with high-intensity exercise conditions. It was shown that hypergravity promoted mass transfer in the LCS to the deep lacunae, and the number of particles in lacunae increased more significantly from normal gravity to hypergravity the further away from the Haversian canal. The microgravity environment inhibited particles transport in the LCS to deep lacunae. Under normal gravity and microgravity, the number of particles in lacunae increased greatly when doing high-intensity exercise compared to stationary standing. This paper presents the first simulation of mass transfer within the LCS with different gravity fields combined with high-intensity exercise using the finite element method. The research suggested that hypergravity can greatly promote mass transfer in the LCS to deep lacunae, and microgravity strongly inhibited this mass transfer; high-intensity exercise increased the mass transfer rate in the LCS. This study provided a new strategy to combat and treat microgravity-induced osteoporosis.

Evaluating the Role of Canalicular Morphology and Perilacunar Region Properties on Local Mechanical Environment of Lacunar-Canalicular Network Using Finite Element Modeling

Physiological and pathological processes such as aging, diseases, treatments, and lactation can alter lacunar-canalicular network (LCN) morphology and perilacunar region properties. These modifications can impact the mechanical environment of osteocytes which in turn can influence osteocyte mechanosensitivity and the remodeling process. In this study, we aim to evaluate how the modifications in the canalicular morphology, lacunar density, and the perilacunar region properties influence the local mechanical environment of LCN and the apparent bone properties using three-dimensional finite element (FE) modeling. The simulation results showed that a 50% reduction in perilacunar elastic modulus led to about 7% decrease in apparent elastic modulus of the bone. The increase in canalicular density, length, and diameter did not influence the strain amplification in the models but they increased the amount of highly strained bone around LCN. Change in lacunar density did not influence the strain amplification and the amount of highly strained regions on LCN surfaces. Reduction in perilacunar elastic modulus increased both the strain amplification and the volume of highly strained tissue around and on the surface of LCN. The FE models of LCN in this study can be utilized to quantify the influence of modifications in canalicular morphology, lacunar density, and perilacunar region properties on the apparent bone properties and the local mechanical environment of LCN. Although this is a numerical study with idealized models, it provides important information on how mechanical environment of osteocytes is influenced by the modifications in LCN morphology and perilacunar region properties due to physiological and pathological processes.

Study on mass transfer in the bone lacunar-canalicular system under different gravity fields

INTRODUCTION

The bone lacunar-canalicular system (LCS) is an important microstructural basis for signaling and material transport in bone tissue, guaranteeing normal physiological processes in tissues. Spaceflight astronauts and elderly osteoporosis are related to its function, so it is necessary to reveal the mass transfer laws in bone microstructure under different gravity fields to provide insight for effective clinical treatment.

MATERIALS AND METHODS

Using the natural LCS structure of bovine tibial cortical bone as the object, the mass transfer experiments on cortical bone were conducted by using sodium fluorescein tracer through different frequency pulsating pressure provided by dynamic perfusion loading device and different high G environments provided by high-speed centrifuge to analyze the mass transfer laws under different gravity fields and different pulsating pressures.

RESULTS

The fluorescence intensity of lacunae within the osteon was lower the farther away from the Haversian canal. As the gravity field magnitude increased, the fluorescence intensity within each lacuna enhanced, and the more distant the lacunae from the Haversian canal, the greater the fluorescence intensity enhancement. High-frequency pulsating pressure simulated high-intensity exercise in humans can improve mass transfer efficiency in the LCS.

CONCLUSION

High-intensity exercise may greatly increase solute molecules, nutrients, and signaling molecules in osteocytes and improve the activity of osteocytes. Hypergravity can enhance the transport of solute molecules, nutrients, and signaling molecules in the LCS, especially promoting mass transfer to deep layer lacunae. Conversely, mass transfer to deep layer lacunae may be inhibited under microgravity, causing bone loss and ultimately leading to osteoporosis.

An Improved Method for Assessing Macroscale Diffusion through Osteonal Bone

The characterization of diffusion through biological tissues has played an important role in fundamental medical research and product development. Understanding the diffusion phenomena allows for the identification of new concepts in fundamental science, evolving medical knowledge and improving future standards and protocols. To illustrate, the structure of cortical bone changes upon the onset of osteoporosis, altering the limited porous compartment through which nutrients and essential signaling molecules travel to bone cells. Estrogen hormone replacement therapy (HRT) is one of the gold standard treatments to attempt to mitigate the effects that this structural change exerts in menopausal osteoporosis patients; however, HRT effectiveness is often variable in these patients, likely due to variability in bone structure and physiology, and thus transport rates. Scientists have studied diffusion in cortical bone tissue for decades. Current methodological standards include fluorescence recovery after photobleaching and computed tomography finite element analysis. Both techniques limit areas of tissue to microscale (1-100 µm2) analysis - only examining a few osteocytes within the structure at a time - and adopt assumptions that oversimplify in vivo tissue structure and transport phenomena. As well, the range of diffusion tracers is limited by the sensitivities of the analytical equipment, typically requiring tracer concentrations in the micromolar range. Herein is described a novel device for directly assessing the diffusion coefficient of 3H-estradiol at 37°C in macroscale osteonal bone specimens (1.4 cm2) - assessing a much larger portion of the total tissue than previously reported - while using radioisotope tracers for much higher sensitivity, thus achieving physiologically relevant estradiol concentrations. The current diffusion chamber device represents a cost-effective and validated method to mitigate these shortcomings. The device provides long-term diffusion data through macroscale (greater than 1 mm2) tissue areas, presenting a more physiologically accurate way to assess cortical bone diffusion. The device can assess solute diffusion through other tissues or materials and may easily be scaled up to run multiple diffusion experiments simultaneously.

Structural role of osteocyte lacunae on mechanical properties of bone matrix: A cohesive finite element study

Despite the extensive studies on biological function of osteocytes, there are limited studies that evaluated the structural role of osteocyte lacunae on local mechanical properties of the bone matrix. As a result, the goal of this study was to elucidate the independent contribution of osteocyte lacunae structure on mechanical properties and fracture behavior of the bone matrix uncoupled from its biological effects and bone tissue composition variation. This study combined cohesive finite element modeling with experimental data from a lactation rat model to evaluate the influence of osteocyte lacunar area porosity, density, size, axis ratio, and orientation on the elastic modulus, ultimate strength, and ultimate strain of the bone matrix as well as on local crack formation and propagation. It also performed a parametric study to isolate the influence of a single osteocyte lacunae structural property on the mechanical properties of the bone matrix. The experimental measurements demonstrated statistically significant differences in lacunar size between ovariectomized rats with lactation history and virgin groups (both ovariectomized and intact) and in axis ratio between rats with lactation history and virgins. There were no differences in mechanical properties between virgin and lactation groups as determined by the finite element simulations. However, there were statistically significant linear relationships between the physiological range of osteocyte lacunar area porosity, density, size, and orientation and the elastic modulus and ultimate strength of the bone matrix in virgin and lactation rats. The parametric study also revealed similar but stronger relationships between elastic modulus and ultimate strength and lacunar density, size, and orientation. The simulations also demonstrated that the osteocyte lacunae guided the crack propagation through local stress concentrations. In summary, this study enhanced the limited knowledge on the structural role of osteocyte lacunae on local mechanical properties of the bone matrix. These data are important in gaining a better understanding of the mechanical implications of the local modifications due to osteocytes in the bone matrix.

Osteocytes and Weightlessness

PURPOSE OF REVIEW

Osteocytes are considered to be the cells responsible for mastering the remodeling process that follows the exposure to unloading conditions. Given the invasiveness of bone biopsies in humans, both rodents and in vitro culture systems are largely adopted as models for studies in space missions or in simulated microgravity conditions models on Earth.

RECENT FINDINGS

After a brief recall of the main changes in bone mass and osteoclastic and osteoblastic activities in space-related models, this review focuses on the potential role of osteocytes in directing these changes. The role of the best-known signalling molecules is questioned, in particular in relation to osteocyte apoptosis. The mechanotransduction actors identified in spatial conditions and the problems related to fluid flow and shear stress changes, probably enhanced by the alteration in fluid flow and lack of convection during spaceflight, are recalled and discussed.