Hierarchical Structure and Properties of the Bone at Nano Level

Bone is a highly hierarchical complex structure that consists of organic and mineral components represented by collagen molecules (CM) and hydroxyapatite crystals (HAC), respectively. The nanostructure of bone can significantly affect its mechanical properties. There is a lack of understanding how collagen fibrils (CF) in different orientations may affect the mechanical properties of the bone. The objective of this study is to investigate the effect of interaction, orientation, and hydration on atomic models of the bone composed of collagen helix (CH) and HAC, using molecular dynamics simulations and therefrom bone-related disease origins. The results demonstrate that the mechanical properties of the bone are affected significantly by the orientation of the CF attributed to contact areas at 0° and 90° models. The molecular dynamics simulation illustrated that there is significant difference (p < 0.005) in the ultimate tensile strength and toughness with respect to the orientation of the hydrated and un-hydrated CF. Additionally, the results indicated that having the force in a longitudinal direction (0°) provides more strength compared with the CF in the perpendicular direction (90°). Furthermore, the results show that substituting glycine (GLY) with any other amino acid affects the mechanical properties and strength of the CH, collagen–hydroxyapatite interface, and eventually affects the HAC. Generally, hydration dramatically influences bone tissue elastic properties, and any change in the orientation or any abnormality in the atomic structure of either the CM or the HAC would be the main reason of the fragility in the bone, affecting bone pathology.


Introduction
Bone tissue contains mineral phase (70%) represented by hydroxyapatite (HAC), organic phase (30%) represented by collagen (CO), non-collagenous proteins, bone cells, and water with different percentages [1][2][3]. These contents form the complex structure and provide the material properties of the bone. Investigating the bone at the nano level (<1 µm) considers the architecture of the collagen fibrils alignment. Investigating and understanding damage development in a bone at nano level is as important as the macro and micro levels.
The bone mainly is composed of type I collagen, which as a triple helix protein has the ability to be organized into fibers that provide strength and flexibility to the bone. At the molecular level, the collagen fibril is a triple helix that includes one α2-chain and two α1-chains that are composed of glycine (GLY), proline (PLY) and hydroxyproline (HYP) [4][5][6][7]. Type I collagen can be found in other tissues such as ligaments, tendons, and skin [3]. Osteocalcin, which is another non-collagenous organic protein, has a significant role in new bone mineralization. According to Burr and Allen (Eds.) (2019, p. 127), hydroxyapatite, which is the main component of the mineral phase, is composed of calcium, phosphate, and hydroxide, where the average hydroxyapatite crystal size is 50 × 25 × 2~3 nm [4]. The mineral and collagen orientation significantly affects the mechanical properties of the bone, where bone stiffness primarily depends on the mineral phase, as the mineral density is highly correlated to the strength of the bone [4].
There are different bone diseases that affect the molecular structure of the tissue. Osteogenesis imperfecta (OI) is a genetic disorder that affects bone strength as a result of mutation in the genes (COL1A1/COL1A2) responsible for collagen type I production [8][9][10][11][12][13]. It occurs as a result of the substitution of GLY with valine, arginine, aspartic acid, glutamic acid, and cystine. This substitution affects the mechanical properties and strength of the collagen helix, collagen-hydroxyapatite interface, and eventually affects the mineral part. In general, the abnormality in the atomic structure of either the collagen or the hydroxyapatite is the main reason of the fragility in OI bone tissues.
There are several techniques that are used to characterize the orientation of the collagen fibrils. Some of these techniques are polarized Raman spectroscopy [14,15], polarized Fourier transform infrared spectroscopy [16], polarized second harmonic generation [17,18], small-angle light scattering spectroscopy [19][20][21][22], elastic scattering spectroscopy [23,24], electron transmission diffraction [25,26], electron backscatter diffraction [27], magnetic resonance imaging [28], and microwave method [29]. The main limitations with these techniques are the limited resolution, they are time consuming, costly, and lack the ability to visualize the three-dimensional orientation. Additionally, none of these techniques had the ability to provide a quantitative representation of the bone. In the current study, X-ray based technique is used to observe the bone at the nano level. The orientation of the collagen can be observed by using the X-ray based technique. A molecular dynamics simulation can be used to investigate the mechanical behavior of the bone at the nano level and the collagen helix (CH) and HAC oriented in different modes depending on the orientation observed using the X-ray-based technique.
Experimental and theoretical studies have been performed to understand the atomic structure of the bone [30][31][32][33][34][35][36][37][38][39][40][41]. Experimentally, it is still a challenge to characterize bone atomic and structural integrity at the nano level. Additionally, most studies were focused on investigating CF and fibers under tension [30,31], and investigating the effect of bone structure at nano level on bone failure [32,33]. On the other hand, the theoretical work has been performed mostly to understand the effect of the CM and HAC structure and properties on biomaterials design, such as bone scaffold hierarchical design [34][35][36][37]. A molecular dynamics simulation was exploited to investigate the variation in the material properties of the bone at the molecular level assemblies [38][39][40], and how changing the residue sequence of the collagen helix affects the molecular mechanical properties of the bone [41], and how hydration significantly affects the bone at nano level [42]. To the best of our knowledge, there is a lack in understanding the mechanical properties of the bone at the nano level with respect to CO and HAC orientation. HA surface, hydration, and the chemical environment are important features that need to be investigated. The objective of this study is to elucidate the unique contributions that nano-modeling of the bone provides in the framework of collage orientation and the effect on the bone mechanical properties. The CO-HAC interface is of fundamental importance in investigating bone under healthy and pathological status. In this study, an attempt was made to create a new model for the collagen fibrils of the bone from Swiss Light Source (SLS) data, provided by the Paul Scherrer Institute (PSI) [6], that maps the collagen order and to understand the change in the architecture of the CF and HAC configuration with damage propagation. The ADF 2019.3 modeling suite was used to perform the molecular dynamics simulation.

