Next Article in Journal
Nanocurcumin Release from Self-Cured Acrylic Resins; Effects on Antimicrobial Action and Flexural Strength
Next Article in Special Issue
Determining the Relationship between Mechanical Properties and Quantitative Magnetic Resonance Imaging of Joint Soft Tissues Using Patient-Specific Templates
Previous Article in Journal
Empowering Hand Rehabilitation with AI-Powered Gesture Recognition: A Study of an sEMG-Based System
Previous Article in Special Issue
Finite Element Modeling of Residual Hearing after Cochlear Implant Surgery in Chinchillas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model

1
Department of Mechanical and Aerospace Engineering, “Sapienza” University of Rome, Via Eudossiana, 18, 00184 Rome, Italy
2
Department of Biomedical Sciences, University of Sassari, Viale San Pietro, 43/B, 07100 Sassari, Italy
3
Research Unit of Orthopaedic and Trauma Surgery, Fondazione Policlinico Universitario Campus Bio-Medico, Via Alvaro del Portillo, 200, 00128 Rome, Italy
4
Research Unit of Orthopaedic and Trauma Surgery, Department of Medicine and Surgery, Università Campus Bio-Medico di Roma, Via Alvaro del Portillo, 21, 00128 Rome, Italy
*
Author to whom correspondence should be addressed.
Bioengineering 2023, 10(5), 558; https://doi.org/10.3390/bioengineering10050558
Submission received: 31 March 2023 / Revised: 27 April 2023 / Accepted: 4 May 2023 / Published: 7 May 2023
(This article belongs to the Special Issue Computational Biomechanics)

Abstract

:
Bone tissue is mainly composed at the nanoscale of apatite minerals, collagen molecules and water that form the mineralized collagen fibril (MCF). In this work, we developed a 3D random walk model to investigate the influence of bone nanostructure on water diffusion. We computed 1000 random walk trajectories of water molecules within the MCF geometric model. An important parameter to analyse transport behaviour in porous media is tortuosity, computed as the ratio between the effective path length and the straight-line distance between initial and final points. The diffusion coefficient is determined from the linear fit of the mean squared displacement of water molecules as a function of time. To achieve more insight into the diffusion phenomenon within MCF, we estimated the tortuosity and diffusivity at different quotes in the longitudinal direction of the model. Tortuosity is characterized by increasing values in the longitudinal direction. As expected, the diffusion coefficient decreases as tortuosity increases. Diffusivity outcomes confirm the findings achieved by experimental investigations. The computational model provides insights into the relation between the MCF structure and mass transport behaviour that may contribute to the improvement of bone-mimicking scaffolds.

1. Introduction

Bone is a mineralized connective tissue with a complex arrangement of structures spanning from the nano- to macro-scale. The structural organization of bone is characterized at the nanoscale by basic building blocks composed of type I collagen, apatite minerals, water and a small percentage of non-collagenous proteins and proteoglycans. Due to the complex architecture of bone tissue [1], current information concerning the correlation between the spatial organization of bone main nanocomponents and the mechanical properties, metabolic functions and transport phenomena is still elusive [2]. In fact, analysis of the cause–effect relation between bone nanostructure and mass transport properties is a topic in constant evolution [3,4,5].
The importance of diffusion for cell metabolism and the activity of proteins is well established [5]. In the context of bone nanostructure, water is responsible for the diffusion of nutrients, metabolic exchange and ion transport, and contributes to the mechanism of bone remodelling, adaptation and repair. Water is also involved in the nucleation process of apatite crystals, participating in the stabilization of mineral structure and mediating the interaction between mineral and collagen [6]. Water movement also allows bone tissue to withstand stress with less deformation [7]. The water content in bone, i.e., 10–25 wt% [8], varies with age and disease state, decreasing with increasing mineralization [9]. This can lead to an increase in bone stiffness of up to 15% for dehydrated tissue, characterized by a high mineralization condition of 40 wt% mineral, as estimated by [3].
In bone, the metabolic traffic and interchange of signalling molecules, physiological solutes and fluids are strongly dependent on the transport pathways, comprising, at the smallest hierarchical structural level, the interconnected pores within the apatite–collagen matrix [8]. The scientific literature [8,10,11] agrees that water occupies different porosity levels of bone tissue, which are nested hierarchically inside one another: collagen–apatite (10 nm), canalicular (100 nm), lacunar (up to 8 μm), vascular (50 μm), and the intertrabecular porosity (up to 1 mm).
The diffusion phenomenon within bone tissue has been studied at different length scales. At the microscale, water diffusion has been investigated in cortical bone specimens, by means of nuclear magnetic resonance [12,13]. Although this technique can non-destructively assess the water distribution in bone tissue, its resolution is limited. Subsequently, at the lacunar canalicular level, i.e., up to 8 μm, the role of fluid phase has been analysed by means of experimental and computational methods [6,14,15], especially in relation to its mechanotransduction functions. Nonetheless, a direct link between interstitial fluid flow and bone adaptation needs to be more completely established [13].
At the nanoscale porosity, initial studies have indicated that water is present as bound water. However, in the last decade, studies have introduced the hypothesis that free water may also exist at this level. Marinozzi and co-workers [16,17] performed an experimental analysis of diffusion within a single trabecula, which led to new insights for bone tissue characterization up to the length scale of its principal constituents, i.e., the collagen matrix and apatite mineral.
Within this framework, computational models have been developed as a powerful engineering tool to analyse and extend experimental findings [18]. In the literature, different in silico studies have been developed to mimic the diffusion phenomenon. With regard to the biological field, diffusion has been investigated in soft tissues such as brain [19,20] and muscle [21] tissues, and in hard mineralized tissues [22,23,24]. To date, most studies have considered 2D models, which may lead to limited predictive outcomes [19,20]. Another crucial issue lies in the difficulty of developing a method able to estimate the effects of the local variations in geometry [19].
Overall, the assessment of the diffusion coefficient can be based on empirical correlations with pore structure or on theoretical models with idealized geometry. According to [25] the influence of the structure is taken into account by means of three factors: the void space inside the porous medium, i.e., porosity, the path length that the particle has to travel in order to cross the medium, i.e., tortuosity, and the reduction of the effective flow due to the changes in the cross section of the pores, i.e., constrictivity. As highlighted by [26], the tortuosity is the main parameter that most influences the diffusion coefficient within the MCF. Tortuosity (τ) is a quantitative measure of the reduction in diffusive flux caused by the sinuous path imposed by the obstacles, compared to the straightest path in an unrestricted medium, in the direction of the flow [27]. Overall, it is influenced by the shape, arrangement and inclinations of the obstacles.
Although the existence of tortuosity and constrictivity factors in the assessment of diffusivity is recognized, a general method to extract them from experimental investigations has not been determined. The literature proposes different approaches for estimating these parameters, such as diffusion experiments [16], ultrasonic reflectivity methods [28], NMR measurements [29], analytical models [30,31] and 3D image analysis [32]. Unfortunately, any straightforward relation between the structure of the porous network and the coefficients of the diffusion process is still a topic under investigation. For instance, numerical models have been implemented to obtain a more explicit expression for tortuosity [30,31]. Overall, there are two methods presented in the literature for the calculation of tortuosity. One approach is based on the development of a model that reproduces, in detail, a region of the porous structure; then, different paths are implemented in order to determine how this structure will influence the diffusion phenomenon [31].
A second method avoids the issues related to the implementation of single paths by mimicking the diffusion process itself, e.g., the random walk approach [33]. The random walk (RW) technique represents a common method used for the investigation of mass transport dynamics, and also in biological systems [33,34,35]. The RW method describes the diffusion phenomenon in a simple manner, without the need for more complex operations [36,37], and reveals features that could be difficult to discern with other approaches. It represents, therefore, an appropriate method to investigate bone nanostructure, since experimental investigations at nanoscale are still challenging. Trajectories developed in RW models are defined as unbiased, meaning that the particle is equally likely to move in each direction, and is uncorrelated in direction, i.e., the direction of motion at a given time is independent of the directions at all preceding times [36,38].
This study aims to achieve information about the influence of the collagen–apatite structure on water diffusion at the nanoscale. The analysis is performed considering a 3D random walk model of water diffusion within the MCF, and represents an extension and improvement of a previous study [39]. In fact, we provide here a thorough analysis of the diffusion coefficient and tortuosity factor in the longitudinal direction of the fibril.

