Human cortical bone: the SiNuPrOs model

Many models that have been developed for cortical bone oversimplify much of the architectural and physical complexity. With SiNuPrOs model, a more complete approach is investigated: it is multiscale because it contains five structural levels and multi physic because it takes into account simultaneously structure (with various properties: elasticity, piezoelectricity, porous medium), fluid and mineralization process modelization. The multiscale aspect is modeled by using 18 structural parameters in a specific application of the mathematical theory of homogenization and 10 other physical parameters are necessary for the multi physic aspect. The modelization of collagen as a piezoelectric medium has needed the development of a new behaviour law allowing a better simulation of the effect of a medium considered as evolving during a mineralization process.

Then the main interest of SiNuPrOs deals with the possibility to study, at each level of the cortical architecture, either the elastic properties or the fluid motion or the piezoelectric effects or both of them. All these possibilities constitute a very large work and all this mass of information (fluid aspects, even at the nanoscopic scale, piezoelectric phenomena and simulations) will be presented in several papers. This first one is only devoted to the presentation of this model with an application to the computation of elastic properties at the macroscopic scale.

The computational methods have been packed into software also called SiNuPrOs and allowing a large number of predictive simulations corresponding to various different configurations.     

Collagen's role in the cortical bone's behaviour: a numerical approach

Literature devoted to experimental measurements of the elastic properties of the human cortical bone gives us a relatively wide spectrum of values. This proves that the result depends on the bone itself and on the area from where the sample has been taken. Ultrasonic measurement, which is a fine technology, points out complex maps. The reason of this very strong heterogeneity is not completely explained. The present study is based on a numerical model of the human cortical bone, the SiNuPrOs model and aims to suggest an explanation. If one admits that mineral apposition occurs around collagen fibres, the spatial orientation of these fibres would have an important consequence on the elastic properties of the medium. On the basis of homogenisation theory allowing to compute all the components of the elasticity tensor, this study quantifies the main influence of this architectural orientation and its effect on the anisotropy of the cortical bone.     

Collagen's role in the cortical bone's behaviour: a numerical approach

Literature devoted to experimental measurements of the elastic properties of the human cortical bone gives us a relatively wide spectrum of values. This proves that the result depends on the bone itself and on the area from where the sample has been taken. Ultrasonic measurement, which is a fine technology, points out complex maps. The reason of this very strong heterogeneity is not completely explained. The present study is based on a numerical model of the human cortical bone, the SiNuPrOs model and aims to suggest an explanation. If one admits that mineral apposition occurs around collagen fibres, the spatial orientation of these fibres would have an important consequence on the elastic properties of the medium. On the basis of homogenisation theory allowing to compute all the components of the elasticity tensor, this study quantifies the main influence of this architectural orientation and its effect on the anisotropy of the cortical bone.     

Determination of the heterogeneous anisotropic elastic properties of human femoral bone: From nanoscopic to organ scale

Cortical bone is a multiscale composite material. Its elastic properties are anisotropic and heterogeneous across its cross-section, due to endosteal bone resorption which might affect bone strength. The aim of this paper was to describe a homogenization method leading to the estimation of the variation of the elastic coefficients across the bone cross-section and along the bone longitudinal axis. The method uses the spatial variations of bone porosity and of the degree of mineralization of the bone matrix (DMB) obtained from the analysis of 3-D synchrotron micro-computed tomography images. For all three scales considered (the foam (100 nm), the ultrastructure (5 μm) and the mesoscale (500 μm)), the elastic coefficients were determined using the Eshelby’s inclusion problem. DMB values were used at the scale of the foam. Collagen was introduced at the scale of the ultrastructure and bone porosity was introduced at the mesoscale. The pores were considered as parallel cylinders oriented along the bone axis. Each elastic coefficient was computed for different regions of interest, allowing an estimation of its variations across the bone cross-section and along the bone longitudinal axis. The method was applied to a human femoral neck bone specimen, which is a site of osteoporotic fracture. The computed elastic coefficients for cortical bone were in good agreement with experimental results, but some discrepancies were obtained in the endosteal part (trabecular bone). These results highlight the importance of accounting for the heterogeneity of cortical bone properties across bone cross-section and along bone longitudinal axis.     

