Predicting cortical bone adaptation to axial loading in the mouse tibia

Department of Bioengineering, Imperial College London, London, UK
Department of Comparative Biomedical Sciences, The Royal Veterinary College, London, UK
Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA, USA
DOI
10.7287/peerj.preprints.1140v1
Subject Areas
Bioengineering, Computational Biology, Orthopedics
Keywords
Bonemechanobiology, Bone, mechanobiology, fluid-flow, mouse, tibia, FEA
Copyright
© 2015 Pereira et al.
Licence
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ PrePrints) and either DOI or URL of the article must be cited.
Cite this article
Pereira AF, Javaheri B, Pitsillides A, Shefelbine S. 2015. Predicting cortical bone adaptation to axial loading in the mouse tibia. PeerJ PrePrints 3:e1140v1

Abstract

The development of predictive mathematical models can contribute to a deeper understanding of the specific stages of bone mechanobiology and the process by which bone adapts to mechanical forces. The objective of this work was to predict, with spatial accuracy, cortical bone adaptation to mechanical load, in order to better understand the mechanical cues that might be driving adaptation. The axial tibial loading model was used to trigger cortical bone adaptation in C57BL/6 mice and provide relevant biological and biomechanical information. A method for mapping cortical thickness in the mouse tibia diaphysis was developed, allowing for a thorough spatial description of where bone adaptation occurs. Poroelastic finite-element (FE) models were used to determine the structural response of the tibia upon axial loading and interstitial fluid velocity as the mechanical stimulus. FE models were coupled with mechanobiological governing equations, which accounted for non-static loads and assumed that bone responds instantly to local mechanical cues in an on-off manner.

Author Comment

This is the 2nd version of a submission to Journal of The Royal Society Interface.

Supplemental Information

Figure S1 and Figure S2

Figure S1: Vector field visualisation of peak interstitial fluid velocity during load (a) and unload (b). Cross sections at Z = 0:5 with contour plots of pore pressure, p, and unit length fluid flow vectors, V, represented by the red arrows (right). Regions in compression experience an influx of fluid during loading, while the reverse happens in regions in tension. Figure S2: Boxplot representation of tibial changes in second moment of area about the minor axis, as a function of the normalised diaphyseal length Z.

DOI: 10.7287/peerj.preprints.1140v1/supp-1