Thickness of hydroxyapatite nanocrystal controls mechanical properties of the collagen-hydroxyapatite interface

Hierarchical Collagen/Apatite Co-assembly for Injection of Mineralized Fibrillar Tissue Analogues

Mineralized biological tissues rich in type I collagen (e.g., bone and dentin) exhibit complex anisotropic suprafibrillar organizations in which the organic and inorganic moieties are intimately coassembled over several length scales. Above a critical size, a defect in such tissue cannot be self-repaired. Biomimetic materials with a composition and microstructure similar to that of bone have been shown to favorably influence bone regeneration. This highlights the value of developing a similar formulation in an injectable form to enable minimally invasive techniques. Here, we report on the fabrication and application potential of an injectable collagen/CHA (carbonated hydroxyapatite) cell-free hydrogel. The organic part consists of spray-dried nondenatured and dense collagen microparticles, while the inorganic part consists of biomimetic apatite mineral. By mixing both powders at desired tissue-like ratios with an aqueous solvent in one step, spontaneous co-self-assembly occurs, leading to the formation of a mineralized matrix with suprafibrillar tissue-like features thanks to the induced liquid crystalline properties of collagen on one hand and apatite on the other hand. When injected into soft tissue, the mineralized collagen hydrogel free of chemical cross-linking agents exhibits suitable cohesion and is biocompatible. Preliminary in vitro tests in a tooth cavity model show its integration onto dentin with a biomimetic interface. Based on the results, this versatile injectable mineralized collagen hydrogel shows promising potential as a biomaterial for bone tissue repair and mineralized tissue-like ink for bioprinting applications.

On the glow of cremated remains: long-lived green photo-luminescence of heat-treated human bones

The long-lived green luminescence of human bone (that has been heated to 600 °C for a short duration) is attributed to a carbon quantum dot material (derived from collagen) encapsulated and protected by an inorganic matrix (derived from bone apatite) and is more intense in dense rigid and crystalline parts of (healthy) human bones. The strong collagen-apatite interaction results (upon decomposition) in a protective inorganic environment of the luminescent centers allowing long-lived triplet-based emission of a carbon (quantum) dot-like material at room temperature, as well as resilience against oxidation between 550 and 650 °C. The graphitic black phase (obtained upon heating around 400 °C) is a precursor to the luminescent carbon-based material, that is strongly interacting with the crystalline inorganic matrix. Human bone samples that have been heated to 600 °C were subjected to steady-state and time-resolved spectroscopy. Excitation-emission matrix (EEM) luminescence spectroscopy revealed a broad range of excitation and emission wavelengths, indicating a heterogeneous system with a broad density of emissive states. The effect of low temperature on the heat-treated bone was studied with Cryogenic Steady State Luminescence Spectroscopy. Cooling the bone to 80 K leads to a slight increase in total emission intensity as well as an intensity increase towards to red part of the spectrum, incompatible with a defect state model displaying luminescent charge recombination in the inorganic matrix. Time-resolved spectroscopy with an Optical Multichannel Analyzer (OMA) and Time Correlated Single Photon Counting (TCSPC) of these samples showed that the decay could be fitted with a multi-exponential decay model as well as with second-order decay kinetics. Confocal Microscopy revealed distinct (plywood type) structures in the bone and high intensity-fast decay areas as well as a spatially heterogeneous distribution of green and (fewer) red emissive species. The use of the ATTO 565 dye aided in bone-structure visualization by chemical adsorption. Conceptually our data interpretation corresponds to previous reports from the material science field on luminescent powders.

Curved mineral platelets in bone

Bone is a composite material principally made up of a mineral phase (apatite) and collagen fibrils. The mineral component of bone occurs in the form of polycrystalline platelets 2-6 nm in thickness. These platelets are packed and probably glued together in stacks of two or more, ranging up to >30 platelets. Here we show that most of these stacks are curved flat sheets whose cylindrical axes are oriented parallel to the long axes of collagen fibrils. Consequently, the curvature of the platelets is not detectable in TEM sections cut parallel to the collagen fibril axes. The radius of curvature around these axes ranges from about 25 nm (the average radius of the collagen fibrils) to 100's of nm. The shapes of these curved forms contribute to the compressive strength of bone. STATEMENT OF SIGNIFICANCE: Bone, the material of which bones are made, is mainly composed of a protein, collagen, and the mineral apatite (calcium phosphate). The crystals have long been known to be flat plates about 5 nanometers (nm) thick. Here we show that the crystals are bound together in curved platelets with a radius of curvature between 25 and several hundred nm, which weave between fibrils of collagen. Some platelets wrap tightly around fibrils. The platelets form stacks of from two to up to 30. The crystals in the platelets are all oriented parallel to the cylindrical fibrils even though most crystals are not in contact with collagen. These curved structures provide greater strength to bone.

Multiscale and multidisciplinary analysis of aging processes in bone

The world population is increasingly aging, deeply affecting our society by challenging our healthcare systems and presenting an economic burden, thus turning the spotlight on aging-related diseases: exempli gratia, osteoporosis, a silent disease until you suddenly break a bone. The increase in bone fracture risk with age is generally associated with a loss of bone mass and an alteration in the skeletal architecture. However, such changes cannot fully explain increased fragility with age. To successfully tackle age-related bone diseases, it is paramount to comprehensively understand the fundamental mechanisms responsible for tissue degeneration. Aging mechanisms persist at multiple length scales within the complex hierarchical bone structure, raising the need for a multiscale and multidisciplinary approach to resolve them. This paper aims to provide an overarching analysis of aging processes in bone and to review the most prominent outcomes of bone aging. A systematic description of different length scales, highlighting the corresponding techniques adopted at each scale and motivating the need for combining diverse techniques, is provided to get a comprehensive description of the multi-physics phenomena involved.

Exploring the hierarchical structure of lamellar bone and its impact on fracture behaviour: A computational study using a phase field damage model

Bone is a naturally occurring composite material composed of a stiff mineral phase and a compliant organic matrix of collagen and non-collagenous proteins (NCP). While diverse mineral morphologies such as platelets and grains have been documented, the precise role of individual constituents, and their morphology, remains poorly understood. To understand the role of constituent morphology on the fracture behaviour of lamellar bone, a damage based representative volume element (RVE) was developed, which considered various mineral morphologies and mineralised collagen fibril (MCF) configurations. This model framework incorporated a novel phase-field damage model to predict the onset and evolution of damage at mineral-mineral and mineral-MCF interfaces. It was found that platelet-based mineral morphologies had superior mechanical performance over their granular counterparts, owing to their higher load-bearing capacity, resulting from a higher aspect ratio. It was also found that MCFs had a remarkable capacity for energy dissipation under axial loading, with these fibrillar structures acting as barriers to crack propagation, thereby enhancing overall elongation and toughness. Interestingly, the presence of extrafibrillar platelet-based minerals also provided an additional toughening through a similar mechanism, whereby these structures also inhibited crack propagation. These findings demonstrate that the two primary constituent materials of lamellar bone play a key role in its toughening behaviour, with combined effect by both mineral and MCFs to inhibit crack propagation at this scale. These results have provided novel insight into the fracture behaviour of lamellar bone, enhancing our understanding of microstructure-property relationships at the sub-tissue level.

In-silico simulation of nanoindentation on bone using a 2D cohesive finite element model

