Study of the combined effects of PTH treatment and mechanical loading in postmenopausal osteoporosis using a new mechanistic PK-PD model

Biomechanics and Modeling in Mechanobiology - One of only a few approved and available anabolic treatments for severe osteoporosis is daily injections of PTH (1-34). This drug has a specific dual...     

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.     

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.     

The effects of anatomical errors on shoulder kinematics computed using multi-body models

Joint motion calculated using multi-body models and inverse kinematics presents many advantages over direct marker-based calculations. However, the sensitivity of the computed kinematics is known to be partly caused by the model and could also be influenced by the participants’ anthropometry and sex. This study aimed to compare kinematics computed from an anatomical shoulder model based on medical images against a scaled-generic model and quantify the effects of anatomical errors and participants’ anthropometry on the calculated joint angles. Twelve participants have had planar shoulder movements experimentally captured in a motion lab, and their shoulder anatomy imaged using an MRI scanner. A shoulder multi-body dynamics model was developed for each participant, using both an image-based approach and a scaled-generic approach. Inverse kinematics have been performed using the two different modelling procedures and the three different experimental motions. Results have been compared using Bland–Altman analysis of agreement and further analysed using multi-linear regressions. Kinematics computed via an anatomical and a scaled-generic shoulder models differed in average from 3.2 to 5.4 degrees depending on the task. The MRI-based model presented smaller limits of agreement to direct kinematics than the scaled-generic model. Finally, the regression model predictors, including anatomical errors, sex, and BMI of the participant, explained from 41 to 80% of the kinematic variability between model types with respect to the task. This study highlighted the consequences of modelling precision, quantified the effects of anatomical errors on the shoulder kinematics, and showed that participants' anthropometry and sex could indirectly affect kinematic outcomes.

     

Mechanobiological osteocyte feedback drives mechanostat regulation of bone in a multiscale computational model

Biomechanics and Modeling in Mechanobiology - Significant progress has been made to identify the cells and signaling molecules involved in the mechanobiological regulation of bone remodeling. It is...     

A fully coupled poroelastic reactive-transport model of cartilage

Cartilage maintains its integrity in a hostile mechanical environment. This task is made more difficult because cartilage has no blood supply, and so nutrients and growth factors need to be transported greater distances than normal to reach cells several millimetres from the cartilage surface. The chondrocytes embedded within the extracellular matrix (ECM) are essential for maintaining the mechanical integrity of the ECM, through a balance of degradation and synthesis of collagen and proteoglycans. A chondrocyte senses various chemical and mechanical signals in its local microenvironment, responding by appropriate adaption of the local ECM. Clearly a 'systems understanding' of cartilage behaviour is of critical importance in developing an integrated understanding of both normal and abnormal physiology of cartilage. In a series of papers, we have developed a reactive-transport porous-media model to investigate the coupled processes of growth factor transport, mechanical deformation and fluid flow, and in this paper, we extend the model to include biosynthesis and degradation of matrix molecules. The model is validated using three independent experimental data sets, it being found that a single set of parameters described the experimental results remarkably well. The model is then employed to make predictions about changes in proteoglycan content under a variety of conditions. This model may prove useful in predicting the behaviour of tissue engineering constructs, or predicting the outcome of repair processes in cartilage.     

Bone refilling in cortical basic multicellular units: insights into tetracycline double labelling from a computational model

Bone remodelling is carried out by ‘bone multicellular units’ ( $\text{ BMU }$ s) in which active osteoclasts and active osteoblasts are spatially and temporally coupled. The refilling of new bone by osteoblasts towards the back of the $\text{ BMU }$ occurs at a rate that depends both on the number of osteoblasts and on their secretory activity. In cortical bone, a linear phenomenological relationship between matrix apposition rate and $\text{ BMU }$ cavity radius is found experimentally. How this relationship emerges from the combination of complex, nonlinear regulations of osteoblast number and secretory activity is unknown. Here, we extend our previous mathematical model of cell development within a single cortical $\text{ BMU }$ to investigate how osteoblast number and osteoblast secretory activity vary along the $\text{ BMU }$ ’s closing cone. The mathematical model is based on biochemical coupling between osteoclasts and osteoblasts of various maturity and includes the differentiation of osteoblasts into osteocytes and bone lining cells, as well as the influence of $\text{ BMU }$ cavity shrinkage on osteoblast development and activity. Matrix apposition rates predicted by the model are compared with data from tetracycline double labelling experiments. We find that the linear phenomenological relationship observed in these experiments between matrix apposition rate and $\text{ BMU }$ cavity radius holds for most of the refilling phase simulated by our model, but not near the start and end of refilling. This suggests that at a particular bone site undergoing remodelling, bone formation starts and ends rapidly, supporting the hypothesis that osteoblasts behave synchronously. Our model also suggests that part of the observed cross-sectional variability in tetracycline data may be due to different bone sites being refilled by $\text{ BMU }$ s at different stages of their lifetime. The different stages of a $\text{ BMU }$ ’s lifetime (such as initiation stage, progression stage, and termination stage) depend on whether the cell populations within the $\text{ BMU }$ are still developing or have reached a quasi-steady state whilst travelling through bone. We find that due to their longer lifespan, active osteoblasts reach a quasi-steady distribution more slowly than active osteoclasts. We suggest that this fact may locally enlarge the Haversian canal diameter (due to a local lack of osteoblasts compared to osteoclasts) near the $\text{ BMU }$ ’s point of origin.     

