Stefan Scheiner

Micromechanics of bone tissue-engineering scaffolds, based on resolution error-cleared computer tomography.

Synchrotron radiation micro-computed tomography (SRmuCT) revealed the microstructure of a CEL2 glass-ceramic scaffold with macropores of several hundred microns characteristic length, in terms of the voxel-by-voxel 3D distribution of the attenuation coefficients throughout the scanned space. The probability density function of all attenuation coefficients related to the macroporous space inside the scaffold gives access to the tomograph-specific machine error included in the SRmuCT measurements (also referred to as instrumental resolution function). After Lorentz function-based clearing of the measured CT data from the systematic resolution error, the voxel-specific attenuation information of the voxels representing the solid skeleton is translated into the composition of the material inside one voxel, in terms of the nanoporosity embedded in a dense CEL2 glass-ceramic matrix. Based on voxel-invariant elastic properties of dense CEL2 glass-ceramic, continuum micromechanics allows for translation of the voxel-specific nanoporosity into voxel-specific elastic properties. They serve as input for Finite Element analyses of the scaffold structure. Youngs modulus of a specific CT-scanned macroporous scaffold sample, predicted from a Finite Element simulation of a uniaxial compression test, agrees well with the experimental value obtained from an ultrasonic test on the same sample. This highlights the satisfactory predictive capabilities of the presented approach.

Layered water in crystal interfaces as source for bone viscoelasticity: arguments from a multiscale approach

Extracellular bone material can be characterised as a nanocomposite where, in a liquid environment, nanometre-sized hydroxyapatite crystals precipitate within as well as between long fibre-like collagen fibrils (with diameters in the 100 nm range), as evidenced from neutron diffraction and transmission electron microscopy. Accordingly, these crystals are referred to as ‘interfibrillar mineral’ and ‘extrafibrillar mineral’, respectively. From a topological viewpoint, it is probable that the mineralisations start on the surfaces of the collagen fibrils (‘mineral-encrusted fibrils’), from where the crystals grow both into the fibril and into the extrafibrillar space. Since the mineral concentration depends on the pore spaces within the fibrils and between the fibrils (there is more space between them), the majority of the crystals (but clearly not all of them) typically lie in the extrafibrillar space. There, larger crystal agglomerations or clusters, spanning tens to hundreds of nanometers, develop in the course of mineralisation, and the micromechanics community has identified the pivotal role, which this extrafibrillar mineral plays for tissue elasticity. In such extrafibrillar crystal agglomerates, single crystals are stuck together, their surfaces being covered with very thin water layers. Recently, the latter have caught our interest regarding strength properties (Fritsch et al. 2009 J Theor Biol. 260(2): 230–252) – we have identified these water layers as weak interfaces in the extrafibrillar mineral of bone. Rate-independent gliding effects of crystals along the aforementioned interfaces, once an elastic threshold is surpassed, can be related to overall elastoplastic material behaviour of the hierarchical material ‘bone’. Extending this idea, the present paper is devoted to viscous gliding along these interfaces, expressing itself, at the macroscale, in the well-known experimentally evidenced phenomenon of bone viscoelasticity. In this context, a multiscale homogenisation scheme is extended to viscoelasticity, mineral-cluster-specific creep parameters are identified from three-point bending tests on hydrated bone samples, and the model is validated by statistically and physically independent experiments on partially dried samples. We expect this model to be relevant when it comes to prediction of time-dependent phenomena, e.g. in the context of bone remodelling.

The influence of bone surface availability in bone remodelling—A mathematical model including coupled geometrical and biomechanical regulations of bone cells