Materials and Method
SLS data was used to map the collagen order; then a computational three-dimensional modeling was performed using MIMICS 22.0 (Materialise-Software Mimics Innovation Suite, Leuven, Belgium) program. A molecular dynamics 3D simulation was performed to investigate the effect of the interaction, orientation, and hydration of the atomic models of the bone composed of CH and HAC. Finally, numerical analysis was performed to investigate the impact of various factors on the mechanical properties of the atomic models of the bone composed of collagen and hydroxyapatite. The statistical analysis was performed using SAS software (SAS Institute Inc. 2015. SAS ® 12 JSL Syntax. SAS Institute Inc., Cary, NC, USA) and JMP (SAS Institute Inc., Cary, NC, USA). One-way ANOVA test was applied in the analysis. The three factors that were considered in the analysis were the collagen fibril orientation, mineral surface, and hydration. Modulus of elasticity (E) was calculated using MATLAB R2021b (The MathWorks Inc., Natick, MA, USA).

Computational 3D Modeling at Nano Level
SLS technology uses three-dimensional scanning small angle X-ray scattering, as shown in Figure 1. The images that were used were for a small bone sample (1 × 1 × 2.5 mm 3 ) extracted from a T12 human vertebra of a 73 year old man, shown in Figure 2. The ImageJ program was used to calculate the frequency of different orientation angles, and those were (0 • , 20 • , 30 • , 45 • , 60 • and 90 • ). Figure 3 illustrates the defining of the material properties of the bone at nano level depending on the orientation of CF. Highly-oriented vs. weaklyoriented fibrils provide anisotropic or isotropic material properties, respectively. On the other hand, previous studies proposed isotropic material properties for the bone at the micro and nano levels [43,44]. The atomic models of CH, HAC, and bone comprised of CH and HAC were modeled. models of the bone composed of collagen and hydroxyapatite. The statist performed using SAS software (SAS Institute Inc. 2015. SAS ® 12 JSL Synt Inc., Cary, NC, USA) and JMP (SAS Institute Inc., Cary, NC, USA). One-w was applied in the analysis. The three factors that were considered in t the collagen fibril orientation, mineral surface, and hydration. Modulu was calculated using MATLAB R2021b (The MathWorks Inc., Natick, M

Computational 3D Modeling at Nano Level
SLS technology uses three-dimensional scanning small angle X-r shown in Figure 1. The images that were used were for a small bone sa mm 3 ) extracted from a T12 human vertebra of a 73 year old man, shown ImageJ program was used to calculate the frequency of different orienta those were (0°, 20°, 30°, 45°, 60° and 90°). Figure 3 illustrates the definin properties of the bone at nano level depending on the orientation of CF. vs. weakly-oriented fibrils provide anisotropic or isotropic material pr tively. On the other hand, previous studies proposed isotropic material p bone at the micro and nano levels [43,44]. The atomic models of CH, HAC prised of CH and HAC were modeled.  The three dimensional scanning small-angle X-ray scattering information [45].

