Anisotropic elastic properties of human femoral cortical bone and relationships with composition and microstructure in elderly

3D quantification of the lacunocanalicular network on human femoral diaphysis through synchrotron radiation-based nanoCT

Osteocytes are the major actors in bone mechanobiology. Within bone matrix, they are trapped close together in a submicrometric interconnected network: the lacunocanalicular network (LCN). The interstitial fluid circulating within the LCN transmits the mechanical information to the osteocytes that convert it into a biochemical signal. Understanding the interstitial fluid dynamics is necessary to better understand the bone mechanobiology. Due to the submicrometric dimensions of the LCN, making it difficult to experimentally investigate fluid dynamics, numerical models appear as a relevant tool for such investigation. To develop such models, there is a need for geometrical and morphological data on the human LCN. This study aims at providing morphological data on the human LCN from measurement of 27 human femoral diaphysis bone samples using synchrotron radiation nano-computed tomography with an isotropic voxel size of 100 nm. Except from the canalicular diameter, the canalicular morphological parameters presented a high variability within one sample. Some differences in terms of both lacunar and canalicular morphology were observed between the male and female populations. But it has to be highlighted that all the canaliculi cannot be detected with a voxel size of 100 nm. Hence, in the current study, only a specific population of large canaliculi that could be characterize. Still, to the authors knowledge, this is the first time such a data set was introduced to the community. Further processing will be achieved in order to provide new insight on the LCN permeability.

Biomimetic cryogel promotes the repair of osteoporotic bone defects through altering the ROS niche via down-regulating the ROMO1

Osteoporosis is a systemic bone disease that is prone to fractures due to decreased bone density and bone quality, and delayed union or nonunion often occurs in osteoporotic fractures. Therefore, it is particularly important to develop tissue engineering materials to promote osteoporotic fracture healing. In this study, a series of biomimetic cryogels prepared from the decellularized extracellular matrix (dECM), methacrylate gelatin (GelMA), and carboxymethyl chitosan (CMCS) via unidirectional freezing, photo- and genipin crosslinking were applied for the regeneration of osteoporotic fractures. Specifically, dECM extracted from normal or osteoporotic rats was applied for the preparation of the cryogels, named as GC-Normal dECM or GC-OVX dECM, respectively. It was verified that the GC-Normal dECM demonstrated superior performance in promoting the proliferation of BMSCs isolated from osteoporotic rats (OVX-BMSCs), and the differentiation of OVX-BMSCs into osteoblasts both in vitro and in vivo. RNA sequencing and further verifications confirmed that GC-Normal dECM cryogel could scavenge the intracellular reactive oxygen species (ROS) in OVX-BMSCs to accelerate the regeneration of osteoporotic fracture by down-regulating the reactive oxygen species modulator 1 (Romo1). The results indicated that by regulating the ROS niche of OVX-BMSCs, biomimetic the GC-Normal dECM cryogel was expected to be a clinical candidate for repairing osteoporotic bone defects.

Effects of different running intensities on the micro-level failure strain of rat femoral cortical bone structures: a finite element investigation

BACKGROUND

Running with the appropriate intensity may produce a positive influence on the mechanical properties of cortical bone structure. However, few studies have discussed the effects of different running intensities on the mechanical properties at different levels, especially at the micro-level, because the micromechanical parameters are difficult to measure experimentally.

METHODS

An approach that combines finite element analysis and experimental data was proposed to predict a micromechanical parameter in the rat femoral cortical bone structure, namely, the micro-level failure strain. Based on the previous three-point bending experimental information, fracture simulations were performed on the femur finite element models to predict their failure process under the same bending load, and the micro-level failure strains in tension and compression of these models were back-calculated by fitting the experimental load-displacement curves. Then, the effects of different running intensities on the micro-level failure strain of rat femoral cortical bone structure were investigated.

RESULTS

The micro-level failure strains of the cortical bone structures expressed statistical variations under different running intensities, which indicated that different mechanical stimuli of running had significant influences on the micromechanical properties. The greatest failure strain occurred in the cortical bone structure under low-intensity running, and the lowest failure strain occurred in the structure under high-intensity running.