This study proposed and validated a 2D finite element (FE) model for conducting in-silico simulations of in-situ nanoindentation tests on mineralized collagen fibrils (MCF) and the extrafibrillar matrix (EFM) within human cortical bone. Initially, a multiscale cohesive FE model was developed by adapting a previous model of bone lamellae, encompassing both MCF and EFM. Subsequently, nanoindentation tests were simulated in-silico using this model, and the resulting predictions were compared to AFM nanoindentation test data to verify the model's accuracy. The FE model accurately predicted nanoindentation results under wet conditions, closely aligning with outcomes obtained from AFM nanoindentation tests. Specifically, it successfully mirrored the traction/separation curve, nanoindentation modulus, plastic energy dissipation, and plastic energy ratio obtained from AFM nanoindentation tests. Additionally, this in-silico model demonstrated its ability to capture alterations in nanoindentation properties caused by the removal of bound water, by considering corresponding changes in mechanical properties of the collagen phase and the interfaces among bone constituents. Notably, significant changes in the elastic modulus and plastic energy dissipation were observed in both MCF and EFM compartments of bone, consistent with observations in AFM nanoindentation tests. These findings indicate that the proposed in-silico model effectively captures the influence of ultrastructural changes on bone's mechanical properties at sub-lamellar levels. Presently, no experimental methods exist to conduct parametric studies elucidating the ultrastructural origins of bone tissue fragility. The introduction of this in-silico model presents an invaluable tool to bridge this knowledge gap in the future.

An Investigation into the Relative Efficacy of High-Velocity Air-Fuel-Sprayed Hydroxyapatite Implants Based on the Crystallinity Index, Residual Stress, Wear, and In-Flight Powder Particle Behavior

Due to its resemblance to the bone, hydroxyapatite (HA) has been widely used for bioactive surface modification of orthopedic implants. However, it undergoes significant thermal decomposition and phase transformations at a high operating temperature, leading to premature implant failure. This investigation uses high-velocity air-fuel (HVAF) spray, an emerging low-temperature thermal spray technique, to deposit HA over the Ti-6Al-4V substrate. Coating characteristics, such as the crystallinity index and phase analysis, were measured using X-ray diffraction, Raman analysis, and Fourier transform infrared spectroscopy, residual stress using the sin2ψ method, and tribological performance by a fretting wear test. The coating retained an over 90% crystallinity index, a crystallite size of 41.04 nm, a compressive residual stress of -229 ± 34.5 MPa, and a wear rate of 1.532 × 10-3 mm3 N-1 m-1. Computational in-flight particle traits of HA particles (5 to 60 μm) were analyzed using computational fluid dynamics; it showed that 90% of particles were deposited at a 700 to 1000 m/s velocity and a 900 to 1450 K temperature with a 2.1 ms mean residence time. In-flight particle oxidation was minimized, and particle impact deformation was maximized, which caused severe plastic deformation, forming crystalline, compressive residual stressed coatings. The thermal decomposition model of low-temperature HVAF-sprayed HA particles helped to understand the implants' crystallinity index, residual stress, and tribological characteristics. Hence, this experimental and computational analysis shows that the HVAF process can be a promising candidate for biomedical applications for having strong and durable implants.

Molecular mechanics and failure mechanisms in B. mori Silk Fibroin-hydroxyapatite composite interfaces: Effect of crystal thickness and surface characteristics

Bombyx mori Silk Fibroin-hydroxyapatite (B. mori SF-HA) bio-nanocomposite is a prospective biomaterial for tissue engineered graft for bone repair. Here, B. mori SF is primarily a soft and tough organic phase, and HA is a hard and stiff mineral phase. In biomaterial design, an understanding about the nanoscale mechanics of SF-HA interface, such as interfacial interaction and interface debonding mechanisms between the two phases is essential for obtaining required functionality. To investigate such nanoscale behavior, molecular dynamics method is a preferred approach. Present study focuses on understanding of the interface debonding mechanisms at SF-HA interface in B. mori SF-HA bio-nanocomposite at nanometer length scale. For this purpose, nanoscale atomistic models of SF-HA interface are also developed based on the HA crystal size and HA surface type (Ca2+ dominated and OH- dominated) in contact with SF. Mechanical behavior analysis of these SF-HA interface models under pull-out type test were performed using Molecular Dynamics (MD) simulations. Surface pull-off strength values in the range of 0.4-0.8 GPa were obtained for SF-HA interface models, for different HA crystal thicknesses, wherein, the pull-off strength values are found to increase with increase in HA thicknesses. Analyses show that deformation mechanisms in SF-HA interface deformation, is a combination of shear deformation in SF phase followed by disintegration of SF phase from HA block. Furthermore, higher rupture force values were obtained for SF-HA interface with Ca2+ dominated HA surface in contact with SF phase, indicating that SF protein has a higher affinity for Ca2+ dominated surface of HA phase. Current work contributes in developing an understanding of mechanistic interactions between organic and inorganic phases in B. mori SF-HA composite nanostructure.

Fatigue-driven compliance increase and collagen unravelling in mechanically tested anterior cruciate ligament

Approximately 300,000 anterior cruciate ligament (ACL) tears occur annually in the United States, half of which lead to the onset of knee osteoarthritis within 10 years of injury. Repetitive loading is known to result in fatigue damage of both ligament and tendon in the form of collagen unravelling, which can lead to structural failure. However, the relationship between tissue's structural, compositional, and mechanical changes are poorly understood. Herein we show that repetitive submaximal loading of cadaver knees causes an increase in co-localised induction of collagen unravelling and tissue compliance, especially in regions of greater mineralisation at the ACL femoral enthesis. Upon 100 cycles of 4× bodyweight knee loading, the ACL exhibited greater unravelled collagen in highly mineralized regions across varying levels of stiffness domains as compared to unloaded controls. A decrease in the total area of the most rigid domain, and an increase in the total area of the most compliant domain was also found. The results highlight fatigue-driven changes in both protein structure and mechanics in the more mineralized regions of the ACL enthesis, a known site of clinical ACL failure. The results provide a starting point for designing studies to limit ligament overuse injury.

Effect of Pore Defects on Uniaxial Mechanical Properties of Bulk Hexagonal Hydroxyapatite: A Molecular Dynamics Study

Hydroxyapatite (HAP) is a calcium apatite bioceramic used in various naturally-derived and synthetic forms for bone repair and regeneration. While useful for the regrowth of osseus tissue, the poor load-bearing capacity of this material relative to other biomaterials is worsened by the propensity for pore formation during the synthetic processing of scaffolds, blocks, and granules. Here we use molecular dynamics (MD) simulations to improve the current understanding of the defect-altered uniaxial mechanical response in hexagonal HAP single crystals relative to defect-free structures. The inclusion of a central spherical pore within a repeated lattice was found to reduce both the failure stress and failure strain in uniaxial tension and compression, with up to a 30% reduction in maximum stress at the point of failure compared to a perfect crystalline structure observed when a 30 Å diameter pore was included. The Z axis ([0 0 0 1] crystalline direction) was found to be the least susceptible to pore defects in tension but the most sensitive to pore inclusion in compression. The deformation mechanisms are discussed to explain the observed mechanical responses, for which charge imbalances and geometric stress concentration factor effects caused by pore inclusion play a significant role.

In silico study on probing atomistic insights into structural stability and tensile properties of Fe-doped hydroxyapatite single crystals

Hydroxyapatite (HA, Ca10PO4(OH)2) is a widely explored material in the experimental domain of biomaterials science, because of its resemblance with natural bone minerals. Specifically, in the bioceramic community, HA doped with multivalent cations (e.g., Mg2+, Fe2+, Sr2+, etc.) has been extensively investigated in the last few decades. Experimental research largely established the critical role of dopant content on mechanical and biocompatibility properties. The plethora of experimental measurements of mechanical response on doped HA is based on compression or indentation testing of polycrystalline materials. Such measurements, and more importantly the computational predictions of mechanical properties of single crystalline (doped) HA are scarce. On that premise, the present study aims to build atomistic models of Fe2+-doped HA with varying Fe content (10, 20, 30, and 40 mol%) and to explore their uniaxial tensile response, by means of molecular dynamics (MD) simulation. In the equilibrated unit cell structures, Ca(1) sites were found to be energetically favourable for Fe2+ substitution. The local distribution of Fe2+ ions significantly affects the atomic partial charge distribution and chemical symmetry surrounding the functional groups, and such signatures are found in the MD analyzed IR spectra. The significant decrease in the intensity of the IR bands found in the Fe-doped HA together with band splitting, because of the symmetry changes in the crystal structure. Another important objective of this work is to computationally predict the mechanical response of doped HA in their single crystal format. An interesting observation is that the elastic anisotropy of undoped HA was not compromised with Fe-doping. Tensile strength (TS) is systematically reduced in doped HA with Fe2+ dopant content and a decrease in TS with temperature can be attributed to the increased thermal agitation of atoms at elevated temperatures. The physics of the tensile response was rationalized in terms of the strain dependent changes in covalent/ionic bond framework (Ca-P distance, P-O bond strain, O-P-O angular strain, O-H bond distance). Further, the dynamic changes in covalent bond network were energetically analyzed by calculating the changes in O-H and P-O bond vibrational energy. Summarizing, the current work establishes our foundational understanding of the atomistic phenomena involved in the structural stability and tensile response of Fe-doped HA single crystals.

