High-resolution image-based simulation reveals membrane strain concentration on osteocyte processes caused by tethering elements

Osteocytes are vital for regulating bone remodeling by sensing the flow-induced mechanical stimuli applied to their cell processes. In this mechanosensing mechanism, tethering elements (TEs) connecting the osteocyte process with the canalicular wall potentially amplify the strain on the osteocyte processes. The ultrastructure of the osteocyte processes and canaliculi can be visualized at a nanometer scale using high-resolution imaging via ultra-high voltage electron microscopy (UHVEM). Moreover, the irregular shapes of the osteocyte processes and the canaliculi, including the TEs in the canalicular space, should considerably influence the mechanical stimuli applied to the osteocytes. This study aims to characterize the roles of the ultrastructure of osteocyte processes and canaliculi in the mechanism of osteocyte mechanosensing. Thus, we constructed a high-resolution image-based model of an osteocyte process and a canaliculus using UHVEM tomography and investigated the distribution and magnitude of flow-induced local strain on the osteocyte process by performing fluid–structure interaction simulation. The analysis results reveal that local strain concentration in the osteocyte process was induced by a small number of TEs with high tension, which were inclined depending on the irregular shapes of osteocyte processes and canaliculi. Therefore, this study could provide meaningful insights into the effect of ultrastructure of osteocyte processes and canaliculi on the osteocyte mechanosensing mechanism.