A new numerical concept for modeling hydroxyapatite in human cortical bone

This research presents a new modelling procedure which allows the computation of the physical properties of the human cortical bone, considered as a strongly heterogeneous medium consisting of bony architecture and the physical properties of the two basic components: the collagen and the hydroxyapatite (Hap). The numerical simulations are based on the homogenisation theory, however, since the size of the Hap crystals are small compared to the size of a collagen stick, a new entity (the elementary volume of mineral content (EVMC)) is defined at the nanoscopic scale. This model permits the testing of all the possible structural configurations that may be present and suggests that the anisotropy of the bone is not only induced by the haversian structure but by the properties of the Hap crystals and their special organisation.     

A new numerical concept for modeling hydroxyapatite in human cortical bone

This research presents a new modelling procedure which allows the computation of the physical properties of the human cortical bone, considered as a strongly heterogeneous medium consisting of bony architecture and the physical properties of the two basic components: the collagen and the hydroxyapatite (Hap). The numerical simulations are based on the homogenisation theory, however, since the size of the Hap crystals are small compared to the size of a collagen stick, a new entity (the elementary volume of mineral content (EVMC)) is defined at the nanoscopic scale. This model permits the testing of all the possible structural configurations that may be present and suggests that the anisotropy of the bone is not only induced by the haversian structure but by the properties of the Hap crystals and their special organisation.     

Auto-adaptative network separator of sources]

This article defines a new convergence criterion for the modelization of networks with modifiable synapses, as presented by Jutten, Herault, Ans ([1], [2]). The network adaptation law, empirically determined in [1], [2], is reconsidered: a theoretical law concerning the evolution of system coefficients and outputs is demonstrated and validated on simulations. The system always converges towards a unique attractor, as predicted by the theoretical convergence criterion. This system evolution towards a particular stable state enhances the interesting physical aspect of this modelization.     

'Universal' microstructural patterns in cortical and trabecular, extracellular and extravascular bone materials: micromechanics-based prediction of anisotropic elasticity

Bone materials are characterized by an astonishing variability and diversity. Still, because of 'architectural constraints' due to once chosen material constituents and their physical interaction, the fundamental hierarchical organization or basic building plans of bone materials remain largely unchanged during biological evolution. Such universal patterns of microstructural organization govern the mechanical interaction of the elementary components of bone (hydroxyapatite, collagen, water; with directly measurable tissue-independent elastic properties), which are here quantified through a multiscale homogenization scheme delivering effective elastic properties of bone materials: at a scale of 10nm, long cylindrical collagen molecules, attached to each other at their ends by approximately 1.5nm long crosslinks and hosting intermolecular water inbetween, form a contiguous matrix called wet collagen. At a scale of several hundred nanometers, wet collagen and mineral crystal agglomerations interpenetrate each other, forming the mineralized fibril. At a scale of 5-10microm, the extracellular solid bone matrix is represented as collagen fibril inclusions embedded in a foam of largely disordered (extrafibrillar) mineral crystals. At a scale above the ultrastructure, where lacunae are embedded in extracellular bone matrix, the extravascular bone material is observed. Model estimates predicted from tissue-specific composition data gained from a multitude of chemical and physical tests agree remarkably well with corresponding acoustic stiffness experiments across a variety of cortical and trabecular, extracellular and extravascular materials. Besides from reconciling the well-documented, seemingly opposed concepts of 'mineral-reinforced collagen matrix' and 'collagen-reinforced mineral matrix' for bone ultrastructure, this approach opens new possibilities in the exploitation of computer tomographic data for nano-to-macro mechanics of bone organs.     

Nonlinear hierarchical multiscale modeling of cortical bone considering its nanoscale microstructure