Treatment of postmenopausal osteoporosis patients with teriparatide for 24 months reverts forming bone quality indices to premenopausal healthy control values

Postmenopausal osteoporosis (PMOP) therapies are frequently evaluated by bone mineral density (BMD) gains against patients receiving placebo (calcium and vitamin D supplementation, a mild bone turnover-suppressing intervention), which is not equivalent to either healthy or treatment-naive PMOP. The aim of the present observational study was to assess the effects of TPTD treatment in PMOP (20 μg, once daily) at 6 (TPTD 6m; n = 28, age 65 ± 7.3 years), and 24 (TPTD 24m; n = 32, age 67.4 ± 6.15 years) months on bone quality indices at actively forming trabecular surfaces (with fluorescent double labels). Data from the TPTD-treated PMOP patients were compared with those in healthy adult premenopausal women (HC; n = 62, age 40.5 ± 10.6 years), and PMOP receiving placebo (PMOP-PLC; n = 94, age 70.6 ± 4.5 years). Iliac crest biopsies were analyzed by Raman microspectroscopy at three distinct tissue ages: mid-distance between the second label and the bone surface, mid-distance between the two labels, and 1 μm behind the first label. Mineral to matrix ratio (MM), mineral maturity/crystallinity (MMC), tissue water (TW), glycosaminoglycan (GAGs), and pyridinoline (Pyd) content were determined. Outcomes were compared by ANCOVA with subject age and tissue age as covariates, and health status as a fixed factor, followed by Sidak's post-hoc testing (significance assigned to p < 0.05). Both TPTD groups increased MM compared to PMOP-PLC. While TPTD 6m had values similar to HC, TPTD 24m had higher values compared to either HC or TPTD 6m. Both TPTD groups had lower MMC values compared to PMOP-PLC and similar to HC. TPTD 6m patients had higher TW content compared to HC, while TPTD 24m had values similar to HC and lower than either PMOP-PLC or TPTD 6m. Both TPTD groups had lower GAG content compared to HC group, while TPTD 6m had higher values compared to PMOP-PLC. Finally, TPTD 6m patients had higher Pyd content compared to HC and lower compared to PMOP-PLC, while TPTD 24m had lower values compared to PMOP-PLC and TPTD 6m, and similar to HC group. The results of the present study indicate that effects of TPTD on forming trabecular bone quality indices depend on treatment duration. At the recommended length of 24 m, TPTD restores bone mineral and organic matrix quality indices (MMC, TW, Pyd content) to premenopausal healthy (HC) levels.

Mineral Crystal Thickness in Calcified Cartilage and Subchondral Bone in Healthy and Osteoarthritic Human Knees

Osteoarthritis (OA) is the most common joint disease, where articular cartilage degradation is often accompanied with sclerosis of the subchondral bone. However, the association between OA and tissue mineralization at the nanostructural level is currently not understood. In particular, it is technically challenging to study calcified cartilage, where relevant but poorly understood pathological processes such as tidemark multiplication and advancement occur. Here, we used state-of-the-art microfocus small-angle X-ray scattering with a 5-μm spatial resolution to determine the size and organization of the mineral crystals at the nanostructural level in human subchondral bone and calcified cartilage. Specimens with a wide spectrum of OA severities were acquired from both medial and lateral compartments of medial compartment knee OA patients (n = 15) and cadaver knees (n = 10). Opposing the common notion, we found that calcified cartilage has thicker and more mutually aligned mineral crystals than adjoining bone. In addition, we, for the first time, identified a well-defined layer of calcified cartilage associated with pathological tidemark multiplication, containing 0.32 nm thicker crystals compared to the rest of calcified cartilage. Finally, we found 0.2 nm thicker mineral crystals in both tissues of the lateral compartment in OA compared with healthy knees, indicating a loading-related disease process because the lateral compartment is typically less loaded in medial compartment knee OA. In summary, we report novel changes in mineral crystal thickness during OA. Our data suggest that unloading in the knee might be involved with the growth of mineral crystals, which is especially evident in the calcified cartilage. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Hydroxyapatite Use in Spine Surgery—Molecular and Clinical Aspect

Hydroxyapatite possesses desirable properties as a scaffold in tissue engineering: it is biocompatible at a site of implantation, and it is degradable to non-toxic products. Moreover, its porosity enables infiltration of cells, nutrients and waste products. The outcome of hydroxyapatite implantation highly depends on the extent of the host immune response. Authors emphasise major roles of the chemical, morphological and physical properties of the surface of biomaterial used. A number of techniques have been applied to transform the theoretical osteoconductive features of HAp into spinal fusion systems—from integration of HAp with autograft to synthetic intervertebral implants. The most popular uses of HAp in spine surgery include implants (ACDF), bone grafts in posterolateral lumbar fusion and transpedicular screws coating. In the past, autologous bone graft has been used as an intervertebral cage in ACDF. Due to the morbidity related to autograft harvesting from the iliac bone, a synthetic cage with osteoconductive material such as hydroxyapatite seems to be a good alternative. Regarding posterolateral lumbar fusion, it requires the graft to induce new bone growth and reinforce fusion between the vertebrae. Hydroxyapatite formulations have shown good results in that field. Moreover, the HAp coating has proven to be an efficient method of increasing screw fixation strength. It can decrease the risk of complications such as screw loosening after pedicle screw fixation in osteoporotic patients. The purpose of this literature review is to describe in vivo reaction to HAp implants and to summarise its current application in spine surgery.

c-Axis-Oriented Platelets of Crystalline Hydroxyapatite in Biomimetic Intrafibrillar Mineralization of Polydopamine-Functionalized Collagen Type I

Mineralized collagen fibrils are important basic building blocks of calcified tissues, such as bone and dentin. Polydopamine (PDA) can introduce functional groups, i.e., hydroxyl and amine groups, on the surfaces of type I collagen (Col-I) as possible nucleation sites of calcium phosphate (CaP) crystallization. Molecular bindings in between PDA and Col-I fibrils (Col-PDA) have been found to significantly reduce the interfacial energy. The wetting effect, mainly hydrophilicity due to the functional groups, escalates the degree of mineralization. The assembly of Col-I molecules into fibrils was initiated at the designated number of collagenous molecules and PDA. In contrast to the infiltration of amorphous calcium phosphate (ACP) precursors into the Col-I matrix by polyaspartic acid (pAsp), this collagen assembly process allows nucleation and ACP to exist in advance by PDA in the intrafibrillar matrix. PDA bound to specific sites, i.e., gap and overlap zones, by the regular arrangement of Col-I fibrils enhanced ACP nucleation and thus mineralization. As a result, the c-axis-oriented platelets of crystalline hydroxyapatite in the Col-I fibril matrix were observed in the enhanced mineralization through PDA functionalization.

Non‐Obese MKR Mouse Model of Type 2 Diabetes Reveals Skeletal Alterations in Mineralization and Material Properties

