Prolonging disuse in aged mice amplifies cortical but not trabecular bones’ response to mechanical loading

Objective: Short-term neurectomy-induced disuse (SN) has been shown to restore load responses in aged mice. We examined whether this restoration was further enhanced in both cortical and trabecular bone by simply extending the SN. Methods: Following load: strain calibration, tibiae in female C57BL/J6 mice at 8, 14 and 20 weeks and 18 months (n=8/group) were loaded and bone changes measured. Effects of long-term SN examined in twenty-six 18 months-old mice, neurectomised for 5 or 100 days with/without subsequent loading. Cortical and trabecular responses were measured histomorphometrically or by micro-computed tomography. Results: Loading increased new cortical bone formation, elevating cross-sectional area in 8, 14 and 20 week-old (p <0.05), but not 18 month-old aged mice. Histomorphometry showed that short-term SN reinstated load-responses in aged mice, with significant 33% and 117% increases in bone accrual at 47% and 37%, but not 27% of tibia length. Cortical responses to loading was heightened and widespread, now evident at all locations, following prolonged SN (108, 167 and 98% at 47, 37 and 27% of tibial length, respectively). In contrast, loading failed to modify trabecular bone mass or architecture. Conclusions: Mechanoadaptation become deficient with ageing and prolonging disuse amplifies this response in cortical but not trabecular bone

Structure-mechanics relationships of collagen fibrils in the Osteogenesis Imperfecta Mouse model