We have used a hierarchical multiscale modeling scheme for the analysis of cortical bone considering it as a nanocomposite. This scheme consists of definition of two boundary value problems, one for macroscale, and another for microscale. The coupling between these scales is done by using the homogenization technique. At every material point in which the constitutive model is needed, a microscale boundary value problem is defined using a macroscopic kinematical quantity and solved. Using the described scheme, we have studied elastic properties of cortical bone considering its nanoscale microstructural constituents with various mineral volume fractions. Since the microstructure of bone consists of mineral platelet with nanometer size embedded in a protein matrix, it is similar to the microstructure of soft matrix nanocomposites reinforced with hard nanostructures. Considering a representative volume element (RVE) of the microstructure of bone as the microscale problem in our hierarchical multiscale modeling scheme, the global behavior of bone is obtained under various macroscopic loading conditions. This scheme may be suitable for modeling arbitrary bone geometries subjected to a variety of loading conditions. Using the presented method, mechanical properties of cortical bone including elastic moduli and Poisson's ratios in two major directions and shear modulus is obtained for different mineral volume fractions.     

Hyperlipidemia affects multiscale structure and strength of murine femur

To improve bone strength prediction beyond limitations of assessment founded solely on the bone mineral component, we investigated the effect of hyperlipidemia, present in more than 40% of osteoporotic patients, on multiscale structure of murine bone. Our overarching purpose is to estimate bone strength accurately, to facilitate mitigating fracture morbidity and mortality in patients. Because (i) orientation of collagen type I affects, independently of degree of mineralization, cortical bone?s micro-structural strength; and, (ii) hyperlipidemia affects collagen orientation and μCT volumetric tissue mineral density (vTMD) in murine cortical bone, we have constructed the first multiscale finite element (mFE), mouse-specific femoral model to study the effect of collagen orientation and vTMD on strength in Ldlr^−/−, a mouse model of hyperlipidemia, and its control wild type, on either high fat diet or normal diet. Each µCT scan-based mFE model included either element-specific elastic orthotropic properties calculated from collagen orientation and vTMD (collagen-density model) by experimentally validated formulation, or usual element-specific elastic isotropic material properties dependent on vTMD-only (density-only model). We found that collagen orientation, assessed by circularly polarized light and confocal microscopies, and vTMD, differed among groups and that microindentation results strongly correlate with elastic modulus of collagen-density models (r²=0.85, p=10⁻⁵). Collagen-density models yielded (1) larger strains, and therefore lower strength, in simulations of 3-point bending and physiological loading; and (2) higher correlation between mFE-predicted strength and 3-point bending experimental strength, than density-only models. This novel method supports ongoing translational research to achieve the as yet elusive goal of accurate bone strength prediction.     

Apparent Young's modulus of human radius using inverse finite-element method

The ability to assess the elastic and failure properties of cortical bone at the radial diaphysis has a clinical importance. A new generation of quantitative ultrasound (QUS) devices and peripheral quantitative computed tomography (p-QCT) has been developed to assess non-invasively bone material and structural properties at the distal radius. This anatomical site is characterized by a thin cortical thickness that complicates traditional mechanical testing methods on specimens. Until now, mechanical properties of cortical bone at distal radius (e.g., elastic modulus, yield stress and strain) remain rarely studied probably due to experimental difficulties. The present study introduces an inverse finite-element method strategy to measure the elastic modulus and yield properties of human cortical specimens of the radial diaphysis. Twenty millimeter-thick portions of diaphysis were cut from 40 human radii (ages 45-90) for biomechanical test. Subsequently the same portion was modeled in order to obtain a specimen-specific three dimensional finite-element model (3D-FEM). Longitudinal elastic constants at the apparent level and stress characterizations were performed by coupling mechanical parameters with isotropic linear-elastic simulations. The results indicated that the mean apparent Young's modulus for radial cortical bone was 16 GPa (SD 1.8) and the yield stress was 153 MPa (SD 33). Breaking load was 12,946 N (SD 3644), cortical thickness 2.9 mm (SD 0.6), structural effective strain at the yield (epsilon(y)=0.0097) and failure (epsilon(u)=0.0154) load were also calculated. The 3D-FEM strategy described here may help to investigate bone mechanical properties when some difficulties arise from machining mechanical sample.     

Influence of physical activity on tibial bone material properties in laying hens