Obesity is a common comorbidity of type 2 diabetes (T2D). Therefore, increased risk of fragility fractures in T2D is often confounded by the effects of obesity. This study was conducted to elucidate the mechanistic basis by which T2D alone leads to skeletal fragility. We hypothesized that obesity independent T2D would deteriorate bone's material quality by accumulating defects in the mineral matrix and undesired modifications in its organic matrix associated with increased oxidative stress and hyperglycemia. To test this hypothesis, we used 15‐week‐old male non‐obese mice with engineered muscle creatine kinase promoter/human dominant negative insulin growth factor 1 (IGF‐I) receptor (MKR) and FVB/N wild‐type (WT) controls (n = 12/group). MKR mice exhibit reduced insulin production and loss of glycemic control leading to diabetic hyperglycemia, verified by fasting blood glucose measurements (>250 mg/dL), without an increase in body weight. MKR mice showed a significant decrease in femoral radial geometry (cortical area, moment of inertia, cortical thickness, endosteal diameter, and periosteal diameter). Bone mineral density (BMD), as assessed by micro–computed tomography (μCT), remained unchanged; however, the quality of bone mineral was altered. In contrast to controls, MKR mice had significantly increased hydroxyapatite crystal thickness, measured by small‐angle X‐ray scattering, and elongated c‐axis length of the crystals evaluated by confocal Raman spectroscopy. There was an increase in changes in the organic matrix of MKR mice, associated with enhanced glycoxidation (carboxymethyl‐lysine [CML] and pentosidine) and overall glycation (fluorescent advanced glycation end products), both of which were associated with various measures of bone fragility. Moreover, increased CML formation positively correlated with elongated mineral crystal length, supporting the role of this negatively charged side chain to attract calcium ions, promote growth of hydroxyapatite, and build a physical link between mineral and collagen. Collectively, our results show, for the first time, changes in bone matrix in a non‐obese T2D model in which skeletal fragility is attributable to alterations in the mineral quality and undesired organic matrix modifications. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Poor bone matrix quality: What can be done about it?

PURPOSE OF THE REVIEW

Bone's ability to withstand load resisting fracture and adapting to it highly depends on the quality of its matrix and its regulators. This review focuses on the contribution of bone quality to fracture resistance and possible therapeutic targets for skeletal fragility in aging and disease.

RECENT FINDINGS

The highly organized, hierarchical composite structure of bone extracellular matrix together with its (re)modeling mechanisms and microdamage dynamics determines its stiffness, strength, and toughness. Aging and disease affect the biological processes regulating bone quality, thus resulting in defective extracellular matrix and bone fragility. Targeted therapies are being developed to restore bone's mechanical integrity. However, their current limitations include low tissue selectivity and adverse side effects. Biological and mechanical insights into the mechanisms controlling bone quality, together with advances in drug delivery and studies in animal models, will accelerate the development and translation to clinical application of effective targeted-therapeutics for bone fragility.

Molecular origin of viscoelasticity in mineralized collagen fibrils

Bone is mineralized tissue constituting the skeletal system, supporting and protecting the body's organs and tissues. In addition to such fundamental mechanical functions, bone also plays a remarkable role in sound conduction. From a mechanical standpoint, bone is a composite material consisting of minerals and collagen arranged in multiple hierarchical structures, with a complex anisotropic viscoelastic response, capable of transmitting and dissipating energy. At the molecular level, mineralized collagen fibrils are the basic building blocks of bone tissue, and hence, understanding bone properties down to fundamental tissue structures enables better identification of the mechanisms of structural failures and damage. While efforts have focused on the study of micro- and macro-scale viscoelasticity related to bone damage and healing based on creep, mineralized collagen has not been explored at the molecular level. We report a study that aims at systematically exploring the viscoelasticity of collagenous fibrils with different mineralization levels. We investigate the dynamic mechanical response upon cyclic and impulsive loads to observe the viscoelastic phenomena from either shear or extensional strains via molecular dynamics. We perform a sensitivity analysis with several key benchmarks: intrafibrillar mineralization percentage, hydration state, and external load amplitude. Our results show an increase of the dynamic moduli with an increase of the mineral percentage, pronounced at low strains. When intrafibrillar water is present, the material softens the elastic component, but considerably increases its viscosity, especially at high frequencies. This behavior is confirmed from the material response upon impulsive loads, in which water drastically reduces the relaxation times throughout the input velocity range by one order of magnitude, with respect to the dehydrated counterparts. We find that, upon transient loads, water has a major impact on the mechanics of mineralized fibrillar collagen, being able to improve the capability of the tissue to passively and effectively dissipate energy, especially after fast and high-amplitude external loads. Our study provides knowledge of bone mechanics in relation to pathologies deriving from dehydration or traumas. Moreover, these findings show the potential for being used in designing new bioinspired materials not limited to tissue engineering applications, in which passive mechanisms for dissipating energy can prevent structural failures.

Interfacial bonding between mineral platelets in bone and its effect on mechanical properties of bone

Bone is a composite material consisting principally of apatite mineral, collagen fibrils, non-collagenous proteins, and other organic species. Recent electron microscopy studies have shown that the mineral in bone occurs as stacks of thin polycrystalline sheets ("mineral lamellae," MLs) which surround and lie between the collagen fibrils. We focus on the effect of the interface between these mineral lamellae on the mechanical properties of bone. Previous studies on bone treated with sodium hypochlorite (NaClO) to remove all organic material showed a greatly weakened mineral framework. Here, we treated femoral cortical bone with ethylenediamine (EDA), which only removes collagen, to study the effect of its removal on bone properties. We tested the degree of completion of the treatment by Raman spectroscopy and thermogravimetric analysis. When only collagen is removed, a continuous mineral structure remains and is less weakened than by NaClO treatment. Transmission electron microscopy study of finely ground particles of the EDA treated bone shows that stacks of MLs remain joined, whereas in NaClO treated bone, only isolated crystals are present. Thus, we infer that the MLs in bone are held together in stacks by an organic glue, which is destroyed by NaClO, but which survives the EDA treatment. We show that this glue may contribute to the stiffness, strength, and energy absorption of bone. Further studies are needed to discover the chemical nature of this glue. This study provides a starting point for such investigations.

Study of nanoscale deformation mechanisms in bulk hexagonal hydroxyapatite under uniaxial loading using molecular dynamics

Hydroxyapatite (HAP) is a natural bioceramic which is currently used in scaffolds and coatings for the regrowth of osseous tissue but offers poor load-bearing capacity compared to other biomaterials. The deformation mechanisms responsible for the mechanical behavior of HAP are not well understood, although the advent of multiscale modeling offers the promise of improvements in many materials through computational materials science. This work utilizes molecular dynamics to study the nanoscale deformation mechanisms of HAP in uniaxial tension and compression. It was found that deformation mechanisms vary with loading direction in tension and compression leading to significant compression/tension asymmetry and crystal anisotropy. Bond orientation and geometry relative to the loading direction was found to be an indicator of whether a specific bond was involved in the deformation of HAP in each loading case. Tensile failure mechanisms were attributed to stretching and failure in loading case-specific ionic bond groups. The compressive failure mechanisms were attributed to coulombic repulsion in each case, although loading case-specific bond group rotation and displacement were found to affect specific failure modes. The elastic modulus was the highest for both tension and compression along the Z direction (i.e. normal to the basal plane), followed by Y and X.

Mechanics of Mineralized Collagen Fibrils upon Transient Loads

Collagen is a key structural protein in the human body, which undergoes mineralization during the formation of hard tissues. Earlier studies have described the mechanical behavior of bone at different scales, highlighting material features across hierarchical structures. Here we present a study that aims to understand the mechanical properties of mineralized collagen fibrils upon tensile/compressive transient loads, investigating how the kinetic energy propagates and it is dissipated at the molecular scale, thus filling a gap of knowledge in this area. These specific features are the mechanisms that nature has developed to passively dissipate stress and prevent structural failures. In addition to the mechanical properties of the mineralized fibrils, we observe distinct nanomechanical behaviors for the two regions (i.e., overlap and gap) of the D-period to highlight the effect of the mineralization. We notice decreasing trends for both wave speeds and Young's moduli over input velocity with a marked strengthening effect in the gap region due to the accumulation of the hydroxyapatite. In contrast, the dissipative behavior is not affected by either loading conditions or the mineral percentage, showing a stronger damping effect upon faster inputs compatible to the bone behavior at the macroscale. Our results offer insights into the dissipative behavior of mineralized collagen composites to design and characterize bioinspired composites for replacement devices (e.g., prostheses for sound transmission or conduction) or optimized structures able to bear transient loads, for example, impact, fatigue, in structural applications.