CONCLUSIONS

Moderate and low-intensity running were effective in enhancing the micromechanical properties, whereas high-intensity running led to the weakening of the micromechanical properties of cortical bone. Based on these, the changing trends in the micromechanical properties were exhibited, and the effects of different running intensities on the fracture performance of rat cortical bone structures could be discussed in combination with the known mechanical parameters at the macro- and nano-levels, which provided the theoretical basis for reducing fracture incidence through running exercise.

Tissue-Engineered Nanomaterials Play Diverse Roles in Bone Injury Repair

Nanomaterials with bone-mimicking characteristics and easily internalized by the cell could create suitable microenvironments in which to regulate the therapeutic effects of bone regeneration. This review provides an overview of the current state-of-the-art research in developing and using nanomaterials for better bone injury repair. First, an overview of the hierarchical architecture from the macroscale to the nanoscale of natural bone is presented, as these bone tissue microstructures and compositions are the basis for constructing bone substitutes. Next, urgent clinical issues associated with bone injury that require resolution and the potential of nanomaterials to overcome them are discussed. Finally, nanomaterials are classified as inorganic or organic based on their chemical properties. Their basic characteristics and the results of related bone engineering studies are described. This review describes theoretical and technical bases for the development of innovative methods for repairing damaged bone and should inspire therapeutic strategies with potential for clinical applications.

The influence of intra-cortical microstructure on the contrast in ultrasound images of the cortex of long bones: A 2D simulation study

Decreased thickness of the bone cortex due to bone loss in the course of ageing and osteoporosis is associated with reduced bone strength. Cortical thickness measurement from ultrasound images was recently demonstrated in young adults. This requires the identification of both the outer (periosteum) and inner (endosteum) surfaces of the bone cortex. However, with bone loss, the cortical porosity and the size of the vascular pores increase resulting in enhanced ultrasound scattering which may prevent the detection of the endosteum. The aim of this work was to study the influence of cortical bone microstructure variables, such as porosity and pore size, on the contrast of the endosteum in ultrasound images. We wanted to estimate the range of these variables for which ultrasound imaging of the endosteum is feasible. We generated synthetic data using a two-dimensional time-domain code to simulate the propagation of elastodynamic waves. A synthetic aperture imaging sequence with an array transducer operating at a center frequency of 2.5 MHz was used. The numerical simulations were conducted for 105 cortical microstructures obtained from high resolution X-ray computed tomography images of ex vivo bone samples with a porosity ranging from 2% to 24 %. Images were reconstructed using a delay-and-sum (DAS) algorithm with optimized f-number, correction of refraction at the periosteum, and sample-specific wave-speed. We observed a range variation of 18 dB of endosteum contrast in our data set depending on the bone microstructure. We found that as porosity increases, speckle intensity inside the bone cortex increases whereas the intensity of the signal from the endosteum decreases. Also, a microstructure with large pores (diameter >250 μm) was associated with poor endosteum visibility, compared with a microstructure with equal porosity but a more narrow distribution of pore sizes. These findings suggest that ultrasound imaging of the bone cortex with a probe operating at a central frequency of 2.5 MHz using refraction-corrected DAS is capable of detecting the endosteum of a cortex with moderate porosity (less than about 10%) if the largest pores remain smaller than about 200 μm.

Using ultrasonic attenuation in cortical bone to infer distributions on pore size

In this work we infer the underlying distribution on pore radius in human cortical bone samples using ultrasonic attenuation data. We first discuss how to formulate polydisperse attenuation models using a probabilistic approach and the Waterman Truell model for scattering attenuation. We then compare the Independent Scattering Approximation and the higher-order Waterman Truell models' forward predictions for total attenuation in polydisperse samples. Following this, we formulate an inverse problem under the Prohorov Metric Framework coupled with variational regularization to stabilize this inverse problem. We then use experimental attenuation data taken from human cadaver samples and solve inverse problems resulting in nonparametric estimates of the probability density function on pore radius. We compare these estimates to the "true" microstructure of the bone samples determined via microCT imaging. We find that our methodology allows us to reliably estimate the underlying microstructure of the bone from attenuation data.