2. Materials and Methods

A 3-D model of the embedded apatite crystals within the collagen fibril has been developed according to the model proposed by Petruska et al. [40] and Jäger et al. [41]. The dimensions of the MCF are still debated in the literature. In this study, we considered a MCF diameter of 200 nm, following [42], while the MCF length is approximately 1000 nm [43].
Type I collagen is composed of a triple helical polypeptide chain molecule referred to as tropocollagen molecule. The collagen molecules are cylindrically shaped with an approximate diameter of 1.23 nm and a length of roughly 300 nm [44].
Apatite minerals are described mainly as platelet-like-shaped, with a low dispersion in thickness, i.e., 2–5 nm, but a wider spread in length (50–170 nm) and width (5–90 nm) [45,46]. The dimensions of apatite platelets have been obtained with random extractions from Gaussian probability distribution functions in the ranges indicated in the literature [45,46]. In Table 1, we report the average dimensions of the apatite crystals considered in this model. In the width and thickness directions, the distances between neighbouring apatite platelets, i.e., aW and aT, respectively, are of the same order as crystal thickness, as achieved from experimental TEM observations [47]. In the longitudinal direction, the mineral platelets are arranged in a staggered pattern, with a periodicity of roughly 67 nm [40,41]. In the length direction, the distance between two consecutive apatite crystals, i.e., aL, is determined according to Equation (1) [41]:
+ a L 2 = D p e r i o d
where ℓ denotes the length of the apatite mineral, and Dperiod is the length of the axial period, i.e., Dperiod = 67 nm. Vercher-Martinez and co-workers [48] suggested parallel layers in the equatorial plane, i.e., the normal plane to the longitudinal direction.
The geometrical model represents a MCF with a mineralization degree of 32%, in agreement with previous studies [41,49]. The RW model is implemented on an MCF configuration achieved after roughly 6·106 random inclinations and rotations of the platelets. We performed random displacements and rotations for each apatite platelet in order to better mimic the structure of bone tissue at the nanoscale. We started from a perfectly aligned configuration of apatite platelets and performed random perturbations to achieve inclinations for each apatite crystal, in agreement with the range available in the literature [50,51], i.e., within ±20 degrees. The method is described in detail in a previous study [42].
The RW model mimics the trajectory of water particles, taking into account the arrangement of the collagen–apatite structure (Figure 1). The implementation of the RW model follows a previous study by Bini and colleagues [39]. To adapt the mathematical model of RW to the spatial configuration of collagen–apatite matrix, we have assumed that apatite platelets are impermeable, and we introduced a reflection condition concerning collagen fibrils, as implemented in previous studies [33]. We also hypothesise that particle–particle interactions are negligible, since no saturation phenomenon occurs during the initial period of the diffusion process [34,37]. Moreover, in the RW model, the time interval Δt is set to 2·10−10 s, so as to ensure that the length of each displacement is smaller than the minimum dimensions of the structural obstacles, i.e., the cross-section of tropocollagen molecules [34,39]. In fact, the displacement δj, in each coordinate direction j, i.e., j = W, T, L, is sampled from a Gaussian distribution with zero mean and standard deviation, calculated as follows:
σ = 2 · D 0 · t
where D0 = 2.66·10−9 m2·s−1 is the water diffusivity in a homogeneous medium at 27 °C [52]. This requirement was implemented to avoid the water particles skipping the obstacles, i.e., collagen fibrils and apatite crystals, which would have led to a distorted estimation of the diffusion coefficient.
The RW simulation starts by randomly positioning the water particles in the equatorial plane at the lower extremity of the MCF, i.e., coordinate L = 0 nm. Then, the particle randomly selects the coordinates of the next position. If the particle does not encounter any obstacle, it moves to the new position; otherwise, the particle stays at the current location and new coordinates are extracted randomly. This random walking procedure is repeated until the particle reaches the upper extremity of the MCF. According to the method proposed by [37], the diffusion coefficient is proportional to the slope of the linear fit of the mean square displacement (MSD) of the water particle versus time (τ).
D = M S D ( τ ) 2 · d S · τ
where dS = 3 is the system dimension. Thus, for 3D models, D equals one-sixth of the slope of the linear fit of the particle MSD as a function of τ. The MSD of the single trajectory (MSDi) is calculated as follows [53]:
M S D i n τ = 1 N n · j = 0 N n r j τ + n τ r j τ 2
where r (jτ), and r (jτ + nτ), are the particle coordinates at time (jτ) and (nτ + jτ), respectively, and n is the time lag, i.e., n = 1, 2, …, N − 1.
The least-squared line fit of the single MSD curve (MSDi) is expressed by Equation (5) [54]:
M S D i = β 0 + β 1 · t
where β0, β1 are the intercept and the slope of the linear fit. The diffusion coefficient (Di) of the individual trajectories is obtained from Equations (3) and (5):
D i = 1 6 β 1
In the available literature concerning the 3D random walk method applied to biological tissues, the number of trajectories is extremely variable, including analyses based on 30 [55], 100 [56] or 350 [33] trajectories. Saxton [53] highlighted that the numerical values of the diffusion coefficient are reliable when considering a number of trajectories (NT) greater than 100. In this study, we considered 1000 trajectories.
Tortuosity is an important parameter in the analysis of the diffusion phenomenon. Here, the tortuosity factor is evaluated using the random walk simulations. In this context, tortuosity is a measure of the diffusive nature of the random walk in relation to the complexity of the MCF anisotropic structure. The tortuosity of the single trajectory is defined as the ratio between the effective length of the particle paths and the Euclidean distance between the initial and final points [27]:
T o r t u o s i t y i = E f f e c t i v e   L e n g t h   ( i ) E u c l i d e a n   d i s t a n c e   ( i ) = = j = 0 N 1 r W j τ + τ r W ( j τ ) 2 + r T j τ + τ r T ( j τ ) 2 + r L j τ + τ r L ( j τ ) 2 r W e n d r W s t a r t 2 + r T e n d r T s t a r t 2 + r L e n d r L s t a r t 2
where rW (jτ + τ), rT (jτ + τ), rL (jτ + τ) are the coordinates in direction W, T, L of the particle position at time (jτ + τ), and rW (jτ), rT (jτ), rL (jτ) of the particle position at time (jτ). The terms rW start, rT start rL start represent the coordinates of the initial point of the trajectory and rW end, rT end rL end are the coordinates of the final point of the path.
Since the MCF structure is characterized by high complexity, in this study, we aim to provide an analysis of the variation of tortuosity and diffusivity along the MCF. We considered 7 reference quotes in the L direction, with a step equal to the length of the MCF building block, i.e., 134 nm, in accordance with [40,41].
In order to analyse the influence of the structure, we considered sub-trajectories of the random paths that span the MCF, each starting from the inferior extremity, i.e., L = 0 nm, and reaching the equatorial planes defined by the reference quotes considered in the L direction. For each sub-trajectory, we calculated the diffusion coefficient and the tortuosity by means of Equation (6) and Equation (7), respectively.
Moreover, for each reference quote, we estimated the average tortuosity and average diffusivity from all sub-trajectories, found using the following equations:
T o r t u o s i t y = 1 N s T i = 1 N s T T o r t u o s i t y i       = 1 N s T i = 1 N s T j = 0 N 1 r W j τ + τ r W ( j τ ) 2 + r T j τ + τ r T ( j τ ) 2 + r L j τ + τ r L ( j τ ) 2 r W e n d r W s t a r t 2 + r T e n d r T s t a r t 2 + r L e n d r L s t a r t 2
D = 1 N s T i = 1 N s T D i
where NsT is the number of sub-trajectories.