An embedded element based 2D finite element model for the strength prediction of mineralized collagen fibril using Monte-Carlo type of simulations

A computationally efficient statistical model for the prediction of the strength of mineralized collagen fibril (a basic building block of bone) is presented by taking into account the uncertainties associated with the geometrical and material parameters of collagen and mineral phases. The mineral plates have been considered as one-dimensional bar elements embedded in the two-dimensional plane stress collagen matrix. The mineral phase is considered as linear elastic and a hyperelastic material model is adopted for the collagen phase. Further, the crack initiation and propagation in the collagen phase have been modeled using a damage plasticity approach. Different realizations of the arrangement of mineral plates have been generated to account for the associated geometrical uncertainties using an in-house MATLAB® code. Monte-Carlo type simulations have been performed on the different realizations of mineralized collagen fibril to predict its characteristic stress-strain response under tensile load. The characteristic strength of 3.64 GPa is obtained for mineralized collagen fibril using Weibull's analysis which is found to be in agreement with the molecular dynamics simulation data and numerical studies reported in the past. A parameter sensitivity analysis concluded that mineral modulus has a significant effect on the overall tangent modulus of mineralized collagen fibril in large strain regime.

Polymer-Controlled Growth and Wrapping of Perovskite Single Crystals Leading to Better Device Stability and Performance

A commodity-scale polymer is used for controlling the nucleation and growth of single crystals of organolead halide perovskite. The polymer [polystyrene (PS)] cross-links and strongly interacts with PbI₂ and MAI (MAPbI₃ perovskite precursors) resulting in the control of the crystallization process. The PS concentration modulates the nucleation time, crystal size, and the number of perovskite single crystals. In addition, the PS-based MAPbI₃ crystals show an enhanced performance as well as improved thermal and environmental stability. Specifically, the PS-MAPbI₃ crystals show 3 times higher photocurrent than plain MAPbI₃ crystals and maintain a stable structure for more than 50 days (1200 h) under continuous 0.1 sun illumination in the air with a relative humidity of 40-45%. The improved performance and stability are attributed to the direct interaction between the PS and perovskite, which greatly reduces the ion migration, defect traps, and charge recombination and improves the carrier mobility and lifetime.

Comparative Study on Factors Governing Binding Mechanisms in Polylactic Acid-Hydroxyapatite and Polyethylene-Hydroxyapatite Systems via Molecular Dynamics Simulations

Binding mechanisms in polylactic acid-hydroxyapatite (PLA-HAp) and polyethylene-hydroxyapatite (PE-HAp) systems are comparatively elucidated on HAp (110) surfaces in unprecedented detail using molecular dynamics simulations conducted with the systematically varying number of monomers (N) between 10 and 400 at 310 K (NVT). Although PE seems to gradually cover the HAp surface more effectively compared to PLA, evident from the corresponding radius of gyration and occupied area values, the interface density and total binding energy in PLA-HAp systems is higher compared to those of PE-HAp systems. It is shown that a linear relationship between the binding energy and the surface area occupied by the monomer exists, consistent with our finding that binding energy converges to a limiting value with respect to monomer size on a constant surface area. The major constituent of the total binding energy is, rather surprisingly, shown to be the energy change in the bulk structure in HAp upon interaction; the next most important contributor is found to be the energy corresponding to surface-polymer interactions. The interplay between mainly these two contributors, acting in different fashions in two systems investigated here, seems to control the total binding energies. Increasing monomer size N initially results in enhanced densification of the interface in the HAp-PLA system up until N ≈ 200 with the positioning of mainly ═O units of PLA onto the HAp surface, consistent with the increasing Ca-O coordination numbers. Further increases in PLA size (N > 200) result in decreasing intensities of the peaks in the concentration profile consistent with the decreasing surface-polymer interaction energies while increased stabilization of the energy of the bulk is pronounced in this region. On the other hand, increasing N leads to a constantly increasing concentration at the interface in PE-HAp systems; -H atoms of the PE chain are positioned closer to the HAp surface than are -C atoms. These changes are coupled with increasing surface-polymer interaction energies in PE-HAp complexes, while slight destabilization in the energy of the bulk is observed for N > 100. A detailed examination of binding mechanisms in these technologically important systems as presented here is essential in material discovery; this valuable information, that will not be available from experiments can be attained through molecular simulations. The current study, to the best of our knowledge, comprises one of the first steps in achieving this goal for PLA/PE-HAp systems.

Microstructure, mineral and mechanical properties of teleost intermuscular bones

There is an increasing interest in understanding teleost bone biomechanics in several scientific communities, for instance as interesting biomaterials with specific structure-function relationships. Intermuscular bones of teleost fish have previously been described to play a role in the mechanical force transmission between muscle and bone, but their biomechanical properties are not yet fully described. Here, we have investigated intermuscular bones (IBs) of the North Atlantic Herring with regard to their structure and micro-architecture, mineral-related properties, and micro-mechanical tensile properties. A total of 115 IBs from 18 fish were investigated. One cohort of IBs, containing 20 bones from 2 smaller fish and 23 bones of 3 larger fish, was used for mechanical testing, wide-angle X-ray scattering, and scanning electron microscopy. Another cohort, containing 36 bones from 7 smaller fish and 36 bones from 6 larger fish, was used for microCT. Results show some astonishing properties of the IBs: (i) IBs present higher ductility, lower Young's modulus but similar strength and TMD (Tissue Mineral Density) compared to mammalian bone, and (ii) IBs from small fish were 49% higher in Young's modulus than fish bones from larger fish while their TMD was not statistically different and crystal length was 8% higher in large fish bones. Our results revealed that teleost IB presents a hybrid nature of soft and hard tissue that differs from other bone types, which might be associated with their evolution from mineralized tendons. This study provides new data regarding teleost fish bone biomechanical and micro-structural properties.

On the Origins of Fracture Toughness in Advanced Teleosts: How the Swordfish Sword's Bone Structure and Composition Allow for Slashing under Water to Kill or Stun Prey

The osseous sword of a swordfish (Xiphias gladius) is specialized to incapacitate prey with stunning blows. Considering the sword's growth and maturation pattern, aging from the sword's base to the tip, while missing a mechanosensitive osteocytic network, an in-depth understanding of its mechanical properties and bone quality is lacking. Microstructural, compositional, and nanomechanical characteristics of the bone along the sword are investigated to reveal structural mechanisms accounting for its exceptional mechanical competence. The degree of mineralization, homogeneity, and particle size increase from the base toward the tip, reflecting aging along its length. Fracture experiments reveal that crack-growth toughness vastly decreases at the highly and homogeneously mineralized tip, suggesting the importance of aging effects. Initiation toughness, however, is unchanged suggesting that aging effects on this hierarchical level are counteracted by constant mineral/fibril interaction. In conclusion, the sword of the swordfish provides an excellent model reflecting base-to-tip-wise aging of bone, as indicated by increasing mineralization and decreasing crack-growth toughness toward the tip. The hierarchical, structural, and compositional changes along the sword reflect peculiar prerequisites needed for resisting high mechanical loads. Further studies on advanced teleosts bone tissue may help to unravel structure-function relationships of heavily loaded skeletons lacking mechanosensing cells.

Analysis of the vibrational and sound spectrum of over 100,000 protein structures and application in sonification