A preliminary study exploring the mechanical properties of normal and Mgp-deficient mouse femurs during early growth

Matrix Gla protein (MGP) is mostly known to be a calcification inhibitor, as its absence leads to ectopic calcification of different tissues such as cartilage or arteries. MGP deficiency also leads to low bone mass and delayed bone growth. In the present contribution, we investigate the effect of MGP deficiency on the structural and material mechanical bone properties by focusing on the elastic response of femurs undergoing three-points bending. To this aim, biomechanical tests are performed on femurs issued from Mgp-deficient mice at 14, 21, 28, and 35 days of postnatal life and compared to healthy control femurs. µCT acquisitions enable to reconstruct bone geometries and are used to construct subject-specific finite element models avoiding some of the reported limitations concerning the use of beam-like assumptions for small bone samples. Our results indicate that MGP deficiency may be associated to differences in both structural and material properties of femurs during early stages of development. MGP deficiency appears to be related to a decrease in bone dimensions, compensated by higher material properties resulting in similar structural bone properties at P35. The search for a unique density-elasticity relationship based on calibrated bone mineral density (BMD) indicates that MGP deficiency may affect bone tissue in several ways, that may not be represented uniquely from the quantification of BMD. Despite of its limitation to elastic response, the present preliminary study reports for the very first time the mechanical skeletal properties of Mgp-deficient mice at early stages of development.

Inverse Problem in Resonant Ultrasound Spectroscopy With Sampling and Optimization: A Comparative Study on Human Cortical Bone

The Bayesian inference with prior knowledge has been proposed recently to solve the inverse problem in resonant ultrasound spectroscopy. It allows inferring the elastic properties of high damping materials, such as cortical bone with less dependence on the initial guessed values. In this method, the estimation of the stiffness coefficients is expressed as a probabilistic solution to the inverse problem, which can be achieved by sampling or optimization methods. However, the detailed performance comparison of these two strategies applied to high damping materials has not been fully studied. In this work, the full stiffness tensor of 52 transversely isotropic cortical bone specimens was obtained using Markov chain Monte Carlo (MCMC) sampling and particle swarm optimization (PSO), respectively. Results showed that the local probability distributions of stiffness coefficients estimated by the two methods are consistent. Compared with MCMC, the average calculation speed of PSO is ten times faster [614 s ± 59 s (MCMC) versus 53 s ± 22 s (PSO)]. The mean standard error between theoretical and experimental resonant frequencies was slightly smaller for PSO compared with MCMC. In conclusion, PSO, a global optimization strategy, is suitable to solve the inverse problem for high damping materials.

Ultrasound Scattering in Cortical Bone

Recent advances in imaging of bone microstructure have led to a growing recognition of the role of cortical microstructure in osteoporosis. It is now accepted that the assessment of the microstructure of cortical porosity is essential to assess bone mechanical competence and predict fracture risk. Cortical porosity affects the propagation of ultrasound waves because pores act as ultrasound scatterers. Scattering by the porosity is an opportunity that should be leveraged to extract quantitative information about cortical microstructure. Scattering by the pores affects a number of ultrasound parameters that should be quantified, including attenuation, backscatter coefficient, ultrasound diffusivity, and their frequency dependence. Measuring these ultrasound parameters and developing models that describe their dependence upon parameters of cortical microstructure is the key to solve inverse problems that will allow the quantitative assessment of cortical porosity and ultimately will improve the non-invasive ultrasound-based evaluation of bone mechanical competence and fracture risk. In this chapter, we present recent advances in measuring and modeling those parameters in cortical bone.

The spatial form periosteal-bone complex promotes bone regeneration by coordinating macrophage polarization and osteogenic-angiogenic events