3. Results

In this study, we developed a computational model simulating the 3D random walk of water particles in the regions confined by apatite platelets and tropocollagen molecules, in order to provide a thorough analysis of the diffusion phenomenon along the MCF. The effect of the structural organization of apatite minerals and collagen on the diffusivity is highlighted for different reference quotes in the longitudinal direction. We estimated the tortuosity and the diffusivity for random walk sub-trajectories that reach the aforementioned reference quotes. In Figure 2, we report the diffusion coefficient as a function of the tortuosity, and the corresponding fit for two extreme cases, characterized by the lowest (Figure 2c) and highest (Figure 2b) values of tortuosity. We also calculated the 95% confidence interval for the diffusion coefficient values estimated for the two quotes.
According to its definition (Equation (7)), the tortuosity is characterized by values greater than 1. From Figure 2, it emerges that for the reference quote sited at the lower extremity of the MCF, i.e., L = 134 nm, the tortuosity assumes modest values, with a maximum of 30. However, it can be noticed that at this reference quote, the majority of sub-trajectories are characterized by low tortuosity, with values up to 5. Conversely, at the upper position of the MCF, i.e., L = 938 nm, the maximum of tortuosity increases up to six times with respect to the maximum value estimated at the MCF lower extremity, reaching values up to 180. Differing from quote L = 134 nm, at the upper quote L = 938 nm, the majority of sub-trajectories span a wider interval of values, i.e., from 40 to 160.
Moreover, the analysis of the outcomes reported in Figure 2 shows that the diffusion coefficient is inversely proportional to the tortuosity. For instance, Figure 2c highlights that for lower values of tortuosity, the diffusivity is characterized by a wider interval of values, spanning from roughly 4.5·10−10 (m2·s−1) to 2.5·10−9 (m2·s−1). Specifically, a significant variation in the diffusion coefficient is estimated to correspond with the modest values of tortuosity, i.e., up to 5, while for higher values of tortuosity, the decrease in the diffusivity coefficient is characterized by a smooth slope.
Conversely, when considering the longer sub-trajectories that reach the upper quote and are thus characterized by higher tortuosity (Figure 2b), the diffusivity interval is narrower, i.e., from 1·10−10(m2·s−1) to 1.5·10−9 (m2·s−1). Consequently, at the upper reference quote L = 938 nm, the diminution of the diffusion coefficient for increasing tortuosity can be described as a modest slope.
We further developed the analysis of local diffusion phenomenon by calculating the mean values of diffusivity and tortuosity for each reference quote considered (Figure 3), applying Equation (8) and Equation (9), respectively. Overall, the average tortuosity shows an increasing trend for random walk sub-trajectories of a higher extension. It can be noticed that for sub-trajectories that reach the upper reference quote L = 938 nm, the average tortuosity increases noticeably, roughly ten times with respect to the case characterized by shortest sub-trajectories, i.e., reference quote L = 134 nm. It is worth pointing out that in this study, tortuosity is an indicator of the complexity of the random walk path of the water particle, which is influenced by the spatial arrangement of the apatite platelets and the collagen matrix. As expected, for sub-trajectories that span roughly the whole length of MCF, the average tortuosity is characterized by higher values due to the increased occurrence of obstacles, i.e., mineral platelets and tropocollagen, along the MCF.
Consequently, the average diffusion coefficient is characterized by a decreasing trend with the increasing length of the random walk sub-trajectories. Specifically, the diffusivity diminishes for longer sub-trajectories (by roughly 20% with respect to the value estimated for the reference quote characterized by shortest sub-trajectories, i.e., L = 134 nm).