Bone is a biomaterial undergoing continuous renewal. The renewal process is known as bone remodelling and is operated by bone-resorbing cells (osteoclasts) and bone-forming cells (osteoblasts). An important function of bone remodelling is the repair of microcracks accumulating in the bone matrix due to mechanical loading. Cell–cell communication between cells of the osteoclastic lineage and cells of the osteoblastic lineage is thought to couple resorption and formation so as to preserve bone integrity and achieve homeostatic bone renewal. Both biochemical and biomechanical regulatory mechanisms have been identified in this coupling. Many bone pathologies are associated with an alteration of bone cell interactions and a consequent disruption of bone homeostasis. In osteoporosis, for example, this disruption leads to long-term bone loss and increased fragility, and can ultimately result in fractures. Here we focus on an additional and poorly understood potential regulatory mechanism of bone cells, that involves the morphology of the microstructure of bone. Bone cells can only remove and replace bone at a bone surface. However, the microscopic availability of bone surface depends in turn on the ever-changing bone microstructure. The importance of this geometrical dependence is unknown and difficult to quantify experimentally. Therefore, we develop a sophisticated mathematical model of bone cell interactions that takes into account biochemical, biomechanical and geometrical regulations. We then investigate numerically the influence of bone surface availability in bone remodelling within a representative bone tissue sample. Biochemical regulations included in the model involve signalling molecules of the receptor–activator nuclear factor κB pathway (rank–rankl–opg), macrophage colony-stimulating factor (mcsf), transforming growth factor β(tgfβ), and parathyroid hormone (pth). For the biomechanical regulation of bone cells, a multiscale homogenisation scheme is used to determine the microscopic strains generated at the level of the extravascular matrix hosting the osteocytes by macroscopic loading. The interdependence between the bone cells’ activity, which modifies the bone microstructure, and changes in the microscopic bone surface availability, which in turn influences bone cell development and activity, is implemented using a remarkable experimental relationship between bone specific surface and bone porosity. Our model suggests that geometrical regulation of the activation of new remodelling events could have a significant effect on bone porosity and bone stiffness in osteoporosis. On the other hand, geometrical regulation of late stages of osteoblast and osteoclast differentiation seems less significant. We conclude that the development of osteoporosis is probably accelerated by this geometrical regulation in cortical bone, but probably slowed down in trabecular bone.

A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse

We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.

Poromicromechanics reveals that physiological bone strains induce osteocyte-stimulating lacunar pressure

Mechanical loads which are macroscopically acting onto bony organs, are known to influence the activities of biological cells located in the pore spaces of bone, in particular so the signaling and production processes mediated by osteocytes. The exact mechanisms by which osteocytes are actually able to “feel” the mechanical loading and changes thereof, has been the subject of numerous studies, and, while several hypotheses have been brought forth over time, this topic has remained a matter of debate. Relaxation times reported in a recent experimental study of Gardinier et al. (Bone 46(4):1075–1081, 2010) strongly suggest that the lacunar pores are likely to experience, during typical physiological load cycles, not only fluid transport, but also undrained conditions. The latter entail the buildup of lacunar pore pressures, which we here quantify by means of a thorough multiscale modeling approach. In particular, the proposed model is based on classical poroelasticity theory, and able to account for multiple pore spaces. First, the model reveals distinct nonlinear dependencies of the resulting lacunar (and vascular) pore pressures on the underlying bone composition, highlighting the importance of a rigorous multiscale approach for appropriate computation of the aforementioned pore pressures. Then, the derived equations are evaluated for macroscopic (uniaxial as well as hydrostatic) mechanical loading of physiological magnitude. The resulting model-predicted pore pressures agree very well with the pressures that have been revealed, by means of in vitro studies, to be of adequate magnitude for modulating the responses of biological cells, including osteocytes. This underlines that osteocytes may respond to many types of loading stimuli at the same time, in particular so to fluid flow and hydrostatic pressure.

Mathematical modeling of postmenopausal osteoporosis and its treatment by the anti-catabolic drug denosumab

Denosumab, a fully human monoclonal antibody, has been approved for the treatment of postmenopausal osteoporosis. The therapeutic effect of denosumab rests on its ability to inhibit osteoclast differentiation. Here, we present a computational approach on the basis of coupling a pharmacokinetics model of denosumab with a pharmacodynamics model for quantifying the effect of denosumab on bone remodeling. The pharmacodynamics model comprises an integrated systems biology-continuum micromechanics approach, including a bone cell population model, considering the governing biochemical factors of bone remodeling (including the action of denosumab), and a multiscale micromechanics-based bone mechanics model, for implementing the mechanobiology of bone remodeling in our model. Numerical studies of postmenopausal osteoporosis show that denosumab suppresses osteoclast differentiation, thus strongly curtailing bone resorption. Simulation results also suggest that denosumab may trigger a short-term bone volume gain, which is, however, followed by constant or decreasing bone volume. This evolution is accompanied by a dramatic decrease of the bone turnover rate by more than one order of magnitude. The latter proposes dominant occurrence of secondary mineralization (which is not anymore impeded through cellular activity), leading to higher mineral concentration per bone volume. This explains the overall higher bone mineral density observed in denosumab-related clinical studies. Copyright

Strength increase during ceramic biomaterial-induced bone regeneration: a micromechanical study