Laying hens develop a type of osteoporosis that arises from a loss of structural bone, resulting in high incidence of fractures. In this study, a comparison of bone material properties was made for lines of hens created by divergent selection to have high and low bone strength and housed in either individual cages, with restricted mobility, or in an aviary system, with opportunity for increased mobility. Improvement of bone biomechanics in the high line hens and in aviary housing was mainly due to increased bone mass, thicker cortical bone and more medullary bone. However, bone material properties such as cortical and medullary bone mineral composition and crystallinity as well as collagen maturity did not differ between lines. However, bone material properties of birds from the different type of housing were markedly different. The cortical bone in aviary birds had a lower degree of mineralization and bone mineral was less mature and less organized than in caged birds. These differences can be explained by increased bone turnover rates due to the higher physical activity of aviary birds that stimulates bone formation and bone remodeling. Multivariate statistical analyses shows that both cortical and medullary bone contribute to breaking strengthThe cortical thickness was the single most important contributor while its degree of mineralization and porosity had a smaller contribution. Bone properties had poorer correlations with mechanical properties in cage birds than in aviary birds presumably due to the greater number of structural defects of cortical bone in cage birds.     

A Multiscale Theoretical Investigation of Electric Measurements in Living Bone

This paper presents a theoretical investigation of the multiphysical phenomena that govern cortical bone behaviour. Taking into account the piezoelectricity of the collagen–apatite matrix and the electrokinetics governing the interstitial fluid movement, we adopt a multiscale approach to derive a coupled poroelastic model of cortical tissue. Following how the phenomena propagate from the microscale to the tissue scale, we are able to determine the nature of macroscopically observed electric phenomena in bone.     

An approach to multiscale modelling with graph grammars

Background and Aims

Functional–structural plant models (FSPMs) simulate biological processes at different spatial scales. Methods exist for multiscale data representation and modification, but the advantages of using multiple scales in the dynamic aspects of FSPMs remain unclear. Results from multiscale models in various other areas of science that share fundamental modelling issues with FSPMs suggest that potential advantages do exist, and this study therefore aims to introduce an approach to multiscale modelling in FSPMs.

Methods

A three-part graph data structure and grammar is revisited, and presented with a conceptual framework for multiscale modelling. The framework is used for identifying roles, categorizing and describing scale-to-scale interactions, thus allowing alternative approaches to model development as opposed to correlation-based modelling at a single scale. Reverse information flow (from macro- to micro-scale) is catered for in the framework. The methods are implemented within the programming language XL.

Key Results

Three example models are implemented using the proposed multiscale graph model and framework. The first illustrates the fundamental usage of the graph data structure and grammar, the second uses probabilistic modelling for organs at the fine scale in order to derive crown growth, and the third combines multiscale plant topology with ozone trends and metabolic network simulations in order to model juvenile beech stands under exposure to a toxic trace gas.

Conclusions

The graph data structure supports data representation and grammar operations at multiple scales. The results demonstrate that multiscale modelling is a viable method in FSPM and an alternative to correlation-based modelling. Advantages and disadvantages of multiscale modelling are illustrated by comparisons with single-scale implementations, leading to motivations for further research in sensitivity analysis and run-time efficiency for these models.     

Derivation of the mesoscopic elasticity tensor of cortical bone from quantitative impedance images at the micron scale

This paper addresses the relationships between the microscopic properties of bone and its elasticity at the millimetre scale, or mesoscale. A method is proposed to estimate the mesoscale properties of cortical bone based on a spatial distribution of acoustic properties at the microscopic scale obtained with scanning acoustic microscopy. The procedure to compute the mesoscopic stiffness tensor involves (i) the segmentation of the pores to obtain a realistic model of the porosity; (ii) the construction of a field of anisotropic elastic coefficients at the microscopic scale which reflects the heterogeneity of the bone matrix; (iii) finite element computations of mesoscopic homogenized properties. The computed mesoscopic properties compare well with available experimental data. It appears that the tissue anisotropy at the microscopic level has a major effect on the mesoscopic anisotropy and that assuming the pores filled with an incompressible fluid or, alternatively, empty, leads to significantly different mesoscopic properties.     

Analytical methods to determine the effective mesoscopic and macroscopic elastic properties of cortical bone