Bone defects associated with soft tissue injuries are an important cause of deformity that threatens people’s health and quality of life. Although bone substitutes have been extensively explored, effective biomaterials that can coordinate early inflammation regulation and subsequent repair events are still lacking. We prepared a spatial form periosteal bone extracellular matrix (ECM) scaffold, which has advantages in terms of low immunogenicity, good retention of bioactive ingredients, and a natural spatial structure. The periosteal bone ECM scaffold with the relatively low-stiffness periosteum (41.6 ± 3.7 kPa) could inhibit iNOS and IL-1β expression, which might be related to actin-mediated YAP translocation. It also helped to promote CD206 expression with the potential influence of proteins related to immune regulation. Moreover, the scaffold combined the excellent properties of decalcified bone and periosteum, promoted the formation of blood vessels, and good osteogenic differentiation (RUNX2, Col 1α1, ALP, OPN, and OCN), and achieved good repair of a cranial defect in rats. This scaffold, with its natural structural and biological advantages, provides a new idea for bone healing treatment that is aligned with bone physiology.

•

We provided a spatial form periosteal-bone complex.

•

The scaffold preserved major biological components and spatial structure.

•

The periosteum part of the scaffold acted as a physical barrier.

•

The scaffold participated in the transformation of the macrophage phenotype.

•

The scaffold promoted osteogenesis and angiogenesis.

Mapping cortical bone stiffness and mineralization from endosteal to periosteal surfaces of bovine mid-diaphyseal femur

INTRODUCTION

While bone literature abounds with correlations of mechanical stiffness to mineralization, such correlations are reported without relating the findings to specific intracortical locations. This study reports on mapping of stiffness and mineralization distributions in ring-shaped cortical bone samples sliced from mid-diaphyseal bovine femur. Stiffness and mineralization measurements were conducted at points across the intracortical thickness along radial lines emanating from the inner (endosteal) surface to the outer (periosteal) surface. Measurements were taken along approximately 4 mm distance of cortical bone thickness.

MATERIALS AND METHODS

Three experimental techniques were employed: Vickers microhardness (HV), energy-dispersive X-ray (EDX) spectroscopy, and computed tomography (CT). Stiffness values were extracted from the Vickers microhardness tests. Elemental mineralization values (calcium %wt. and phosphorus %wt.) were determined from EDX data. All measurements were repeated on three different femur bones taken from different bovines (collected fresh from butcher).

RESULTS

The study plots stiffness values and elemental mineralization (calcium %wt. and phosphorus %wt.) versus cortical thickness. Both stiffness and Ca %wt. and P %wt. are found to track and to linearly increase when plotted along the radial distance. The stiffness and mineralization trends collected from Vickers and EDX measurements were verified by employing the CT number (Hounsfield units, HU) via CT scans of the same bone samples. Data fitting via statistical methods revealed that all correlations were statistically significant.

CONCLUSION

Starting from endosteal to periosteal surfaces of mid-diaphyseal bovine femur, it was found that stiffness, mineralization, and HU values all exhibit increasing and correlating trends.

Compressive Strength of Iliac Bone ECM Is Not Reduced in Osteogenesis Imperfecta and Increases With Mineralization

Osteogenesis imperfecta (OI) is an inheritable, genetic, and collagen-related disorder leading to an increase in bone fragility, but the origin of its "brittle behavior" is unclear. Because of its complex hierarchical structure, bone behaves differently at various length scales. This study aims to compare mechanical properties of human OI bone with healthy control bone at the extracellular matrix (ECM) level and to quantify the influence of the degree of mineralization. Degree of mineralization and mechanical properties were analyzed under dry conditions in 12 fixed and embedded transiliac crest biopsies (control n = 6, OI type I n = 3, OI type IV n = 2, and OI type III n = 1). Mean degree of mineralization was measured by microcomputed tomography at the biopsy level and the mineral-to-matrix ratio was assessed by Raman spectroscopy at the ECM level. Both methods revealed that the degree of mineralization is higher for OI bone compared with healthy control. Micropillar compression is a novel technique for quantifying post-yield properties of bone at the ECM level. Micropillars (d = 5 μm, h = 10 μm) were fabricated using focused ion beam milling and quasi-statically compressed to capture key post-yield properties such as ultimate strength. The qualitative inspection of the stress-strain curves showed that both OI and healthy control bone have a ductile response at the ECM level. The quantitative results showed that compressive strength is not reduced in OI bone and is increasing with OI severity. Nanoindentation measurements revealed that OI bone tends to have a higher Young's modulus, hardness, and dissipated energy compared with healthy bone. Micropillar strength and indentation modulus increased linearly and significantly (p < .0001) with mineral-to-matrix ratio. In conclusion, this study indicates that compressive mechanical properties of dry OI bone at the iliac crest are not inferior to healthy control at the ECM level and increase with mineralization. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Maximum effect of the heterogeneity of tissue mineralization on the effective cortical bone elastic properties