Bone tissue engineering materials must blend in the targeted physiological environment, in terms of both the materials’ biocompatibility and mechanical properties. As for the latter, a well-adjusted stiffness ensures that the biomaterial’s deformation behavior fits well to the deformation behavior of the surrounding biological tissue, whereas an appropriate strength provides sufficient load-carrying capacity of the biomaterial. Here, a mathematical modeling approach for estimating the macroscopic load that initiates failure of a hierarchically organized, granular, hydroxyapatite-based biomaterial is presented. For this purpose, a micromechanics model is developed for downscaling macroscopically prescribed stress (or strain) states to the level of the needle-shaped hydroxyapatite crystals. Presuming that the biomaterial fails due to the quasi-brittle failure of the most unfavorably stressed hydroxyapatite needle, the downscaled stress tensors are fed into a suitable, Mohr-Coulomb-type failure criterion, based on which the macroscopic failure load is deduced. The change of the biomaterial’s composition in response to placing it in physiological solution, caused by growth of new bone tissue on the granules’s surfaces, on the one hand, and by resorption of the hydroxyapatite crystals, on the other hand, is taken into account by means of suitable evolution laws. Numerical studies show how the macroscopic load-carrying capacity of the biomaterial is influenced by its design parameters. The presented modeling approach could prove beneficial for the design process of the studied biomaterials (as well as similarly composed biomaterials), particularly in terms of optimizing its mechanical performance.

Two-scale model for electro-diffusive transport through charged porous materials

The electro-diffusive transport behavior of ions through charged porous media is commonly described by phenomenological theories, directly formulated on the observation scale of the macroscopic material and disregarding the underlying physical principles known on the particle scale. In this paper, to remedy this deficit, a new approach is presented, by derivation of a generalized macroscopic mathematical framework valid for both charged and uncharged porous media, on the basis of classical Poisson and Nernst-Planck equations. Due to the distinctive nonlinearity of the sought-after quantities, an iterative, numerical two-scale homogenization scheme is employed for evaluation of the resulting governing equations, giving finally access to the effective macroscopic diffusion coefficients and fixed charge concentration as functions of the electrolyte background concentration and of the surface charge. While this methodology can be applied for a large number of different pore geometries, we restrict ourselves, for the time being, to linear cylindrical pores, to demonstrate the capability of the presented approach. For the chosen setup and boundary conditions, the applied two-scale model reveals deviations of the effective macroscopic diffusion coefficients from corresponding self-diffusion coefficients amounting to up to ≈ 12%. These results corroborate the relevance of the presented approach.

Self-Consistent Channel Approach for Upscaling Chloride Diffusivity in Cement Pastes

Chloride ingress into concrete is a major cause for material degradation, such as cracking due to corrosion-induced steel reinforcement expansion. Corresponding transport processes encompass diffusion, convection, and migration, and their mathematical quantification as a function of the concrete composition remains an unrevealed enigma. Approaching the problem step by step, we here concentrate on the diffusivity of cement paste, and how it follows from the microstructural features of the material and from the chloride diffusivity in the capillary pore spaces. For this purpose, we employ advanced self-consistent homogenization theory as recently used for permeability upscaling, based on the resolution of the pore space as pore channels being oriented in all space directions, resulting in a quite compact analytical relation between porosity, pore diffusivity, and the overall diffusivity of the cement paste. This relation is supported by experiments and reconfirms the pivotal role that layered water most probably plays for the reduction of the pore diffusivity, with respect to the diffusivity found under the chemical condition of a bulk solution.

Multiscale Mathematical Modeling in Dental Tissue Engineering: Toward Computer-Aided Design of a Regenerative System Based on Hydroxyapatite Granules, Focussing on Early and Mid-Term Stiffness Recovery

We here explore for the very first time how an advanced multiscale mathematical modeling approach may support the design of a provenly successful tissue engineering concept for mandibular bone. The latter employs double-porous, potentially cracked, single millimeter-sized granules packed into an overall conglomerate-type scaffold material, which is then gradually penetrated and partially replaced by newly grown bone tissue. During this process, the newly developing scaffold-bone compound needs to attain the stiffness of mandibular bone under normal physiological conditions. In this context, the question arises how the compound stiffness is driven by the key design parameters of the tissue engineering system: macroporosity, crack density, as well as scaffold resorption/bone formation rates. We here tackle this question by combining the latest state-of-the-art mathematical modeling techniques in the field of multiscale micromechanics, into an unprecedented suite of highly efficient, semi-analytically defined computation steps resolving several levels of hierarchical organization, from the millimeter- down to the nanometer-scale. This includes several types of homogenization schemes, namely such for porous polycrystals with elongated solid elements, for cracked matrix-inclusion composites, as well as for assemblies of coated spherical compounds. Together with the experimentally known stiffnesses of hydroxyapatite crystals and mandibular bone tissue, the new mathematical model suggests that early stiffness recovery (i.e., within several weeks) requires total avoidance of microcracks in the hydroxyapatite scaffolds, while mid-term stiffness recovery (i.e., within several months) is additionally promoted by provision of small granule sizes, in combination with high bone formation and low scaffold resorption rates.