We compare theoretical predictions of the effective elastic moduli of cortical bone at both the meso- and macroscales. We consider the efficacy of three alternative approaches: the method of asymptotic homogenization, the Mori-Tanaka scheme and the Hashin-Rosen bounds. The methods concur for specific engineering moduli such as the axial Young's modulus but can vary for others. In a past study, the effect of porosity alone on mesoscopic properties of cortical bone was considered, taking the matrix to be isotropic. Here, we consider the additional influence of the transverse isotropy of the matrix. We make the point that micromechanical approaches can be used in two alternative ways to predict either the macroscopic (size of cortical bone sample) or mesoscopic (in between micro- and macroscales) effective moduli, depending upon the choice of representative volume element size. It is widely accepted that the mesoscale behaviour is an important aspect of the mechanical behaviour of bone but models incorporating its effect have started to appear only relatively recently. Before this only macroscopic behaviour was addressed. Comparisons are drawn with experimental data and simulations from the literature for macroscale predictions with particularly good agreement in the case of dry bone. Finally, we show how predictions of the effective mesoscopic elastic moduli can be made which retain dependence on the well-known porosity gradient across the thickness of cortical bone.     

Spatial distribution of tissue level properties in a human femoral cortical bone

The mechanical properties of cortical bone are determined by a combination bone tissue composition, and structure at several hierarchical length scales. In this study the spatial distribution of tissue level properties within a human femoral shaft has been investigated. Cylindrically shaped samples (diameter: 4.4mm, N=56) were prepared from cortical regions along the entire length (20-85% of the total femur length), and around the periphery (anterior, medial, posterior and lateral quadrants). The samples were analyzed using scanning acoustic microscopy (SAM) at 50MHz and synchrotron radiation micro computed tomography (SRμCT). For all samples the average cortical porosity (Ct.Po), tissue elastic coefficients (c(ij)) and the average tissue degree of mineralization (DMB) were determined. The smallest coefficient of variation was observed for DMB (1.8%), followed by BV/TV (5.4%), c(ij) (8.2-45.5%), and Ct.Po (47.5%). Different variations with respect to the anatomical position were found for DMB, Ct.Po and c(ij). These data address the anatomical variations in anisotropic elastic properties and link them to tissue mineralization and porosity, which are important input parameters for numerical multi-scale bone models.     

Impact of the porous microstructure on the overall elastic properties of the osteonal cortical bone

Mechanical properties of bones are largely determined by their microstructure. The latter comprises a large number of diverse pores. The present paper analyzes a connection between structure of the porous space of the osteonal cortical bone and bone's overall anisotropic elastic moduli. The analysis is based on recent developments in the theory of porous materials that predict the anisotropic effective moduli of porous solids in terms of pores' shapes, orientations and densities. Bone's microstructure is modeled using available micrographs. The calculated anisotropic elastic constants for porous cortical bone are, mostly, in agreement with available experimental data. The influence of each of the pore types on the overall moduli is examined. The results of the analysis can also be used to estimate the extent of mineralization (hydroxyapatite content) if the overall porosity and the effective moduli are known and, vice versa, to estimate porosity from the measured moduli and the extent of mineralization.     

Gating mechanisms of mechanosensitive channels of large conductance, I: a continuum mechanics-based hierarchical framework

A hierarchical simulation framework that integrates information from molecular dynamics (MD) simulations into a continuum model is established to study the mechanical response of mechanosensitive channel of large-conductance (MscL) using the finite element method (FEM). The proposed MD-decorated FEM (MDeFEM) approach is used to explore the detailed gating mechanisms of the MscL in Escherichia coli embedded in a palmitoyloleoylphosphatidylethanolamine lipid bilayer. In Part I of this study, the framework of MDeFEM is established. The transmembrane and cytoplasmic helices are taken to be elastic rods, the loops are modeled as springs, and the lipid bilayer is approximated by a three-layer sheet. The mechanical properties of the continuum components, as well as their interactions, are derived from molecular simulations based on atomic force fields. In addition, analytical closed-form continuum model and elastic network model are established to complement the MDeFEM approach and to capture the most essential features of gating. In Part II of this study, the detailed gating mechanisms of E. coli-MscL under various types of loading are presented and compared with experiments, structural model, and all-atom simulations, as well as the analytical models established in Part I. It is envisioned that such a hierarchical multiscale framework will find great value in the study of a variety of biological processes involving complex mechanical deformations such as muscle contraction and mechanotransduction.