We report a high-throughput method that enables us to automatically compute the vibrational spectra of more than 100,000 proteins available in the Protein Data Bank to date, in a consistent manner. Using this new algorithm we report a comprehensive database of the normal mode frequencies of all known protein structures, which has not been available before. We then use the resulting frequency spectra of the proteins to generate audible sound by overlaying the molecular vibrations and translating them to the audible frequency range using the music theoretic concept of transpositional equivalence. The method, implemented as a Max audio device for use in a digital audio workstation (DAW), provides unparalleled insights into the rich vibrational signatures of protein structures, and offers a new way for creative expression by using it as a new type of musical instrument. This musical instrument is fully defined by the vibrational feature of almost all known protein structures, making it fundamentally different from all the traditional instruments that are limited by the material properties of a few types of conventional engineering materials, such as wood, metals or polymers.

Natural hydrogel in American lobster: A soft armor with high toughness and strength

Homarus americanus, known as American lobster, is fully covered by its exoskeleton composed of rigid cuticles and soft membranes. These soft membranes are mainly located at the joints and abdomen to connect the rigid cuticles and greatly contribute to the agility of the lobster in swimming and preying. Herein, we show that the soft membrane from American lobster is a natural hydrogel (90% water) with exceptionally high toughness (up to 24.98 MJ/m³) and strength (up to 23.36 MPa), and is very insensitive to cracks. By combining experimental measurements and large-scale computational modeling, we demonstrate that the unique multilayered structure in this membrane, achieved through the ordered arrangement of chitin fibers, plays a crucial role in dissipating energy during rupture and making this membrane tough and damage tolerant. The knowledge learned from the soft membrane of natural lobsters sheds light on designing synthetic soft, yet strong and tough materials for reliable usage under extreme mechanical conditions, including a flexible armor that can provide full-body protection without sacrificing limb mobility. STATEMENT OF SIGNIFICANCE: A body armor to provide protection to people who are at risk of being hurt is only enabled by using a material that is tough and strong enough to prevent mechanical penetration. However, most modern body armors sacrifice limb protection to gain mobility, simply because none of the existing armor materials are flexible enough and they all inhibit movement of the arms and legs. Herein, we focus on the mechanics and mesoscopic structure of American lobsters' soft membrane and explore how such a natural flexible armor is designed to integrate flexibility and toughness. The knowledge learned from this study is useful to design a flexible armor for full-body protection under extreme mechanical conditions.

Computational investigation of ultrastructural behavior of bone using a cohesive finite element approach

Bone ultrastructure at sub-lamellar length scale is a key structural unit in bone that bridges nano- and microscale hierarchies of the tissue. Despite its influence on bulk response of bone, the mechanical behavior of bone at ultrastructural level remains poorly understood. To fill this gap, in this study, a two-dimensional cohesive finite element model of bone at sub-lamellar level was proposed and analyzed under tensile and compressive loading conditions. In the model, ultrastructural bone was considered as a composite of mineralized collagen fibrils (MCFs) embedded in an extrafibrillar matrix (EFM) that is comprised of hydroxyapatite (HA) polycrystals bounded via thin organic interfaces of non-collagenous proteins (NCPs). The simulation results indicated that in compression, EFM dictated the pre-yield deformation of the model, then damage was initiated via relative sliding of HA polycrystals along the organic interfaces, and finally shear bands were formed followed by delamination between MCF and EFM and local buckling of MCF. In tension, EFM carried the most of load in pre-yield deformation, and then an array of opening-mode nano-cracks began to form within EFM after yielding, thus gradually transferring the load to MCF until failure, which acted as crack bridging filament. The failure modes, stress-strain curves, and in situ mineral strain of ultrastructural bone predicted by the model were in good agreement with the experimental observations reported in the literature, thus suggesting that this model can provide new insights into sub-microscale mechanical behavior of bone.

The effect of polyaspartate chain length on mediating biomimetic remineralization of collagenous tissues

Formation of hydroxyapatite (HAP) within collagen fibrils, as found in bone, dentine and cementum, is thought to be mediated by proteins rich in aspartate (Asp) and glutamate such as osteopontin and bone sialoprotein, respectively. Indeed polyaspartate (pAsp), a homopolymer analogue of such proteins, has been shown to induce intrafibrillar mineralization of collagen from solutions of calcium and phosphate that are supersaturated with respect to HAP. To elucidate the role of pAsp in mineralization of collagen, we explored the effect of pAsp chain length on in vitro HAP deposition in demineralized mouse periodontal tissue sections. Through characterization of both tissue sections and mineralizing solution, we show that chain length contributes to the effectiveness of pAsp in mediating intrafibrillar mineralization. This function appears to be associated with inhibition of otherwise kinetically favoured crystallization in the bulk solution, which allows for intrafibrillar crystallization, though this does not preclude the possibility of a more active role for pAsp in addition. Inhibition of crystallization in solution by pAsp occurs by slowing the growth of amorphous calcium phosphate and stabilization of this phase, rather than by sequestration of Ca²⁺ ions. These results suggest that the length of Asp-rich sequences of mineralizing proteins may be essential to their function, and could also be useful in optimization of mineralized tissue replacement synthesis.

Structural role of osteocalcin and osteopontin in energy dissipation in bone

Non-collagenous proteins are a vital component of bone matrix. Amongst them, osteocalcin (OC) and osteopontin (OPN) hold special significance due to their intimate interaction with the mineral and collagenous matrix in bone. Both proteins have been associated with microdamage and fracture, but their structural role in energy dissipation is unclear. This study used bone tissue from genetic deficient mice lacking OC and/or OPN and subjected them to a series of creep-fatigue-creep tests. To this end, whole tibiae were loaded in four-point bending to 70% stiffness loss which captured the three characteristic phases of fatigue associated with initiation, propagation, and coalescence of microdamage. Fatigue loading preceded and followed creep tests to determine creep and dampening parameters. Microdamage in the form of linear microcracks and diffuse damage were analyzed by histology. It was shown that OC and OPN were 'activated' following stiffness loss associated with fatigue damage where they facilitated creep and dampening parameters (i.e. increased energy dissipation). More specifically, post-fatigue creep rate and dampening were significantly greater in wild-types (WTs) than genetic deficient mice (p < 0.05). These results were supported by microdamage analysis which showed significant increase in creep-associated diffuse damage formation in WTs compared to genetic deficient groups (p < 0.05). Based on these findings, we propose that during local yield events, OC and OPN rely on ionic interactions of their charged side chains and on hydrogen bonding to dissipate energy in bone.

Biomimetic Bone-like Hydroxyapatite by Mineralization on Supramolecular Porous Fiber Networks

Hydroxyapatite (HA), the main inorganic component of bone tissue, is mineralized with collagen fibril scaffolds during bone formation. Inspired by the process, a self-assembled porous network architecture was designed and synthesized by using the 2-ureido-4[1H]-pyrimidone (UPy) modified glycerol molecule UPy-Gly, which was further utilized as a template for biomimetic mineralization. When incubated in simulated body fluid (SBF), the HA nucleus first formed in the holes of the template by the induction of hydroxyls on the surface, grew along the nanofibers, and fused with the template to fabricate hydroxyapatite composites (UPy-Gly/HA). Transmission electron microscopic observation demonstrates that the mineral clusters are accumulated by lamella-like nano hydroxyapatite and the elasticity modulus measured by atomic force microscopy is about 5.5 GPa, which is quite close to the natural cancellous bone tissue of human both in structure and in mechanical properties. The Cell Counting Kit 8 (CCK-8) assay of UPy-Gly and UPy-Gly/HA shows noncytotoxicity to mouse fibroblast L-929 cells. This bioinspired composite will be a promising material for potential use in bone tissue implantation and regeneration engineering.

Clustering of hydroxyapatite on a super-twisted collagen microfibril under mechanical tension

Atomistic simulation of biomineralization of a super-twisted collagen microfibril reveals that mechanical stimulation facilitates clustering and growth of hydroxyapatite onto collagen.