Molecular Dynamics Simulation (MDS)
In this study, a molecule of CO-HA composite simulation is presented, and a seri of simulations were carried out to investigate the mechanical properties of the bone nano level. Our attempt was to design a realistic representation of bone and expose it uniaxial tensile loading. We investigated the properties of CO-HA interface and analyze the deformation associated with both hydrated and un-hydrated conditions. Our mod aimed to be as accurate as possible a representation of bone tissue at the atomic leve taking into consideration the influence of the nano level on upper micro and macro lev mechanical behavior, in order to understand the complex structure of the bone.

Molecular Dynamics Simulation (MDS)
In this study, a molecule of CO-HA composite simulation is presented, and a series of simulations were carried out to investigate the mechanical properties of the bone at nano level. Our attempt was to design a realistic representation of bone and expose it to uniaxial tensile loading. We investigated the properties of CO-HA interface and analyzed the deformation associated with both hydrated and un-hydrated conditions. Our model aimed to be as accurate as possible a representation of bone tissue at the atomic level, taking into consideration the influence of the nano level on upper micro and macro level mechanical behavior, in order to understand the complex structure of the bone.

Molecular Dynamics Simulation (MDS)
In this study, a molecule of CO-HA composite simulation is presented, and a series of simulations were carried out to investigate the mechanical properties of the bone at nano level. Our attempt was to design a realistic representation of bone and expose it to uniaxial tensile loading. We investigated the properties of CO-HA interface and analyzed the deformation associated with both hydrated and un-hydrated conditions. Our model aimed to be as accurate as possible a representation of bone tissue at the atomic level, taking into consideration the influence of the nano level on upper micro and macro level mechanical behavior, in order to understand the complex structure of the bone.   HA including calcium, phosphate and hydroxide modeled in ADF program (right).

Forcefield
In order to perform the simulation, we developed a forcefield file that includes a previously used HA forcefield [46,47] combined with a CHARMM22 forcefield [48]. The CHARMM forcefield was used previously and has shown it is an adequate representation for the behavior of the bone at the atomic level [49][50][51][52][53]. All the molecular dynamics simulations were performed in an ADF package that has the capabilities to simulate enormous biomolecules. Van der Waals parameters were used to describe the CO-HA interaction. Equation (1) [47] was used to define the potential energy: and

Forcefield
In order to perform the simulation, we developed a forcefield file that includes a previously used HA forcefield [46,47] combined with a CHARMM22 forcefield [48]. The CHARMM forcefield was used previously and has shown it is an adequate representation for the behavior of the bone at the atomic level [49][50][51][52][53]. All the molecular dynamics simulations were performed in an ADF package that has the capabilities to simulate enormous biomolecules. Van der Waals parameters were used to describe the CO-HA interaction. Equation (1) [47] was used to define the potential energy: (2)