4. Discussion

Water plays a fundamental role in the provision of nutrients and the interaction between apatite crystals and protein constituents [6], enhancing the strength of the mineral phase within the MCF. An impairment of the diffusion process may also affect the bone-remodelling process [11], since mass transport contributes actively to the signalling functions that are involved in tissue growth.
Experimental analysis of the diffusion phenomenon at the bone nanoscale is still limited in the literature.
Thus, in this study, we used computational modelling to provide a more detailed analysis of the water diffusion phenomenon within the entire MCF. In fact, computational modelling can simulate the physical, chemical and biological aspects of the real system and study their interactions from a cause–effect approach. Moreover, the model developed here allows us to analyse locally, at different quotes in the longitudinal direction, how the random walk of water particles is influenced by the structural constraints derived from the spatial arrangement of apatite platelets and the collagen matrix.
It should be noted that previous studies [26,50] have investigated diffusivity, mostly considering the building blocks of the MCF. The potentiality of the RW technique with regard to the diffusion coefficient has been observed in a former study [39], in which RW simulations were applied to investigate the influence of the mineralization degree on the diffusion coefficient. In the study of Bini and colleagues [39], the computational model is used to investigate different mineral contents, i.e., low (7%), intermediate (32%), and high (42%) volume fractions. The outcomes of [39] highlight that an increase in the mineral content within the MCF leads to a diminution of the diffusion coefficient. In addition to the results presented in [39], we focused here on the local characterization of the diffusivity within the MCF characterized by 32% mineralization, which represents a physiological value of the MCF mineral content [26,41,42,49]. Moreover, the computational model allowed us to perform an analysis of the tortuosity of water particles’ trajectories, and their local variation within the entire MCF.
The collagen–apatite structure influences both tortuosity and the diffusion coefficient. The diffusivity is significantly lower than the free water diffusion coefficient because of the steric constraints of the apatite mineral configuration and collagen. Overall, the diffusion coefficient is characterized by an order of magnitude of 10−10 (m2·s−1), in good agreement with previous outcomes from experimental [16] and computational investigations [22]. Marinozzi et al. [16] estimated a diffusion coefficient of (4.12 ± 0.8) 10−10 (m2·s−1) from experimental measurement during the swelling of a single human trabeculae. Lemaire et al. [22] investigated using molecular dynamics water flow within the collagen–apatite network and achieved a diffusivity value of 8·10−10 m2·s−1. It is worth pointing out that studies based on the molecular dynamics method take into account a limited region of the MCF due to the high computational cost of this technique. Conversely, the RW model presented here considers the entire MCF in its assessment of water diffusion coefficient. The results are also in good agreement with previous studies that have analysed the diffusion phenomenon in bone tissue from animal samples. For instance, Gul-E-Noor et al. [57] reported a water diffusivity of 6.2 ± 0.7·10−10 m2·s−1 from measurements in goat cortical bone. In the work of Wang et al. [58], the fluorescein diffusion coefficient at the lacunar-canalicular level of mouse bone samples was determined to be 3.3 ± 0.6·10−10 m2·s−1.
Moreover, it is worth pointing out that local analysis of the tortuosity and diffusivity at different quotes in the longitudinal direction of the MCF allows us to achieve a first assessment of the obstructive effect related to apatite platelets and the collagen matrix. As expected, longer sub-trajectories are characterized by higher values of tortuosity and a lower diffusion coefficient, since these parameters are affected by the spatial arrangement of the MCF main nanocomponents. This study focused on the analysis of the random trajectories developed within the MCF. However, a more detailed analysis and quantification of the obstructive factor that characterizes the nanoconfined structure of the MCF will be the subject of future work.
The outcomes that emerged from this analysis related to the mass transport behaviour within the MCF nanostructure may facilitate the development of appropriate synthetic systems. Biomimetic scaffolds for bone should match the properties of the original tissue, such as osteoinductivity, osteoconductivity, biocompatibility, suitable porosity and pore interconnectivity, in order to allow the diffusion of cells, growth factors, nutrients, and waste products [59]. In addition, the appropriate architecture of the scaffold should be developed in order to meet the mechanical demand of the tissue environment. Application of nanotechnology [60], high-performance biomaterials [61,62,63] and computational models [24,26,39] can contribute to avoiding cell death due to insufficient nutrient transport, inadequate integration of the regenerated tissue with the surrounding native tissue, and mismatching of scaffold properties with respect to the host tissue.
Random walk is a mathematical model commonly used to describe the diffusion of molecules or particles in biological tissues that allows us to investigate the influence of structural hindrance on diffusivity. However, some limitations should be highlighted. We considered an ideal geometry for tropocollagen and apatite crystals. Tropocollagen is modelled with long cylinders, while mineral platelets are presented as regular, thin parallelepipeds. This strategy allowed us to represent the relevant characteristics of bone nanostructure whilst also maintaining a low computational cost. In this work, we assumed that the water particles start their path from the lower extremity of the MCF, i.e., L = 0 nm. A future study will also consider the propagation of water particles from the lateral surface of the MCF in order to mimic immersion in water [16].