Theoretical investigation of the role of the RANK-RANKL-OPG system in bone remodeling

The RANK-RANKL-OPG system is an essential signaling pathway involved in bone cell-cell communication, with ample evidence that modification of the RANK-RANKL-OPG signaling pathway has major effects on bone remodeling. The first focus of this paper is to demonstrate that a theoretical model of bone cell-cell interactions is capable of qualitatively reproducing changes in bone associated with RANK-RANKL-OPG signaling. To do this we consider either biological experiments or bone diseases related to receptor and/or ligand deficiencies, including RANKL over-expression, ablation of OPG production and/or RANK receptor modifications. The second focus is to investigate a wide range of possible therapeutic strategies for re-establishing bone homeostasis for various pathologies of the RANK-RANKL-OPG pathway. These simulations indicate that bone diseases associated with the RANK-RANKL-OPG pathway are very effective in triggering bone resorption compared to bone formation. These results align with Hofbauer's “convergence hypothesis”, which states that catabolic bone diseases most effectively act through the RANK-RANKL-OPG system. Additionally, we demonstrate that severity of catabolic bone diseases strongly depends on how many components of this pathway are affected. Using optimization algorithms and the theoretical model, we identify a variety of successful “virtual therapies” for different disease states using both single and dual therapies.     

The effect of cyclic deformation and solute binding on solute transport in cartilage

Diffusive transport must play an important role in transporting nutrients into cartilage due to its avascular nature. Recent theoretical studies generally support the idea that cyclic loading enhances large molecule transport through advection. However, to date, reactive transport, i.e. the effects of solute binding, has not yet been taken into consideration in cyclically deformed cartilage. In the present study, we develop a reactive transport model to describe the potential role of binding of solute within cyclically deformed cartilage. Our results show that binding does have a significant effect on transport, particularly for the low IGF-I concentrations typical of synovial fluid. A dynamic loading regime of high strain magnitudes (up to 10%) in combination with high frequencies (e.g. 1 Hz) was seen to produce the most dramatic results with enhanced total uptake ratio as high as 25% averaged over the first 5h of cyclic loading.     

Mathematical modelling of bone adaptation of the metacarpal subchondral bone in racehorses

In Thoroughbred racehorses, fractures of the distal limb are commonly catastrophic. Most of these fractures occur due to the accumulation of fatigue damage from repetitive loading, as evidenced by microdamage at the predilection sites for fracture. Adaptation of the bone in response to training loads is important for fatigue resistance. In order to better understand the mechanism of subchondral bone adaptation to its loading environment, we utilised a square root function defining the relationship between bone volume fraction \((f_{BM} )\) and specific surface \((S_v )\) of the subchondral bone of the lateral condyles of the third metacarpal bone (MCIII) of the racehorse, and using this equation, developed a mathematical model of subchondral bone that adapts to loading conditions observed in vivo. The model is expressed as an ordinary differential equation incorporating a formation rate that is dependent on strain energy density. The loading conditions applied to a selected subchondral region, i.e. volume of interest, were estimated based on joint contact forces sustained by racehorses in training. For each of the initial conditions of \(f_{BM} \) we found no difference between subsequent homoeostatic \(f_{BM} \) at any given loading condition, but the time to reach equilibrium differed by initial \(f_{BM} \) and loading condition. We found that the observed values for \(f_{BM} \) from the mathematical model output were a good approximation to the existing data for racehorses in training or at rest. This model provides the basis for understanding the effect of changes to training strategies that may reduce the risk of racehorse injury.