Forcefield
In order to perform the simulation, we developed a forcefield file that includes a previously used HA forcefield [46,47] combined with a CHARMM22 forcefield [48]. The CHARMM forcefield was used previously and has shown it is an adequate representation for the behavior of the bone at the atomic level [49][50][51][52][53]. All the molecular dynamics simulations were performed in an ADF package that has the capabilities to simulate enormous biomolecules. Van der Waals parameters were used to describe the CO-HA interaction. Equation (1) [47] was used to define the potential energy: where and The definitions of all the parameters are shown in Table 1. Each CH is composed of three chains and 1014 residues, while each HA unit cell is composed of forty-four atoms with the lattice parameters shown in Table 2. Different orientations are proposed in the current study. The chain sequence of α 1 and α 2 type 1 collagen is shown in Table 3. These orientations are defined with respect to SLS results of collagen orientation. Each HA crystal has the same dimension and geometry, where its height is 1.6 nm. Additionally, each collagen helix length is 14.9 nm, and 1.5 nm in diameter.  Table 3. The chain sequence of α 1 and α 2 type I collagen.
The Amino Acids Sequence To develop the model, each HA unit cell was duplicated 20, 2, and 4 times in the X, Y, and Z directions, respectively. The literature proposed that interaction with different planes of HA plays a significant role in the morphology of the material [54]. To further understand the CO-HA surface interaction, we investigated the effect of CO vs. OH planes on the total mechanical properties of the bone at nano level and how much difference the interaction plane can cause. The interface of the nano constituents was modeled. Python code was used to extract the stress-strain data for each simulation. Once the stress-strain curve was plotted, the mechanical properties were calculated. The collagen helix and HA were oriented in different modes and the simulation was performed to investigate the effect of orientation on the bone mechanical properties. Bone damage was investigated by introducing bone disease (OI) to the model, where the GLY amino acid in the collagen helix was substituted with valine, arginine, aspartic acid, glutamic acid, and cystine. We also introduced a crack with different geometries. In both methods, the mechanical properties were calculated.

Results and Discussion
The current research focuses on understanding how the change in collagen fibrils orientation can improve or deteriorate the structure of the bone. The MDS revealed that having different orientations of CH and HA significantly (<0.001) affect the mechanical properties of the bone. Figure 6 presents a schematic of the bone structure at each hierarchical level, illustrating the change in the modulus of elasticity.
Bioengineering 2022, 9, x FOR PEER REVIEW 7 of 18 interaction plane can cause. The interface of the nano constituents was modeled. Python code was used to extract the stress-strain data for each simulation. Once the stress-strain curve was plotted, the mechanical properties were calculated. The collagen helix and HA were oriented in different modes and the simulation was performed to investigate the effect of orientation on the bone mechanical properties. Bone damage was investigated by introducing bone disease (OI) to the model, where the GLY amino acid in the collagen helix was substituted with valine, arginine, aspartic acid, glutamic acid, and cystine. We also introduced a crack with different geometries. In both methods, the mechanical properties were calculated.