5. Conclusions

The properties of bone are closely related to its hierarchical structure. The interplay between the main constituents of bone at the nanoscale affects its structural characteristics and its mass transport properties. Information concerning the behaviour and interactions of water at the nanostructure level of bone tissue is limited in the literature. In this study, we analysed tortuosity and diffusivity along the single MCF in order to predict the local variation of these parameters. The agreement between the numerical values of the diffusion coefficient determined by the RW model and the literature confirms the suitability of the geometric dimensions and spatial arrangement considered for apatite platelets and collagen. In addition, the computational model provides insights into the MCF structure and the cause–effect relationship with mass transport properties, thereby enhancing our understanding of bone mineralization.
The outcomes presented here may contribute to improving the design of smart structural nanomaterials and biomimetic scaffolds. Nowadays, the optimization of bioscaffolds mimicking the hierarchical structure of bone represents a well-known challenge. The RW model can be applied in combination with experimental investigations, and could provide new insights for more accurate tissue characterization, which could be of use in controlling the mechanical and biological performance of bioscaffolds.

Author Contributions

Conceptualization, F.B. and F.M.; methodology, F.B., A.P. and F.M.; software, A.P.; validation, F.B. and A.P.; formal analysis, F.B. and A.P.; investigation, F.B. and A.P.; resources, F.B. and F.M.; data curation, F.B. and A.P.; writing—original draft preparation, F.B. and A.P.; writing—review and editing, F.B., A.P., A.M. and F.M.; visualization, F.B. and A.P.; supervision, F.B., A.M. and F.M.; project administration, F.B.; funding acquisition, not applicable. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Niu, Y.; Du, T.; Liu, Y. Biomechanical Characteristics and Analysis Approaches of Bone and Bone Substitute Materials. J. Funct. Biomater. 2023, 14, 212. [Google Scholar] [CrossRef]
  2. Liebi, M.; Georgiadis, M.; Menzel, A.; Schneider, P.; Kohlbrecher, J.; Bunk, O.; Guizar-Sicairos, M. Nanostructure Surveys of Macroscopic Specimens by Small-Angle Scattering Tensor Tomography. Nature 2015, 527, 349–352. [Google Scholar] [CrossRef]
  3. Fielder, M.; Nair, A.K. Effects of Hydration and Mineralization on the Deformation Mechanisms of Collagen Fibrils in Bone at the Nanoscale. Biomech. Model. Mechanobiol. 2019, 18, 57–68. [Google Scholar] [CrossRef]
  4. Chen, X.; Qian, T.; Hang, F.; Chen, X. Water Promotes the Formation of Fibril Bridging in Antler Bone Illuminated by in Situ AFM Testing. J. Mech. Behav. Biomed. Mater. 2021, 120, 104580. [Google Scholar] [CrossRef] [PubMed]
  5. Surowiec, R.K.; Allen, M.R.; Wallace, J.M. Bone Hydration: How We Can Evaluate It, What Can It Tell Us, and Is It an Effective Therapeutic Target? Bone Rep. 2022, 16, 101161. [Google Scholar] [CrossRef]
  6. Wang, Y.; Von Euw, S.; Fernandes, F.M.; Cassaignon, S.; Selmane, M.; Laurent, G.; Pehau-Arnaudet, G.; Coelho, C.; Bonhomme-Coury, L.; Giraud-Guille, M.M.; et al. Water-Mediated Structuring of Bone Apatite. Nat. Mater. 2013, 12, 1144–1153. [Google Scholar] [CrossRef] [PubMed]
  7. Rai, R.K.; Sinha, N. Dehydration-Induced Structural Changes in the Collagen-Hydroxyapatite Interface in Bone by High-Resolution Solid-State NMR Spectroscopy. J. Phys. Chem. C 2011, 115, 14219–14227. [Google Scholar] [CrossRef]
  8. Cowin, S.C. Bone Mechanics Handbook, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2001; ISBN 9780429125447. [Google Scholar]
  9. Burr, D.B. Changes in Bone Matrix Properties with Aging. Bone 2019, 120, 85–93. [Google Scholar] [CrossRef]
  10. Maghsoudi-Ganjeh, M.; Wang, X.; Zeng, X. Computational Investigation of the Effect of Water on the Nanomechanical Behavior of Bone. J. Mech. Behav. Biomed. Mater. 2020, 101, 103454. [Google Scholar] [CrossRef]
  11. Fritton, S.P.; Weinbaum, S. Fluid and Solute Transport in Bone: Flow-Induced Mechanotransduction. Annu. Rev. Fluid Mech. 2009, 41, 347–374. [Google Scholar] [CrossRef]
  12. Fernández-Seara, M.A.; Wehrli, S.L.; Wehrli, F.W. Diffusion of Exchangeable Water in Cortical Bone Studied by Nuclear Magnetic Resonance. Biophys. J. 2002, 82, 522–529. [Google Scholar] [CrossRef] [PubMed]
  13. Ni, Q.; Nyman, J.S.; Wang, X.; De Los Santos, A.; Nicolella, D.P. Assessment of Water Distribution Changes in Human Cortical Bone by Nuclear Magnetic Resonance. Meas. Sci. Technol. 2007, 18, 715–723. [Google Scholar] [CrossRef]
  14. Samuel, J.; Park, J.; Almer, J.; Wang, X. Effect of Water on Nanomechanics of Bone Is Different between Tension and Compression. J. Mech. Behav. Biomed. Mater. 2016, 57, 128–138. [Google Scholar] [CrossRef] [PubMed]
  15. Van Tol, A.F.; Schemenz, V.; Wagermaier, W.; Roschger, A.; Razi, H.; Vitienes, I.; Fratzl, P.; Willie, B.M.; Weinkamer, R. The Mechanoresponse of Bone Is Closely Related to the Osteocyte Lacunocanalicular Network Architecture. Proc. Natl. Acad. Sci. USA 2020, 117, 32251–32259. [Google Scholar] [CrossRef] [PubMed]
  16. Marinozzi, F.; Bini, F.; Marinozzi, A. Water Uptake and Swelling in Single Trabeculæ from Human Femur Head. Biomatter 2014, 4, e28237. [Google Scholar] [CrossRef]
  17. Marinozzi, F.; Bini, F.; Quintino, A.; Corcione, M.; Marinozzi, A. Experimental Study of Diffusion Coefficients of Water through the Collagen: Apatite Porosity in Human Trabecular Bone Tissue. Biomed Res. Int. 2014, 2014, 796519. [Google Scholar] [CrossRef]
  18. Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. A Nanostructured Look of Collagen Apatite Porosity into Human Mineralized Collagen Fibril. Biocell 2022, 46, 2225–2229. [Google Scholar] [CrossRef]
  19. Hrabe, J.; Hrabětová, S.; Segeth, K. A Model of Effective Diffusion and Tortuosity in the Extracellular Space of the Brain. Biophys. J. 2004, 87, 1606–1617. [Google Scholar] [CrossRef]
  20. Nicholson, C. Diffusion and Related Transport Mechanisms in Brain Tissue. Rep. Prog. Phys. 2001, 64, 815–884. [Google Scholar] [CrossRef]
  21. Carneiro, I.; Carvalho, S.; Henrique, R.; Oliveira, L.M.; Tuchin, V.V. A Robust Ex Vivo Method to Evaluate the Diffusion Properties of Agents in Biological Tissues. J. Biophotonics 2019, 12, e201800333. [Google Scholar] [CrossRef]
  22. Lemaire, T.; Pham, T.T.; Capiez-Lernout, E.; de Leeuw, N.H.; Naili, S. Water in Hydroxyapatite Nanopores: Possible Implications for Interstitial Bone Fluid Flow. J. Biomech. 2015, 48, 3066–3071. [Google Scholar] [CrossRef] [PubMed]
  23. Bini, F.; Pica, A.; Novelli, S.; Marinozzi, A.; Marinozzi, F.; Fabiano Bini, Andrada Pica, Simone Novelli, Andrea Marinozzi, F. 3D-FEM Modeling of Iso-Concentration Maps in Single Trabecula from Human Femur Head. In VipIMAGE 2019. VipIMAGE 2019. Lecture Notes in Computational Vision and Biomechanics; Tavares, J.M.R.S., Jorge, R.M.N., Eds.; Springer: Berlin/Heidelberg, Germany, 2019; Volume 34, ISBN 978-3-030-32039-3. [Google Scholar]
  24. Bini, F.; Pica, A.; Novelli, S.; Pecci, R.; Bedini, R.; Marinozzi, A.; Marinozzi, F. 3D FEM Model to Simulate Brownian Motion inside Trabecular Tissue from Human Femoral Head. Comput. Methods Biomech. Biomed. Eng. Imaging Vis. 2022, 10, 500–507. [Google Scholar] [CrossRef]
  25. van Brakel, J.; Heertjes, P.M. Analysis of Diffusion in Macroporous Media in Terms of a Porosity, a Tortuosity and a Constrictivity Factor. Int. J. Heat Mass Transf. 1974, 17, 1093–1103. [Google Scholar] [CrossRef]
  26. Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. A 3D Model of the Effect of Tortuosity and Constrictivity on the Diffusion in Mineralized Collagen Fibril. Sci. Rep. 2019, 9, 2658. [Google Scholar] [CrossRef] [PubMed]
  27. Ghanbarian, B.; Hunt, A.G.; Ewing, R.P.; Sahimi, M. Tortuosity in Porous Media: A Critical Review. Soil Sci. Soc. Am. J. 2013, 77, 1461–1477. [Google Scholar] [CrossRef]
  28. Moussatov, A.; Ayrault, C.; Castagnède, B. Porous Material Characterization-Ultrasonic Method for Estimation of Tortuosity and Characteristic Length Using a Barometric Chamber. Ultrasonics 2001, 39, 195–202. [Google Scholar] [CrossRef]
  29. Rottreau, T.J.; Parlett, C.M.A.; Lee, A.F.; Evans, R. Diffusion NMR Characterization of Catalytic Silica Supports: A Tortuous Path. J. Phys. Chem. C 2017, 121, 16250–16256. [Google Scholar] [CrossRef]
  30. Matyka, M.; Khalili, A.; Koza, Z. Tortuosity-Porosity Relation in Porous Media Flow. Phys. Rev. E 2008, 78, 026306. [Google Scholar] [CrossRef]
  31. Yun, M.; Yu, B.; Xu, P.; Wu, J. Geometrical Models for Tortuosity of Streamlines in Three-Dimensional Porous Media. Can. J. Chem. Eng. 2006, 84, 301–309. [Google Scholar] [CrossRef]
  32. Wu, Y.S.; van Vliet, L.J.; Frijlink, H.W.; van der Voort Maarschalk, K. The Determination of Relative Path Length as a Measure for Tortuosity in Compacts Using Image Analysis. Eur. J. Pharm. Sci. 2006, 28, 433–440. [Google Scholar] [CrossRef]
  33. Jin, S.; Verkman, A.S. Single Particle Tracking of Complex Diffusion in Membranes: Simulation and Detection of Barrier, Raft, and Interaction Phenomena. J. Phys. Chem. B 2007, 111, 3625–3632. [Google Scholar] [CrossRef] [PubMed]
  34. Momot, K.I. Diffusion Tensor of Water in Model Articular Cartilage. Eur. Biophys. J. 2011, 40, 81–91. [Google Scholar] [CrossRef] [PubMed]
  35. Stylianopoulos, T.; Poh, M.Z.; Insin, N.; Bawendi, M.G.; Fukumura, D.; Munn, L.L.; Jain, R.K. Diffusion of Particles in the Extracellular Matrix: The Effect of Repulsive Electrostatic Interactions. Biophys. J. 2010, 99, 1342–1349. [Google Scholar] [CrossRef] [PubMed]
  36. Berg, H.C. Random Walks in Biology; Princeton University Press: Princeton, NJ, USA, 1983. [Google Scholar]
  37. Einstein, A.; Furth, R. Investigations on the Theory of the Brownian Movement; Dover Publications: New York, NY, USA, 1956. [Google Scholar]
  38. Codling, E.A.; Plank, M.J.; Benhamou, S. Random Walk Models in Biology. J. R. Soc. Interface 2008, 5, 813–834. [Google Scholar] [CrossRef]
  39. Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. 3D Random Walk Model of Diffusion in Human Hypo- and Hyper- Mineralized Collagen Fibrils. J. Biomech. 2021, 125, 110586. [Google Scholar] [CrossRef] [PubMed]
  40. Petruska, J.A.; Hodge, A.J. A Subunit Model for the Tropocollagen Macromolecule. Proc. Natl. Acad. Sci. USA 1964, 51, 871–876. [Google Scholar] [CrossRef]
  41. Jäger, I.; Fratzl, P. Mineralized Collagen Fibrils: A Mechanical Model with a Staggered Arrangement of Mineral Particles. Biophys. J. 2000, 79, 1737–1746. [Google Scholar] [CrossRef]
  42. Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. Percolation Networks inside 3D Model of the Mineralized Collagen Fibril. Sci. Rep. 2021, 11, 11398. [Google Scholar] [CrossRef]
  43. Buehler, M.J. Nature Designs Tough Collagen: Explaining the Nanostructure of Collagen Fibrils. Proc. Natl. Acad. Sci. USA 2006, 103, 12285–12290. [Google Scholar] [CrossRef]
  44. Lees, S. Considerations Regarding the Structure of the Mammalian Mineralized Osteoid from Viewpoint of the Generalized Packing Model. Connect. Tissue Res. 1987, 16, 281–303. [Google Scholar] [CrossRef]
  45. Fratzl, P.; Gupta, H.S.; Paschalis, E.P.; Roschger, P. Structure and Mechanical Quality of the Collagen-Mineral Nano-Composite in Bone. J. Mater. Chem. 2004, 14, 2115–2123. [Google Scholar] [CrossRef]
  46. Olszta, M.J.; Cheng, X.; Jee, S.S.; Kumar, R.; Kim, Y.Y.; Kaufman, M.J.; Douglas, E.P.; Gower, L.B. Bone Structure and Formation: A New Perspective. Mater. Sci. Eng. R Rep. 2007, 58, 77–116. [Google Scholar] [CrossRef]
  47. Landis, W.J.; Song, M.J.; Leith, A.; McEwen, L.; McEwen, B.F. Mineral and Organic Matrix Interaction in Normally Calcifying Tendon Visualized in Three Dimensions by High-Voltage Electron Microscopic Tomography and Graphic Image Reconstruction. J. Struct. Biol. 1993, 110, 39–54. [Google Scholar] [CrossRef] [PubMed]
  48. Vercher-Martínez, A.; Giner, E.; Arango, C.; Javier Fuenmayor, F. Influence of the Mineral Staggering on the Elastic Properties of the Mineralized Collagen Fibril in Lamellar Bone. J. Mech. Behav. Biomed. Mater. 2015, 42, 243–256. [Google Scholar] [CrossRef]
  49. Depalle, B.; Qin, Z.; Shefelbine, S.J.; Buehler, M.J. Large Deformation Mechanisms, Plasticity, and Failure of an Individual Collagen Fibril with Different Mineral Content. J. Bone Miner. Res. 2016, 31, 380–390. [Google Scholar] [CrossRef]
  50. Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. 3D Diffusion Model within the Collagen Apatite Porosity: An Insight to the Nanostructure of Human Trabecular Bone. PLoS ONE 2017, 12, e0189041. [Google Scholar] [CrossRef]
  51. Xu, Y.F.; Nudelman, F.; Eren, E.D.; Wirix, M.J.M.; Cantaert, B.; Nijhuis, W.H.; Hermida-Merino, D.; Portale, G.; Bomans, P.H.H.; Ottmann, C.; et al. Intermolecular Channels Direct Crystal Orientation in Mineralized Collagen. Nat. Commun. 2020, 11, 5068. [Google Scholar] [CrossRef]
  52. Chiavazzo, E.; Fasano, M.; Asinari, P.; Decuzzi, P. Scaling Behaviour for the Water Transport in Nanoconfined Geometries. Nat. Commun. 2014, 5, 3565. [Google Scholar] [CrossRef]
  53. Saxton, M.J. Single-Particle Tracking: The Distribution of Diffusion Coefficients. Biophys. J. 1997, 72, 1744–1753. [Google Scholar] [CrossRef]
  54. Weimann, L.; Ganzinger, K.A.; McColl, J.; Irvine, K.L.; Davis, S.J.; Gay, N.J.; Bryant, C.E.; Klenerman, D. A Quantitative Comparison of Single-Dye Tracking Analysis Tools Using Monte Carlo Simulations. PLoS ONE 2013, 8, e64287. [Google Scholar] [CrossRef]
  55. Monnier, N.; Guo, S.M.; Mori, M.; He, J.; Lénárt, P.; Bathe, M. Bayesian Approach to MSD-Based Analysis of Particle Motion in Live Cells. Biophys. J. 2012, 103, 616–626. [Google Scholar] [CrossRef] [PubMed]
  56. Huang, F.; Watson, E.; Dempsey, C.; Suh, J. Real-Time Particle Tracking for Studying Intracellular Trafficking of Pharmaceutical Nanocarriers. Methods Mol. Biol. 2013, 991, 211–223. [Google Scholar] [CrossRef] [PubMed]
  57. Gul-E-Noor, F.; Singh, C.; Papaioannou, A.; Sinha, N.; Boutis, G.S. Behavior of Water in Collagen and Hydroxyapatite Sites of Cortical Bone: Fracture, Mechanical Wear, and Load Bearing Studies. J. Phys. Chem. C 2015, 119, 21528–21537. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, L.; Wang, Y.; Han, Y.; Henderson, S.C.; Majeska, R.J.; Weinbaum, S.; Schaffler, M.B. In Situ Measurement of Solute Transport in the Bone Lacunar-Canalicular System. Proc. Natl. Acad. Sci. USA 2005, 102, 11911–11916. [Google Scholar] [CrossRef]
  59. Oosterlaken, B.M.; Vena, M.P.; de With, G. In Vitro Mineralization of Collagen. Adv. Mater. 2021, 33, 2004418. [Google Scholar] [CrossRef]
  60. Araneo, R.; Rinaldi, A.; Notargiacomo, A.; Bini, F.; Marinozzi, F.; Pea, M.; Lovat, G.; Celozzi, S. Effect of the Scaling of the Mechanical Properties on the Performances of ZnO Piezo-Semiconductive Nanowires. AIP Conf. Proc. 2014, 1603, 14–22. [Google Scholar] [CrossRef]
  61. Othman, N.S.; Jaafar, M.S.; Rahman, A.A.; Othman, E.S.; Rozlan, A.A. Ultrasound Speed of Polymer Gel Mimicked Human Soft Tissue within Three Weeks. Int. J. Biosci. Biochem. Bioinform. 2011, 1, 223–225. [Google Scholar] [CrossRef]
  62. Blyakhman, F.A.; Makarova, E.B.; Fadeyev, F.A.; Lugovets, D.V.; Safronov, A.P.; Shabadrov, P.A.; Shklyar, T.F.; Melnikov, G.Y.; Orue, I.; Kurlyandskaya, G.V. The Contribution of Magnetic Nanoparticles to Ferrogel Biophysical Properties. Nanomaterials 2019, 9, 232. [Google Scholar] [CrossRef]
  63. Li, Y.; Huang, G.; Zhang, X.; Li, B.; Chen, Y.; Lu, T.; Lu, T.J.; Xu, F. Magnetic Hydrogels and Their Potential Biomedical Applications. Adv. Funct. Mater. 2013, 23, 660–672. [Google Scholar] [CrossRef]