The mineralization level is heterogeneous in cortical bone extracellular matrix as a consequence of remodeling. Models of the effective elastic properties at the millimeter scale have been developed based on idealizations of the vascular pore network and matrix properties. Some popular models do not take into account the heterogeneity of the matrix. However, the errors on the predicted elasticity when the difference in elastic properties between osteonal and interstitial tissues is not modeled have not been quantified. This work provides an estimation of the maximum error. We compare the effective elasticity of a representative volume element (RVE) assuming (1) different elastic properties in osteonal and interstitial tissues vs. (2) average matrix properties. In order to account for the variability of bone microstructure, we use a collection of high resolution images of the pore network to build RVEs. In each RVE we assumed a constant osteonal wall thickness and we artificially varied this thickness between 35 and 140 [Formula: see text]m to create RVEs with different amounts of osteonal tissue. The homogenization problem was solved with a fast Fourier transform (FFT)-based numerical scheme. We found that the error depends on pore volume fraction and varies on average from 1 to [Formula: see text] depending on the assumed diameter of the osteons. The results suggest that matrix heterogeneity may be disregarded in cortical bone models in most practical cases.

Bulk Wave Velocities in Cortical Bone Reflect Porosity and Compression Strength

The goal of this study was to evaluate whether ultrasonic velocities in cortical bone can be considered as a proxy for mechanical quality of cortical bone tissue reflected by porosity and compression strength. Micro-computed tomography, compression mechanical testing and resonant ultrasound spectroscopy were used to assess, respectively, porosity, strength and velocity of bulk waves of both shear and longitudinal polarisations propagating along and perpendicular to osteons, in 92 cortical bone specimens from tibia and femur of elderly human donors. All velocities were significantly associated with strength (r = 0.65-0.83) and porosity (r = -0.64 to -0.77). Roughly, according to linear regression models, a decrease in velocity of 100 m/s corresponded to a loss of 20 MPa in strength (which is approximately 10% of the largest strength value) and to an increase in porosity of 5%. These results provide a rationale for the in vivo measurement of one or several velocities for the diagnosis of bone fragility.

Cortical bone viscoelastic damping assessed with resonant ultrasound spectroscopy reflects porosity and mineral content

Viscoelasticity is an essential property of bone related to fragility, which is altered in aging and bone disease. Bone viscoelastic behavior is attributed to several mechanisms involving collagen and mineral properties, porosities, and bone hierarchical tissue organization. We aimed to assess the relationships between cortical bone viscoelastic damping measured with Resonant Ultrasound Spectroscopy (RUS), microstructural and compositional characteristics. We measured 52 bone specimens from the femur of 26 elderly human donors. RUS provided a shear damping coefficient at a frequency of the order of 150 kHz. The characteristics of the structure of the vascular pore network and tissue mineral density were measured using synchrotron radiation high-resolution computed tomography (SR-μCT). Fourier transformed infrared microspectroscopy (FTIRM) was used to quantify mineral-to-organic phase ratio, mineral maturity, crystallinity, and collagen maturity. Cross-links were quantified from biochemistry. Viscoelastic damping was found to increase with vascular porosity (r=0.68), to decrease with the degree of mineralization of the extravascular matrix (r=-0.68), and was marginally affected by collagen. We built a multilinear model suggesting that when porosity is controlled, the variation of mineral content explains a small additional part of the variability of damping. The work supports the consideration of viscoelasticity measurement as a potential biomarker of fragility and provides a documentation of bone viscoelastic behavior and its determinants in a frequency range rarely investigated.