Results and Discussion
The current research focuses on understanding how the change in collagen fibrils orientation can improve or deteriorate the structure of the bone. The MDS revealed that having different orientations of CH and HA significantly (<0.001) affect the mechanical properties of the bone. Figure 6 presents a schematic of the bone structure at each hierarchical level, illustrating the change in the modulus of elasticity. The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength (UTS) and E with respect to the orientation of the hydrated and un-hydrated CF, as shown in Figure 7. The results show that having the force in a longitudinal direction (0°) provides more strength and toughness compared with the CF in the perpendicular direction (90°). The effect of hydration was observed, and the results showed that the stiffness of the CF increased in un-hydrated models. This means that the loss of water increases the rigidity of the CF and eventually makes the bone more susceptible to fracture. Figure 8 shows the molecular structure of the CF at the initial and failure stages. In general, the simulation illustrated that hydration (tightly bound water) decreases the ultimate tensile strength (UTS), toughness, and E of the CF. The MDS points to a loss of water in the collagen causing an increase in the UTS of bone. The simulation illustrated that there is a significant difference (p < 0.005) in UTS with respect to the orientation of the hydrated and un-hydrated collagen fibrils, as shown in Figure 9. However, there was no significant difference (p > 0.05) in the results with respect to the calcium vs. hydroxide mineral surfaces, as shown in Figure 10. The simulation demonstrated that the modulus Figure 6. The hierarchical structure of the bone starts at the femoral bone macrolevel [55] and ends with the nanostructure and molecular level.
The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength (UTS) and E with respect to the orientation of the hydrated and un-hydrated CF, as shown in Figure 7. The results show that having the force in a longitudinal direction (0 • ) provides more strength and toughness compared with the CF in the perpendicular direction (90 • ). The effect of hydration was observed, and the results showed that the stiffness of the CF increased in un-hydrated models. This means that the loss of water increases the rigidity of the CF and eventually makes the bone more susceptible to fracture. Figure 8 shows the molecular structure of the CF at the initial and failure stages. In general, the simulation illustrated that hydration (tightly bound water) decreases the ultimate tensile strength (UTS), toughness, and E of the CF. The MDS points to a loss of water in the collagen causing an increase in the UTS of bone. The simulation illustrated that there is a significant difference (p < 0.005) in UTS with respect to the orientation of the hydrated and un-hydrated collagen fibrils, as shown in Figure 9. However, there was no significant difference (p > 0.05) in the results with respect to the calcium vs. hydroxide mineral surfaces, as shown in Figure 10. The simulation demonstrated that the modulus of elasticity (E) calculated from the slope of a linear fit of the stress-strain curves ( Table 4). The results showed that E was the was higher in un-hydrated simulations (2.29 ± 0.78 GPa) when compared with the hydrated simulations (1.99 ± 0.74 GPa). The results indicated that E was the highest in the longitudinal direction (3.53 GPa) when compared with the CF in the perpendicular direction (1.54 GPa). However, having different mineral surfaces showed no significant difference in E (p > 0.05), as shown in Figures 11 and 12. The constitutive equations of modulus of elasticity with respect to collagen fibril orientation is as follows: with OH mineral surface: For              The MDS with proposed OI bone disease revealed that having different amino acids in the collagen helix significantly (<0.001) affects the mechanical properties of the bone. The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength and toughness with respect to different amino acids. The results show that substituting GLY with any other amino acid affects the mechanical properties and strength of the CH, collagen-hydroxyapatite interface, and eventually affects the HAC, as shown in Figures 13 and 14. In general, the abnormality in the atomic structure of either the collagen or the hydroxyapatite is the main reason of the fragility in OI bone tissues. The MDS with proposed OI bone disease revealed that having different amino acids in the collagen helix significantly (<0.001) affects the mechanical properties of the bone. The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength and toughness with respect to different amino acids. The results show that substituting GLY with any other amino acid affects the mechanical properties and strength of the CH, collagen-hydroxyapatite interface, and eventually affects the HAC, as shown in Figures 13 and 14. In general, the abnormality in the atomic structure of either the collagen or the hydroxyapatite is the main reason of the fragility in OI bone tissues.

Validation with Experimental Work
A comparison between the MDS results and previous experimental results was evaluated. Gupta et al. [56] investigated 36 bone samples (29 wet samples and 7 dry samples). The dimens(13.9 ± 3.4 GPa)ion of each sample was 50 µm × 150 µm × 3 mm, and a tensile test was performed on all the models. The comparison shows agreement between the current MDS results and the experimental testing, as the dry samples showed higher modulus than the wet samples (11.5 ± 3.7 GPa). Note that the MDS show relatively good correlation with the results obtained by Liu, Y. et al. [57], where they investigated 12 samples with 305 ± 79 nm diameter and 12.7 9 ± 8.1 µm length, as shown in Figures 15 and 16. The simulation was compared with Yamamoto, N. [58] experimental results. In this study, ten samples with 410 ± 60 nm were investigated. As can be seen from Figure 17, there was an agreement between the simulation and the previous experiment results in terms of ultimate tensile strength.