Figure 1. (a) 3D representation of the mineralized collagen fibril (MCF) and a typical trajectory of water molecule. The MCF is characterized by 32% of mineral volume fraction, and it is composed of tropocollagen (light blue) and plate-like apatite crystals (grey). The depicted configuration is obtained after roughly overall 6·106 moves and rotations. (b) Enlarged view of a segment of the 3D random walk trajectory to highlight the path confinement within the channels created between apatite platelets and tropocollagen.
Figure 1. (a) 3D representation of the mineralized collagen fibril (MCF) and a typical trajectory of water molecule. The MCF is characterized by 32% of mineral volume fraction, and it is composed of tropocollagen (light blue) and plate-like apatite crystals (grey). The depicted configuration is obtained after roughly overall 6·106 moves and rotations. (b) Enlarged view of a segment of the 3D random walk trajectory to highlight the path confinement within the channels created between apatite platelets and tropocollagen.
Bioengineering 10 00558 g001
Figure 2. Diffusion coefficient achieved from MSD data analysis versus tortuosity for two reference quotes (a) that are characterized by minimum (c) and maximum (b) values of tortuosity. The light grey bands in (b,c) represent the confidence interval at 95%.
Figure 2. Diffusion coefficient achieved from MSD data analysis versus tortuosity for two reference quotes (a) that are characterized by minimum (c) and maximum (b) values of tortuosity. The light grey bands in (b,c) represent the confidence interval at 95%.
Bioengineering 10 00558 g002
Figure 3. (a) Equatorial planes depicted correspondingto each reference plane considered in the longitudinal direction of the MCF. Average tortuosity (b) and average diffusion coefficient (c) are represented as a function of the reference quote in the longitudinal (L) direction. The light grey bands in (b,c) represent the confidence interval at 95%.
Figure 3. (a) Equatorial planes depicted correspondingto each reference plane considered in the longitudinal direction of the MCF. Average tortuosity (b) and average diffusion coefficient (c) are represented as a function of the reference quote in the longitudinal (L) direction. The light grey bands in (b,c) represent the confidence interval at 95%.
Bioengineering 10 00558 g003
Table 1. Average dimensions of apatite platelets.
Table 1. Average dimensions of apatite platelets.
ParameterMean Value
Width41.80 nm
Thickness3.55 nm
Length94.51 nm
aW 13.19 nm
aT 2.35 nm
aL 39.49 nm
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bini, F.; Pica, A.; Marinozzi, A.; Marinozzi, F. 3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model. Bioengineering 2023, 10, 558. https://doi.org/10.3390/bioengineering10050558

AMA Style

Bini F, Pica A, Marinozzi A, Marinozzi F. 3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model. Bioengineering. 2023; 10(5):558. https://doi.org/10.3390/bioengineering10050558

Chicago/Turabian Style

Bini, Fabiano, Andrada Pica, Andrea Marinozzi, and Franco Marinozzi. 2023. "3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model" Bioengineering 10, no. 5: 558. https://doi.org/10.3390/bioengineering10050558

APA Style

Bini, F., Pica, A., Marinozzi, A., & Marinozzi, F. (2023). 3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model. Bioengineering, 10(5), 558. https://doi.org/10.3390/bioengineering10050558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop