Abstract
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.