Validation with Experimental Work
A comparison between the MDS results and previous experimental results was evaluated. Gupta et al. [56] investigated 36 bone samples (29 wet samples and 7 dry samples). The dimension of each sample was 50 µm × 150 µm × 3 mm, and a tensile test was performed on all the models. The comparison shows agreement between the current MDS results and the experimental testing, as the dry samples showed higher modulus (13.9 ± 3.4 GPa) than the wet samples (11.5 ± 3.7 GPa). Note that the MDS show relatively good correlation with the results obtained by Liu, Y. et al. [57], where they investigated 12 samples with 305 ± 79 nm diameter and 12.7 9 ± 8.1 µm length, as shown in Figures 15 and  16. The simulation was compared with Yamamoto, N. [58] experimental results. In this study, ten samples with 410 ± 60 nm were investigated. As can be seen from Figure 17, there was an agreement between the simulation and the previous experiment results in terms of ultimate tensile strength.     The numerical analysis was performed to understand the impact of various factors on the mechanical properties of the atomic models of the bone composed of collagen and hydroxyapatite. The three factors that were considered in the analysis were the collagen fibril orientation, mineral surface, and hydration. The F-test from the ANOVA indicates that the p-value < 0.05 and it is small and close to zero. This means that there is evidence to reject the claim that the ultimate tensile strength is similar for all the factors and at least one of them is different. The analysis showed that the collagen fibril orientations have a significant effect (p < 0.0001) on the ultimate tensile strength of the atomic models of the bone, and the LSMeans (least square means) differences Student's t-test connecting letter report showed that the ultimate tensile strength is significantly different in terms of orientation. Hydration showed a significant effect as well (p < 0.0001), where un-hydrated models ultimate strength was significantly higher than hydrated models. The LSMeans differences Student's t-test connecting letter report showed no significant difference between Ca and OH mineral surfaces. In general, the statistical analysis showed the collagen fibril orientation and the hydration have significant effects on the ultimate tensile strength, but the mineral surface does not have any significant effect. The numerical analysis was performed to understand the impact of various factors on the mechanical properties of the atomic models of the bone composed of collagen and hydroxyapatite. The three factors that were considered in the analysis were the collagen fibril orientation, mineral surface, and hydration. The F-test from the ANOVA indicates that the p-value < 0.05 and it is small and close to zero. This means that there is evidence to reject the claim that the ultimate tensile strength is similar for all the factors and at least one of them is different. The analysis showed that the collagen fibril orientations have a significant effect (p < 0.0001) on the ultimate tensile strength of the atomic models of the bone, and the LSMeans (least square means) differences Student's t-test connecting letter report showed that the ultimate tensile strength is significantly different in terms of orientation. Hydration showed a significant effect as well (p < 0.0001), where un-hydrated models ultimate strength was significantly higher than hydrated models. The LSMeans differences Student's t-test connecting letter report showed no significant difference between Ca and OH mineral surfaces. In general, the statistical analysis showed the collagen fibril orientation and the hydration have significant effects on the ultimate tensile strength, but the mineral surface does not have any significant effect.  Table 5 shows the prediction equations of stress vs. strain for each collagen fibril orientation. Figure 18 shows the bivariate fit of ultimate tensile strength by orientation, where the ultimate tensile strength is the highest for 0 • and decreases as the orientation angle increases until reaching the minimum at 90 • .

Prediction Equation
where is the stress in un-hydrated environment, is the strain, and θ is the orientation angle, with SSE: 1.473 × 10 4 , R-square: 0.9324, Adjusted R-square: 0.9267, and RMSE: 8.074. This regression model can be used to predict the ultimate tensile strength of the CF with respect to different orientations. The sensitivity analysis demonstrated the A regression model was developed for CF model ultimate tensile strength with respect to strain and the orientation angle using MATLAB.
mineral surface vs. the strain, and the orientation of CF obtained from the constitutive equations for each CF orientation shown in Table 5. The analysis showed the stresses were higher when CF was in the parallel direction (0°) with maximum stress (125.12 MPa). On the other hand, the stresses were the least when the CF was in the perpendicular direction (90°) with maximum stress (60.03 MPa). In general, the sensitivity analysis ( Figure 19) illustrates that the orientation has a significant effect on stress vs. strain distribution where the stress decreases as the angle increases to reach the least at 90°. Figure 19. The sensitivity plot demonstrating the relationship between the maximum stress of the un-hydrated CO-HAC composite with OH mineral surface vs. the strain and the orientation of CF. The analysis illustrates that the orientation has a significant effect on stress vs. strain distribution where the stress decreases as the angle increases to reach the least at 90°.