The collagen molecule, which is the building block of collagen fibrils, is a triple helix of two α1(I) chains and one α2(I) chain. However, in the severe mouse model of osteogenesis imperfecta (OIM), deletion of the COL1A2 gene results in the substitution of the α2(I) chain by one α1(I) chain. As this substitution severely impairs the structure and mechanics of collagen-rich tissues at the tissue and organ level, the main aim of this study was to investigate how the structure and mechanics are altered in OIM collagen fibrils. Comparing results from atomic force microscopy imaging and cantilever-based nanoindentation on collagen fibrils from OIM and wild-type (WT) animals, we found a 33% lower indentation modulus in OIM when air-dried (bound water present) and an almost fivefold higher indentation modulus in OIM collagen fibrils when fully hydrated (bound and unbound water present) in phosphate-buffered saline solution (PBS) compared with WT collagen fibrils. These mechanical changes were accompanied by an impaired swelling upon hydration within PBS. Our experimental and atomistic simulation results show how the structure and mechanics are altered at the individual collagen fibril level as a result of collagen gene mutation in OIM. We envisage that the combination of experimental and modelling approaches could allow mechanical phenotyping at the collagen fibril level of virtually any alteration of collagen structure or chemistry.United States. Dept. of Defense. Presidential Early Career Award for Scientists and EngineersNational Science Foundation (U.S.) (CAREER Award

Predicting cortical bone adaptation to axial loading in the mouse tibia

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. The presented formulation was able to simulate the areas of adaptation and accurately reproduce the distributions of cortical thickening observed in the experimental data with a statistically significant positive correlation (Kendall's τ rank coefficient τ = 0.51, p < 0.001). This work demonstrates that computational models can spatially predict cortical bone mechanoadaptation to a time variant stimulus. Such models could be used in the design of more efficient loading protocols and drug therapies that target the relevant physiological mechanisms

Limb bone scaling in hopping diprotodonts and quadrupedal artiodactyls

Bone adaptation is modulated by the timing, direction, rate, and magnitude of mechanical loads. To investigate whether frequent slow, or infrequent fast, gaits could dominate bone adaptation to load, we compared scaling of the limb bones from two mammalian herbivore clades that use radically different high-speed gaits, bipedal hopping and quadrupedal galloping. Forelimb and hindlimb bones were collected from 20 artiodactyl and 15 diprotodont species (body mass M 1.05 - 1536 kg) and scanned in clinical computed tomography or X-ray microtomography. Second moment of area (Imax) and bone length (l) were measured. Scaling relations (y = axb) were calculated for l vs M for each bone and for Imax vs M and Imax vs l for every 5% of length. Imax vs M scaling relationships were broadly similar between clades despite the diprotodont forelimb being nearly unloaded, and the hindlimb highly loaded, during bipedal hopping. Imax vs l and l vs M scaling were related to locomotor and behavioural specialisations. Low-intensity loads may be sufficient to maintain bone mass across a wide range of species. Occasional high-intensity gaits might not break through the load sensitivity saturation engendered by frequent low-intensity gaits

Limb bone scaling in hopping macropods and quadrupedal artiodactyls

Bone adaptation is modulated by the timing, direction, rate and magnitude of mechanical loads. To investigate whether frequent slow, or infrequent fast, gaits could dominate bone adaptation to load, we compared scaling of the limb bones from two mammalian herbivore clades that use radically different high-speed gaits, bipedal hopping (suborder Macropodiformes; kangaroos and kin) and quadrupedal galloping (order Artiodactyla; goats, deer and kin). Forelimb and hindlimb bones were collected from 20 artiodactyl and 15 macropod species (body mass M 1.05–1536 kg) and scanned in computed tomography or X-ray microtomography. Second moment of area (Imax) and bone length (l) were measured. Scaling relations (y = axb) were calculated for l versus M for each bone and for Imax versus M and Imax versus l for every 5% of length. Imax versus M scaling relationships were broadly similar between clades despite the macropod forelimb being nearly unloaded, and the hindlimb highly loaded, during bipedal hopping. Imax versus l and l versus M scaling were related to locomotor and behavioural specializations. Low-intensity loads may be sufficient to maintain bone mass across a wide range of species. Occasional high-intensity gaits might not break through the load sensitivity saturation engendered by frequent low-intensity gaits

Studies of chain substitution caused sub-fibril level differences in stiffness and ultrastructure of wildtype and oim/oim collagen fibers using multifrequency-AFM and molecular modeling.

Molecular alteration in type I collagen, i.e., substituting the α2 chain with α1 chain in tropocollagen molecule, can cause osteogenesis imperfecta (OI), a brittle bone disease, which can be represented by a mouse model (oim/oim). In this work, we use dual-frequency Atomic Force Microscopy (AFM) and incorporated with molecular modeling to quantify the ultrastructure and stiffness of the individual native collagen fibers from wildtype (+/+) and oim/oim diseased mice humeri. Our work presents direct experimental evidences that the +/+ fibers have highly organized and compact ultrastructure and corresponding ordered stiffness distribution. In contrast, oim/oim fibers have ordered but loosely packed ultrastructure with uncorrelated stiffness distribution, as well as local defects. The molecular model also demonstrates the structural and molecular packing differences between +/+ and oim/oim collagens. The molecular mutation significantly altered sub-fibril structure and mechanical property of collagen fibers. This study can give the new insight for the mechanisms and treatment of the brittle bone disease

Altered lacunar and vascular porosity in osteogenesis imperfecta mouse bone as revealed by synchrotron tomography contributes to bone fragility

Osteogenesis imperfecta (brittle bone disease) is caused by mutations in the collagen genes and results in skeletal fragility. Changes in bone porosity at the tissue level indicate changes in bone metabolism and alter bone mechanical integrity. We investigated the cortical bone tissue porosity of a mouse model of the disease, oim, in comparison to a wild type (WT-C57BL/6), and examined the influence of canal architecture on bone mechanical performance.High-resolution 3D representations of the posterior tibial and the lateral humeral mid-diaphysis of the bones were acquired for both mouse groups using synchrotron radiation-based computed tomography at a nominal resolution of 700 nm. Volumetric morphometric indices were determined for cortical bone, canal network and osteocyte lacunae. The influence of canal porosity architecture on bone mechanics was investigated using microarchitectural finite element (?FE) models of the cortical bone. Bright-field microscopy of stained sections was used to determine if canals were vascular.Although total cortical porosity was comparable between oim and WT bone, oim bone had more numerous and more branched canals (p &lt; 0.001), and more osteocyte lacunae per unit volume compared to WT (p &lt; 0.001). Lacunae in oim were more spherical in shape compared to the ellipsoidal WT lacunae (p &lt; 0.001). Histology revealed blood vessels in all WT and oim canals. ?FE models of cortical bone revealed that small and branched canals, typical of oim bone, increase the risk of bone failure. These results portray a state of compromised bone quality in oim bone at the tissue level, which contributes to its deficient mechanical properties

Spatial relationship between bone formation and mechanical stimulus within cortical bone: Combining 3D fluorochrome mapping and poroelastic finite element modelling

Bone is a dynamic tissue and adapts its architecture in response to biological and mechanical factors. Here we investigate how cortical bone formation is spatially controlled by the local mechanical environment in the murine tibia axial loading model (C57BL/6). We obtained 3D locations of new bone formation by performing ‘slice and view’3D fluorochrome mapping of the entire bone and compared these sites with the regions of high fluid velocity or strain energy density estimated using a finite element model, validated with ex-vivo bone surface strain map acquired ex-vivo using digital image correlation. For the comparison, 2D maps of the average bone formation and peak mechanical stimulus on the tibial endosteal and periosteal surface across the entire cortical surface were created. Results showed that bone formed on the periosteal and endosteal surface in regions of high fluid flow. Peak strain energy density predicted only the formation of bone periosteally. Understanding how the mechanical stimuli spatially relates with regions of cortical bone formation in response to loading will eventually guide loading regime therapies to maintain or restore bone mass in specific sites in skeletal pathologies

A multi-scale modelling framework combining musculoskeletal rigid-body simulations with adaptive finite element analyses, to evaluate the impact of femoral geometry on hip joint contact forces and femoral bone growth

Multi-scale simulations, combining muscle and joint contact force (JCF) from musculoskeletal simulations with adaptive mechanobiological finite element analysis, allow to estimate musculoskeletal loading and predict femoral growth in children. Generic linearly scaled musculoskeletal models are commonly used. This approach, however, neglects subject- and age-specific musculoskeletal geometry, e.g. femoral neck-shaft angle (NSA) and anteversion angle (AVA). This study aimed to evaluate the impact of proximal femoral geometry, i.e. altered NSA and AVA, on hip JCF and femoral growth simulations. Musculoskeletal models with NSA ranging from 120° to 150° and AVA ranging from 20° to 50° were created and used to calculate muscle and hip JCF based on the gait analysis data of a typically developing child. A finite element model of a paediatric femur was created from magnetic resonance images. The finite element model was morphed to the geometries of the different musculoskeletal models and used for mechanobiological finite element analysis to predict femoral growth trends. Our findings showed that hip JCF increase with increasing NSA and AVA. Furthermore, the orientation of the hip JCF followed the orientation of the femoral neck axis. Consequently, the osteogenic index, which is a function of cartilage stresses and defines the growth rate, barely changed with altered NSA and AVA. Nevertheless, growth predictions were sensitive to the femoral geometry due to changes in the predicted growth directions. Altered NSA had a bigger impact on the growth results than altered AVA. Growth simulations based on mechanobiological principles were in agreement with reported changes in paediatric populations

IPCP: Immersive Parallel Coordinates Plots for Engineering Design Processes

Computational engineering design methods and tools are common practice in modern industry. Such approaches are integral in enabling designers to efficiently explore larger and more complex design spaces. However, at the same time, computational engineering design methods tend to dramatically increase the number of candidate solutions that decision-makers must interpret in order to make appropriate choices within a set of solutions. Since all candidate solutions can be represented in digital form together with their assessment criteria, evaluated according to some sort of simulation model, a natural way to explore and understand the complexities of the design problem is to visualize their multidimensional nature. The task now involves the discovery of patterns and trends within the multidimensional design space. In this work, we aim to enhance the design decision-making process by embedding visual analytics into an immersive virtual reality environment. To this end, we present a system called IPCP: immersive parallel coordinates plots. IPCP combines the well-established parallel coordinates visualization technique for high-dimensional data with immersive virtual reality. We propose this approach in order to exploit and discover efficient means to use new technology within a conventional decision-making process. The aim is to provide benefits by enhancing visualizations of 3D geometry and other physical quantities with scientific information. We present the design of this system, which allows the representation and exploration of multidimensional scientific datasets. A qualitative evaluation with two surrogate expert users, knowledgeable in multidimensional data analysis, demonstrate that the system can be used successfully to detect both known and previously unknown patterns in a real-world test dataset, producing an early indicative validation of its suitability for decision support in engineering design processes.Cambridge European and Trinity Hall; Engineering and Physical Sciences Research Council (EPSRC-1788814