An improved interfacial bonding model for material interface modeling

An improved interfacial bonding model was proposed from potential function point of view to investigate interfacial interactions in polycrystalline materials. It characterizes both attractive and repulsive interfacial interactions and can be applied to model different material interfaces. The path dependence of work-of-separation study indicates that the transformation of separation work is smooth in normal and tangential direction and the proposed model guarantees the consistency of the cohesive constitutive model. The improved interfacial bonding model was verified through a simple compression test in a standard hexagonal structure. The error between analytical solutions and numerical results from the proposed model is reasonable in linear elastic region. Ultimately, we investigated the mechanical behavior of extrafibrillar matrix in bone and the simulation results agreed well with experimental observations of bone fracture.

Nanocrystalline hydroxyapatite intervertebral cages induce fusion after anterior cervical discectomy and may be a safe alternative to PEEK or carbon fiber intervertebral cages

PURPOSE

Nanocrystalline hydroxyapatite (nHA) cages have emerged as a new alternative to carbon fiber or polyether ether ketone (PEEK) devices to promote intervertebral fusion. No evidence has been published to date regarding rates of fusion for these devices after anterior cervical discectomy and fusion (ACDF).

METHODS

Eight patients underwent one- or two-level ACDF with nHA intervertebral cages (Nanoss^®-Cervical, Pioneer^® Surgical Technology, Inc., Marquette, MI). Radiographs, neck disability index (NDI) and visual analog scores (VAS) for pain were taken preoperatively and at a minimum of 19 months postoperatively.

RESULTS

At an average follow-up of 21 months, all eight patients (100%) achieved fusion as assessed by plain radiographs. Reduction in preoperative symptomology was comparable to previously published data with a mean reduction of neck VAS of 3, arm VAS of 6 and NDI reduced by 27%. Radiographs showed clear evidence of bridging bone.

CONCLUSIONS

This series provides evidence that nHA intervertebral cages can successfully promote fusion after ACDF and may provide an alternative to carbon fiber and PEEK cages.

Isoexergonic Conformations of Surface-Bound Citrate Regulated Bioinspired Apatite Nanocrystal Growth

The superior biomechanical properties of bone and dentin are dictated, in part, by the unique plate-like morphology of hydroxyapatite (HAP) nanocrysals within a hierarchically assembled collagen matrix. Understanding the mechanism of crystal growth and thus morphology is important to the rational design of bioinspired apatite nanocrystals for orthopedic and dental applications. Citrate has long been proposed to modulate apatite crystal growth, but major questions exist regarding the HAP-bound citrate conformations and the identities of the interacting functional groups and HAP surface sites. Here, we conducted a comprehensive investigation of the mechanism from the angstrom to submicrometer scale by detailed correlation of the results of high-level metadynamics simulations, employing force-fields benchmarked to experiment and density functional theory calculations with the results of high resolution transmission electron microscopy, nuclear magnetic resonance spectroscopy, solution analysis, and thermogravimetric analysis. Crystal morphology changed from needle- to plate-like with increasing citrate concentration. Citrate adsorbed more strongly on the HAP (100) face than on the (001) face, thus resulting in preferential growth in the [001] direction and the plate-like morphology. Two very different bound conformations were obtained, involving interactions of either one or both terminal carboxyl groups with three or five surface calcium ions, respectively, and a hydrogen bond between the citrate hydroxyl and the HAP surface. Remarkably, despite fewer interaction sites in the single bound carboxyl conformation, the structures were isoexergonic, so both exist at equilibrium. Identification of the former conformation is significant because it allows a greater adsorption density than is traditionally assumed and can help explain concentration-dependence of citrate in modulating crystal morphology. These unique results were enabled first by the application of advanced metadynamics, a technique necessary for the accurate simulation of ionic materials but which is rarely employed in the biomaterials and biomineralization fields and second by the detailed correlation of computational, spectroscopic, and analytical results.

Hierarchical Structures of Bone and Bioinspired Bone Tissue Engineering

Bone, as a mineralized composite of inorganic (mostly carbonated hydroxyapatite) and organic (mainly type I collagen) phases, possesses a unique combination of remarkable strength and toughness. Its excellent mechanical properties are related to its hierarchical structures and precise organization of the inorganic and organic phases at the nanoscale: Nanometer-sized hydroxyapatite crystals periodically deposit within the gap zones of collagen fibrils during bone biomineralization process. This hierarchical arrangement produces nanomechanical heterogeneities, which enable a mechanism for high energy dissipation and resistance to fracture. The excellent mechanical properties integrated with the hierarchical nanostructure of bone have inspired chemists and material scientists to develop biomimetic strategies for artificial bone grafts in tissue engineering (TE). This critical review provides a broad overview of the current mechanisms involved in bone biomineralization, and the relationship between bone hierarchical structures and the deformation mechanism. Our goal in this review is to inspire the application of these principles toward bone TE.

Contribution of extrafibrillar matrix to the mechanical behavior of bone using a novel cohesive finite element model

The mechanical behavior of bone is determined at all hierarchical levels, including lamellae (the basic building block of bone) that are comprised of mineralized collagen fibrils and extrafibrillar matrix. The mechanical behavior of mineralized collagen fibrils has been investigated intensively using both experimental and computational approaches. Yet, the contribution of the extrafibrillar matrix to bone mechanical properties is poorly documented. In this study, we intended to address this issue using a novel cohesive finite element (FE) model, in conjunction with the experimental observations reported in the literature. In the FE model, the extrafibrillar matrix was considered as a nanocomposite of hydroxyapatite (HA) crystals bounded through a thin organic interface modeled as a cohesive interfacial zone. The parameters required by the cohesive FE model were defined based on the experimental data reported in the literature. This hybrid nanocomposite model was tested in two loading modes (i.e. tension and compression) and under two hydration conditions (i.e. wet and dry). The simulation results indicated that (1) the failure modes of the extrafibrillar matrix predicted using the cohesive FE model were closely coincided with those experimentally observed in tension and compression tests; (2) the pre-yield deformation (i.e. internal strain) of HA crystals with respect to the applied strain was consistent with that obtained from the synchrotron X-ray scattering measurements irrespective of the loading modes and hydration status; and (3) the mechanical behavior of the extrafibrillar matrix was dictated by the properties of the organic interface between the HA crystals. Taken together, we postulate that the extrafibrillar matrix plays a major role in the pre-yield deformation and the failure mode of bone, thus, giving rise to important insights in the ultrastructural origins of bone fragility.

Biomimetic composite scaffolds containing bioceramics and collagen/gelatin for bone tissue engineering - A mini review

Bone is a natural composite material consisting of an organic phase (collagen) and a mineral phase (calcium phosphate, especially hydroxyapatite). The strength of bone is attributed to the apatite, while the collagen fibrils are responsible for the toughness and visco-elasticity. The challenge in bone tissue engineering is to develop such biomimetic composite scaffolds, having a balance between biological and biomechanical properties. This review summarizes the current state of the field by outlining composite scaffolds made of gelatin/collagen in combination with bioactive ceramics for bone tissue engineering application.

Remodeling by fibroblasts alters the rate-dependent mechanical properties of collagen

The ways that fibroblasts remodel their environment are central to wound healing, development of musculoskeletal tissues, and progression of pathologies such as fibrosis. However, the changes that fibroblasts make to the material around them and the mechanical consequences of these changes have proven difficult to quantify, especially in realistic, viscoelastic three-dimensional culture environments, leaving a critical need for quantitative data. Here, we observed the mechanisms and quantified the mechanical effects of fibroblast remodeling in engineered tissue constructs (ETCs) comprised of reconstituted rat tail (type I) collagen and human fibroblast cells. To study the effects of remodeling on tissue mechanics, stress-relaxation tests were performed on ETCs cultured for 24, 48, and 72h. ETCs were treated with deoxycholate and tested again to assess the ECM response. Viscoelastic relaxation spectra were obtained using the generalized Maxwell model. Cells exhibited viscoelastic damping at two finite time constants over which the ECM showed little damping, approximately 0.2s and 10-30s. Different finite time constants in the range of 1-7000s were attributed to ECM relaxation. Cells remodeled the ECM to produce a relaxation time constant on the order of 7000s, and to merge relaxation finite time constants in the 0.5-2s range into a single time content in the 1s range. Results shed light on hierarchical deformation mechanisms in tissues, and on pathologies related to collagen relaxation such as diastolic dysfunction.