Introduction
Osteocytes are the most abundant cells present inside the bone matrix and are vital for bone remodeling as they regulate osteoclastic bone resorption and osteoblastic bone formation in response to a mechanical stimuli (Tatsumi et al. 2007;Nakashima et al. 2011;Xiong et al. 2011). Typically, osteocytes are sensitive to local mechanical stimuli, especially in cell processes (Adachi et al. 2009;Wu et al. 2011). Results of prior research have indicated that a reduction in mechanical stimuli owing to disuse causes osteoporotic conditions in the bone structure (Aguirre et al. 2006;Liu et al. 2020). Therefore, the clarification of mechanosensing mechanisms related to osteocyte processes is necessary for understanding the cellular contribution toward maintaining the physiological bone structure.
In the mechanosensing of osteocytes, local strain on the cell process membrane is considered a representative indicator of mechanical stimuli as it can activate gating of mechanosensitive ion channels such as Piezo1 (McMahon et al. 2008;Li et al. 2019;Diem et al. 2020;Sasaki et al. 2020;Zhou et al. 2020). Although a ~ 10% strain on the bone cell membranes can induce Ca 2+ responses (You et al. 2000), the strain of mineralized bone tissue during normal locomotion does not exceed ~ 0.1% (Fritton et al. 2000). Thus, a 1 3 strain amplification model was proposed in previous studies to infuse the gap in strain magnitude between the cellular and tissue levels Zeng et al. 1994;You et al. 2001;Han et al. 2004;Wang et al. 2007Wang et al. , 2014Scheiner et al. 2016). In context, the interstitial fluid flow generated by external forces in the canalicular space is considered a principal factor in strain amplification Zeng et al. 1994;Scheiner et al. 2016). Furthermore, the pericellular matrix in the canalicular space, such as proteoglycans, can additionally facilitate strain amplification. Specifically, tethering elements (TEs) connecting the osteocyte process to the canalicular wall, which are regarded as the core protein of perlecan, are supposed to provide drag force to the cell processes via interstitial fluid flow (You et al. 2001;Han et al. 2004;Wang et al. 2007Wang et al. , 2014. The impacts of interstitial fluid flow and TEs on the strain of osteocyte processes have been revealed through theoretical analyses (You et al. 2001;Han et al. 2004;Wang et al. 2007), in which the mechanism of strain amplification via fluid flow and TEs was clarified using idealized cylindrical models representing the osteocyte processes and canaliculi. Despite the significance of this mechanism, the irregular shape of osteocyte processes and canaliculi is important for amplifying the strain as it can induce heterogeneous distribution of interstitial fluid velocity (Kamioka et al. 2012) and thus affect the local strain on the cell processes. Furthermore, the deterioration in the canalicular structure owing to aging or bone diseases (Jast and Jasiuk 2013;Lai et al. 2015; Tiede-Lewis and Dallas 2019) can alter the flow-induced strain of osteocytes and influence subsequent bone remodeling. Therefore, the contribution of the irregular shape of osteocyte processes and canaliculi toward influencing the local osteocyte strain is required to be investigated.
Computational studies utilizing image-based models are useful for quantitatively evaluating the mechanical stimuli required for osteocyte processes. In previous studies, the lacuno-canalicular model constructed from confocal laser scanning images was used to investigate the distribution of interstitial fluid velocity and osteocyte strain (Verbruggen et al. 2012(Verbruggen et al. , 2014(Verbruggen et al. , 2016. Although the confocal laser scanning microscopy cannot appropriately capture the microscopic surface of osteocyte processes and canaliculi, ultrahigh voltage electron microscope (UHVEM) tomography enables observations of their ultrastructure at a nanometer scale (Kamioka et al. 2012;Hosaki-Takamiya et al. 2016). Moreover, computational analysis using such high-resolution image-based models can potentially reveal the influence of the ultrastructure on the local strain of osteocyte processes.
This study aims to characterize the roles of the ultrastructure of osteocyte processes and canaliculi in the osteocyte mechanosensing mechanism via interstitial fluid flow and TEs. Thus, we constructed a high-resolution image-based model of an osteocyte process and a canaliculus using UHVEM tomography (Kamioka et al. 2012). The influence of TEs on the mechanical stimuli for osteocyte processes was investigated by introducing randomly distributed spring models inside the canalicular space of the image-based model. Subsequently, the distribution and magnitude of the flow-induced local strain on the osteocyte process were examined through fluid-structure interaction simulation, the results of which revealed that the heterogeneous inclination of TEs produced by the irregular shape of osteocyte processes and canaliculi is a key factor in activating osteocyte mechanosensing.

Fluid-structure interaction simulation
The deformation of osteocyte processes induced by interstitial fluid flow and TEs was investigated by developing a method for fluid-structure interaction simulation based on previous studies (Takeishi et al. 2019). In addition, interstitial fluid flow through the pericellular matrix in the canaliculus was assumed to follow the Brinkman equation , in which the fluid velocity, u, induced by the pressure gradient, ∇p , and the body force, F, can be described as where μ is the fluid viscosity and k p the permeability of the pericellular matrix. Furthermore, the lattice Boltzmann method (Chen and Doolen 1998) was employed for the numerical analysis of the fluid flow governed by the Brinkman equation.
The osteocyte process membrane was modeled as a thin hyperelastic membrane and discretized into triangular finite elements. The deformation of the hyperelastic membrane was assumed to be governed by Skalak's constitutive law (Skalak et al. 1973), in which the strain energy function can be expressed as where I 1 and I 2 are the first and second strain invariants of the right Cauchy-Green tensor, G s the surface shear elastic modulus, and C the area incompressibility coefficient. The surface Poisson's ratio, s , and Young's modulus, E, of the osteocyte cell membrane can be expressed using the membrane thickness, t, as s = C C+1 and E = 2G s (1+ s) t . In addition, the restoring force, q m , on the membrane node, x m , was calculated from the membrane displacement based on a finite element procedure. Here, all the displacements of the membrane nodes were determined from the advection of the fluid, as described next. A TE was introduced as a linear elastic spring, which connects a point on the membrane of the osteocyte process, x m , to a point on the canalicular wall, x w . The point of attachment of each TE on the canalicular wall, x w , was constant during the simulation. The spring constant, k sp , was determined as k sp = , where E sp , r sp , and l 0 are the Young's modulus, radius, and the natural length of TEs, respectively. Consequently, the restoring force, q t , acting on the membrane node induced by the stretching of the TEs, i.e. the displacement of the membrane node x m , was calculated as Moreover, the immersed boundary method (Peskin 2002) was employed to couple the interstitial fluid flow and cell membrane deformation. The total restoring force acting on the membrane, q = q m + q t , was distributed to the surrounding fluid nodes, x f , and it acted as an external body force, F, expressed as where D(x) is the numerically approximated Dirac delta function given by where Δx is the lattice interval in the fluid domain, which is also described in Sect. 2.2. Conversely, the cell membrane was advected with velocity U obtained by interpolating the surrounding fluid velocity, u, using the approximated Dirac delta function as to satisfy the non-slip condition on the membrane surface. By using the displacement of membrane node x m advected by membrane velocity U, restoring force q m on the membrane node was calculated again.

Reconstruction of image-based models
The three-dimensional image-based model of an osteocyte process and surrounding canaliculus (Fig. 1a) was reconstructed from UHVEM tomographic images (Kamioka et al. 2012), and a Cartesian coordinate system (x, y, z) was set along the longitudinal direction of the osteocyte process with the canaliculus in the z direction. The original resolution of the tomographic images was 1.7 nm/pixel, and the length of the cell process and canaliculus along the z direction was 453.9 nm. Based on these images, the canalicular space was discretized by a cubic lattice with intervals Δx = 5.1 nm (= 3 pixels) for the lattice Boltzmann method. Thereafter, the hyperelastic membrane representing the osteocyte process membrane was discretized into triangular finite elements using MeshLab (Visual Computing Lab, ISTI-CNR, Italy). Typically, the representative length of the individual elements was ~ 4.7 nm. The attachment points of TEs on the osteocyte process were randomly arranged with a number density of 12/37.5 nm −1 along the z direction (Han et al. 2004;You et al. 2004;Wang et al. 2007). Then, the other end of each TE was attached to the nearest point on the canalicular wall. As each TE was inclined under no-flow conditions depending on the irregular shape of the canalicular wall (Fig. 1b), its natural length, l 0 , was set as the length under no-flow condition. In this study, we constructed five models, each of which has a different random configuration of TEs. Furthermore, a mirror symmetric model was connected to the original model with respect to the x-y plane to apply the periodic boundary conditions for the fluid-structure interaction simulation, as depicted in Fig. 1a.
The values of the parameters used in the simulation are listed in Table 1 Sugawara et al. 2008;Kamioka et al. 2012;Takeishi et al. 2014;Wijeratne et al. 2016). Thus, the local strain on the osteocyte process was simulated under reciprocal fluid flow along the positive/ negative z direction (hereinafter called + z/− z flow) with a pressure gradient of |∇p| = 1.0 Pa/nm (Kamioka et al. 2012) along the z direction.

Results
The image-based fluid-structure interaction simulation reflected the flow-induced strain on the osteocyte process (Fig. 2). Although the distributions of flow velocity magnitude under + z and − z flow conditions exhibited slight variations (Fig. 2a), the distributions of TEs with high tension (Fig. 2b) and those of the maximum principal strain on the osteocyte process (Fig. 2c) were significantly dependent on the flow direction. Moreover, the strain concentration was detected in the neighborhood of the TEs with higher tension (Fig. 2b, c). These observations were common in the five models with varying random TE configurations (data not presented).
Quantitative analysis results of the osteocyte process strain revealed that the strain in the vicinity of the TEs was 1.7 times larger than that present farther from the TEs (Fig. 3a). Here, we defined the vicinity of the TEs as the area within 10.2 nm from each TE, which is equivalent to twice of the fluid lattice interval and a quarter of the representative length of TE interval on the cell process membrane. Moreover, this larger local strain was caused by TEs with higher tension (Fig. 3b top). Interestingly, the percentage of TEs with high tension, primarily contributing to the generation of a large local strain, was extremely small (Fig. 3b bottom), indicating that a small number of TEs with high tension are Table 1 Parameter settings for the fluid-structure interaction simulation Sugawara et al. 2008 The TEs with high tension that produce a strain concentration on the osteocyte process were identified by plotting the tension of TEs against their orientation angle under noflow condition from the positive z direction (Fig. 4). As shown in Fig. 4a, TEs with a large value of under the + z flow condition-TEs aligned upstream on the osteocyte process-produced high tension. On the contrary, TEs with a small under the − z flow condition generated high tension (Fig. 4b). Upon focusing on individual TEs, the TE that was aligned upstream in the + z flow ( > ∕2 ) generated a strain concentration on the osteocyte process under the + z flow condition (Fig. 4c), whereas the TE inclined upstream to the − z flow ( < ∕2 ) generated a strain concentration under the − z flow condition (Fig. 4d). These results indicate that the relative angle of the TEs to the flow direction, which can be determined based on the irregular shapes of osteocyte processes and canaliculi, is a crucial factor in producing a strain concentration in the osteocyte process.

Discussion
The influences of the ultrastructure of osteocyte processes and canaliculi on the mechanical stimuli to osteocytes were investigated by evaluating the flow-induced strain on the osteocyte processes through image-based fluid-structure interaction simulation (Fig. 1). Consequently, the simulation results revealed that specific TEs with high tension, which depend on the irregular shapes of osteocyte processes and canaliculi, can generate a local strain concentration in the osteocyte process. The local large strain induced by such TEs can potentially activate mechanical responses of osteocytes by causing structural changes in ion channels on the neighboring cell membrane (Hayakawa et al. 2008). Therefore, the strain concentration on the osteocyte process via TEs is a vital mechanism that is responsible for osteocyte mechanosensing.
Previous studies have theoretically proposed the mechanism of strain amplification in osteocyte mechanosensing, in which the fluid drag on the pericellular matrix induces a large strain on the cell process via TEs (You et al. 2001;Han et al. 2004;Wang et al. 2007). Although this strain amplification mechanism can infuse the gap between cellular and tissue level strains, it was derived using idealized models of osteocyte process and canaliculus, regardless of their irregular shapes (Kamioka et al. 2012). The high-resolution image-based model (Fig. 1) developed using UHVEM tomography (Kamioka et al. 2012;Hosaki-Takamiya et al. 2016) was utilized for conducting fluid-structure interaction simulations. This approach facilitated the evaluations of the distribution and magnitude of local strain on the osteocyte process, considering the heterogeneous inclination of TEs.
The image-based fluid-structure interaction simulation results reveal that a small number of TEs with high tension induced local strain concentrations that were heterogeneously distributed in the osteocyte process. These findings Normalized tension of TEs were commonly observed regardless of the variation in TE configurations, suggesting that TE-induced strain concentration in the osteocyte process is a universal feature associated with osteocyte mechanosensing. In addition, mechanotransductive ion channels on the cell membrane, such as Piezo1, are considered important for osteocyte mechanosensing as they respond to the neighboring membrane strain (McMahon et al. 2008;Li et al. 2019;Diem et al. 2020;Sasaki et al. 2020;Zhou et al. 2020). Thus, the locally concentrated strain on the osteocyte process, which was observed in this study, can cause structural changes in the ion channels as an activation of mechanosensing. Application of local deformation has been demonstrated to induce calcium transients in the vicinity of the stimulated point and cause diffusive wave propagation of the calcium transient to the entire intracellular region (Adachi et al. 2009). Moreover, the specific TEs contributing to the osteocyte strain concentration, which appeared to be inclined toward the upstream interstitial fluid flow in the osteocyte process (Fig. 4), were determined based on the irregular shape of osteocyte processes and canaliculi (Fig. 1b). Therefore, the transformation of the surface structure pertaining to the osteocyte process and canaliculus can considerably influence the distribution and magnitude of the osteocyte strain. This signifies that bone diseases such as osteoporosis that are accompanied by alterations in the lacuna-canalicular structure (Knothe Tate et al. 2004) are closely associated with the abnormal mechanical stimuli transmitted to osteocytes owing to the variation in the distribution of TEs in canaliculi.
One of the limitations of this analysis is that TEs were assumed to be randomly arranged in the canaliculi owing to lack of information regarding their distribution. Thus, future development of imaging technology would enable visualization of the osteocyte-specific distribution of TEs in canaliculi. In context, the image-based fluid-structure interaction simulation considering actual TE distribution will help identify the molecules responsible for osteocyte mechanosensing. This can be achieved by comparing the obtained strain distribution with the location of candidate molecules, which can be detected using techniques such as structured illumination microscopy-based super-resolution microscopy (Kiuchi et al. 2015;Cabahug-Zuckerman et al. 2018). An additional limitation of the present analysis is the restricted size of the analyzed region in relation to the canalicular network scale, which originates from the limitation posed by the electron transmissibility of UHVEM. In a bone tissue, numerous osteocytes form a complex network via their branching processes, which can change the flow velocity inside each canaliculus (van Tol et al. 2020a;. Therefore, the strain concentration regions in the osteocyte network can be clarified with an analysis considering the network structure. In context, focused ion beam-scanning electron microscopy (FIB-SEM) is a possible solution for ior of c TE with > ∕2 (arrow head) and d TE with < ∕2 (arrow head) under no-flow, + z flow, and − z flow conditions. Osteocyte process strain and TEs tension are also depicted. Tension is normalized by its median (0.017 pN) visualizing the network structure at a high resolution (Hashimoto et al. 2017;Hasegawa et al. 2018;Hirashima et al. 2020). Fluid-structure interaction simulation using the FIB-SEM image-based model will enable the evaluation of the impact of the network structure on mechanical stimuli to osteocytes.
In summary, results of the simulation conducted in this study reveal that local strain concentration in the osteocyte process was induced by a small number of TEs with high tension, which are determined by the irregular shapes of osteocyte processes and canaliculi. These results provide meaningful insights into the effect of the ultrastructure of osteocyte processes and canaliculi on osteocyte mechanical stimuli. We anticipate that this study will help in revealing molecular mechanisms that activate osteocyte mechanosensing and understanding the regulatory mechanisms of bone remodeling.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.