Conclusions
The current investigation suggests that the mechanical properties of the bone are affected significantly by the orientation of the CF. Any change in the structure of the collagen-HAC composite would lead to variable bone diseases (OI). These results provide significant insight into the behavior of the collagen-HAC interface and represent a leap in understanding of bone material performance at the nano level. Any change in the orientation would lead to variable bone diseases. The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength, toughness, and modulus of elasticity with respect to the orientation of the hydrated and un-hydrated collagen fibrils. Additionally, the results show that having the force in a longitudinal direction (0°) provides more mechanical properties compared with the collagen fibril in the perpendicular direction (90°). The main reason for that difference is the fact that longitudinal orientation with 0° models had more contact areas than the perpendicular orientation with 90° models. Furthermore, validation showed agreement between the current MDS and the experimental testing. The effect of hydration was observed, and the results showed that the stiffness of the CF increased in un-hydrated models. This means that the loss of water increases the rigidity of the CF and eventually makes the bone more susceptible to fracture. Generally, the water content dramatically influences the elastic properties of the tissues. Additionally, substituting GLY with any other amino acid affects the mechanical properties and strength of the CH, collagen-HAC interface, and eventually has effects on the HAC. Consequently, it has been evidenced that the abnormality in the atomic structure of either the collagen or the hydroxyapatite is the main reason of the fragility in OI bone tissues. Further investigation is needed to understand the effect of the atomic structure and nano level Figure 19. The sensitivity plot demonstrating the relationship between the maximum stress of the un-hydrated CO-HAC composite with OH mineral surface vs. the strain and the orientation of CF. The analysis illustrates that the orientation has a significant effect on stress vs. strain distribution where the stress decreases as the angle increases to reach the least at 90 • .

Conclusions
The current investigation suggests that the mechanical properties of the bone are affected significantly by the orientation of the CF. Any change in the structure of the collagen-HAC composite would lead to variable bone diseases (OI). These results provide significant insight into the behavior of the collagen-HAC interface and represent a leap in understanding of bone material performance at the nano level. Any change in the orientation would lead to variable bone diseases. The MDS illustrated that there is a significant difference (p < 0.005) in the ultimate tensile strength, toughness, and modulus of elasticity with respect to the orientation of the hydrated and un-hydrated collagen fibrils. Additionally, the results show that having the force in a longitudinal direction (0 • ) provides more mechanical properties compared with the collagen fibril in the perpendicular direction (90 • ). The main reason for that difference is the fact that longitudinal orientation with 0 • models had more contact areas than the perpendicular orientation with 90 • models. Furthermore, validation showed agreement between the current MDS and the experimental testing. The effect of hydration was observed, and the results showed that the stiffness of the CF increased in un-hydrated models. This means that the loss of water increases the rigidity of the CF and eventually makes the bone more susceptible to fracture. Generally, the water content dramatically influences the elastic properties of the tissues. Additionally, substituting GLY with any other amino acid affects the mechanical properties and strength of the CH, collagen-HAC interface, and eventually has effects on the HAC. Consequently, it has been evidenced that the abnormality in the atomic structure of either the collagen or the hydroxyapatite is the main reason of the fragility in OI bone tissues. Further investigation is needed to understand the effect of the atomic structure and nano level mechanical properties on higher micro and macro levels of bone tissue. This study represents a step toward a deeper understanding of the structure-mechanical function relationship of the bone.

Data Availability Statement:
The data presented in this study are available upon reasonable request.