STATEMENT OF SIGNIFICANCE

As fibroblasts proliferate within and remodel a tissue, they change the tissue mechanically. Quantifying these changes is critical for understanding wound healing and the development of pathologies such as cardiac fibrosis. Here, we characterize for the first time the spectrum of viscoelastic (rate-dependent) changes arising from the remodeling of reconstituted collagen by fibroblasts. The method also provides estimates of the viscoelastic spectra of fibroblasts within a three-dimensional culture environment. Results are of particular interest because of the ways that fibroblasts alter the mechanical response of collagen at loading frequencies associated with cardiac contraction in humans.

Surface Structure of Hydroxyapatite from Simulated Annealing Molecular Dynamics Simulations

The surface structure of hydroxyapatite (HAP) is crucial for its bioactivity. Using a molecular dynamics simulated annealing method, we studied the structure and its variation with annealing temperature of the HAP (100) surface. In contrast to the commonly used HAP surface model, which is sliced from HAP crystal and then relaxed at 0 K with first-principles or force-field calculations, a new surface structure with gradual changes from ordered inside to disordered on the surface was revealed. The disordering is dependent on the annealing temperature, Tmax. When Tmax increases up to the melting point, which was usually adopted in experiments, the disordering increases, as reflected by its radial distribution functions, structural factors, and atomic coordination numbers. The disordering of annealed structures does not show significant changes when Tmax is above the melting point. The thickness of disordered layers is about 10 Å. The surface energy of the annealed structures at high temperature is significantly less than that of the crystal structure relaxed at room temperature. A three-layer model of interior, middle, and surface was then proposed to describe the surface structure of HAP. The interior layer retains the atomic configurations in crystal. The middle layer has its atoms moved and its groups rotated about their original locations. In the surface layer, the atomic arrangements are totally different from those in crystal. In particular for the hydroxyl groups, they move outward and cover the Ca(2+) ions, leaving holes occupied by the phosphate groups. Our study suggested a new model with disordered surface structures for studying the interaction of HAP-based biomaterials with other molecules.

c-axis preferential orientation of hydroxyapatite accounts for the high wear resistance of the teeth of black carp (Mylopharyngodon piceus)

Biological armors such as mollusk shells have long been recognized and studied for their values in inspiring novel designs of engineering materials with higher toughness and strength. However, no material is invincible and biological armors also have their rivals. In this paper, our attention is focused on the teeth of black carp (Mylopharyngodon piceus) which is a predator of shelled mollusks like snails and mussels. Nanoscratching test on the enameloid, the outermost layer of the teeth, indicates that the natural occlusal surface (OS) has much higher wear resistance compared to the other sections. Subsequent X-ray diffraction analysis reveals that the hydroxyapatite (HAp) crystallites in the vicinity of OS possess c-axis preferential orientation. The superior wear resistance of black carp teeth is attributed to the c-axis preferential orientation of HAp near the OS since the (001) surface of HAp crystal, which is perpendicular to the c-axis, exhibits much better wear resistance compared to the other surfaces as demonstrated by the molecular dynamics simulation. Our results not only shed light on the origin of the good wear resistance exhibited by the black carp teeth but are of great value to the design of engineering materials with better abrasion resistance.

Modelling of bone fracture and strength at different length scales: a review

In this paper, we review analytical and computational models of bone fracture and strength. Bone fracture is a complex phenomenon due to the composite, inhomogeneous and hierarchical structure of bone. First, we briefly summarize the hierarchical structure of bone, spanning from the nanoscale, sub-microscale, microscale, mesoscale to the macroscale, and discuss experimental observations on failure mechanisms in bone at these scales. Then, we highlight representative analytical and computational models of bone fracture and strength at different length scales and discuss the main findings in the context of experiments. We conclude by summarizing the challenges in modelling of bone fracture and strength and list open topics for scientific exploration. Modelling of bone, accounting for different scales, provides new and needed insights into the fracture and strength of bone, which, in turn, can lead to improved diagnostic tools and treatments of bone diseases such as osteoporosis.

Large Deformation Mechanisms, Plasticity, and Failure of an Individual Collagen Fibril With Different Mineral Content

Mineralized collagen fibrils are composed of tropocollagen molecules and mineral crystals derived from hydroxyapatite to form a composite material that combines optimal properties of both constituents and exhibits incredible strength and toughness. Their complex hierarchical structure allows collagen fibrils to sustain large deformation without breaking. In this study, we report a mesoscale model of a single mineralized collagen fibril using a bottom-up approach. By conserving the three-dimensional structure and the entanglement of the molecules, we were able to construct finite-size fibril models that allowed us to explore the deformation mechanisms which govern their mechanical behavior under large deformation. We investigated the tensile behavior of a single collagen fibril with various intrafibrillar mineral content and found that a mineralized collagen fibril can present up to five different deformation mechanisms to dissipate energy. These mechanisms include molecular uncoiling, molecular stretching, mineral/collagen sliding, molecular slippage, and crystal dissociation. By multiplying its sources of energy dissipation and deformation mechanisms, a collagen fibril can reach impressive strength and toughness. Adding mineral into the collagen fibril can increase its strength up to 10 times and its toughness up to 35 times. Combining crosslinks with mineral makes the fibril stiffer but more brittle. We also found that a mineralized fibril reaches its maximum toughness to density and strength to density ratios for a mineral density of around 30%. This result, in good agreement with experimental observations, attests that bone tissue is optimized mechanically to remain lightweight but maintain strength and toughness.

Assessment of the interacting mechanism between Candida rugosa lipases and hydroxyapatite and identification of the hydroxyapatite-binding sequence through proteomics and molecular modelling

Hydroxyapatite (HAP), a calcium-phosphate bioactive ceramic, is actively employed in medical and separation sciences.

A potential mechanism for amino acid-controlled crystal growth of hydroxyapatite

The mineral component of bone, dentin and calcified parts of avian tendon, hydroxyapatite (HAP), has non-stoichiometric composition (idealized as Ca₁₀(PO₄)₆(OH)₂), plate-like morphology and nanometer size. This unique crystal morphology contributes to the physico-chemical and biochemical properties of bone. Thus, understanding the mechanism for the controlled growth of plate-like HAP nanocrystals is significant in the study of bone biomineralization. Previous studies have shown that acidic non-collagenous proteins (ANCPs), which are enriched in the residues of acidic amino acids, may play an important role in HAP crystal growth modulation. In this study, glutamic acid (Glu) and phosphoserine (Ser-OPO₃) were used as model compounds to modify the synthesis of HAP nanocrystals. To identify the mechanisms of amino acids as regulators, X-ray diffraction (XRD), transmission electron microscopy (TEM) and solid state nuclear magnetic resonance (ssNMR) were used. The crystals obtained in the inorganic controls were needle-like, while crystals synthesized in the presence of the amino acids presented a plate-like morphology. The plate-like crystals had a preferred crystal orientation on (300) face, which was lacking in the inorganically grown crystals, indicating preferential adsorption and suppression of growth in specific crystal directions. Ser-OPO₃ was more efficient than Glu in modulating HAP nucleation and crystal growth. Furthermore, NMR revealed interactions between the charged side chain groups in amino acids and the crystal surfaces. These results were successfully explained through our MD simulations for the free energy calculation of amino acid binding on HAP crystal faces. The present study revealed that amino acids may act as effective regulators of HAP morphology without the need to invoke large NCPs in bone biomineralization and in designing bioinspired materials for orthopaedic and dental applications.