Intrinsic properties of osteomalacia bone evaluated by nanoindentation and FTIRM analysis

Osteomalacia is a pathological bone condition consisting in a deficient primary mineralization of the matrix, leading to an accumulation of osteoid tissue and reduced bone mechanical strength. The amounts, properties and organization of bone constituents at tissue level, are known to influence its mechanical properties. It is then important to investigate the relationship between mechanical behavior and tissue composition at this scale in order to provide a better understanding of bone fragility mechanisms associates with this pathology. Our purpose was to analyze the links between ultra-structural properties and the mechanical behavior of this pathological bone tissue (osteomalacia) at tissue level (mineral and osteoid separately, or global). Four bone biopsies were taken from patients with osteomalacia, and subsequently embedded, sectioned, and polished. Then nanoindentation tests were performed to determine local elastic modulus E, contact hardness Hc and true hardness H for both mineralized and organic bone phases and for the global bone. The creep of the bone was also studied using a special indentation procedure in order to assess visco-elasto-plastic (creep) bone behavior. This allowed a detailed study of the rheological models adapted to the bone and to calculate the parameters associated to a Burgers model. Ultra-structural parameters were measured by Fourier Transform InfraRed Microspectroscopy (FTIRM) on the same position as the indents. The use of rheological models confirmed a significant contribution from the organic phase on the viscous character of bone tissue. The elastic E and the elasto-plastic H_c deformation were correlated to both collagen maturity and Mineral/Matrix. The pure plastic deformation H was only correlated to the mineral phase. Our data show that mineral phase greatly affects mechanical variables (moduli and viscosities) and that organic phase (as illustrated in osteoid tissue) may play an important role in the creep behavior of bone. In conclusion, this study brings mechanical and physicochemical values for osteoid and mineral phases.

Micromechanical evaluation of cortical bone using in situ XCT indentation and digital volume correlation

The overall mechanical behaviour of cortical bone is strongly dependant on its microstructure. X-ray computed tomography (XCT) has been widely used to identify the microstructural morphology of cortical tissue (i.e. pore network, Haversian and Volkmann's canals). However, the connection between microstructure and mechanics of cortical bone during plastic deformation is unclear. Hence, the purpose of this study is to provide an in-depth evaluation of the interplay of plastic strain building up in relation to changes in the canal network for cortical bone tissue. In situ step-wise XCT indentation was used to introduce a localised load on the surface of the tissue and digital volume correlation (DVC) was employed to assess the three-dimensional (3D) full-field plastic strain distribution in proximity of the indent. It was observed that regions adjacent to the imprint were under tensile strain, whereas the volume underneath experienced compressive strain. Canal loss and disruption was detected in regions of higher compressive strains exceeding -20000 με and crack formation occurred in specimens where Haversian canals were running parallel to the indentation tip. The results of this study outline the relationship between the micromechanical and structural behaviour of cortical bone during plastic deformation, providing information on cortical tissue fracture pathways.

Assessing the effective elastic properties of the tendon-to-bone insertion: a multiscale modeling approach

The interphase joining tendon to bone plays the crucial role of integrating soft to hard tissues, by effectively transferring stresses across two tissues displaying a mismatch in mechanical properties of nearly two orders of magnitude. The outstanding mechanical properties of this interphase are attributed to its complex hierarchical structure, especially by means of competing gradients in mineral content and collagen fibers organization at different length scales. The goal of this study is to develop a multiscale model to describe how the tendon-to-bone insertion derives its overall mechanical behavior. To this end, the effective anisotropic stiffness tensor of the interphase is predicted by modeling its elastic response at different scales, spanning from the nanostructural to the mesostructural levels, using continuum micromechanics methods. The results obtained at a lower scale serve as inputs for the modeling at a higher scale. The obtained predictions are in good agreement with stochastic finite element simulations and experimental trends reported in literature. Such model has implication for the design of bioinspired bi-materials that display the functionally graded properties of the tendon-to-bone insertion.