Biomechanics of trabecular bone

Development of the CAVEMAN Human Body Model: Validation of Lower Extremity Sub-Injurious Response to Vertical Accelerative Loading

Improving injury prediction accuracy and fidelity for mounted Warfighters has become an area of focus for the U.S. military in response to improvised explosive device (IED) use in both Iraq and Afghanistan. Although the Hybrid III anthropomorphic test device (ATD) has historically been used for crew injury analysis, it is only capable of predicting a few select skeletal injuries. The Computational Anthropomorphic Virtual Experiment Man (CAVEMAN) human body model is being developed to expand the injury analysis capability to both skeletal and soft tissues. The CAVEMAN model is built upon the Zygote 50^th percentile male human CAD model and uses a finite element modeling approach developed for high performance computing (HPC). The lower extremity subset of the CAVEMAN human body model presented herein includes: 28 bones, 26 muscles, 40 ligaments, fascia, cartilage and skin. Sensitivity studies have been conducted with the CAVEMAN lower extremity model to determine the structures critical for load transmission through the leg in the underbody blast (UBB) environment. An evaluation of the CAVEMAN lower extremity biofidelity was also carried out using 14 unique data sets derived by the Warrior Injury Assessment Manikin (WIAMan) program cadaveric lower leg testing. Extension of the CAVEMAN lower extremity model into anatomical tissue failure will provide additional injury prediction capabilities, beyond what is currently achievable using ATDs, to improve occupant survivability analyses within military vehicles.

An investigation of dimensional scaling using cervical spine motion segment finite element models

The paucity of experimental data for validating computational models of different statures underscores the need for appropriate scaling methods so that models can be verified and validated using experimental data. Scaling was investigated using 50th percentile male (M50) and 5th percentile female (F05) cervical spine motion segment (C4-C5) finite element models subject to tension, flexion, and extension loading. Two approaches were undertaken: geometric scaling of the models to investigate size effects (volumetric scaling) and scaling of the force-displacement or moment-angle model results (data scaling). Three sets of scale factors were considered: global (body mass), regional (neck dimensions), and local (segment tissue dimensions). Volumetric scaling of the segment models from M50 to F05, and vice versa, produced correlations that were good or excellent in both tension and flexion (0.825-0.991); however, less agreement was found in extension (0.550-0.569). The reduced correlation in extension was attributed to variations in shape between the models leading to nonlinear effects such as different time to contact for the facet joints and posterior processes. Data scaling of the responses between the M50 and F05 models produced similar trends to volumetric scaling, with marginally greater correlations. Overall, the local tissue level and neck region level scale factors produced better correlations than the traditional global scaling. The scaling methods work well for a given subject, but are limited in applicability between subjects with different morphology, where nonlinear effects may dominate the response.

A novel in silico method to quantify primary stability of screws in trabecular bone

Insufficient primary stability of screws in bone leads to screw loosening and failure. Unlike conventional continuum finite-element models, micro-CT based finite-element analysis (micro-FE) is capable of capturing the patient-specific bone micro-architecture, providing accurate estimates of bone stiffness. However, such in silico models for screws in bone highly overestimate the apparent stiffness. We hypothesized that a more accurate prediction of primary implant stability of screws in bone is possible by considering insertion-related bone damage. We assessed two different screw types and loading scenarios in 20 trabecular bone specimens extracted from 12 cadaveric human femoral heads (N = 5 for each case). In the micro-FE model, we predicted specimen-specific Young's moduli of the peri-implant bone damage region based on morphometric parameters such that the apparent stiffness of each in silico model matched the experimentally measured stiffness of the corresponding in vitro specimen as closely as possible. The standard micro-FE models assuming perfectly intact peri-implant bone overestimated the stiffness by over 330%. The consideration of insertion related damaged peri-implant bone corrected the mean absolute percentage error down to 11.4% for both loading scenarios and screw types. Cross-validation revealed a mean absolute percentage error of 14.2%. We present the validation of a novel micro-FE modeling technique to quantify the apparent stiffness of screws in trabecular bone. While the standard micro-FE model overestimated the bone-implant stiffness, the consideration of insertion-related bone damage was crucial for an accurate stiffness prediction. This approach provides an important step toward more accurate specimen-specific micro-FE models. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2415-2424, 2017.

Biomechanical analysis of the posterior bony column of the lumbar spine

Background

Each part of the rear bone structure can become an anchor point for an attachment device. The objective of this study was to evaluate the stiffness and strength of different parts of the rear lumbar bone structure by axial compression damage experiments.

Methods

Five adult male lumbar bone structures from L2 to L5 were exposed. The superior and inferior articular processes, upper and lower edges of the lamina, and upper and lower edges of the spinous process were observed and isolated and then divided into six groups (n = 10). The specimens were placed between the compaction disc and the load platform in a universal testing machine, which was first preloaded to 5.0 N tension to eliminate water on the surface and then loaded to the specimen curve decline at a constant tension loading rate of 0.01 mm/s, until the specimens had been destroyed.

Results

Significant differences in mechanical properties were found among different parts of the rear lumbar bone structure. Compared with other parts, the lower edge of the lamina has good mechanical properties, which have a high modulus of elasticity; the superior and inferior articular processes have greater ultimate strength, which can withstand greater compressive loads; and the mechanical properties of the spinous process are poor, and it is significantly stiffer and weaker than the lamina and articular processes.

Conclusion

These data can be useful in future spinal biomechanics research leading to better biomechanical compatibility and provide theoretical references for spinal implant materials.

Promotion effect of extracts from plastrum testudinis on alendronate against glucocorticoid-induced osteoporosis in rat spine

Alendronate (ALN) is a key therapeutic used to treat glucocorticoid-induced osteoporosis (GIOP), but may induce severe side effects. We showed earlier that plastrum testudinis extracts (PTE) prevented and treated GIOP in vivo. However, clinically, PTE is seldom used alone. Herein, we reveal the synergistic effect of ALN and PTE can treat GIOP of the rat spine and define the mechanism. Sprague-Dawley rats were randomly assigned to four groups: a vehicle group, a GIOP group, an ALN group, and an ALN+PTE group. Each group was further divided into two experimental phases, including dexamethasone (DXM) intervention and withdrawal. Bone mass, microarchitecture, biomechanics, bone-turnover markers, and histomorphology were evaluated. The mRNA and protein expression levels of CTSK and Runx2 were detemined. We found that ALN+PTE improved bone quantity and quality, bone strength, bone turnover; and mitigated histological damage during glucocorticoid intervention and withdrawal. The therapeutic effect was better than that afforded by ALN alone. ALN+PTE reduced CTSK protein expression, promoted Runx2 mRNA and protein expression to varying extents, and more strongly inhibited bone resorption than did ALN alone. Overall, the synergistic effect mediated by ALN+PTE reversed GIOP during DXM intervention and withdrawal via affecting CTSK and Runx2 expression at mRNA and protein levels.

Functional Adaptation of the Calcaneus in Historical Foot Binding

The normal structure of human feet is optimized for shock dampening during walking and running. Foot binding was a historical practice in China aimed at restricting the growth of female feet for aesthetic reasons. In a bound foot the shock-dampening function normally facilitated by the foot arches is withdrawn, resulting in the foot functioning as a rigid extension of the lower leg. An interesting question inspiring this study regards the nature of adaptation of the heel bone to this nonphysiological function using the parameters of cancellous bone anisotropy and 3D fabric topology and a novel intertrabecular angle (ITA) analysis. We found that the trabecular microarchitecture of the normal heel bone, but not of the bound foot, adapts to function by increased anisotropy and preferred orientation of trabeculae. The anisotropic texture in the normal heel bone consistently follows the physiological stress trajectories. However, in the bound foot heel bone the characteristic anisotropy pattern fails to develop, reflecting the lack of a normal biomechanical input. Moreover, the basic topological blueprint of cancellous bone investigated by the ITA method is nearly invariant in both normal and bound foot. These findings suggest that the anisotropic cancellous bone texture is an acquired characteristic that reflects recurrent loading conditions; conversely, an inadequate biomechanical input precludes the formation of anisotropic texture. This opens a long-sought-after possibility to reconstruct bone function from its form. The conserved topological parameters characterize the generic 3D fabric of cancellous bone, which is to a large extent independent of its adaptation to recurrent loading and perhaps determines the mechanical competence of trabecular bone regardless of its functional adaptation. © 2017 American Society for Bone and Mineral Research.

Next Generation Tissue Engineering of Orthopedic Soft Tissue-to-Bone Interfaces

Soft tissue-to-bone interfaces are complex structures that consist of gradients of extracellular matrix materials, cell phenotypes, and biochemical signals. These interfaces, called entheses for ligaments, tendons, and the meniscus, are crucial to joint function, transferring mechanical loads and stabilizing orthopedic joints. When injuries occur to connected soft tissue, the enthesis must be re-established to restore function, but due to structural complexity, repair has proven challenging. Tissue engineering offers a promising solution for regenerating these tissues. This prospective review discusses methodologies for tissue engineering the enthesis, outlined in three key design inputs: materials processing methods, cellular contributions, and biochemical factors.

Shear stress and von Mises stress distributions in the periphery of an embedded acetabular cup implant during impingement

As literature implies, daily activities of total hip arthroplasty (THA) patients may include movements prone to implant-implant impingement. Thus, high shear stresses may be induced at the acetabular implant-bone interface, increasing the risk of implant loosening. The aim of the current study is to determine whether or not impingement events may pose an actual risk to acetabular periprosthetic bone. An existing experimental workflow was augmented to cover complete three-dimensional strain gage measurement. von Mises and shear stresses were calculated from 1620 measured strain values, collected around a hemispherical cup implant at 2.5 mm interface distance during worst-case impingement loading. A shear stress criterion for acetabular periprosthetic bone was derived from the literature. At the impingement site, magnitudes of von Mises stress amount to 0.57 MPa and tilting shear stress amount to -0.3 MPa at 2.5 mm interface distance. Conclusion can be drawn that worst-case impingement events are unlikely to pose a risk of bone material failure in the periphery around fully integrated cementless acetabular hip implants in otherwise healthy THA patients. As numerical predictions in the literature suggested, it can now be confirmed that impingement moments are unlikely to cause acetabular implant-bone interface fixation failures.

Lower Cervical Spine Motion Segment Computational Model Validation: Kinematic and Kinetic Response for Quasi-Static and Dynamic Loading

Advanced computational human body models (HBM) enabling enhanced safety require verification and validation at different levels or scales. Specifically, the motion segments, which are the building blocks of a detailed neck model, must be validated with representative experimental data to have confidence in segment and, ultimately, full neck model response. In this study, we introduce detailed finite element motion segment models and assess the models for quasi-static and dynamic loading scenarios. Finite element segment models at all levels in the lower human cervical spine were developed from scans of a 26-yr old male subject. Material properties were derived from the in vitro experimental data. The segment models were simulated in quasi-static loading in flexion, extension, lateral bending and axial rotation, and at dynamic rates in flexion and extension in comparison to previous experimental studies and new dynamic experimental data introduced in this study. Single-valued experimental data did not provide adequate information to assess the model biofidelity, while application of traditional corridor methods highlighted that data sets with higher variability could lead to an incorrect conclusion of improved model biofidelity. Data sets with continuous or multiple moment–rotation measurements enabled the use of cross-correlation for an objective assessment of the model and highlighted the importance of assessing all motion segments of the lower cervical spine to evaluate the model biofidelity. The presented new segment models of the lower cervical spine, assessed for range of motion and dynamic/traumatic loading scenarios, provide a foundation to construct a biofidelic model of the spine and neck, which can be used to understand and mitigate injury for improved human safety in the future.

Intraosseous monitoring of drilling in lumbar vertebrae by ultrasound: An experimental feasibility study

The rationale for this project is to evaluate the efficiency of a novel sonographic method for measurements of interosseous distances. The method utilizes a propagating ultrasonic beam through aqueous milieu which is directed as a jet into a drilled tract. We used a plastic model of human L5 vertebra and ex vivo specimen of L5 porcine vertebra and generated 2 mm in diameter tracts in vertebral pedicles. The tracts were created in the "desired" central direction and in the "wrong" medial and lateral directions. The drilled tracts and the residual, up to opposite cortex, distances were measured sonographically and mechanically and compared statistically. We show that "true" mechanical measurements can be predicted from sonographic measurements with correction of 1-3 mm. The correct central route can be distinguished from the wrong misplaced routes. By using the sonographic measurements, a correct direction of drilling in the pedicle of lumbar L5 vertebra can be efficiently monitored.

Time Dependent Behaviour of Trabecular Bone at Multiple Load Levels

The deformation of bone when subjected to loads is not instantaneous but varies with time. To investigate this time-dependent behaviour sixteen bovine trabecular bone specimens were subjected to compressive loading, creep, unloading and recovery at multiple load levels corresponding to apparent strains of 2000-25,000 με. We found that: the time-dependent response of trabecular bone comprises of both recoverable and irrecoverable strains; the strain response is nonlinearly related to applied load levels; and the response is linked to bone volume fraction. Although majority of strain is recovered after the load-creep-unload-recovery cycle some residual strain always exists. The analysis of results indicates that trabecular bone becomes stiffer initially and then experiences stiffness degradation with the increasing load levels. Steady state creep rate was found to be dependent on applied stress level and bone volume fraction with a power law relationship.

Mimetization of the elastic properties of cancellous bone via a parameterized cellular material

Bone tissue mechanical properties and trabecular microarchitecture are the main factors that determine the biomechanical properties of cancellous bone. Artificial cancellous microstructures, typically described by a reduced number of geometrical parameters, can be designed to obtain a mechanical behavior mimicking that of natural bone. In this work, we assess the ability of the parameterized microstructure introduced by Kowalczyk (Comput Methods Biomech Biomed Eng 9:135-147, 2006. doi: 10.1080/10255840600751473 ) to mimic the elastic response of cancellous bone. Artificial microstructures are compared with actual bone samples in terms of elasticity matrices and their symmetry classes. The capability of the parameterized microstructure to combine the dominant isotropic, hexagonal, tetragonal and orthorhombic symmetry classes in the proportions present in the cancellous bone is shown. Based on this finding, two optimization approaches are devised to find the geometrical parameters of the artificial microstructure that better mimics the elastic response of a target natural bone specimen: a Sequential Quadratic Programming algorithm that minimizes the norm of the difference between the elasticity matrices, and a Pattern Search algorithm that minimizes the difference between the symmetry class decompositions. The pattern search approach is found to produce the best results. The performance of the method is demonstrated via analyses for 146 bone samples.

Virtual stress testing of fracture stability in soldiers with severely comminuted tibial fractures

Virtual stress testing (VST) provides a non-invasive estimate of the strength of a healing bone through a biomechanical analysis of a patient's computed tomography (CT) scan. We asked whether VST could improve management of patients who had a tibia fracture treated with external fixation. In a retrospective case-control study of 65 soldier-patients who had tibia fractures treated with an external fixator, we performed VST utilizing CT scans acquired prior to fixator removal. The strength of the healing bone and the amount of tissue damage after application of an overload were computed for various virtual loading cases. Logistic regression identified computed outcomes with the strongest association to clinical events related to nonunion within 2 months after fixator removal. Clinical events (n = 9) were associated with a low tibial strength for compression loading (p < 0.05, AUC = 0.74) or a low proportion of failed cortical bone tissue for torsional loading (p < 0.005, AUC = 0.84). Using post-hoc thresholds of a compressive strength of four times body-weight and a proportional of failed cortical bone tissue of 5%, the test identified all nine patients who failed clinically (100% sensitivity; 40.9% positive predictive value) and over three fourths of those (43 of 56) who progressed to successful healing (76.8% specificity; 100% negative predictive value). In this study, VST identified all patients who progressed to full, uneventful union after fixator removal; thus, we conclude that this new test has the potential to provide a quantitative, objective means of identifying tibia-fracture patients who can safely resume weight bearing. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:805-811, 2017.

Validation of a new multiscale finite element analysis approach at the distal radius

High-resolution peripheral computed tomography is commonly used to evaluate mechanical behavior of the distal radius microstructure using micro-finite element (FE) modeling. However, only a 9mm section is considered and boundary conditions (BCs) are usually simplified (platen-compression), and may not represent physiologic loading. Regardless, these methods are increasingly being used for clinical evaluations. Our goal was to develop and validate a novel multiscale solution that allows for physiologically relevant loading simulations (such as bracing during a fall), and show that mechanical behavior in the distal radius is different under platen BCs. Our approach incorporated bone microstructure together with organ-level radius geometry, by replacing matching continuum regions with micro-FE sections in user-defined regions of interest. Multiscale model predicted strains showed a strong correlation and a significant relationship with measured strains (r=0.836, p<0.001; slope=0.881, intercept=-12.17 µε, p<0.001). Interestingly, platen BC simulated strains were almost 50% lower than measured strains (r=0.835, p<0.001), and strain distributions were clearly different. Our multiscale method demonstrated excellent potential as a computationally efficient alternative for observing true mechanical environment within distal radius microstructure under physiologically accurate loading.

Characterization of synthetic foam structures used to manufacture artificial vertebral trabecular bone

Artificial materials reflecting the mechanical properties of human bone are essential for valid and reliable implant testing and design. They also are of great benefit for realistic simulation of surgical procedures. The objective of this study was therefore to characterize two groups of self-developed synthetic foam structures by static compressive testing and by microcomputed tomography. Two mineral fillers and varying amounts of a blowing agent were used to create different expansion behavior of the synthetic open-cell foams. The resulting compressive and morphometric properties thus differed within and also slightly between both groups. Apart from the structural anisotropy, the compressive and morphometric properties of the synthetic foam materials were shown to mirror the respective characteristics of human vertebral trabecular bone in good approximation. In conclusion, the artificial materials created can be used to manufacture valid synthetic bones for surgical training. Further, they provide novel possibilities for studying the relationship between trabecular bone microstructure and biomechanical properties.

Sex dependent mechanical properties of the human mandibular condyle

The mandibular condyle consists of articular cartilage and subchondral bone that play an important role in bearing loads at the temporomandibular joint (TMJ) during static occlusion and dynamic mastication. The objective of the current study was to examine effects of sex and cartilage on 1) static and dynamic mechanical analysis (DMA) based dynamic energy storage and dissipation for the cartilage-subchondral bone construct of the human mandibular condyle, and 2) their correlations with the tissue mineral density and trabecular morphological parameters of subchondral bone. Cartilage-subchondral bone constructs were obtained from 16 individual human cadavers (9 males, 7 females, 79.00±13.10 years). After scanning with micro-computed tomography, the specimens were subjected to a non-destructive compressive static loading up to 7N and DMA using a cyclic loading profile (-5±2N at 2Hz). After removing the cartilage from the same specimen, the series of loading experiments were repeated. Static stiffness (K) and energy dissipation (W), and dynamic storage (K'), loss (K'') stiffness, and energy dissipation (tan δ) were assessed. Gray values, which are proportional to degree of bone mineralization, and trabecular morphological parameters of the subchondral bone were also measured. After removal of the cartilage, static energy dissipation significantly decreased (p<0.009) but dynamic energy dissipation was not influenced (p>0.064). Many subchondral bone properties were significantly correlated with the overall mechanical behavior of the cartilage-subchondral bone constructs for males (p<0.047) but not females (p>0.054). However, after removal of cartilage from the constructs, all of the significant correlations were no longer found (p>0.057). The current findings indicate that the subchondral bone is responsible for bearing static and dynamic loading in males but not in females. This result indicates that the female condyle may have a mechanically disadvantageous TMJ loading environment.

Does cortical bone thickness in the last sacral vertebra differ among tail types in primates?

OBJECTIVES

The external morphology of the sacrum is demonstrably informative regarding tail type (i.e., tail presence/absence, length, and prehensility) in living and extinct primates. However, little research has focused on the relationship between tail type and internal sacral morphology, a potentially important source of functional information when fossil sacra are incomplete. Here, we determine if cortical bone cross-sectional thickness of the last sacral vertebral body differs among tail types in extant primates and can be used to reconstruct tail types in extinct primates.

MATERIALS AND METHODS

Cortical bone cross-sectional thickness in the last sacral vertebral body was measured from high-resolution CT scans belonging to 20 extant primate species (N = 72) assigned to tail type categories ("tailless," "nonprehensile short-tailed," "nonprehensile long-tailed," and "prehensile-tailed"). The extant dataset was then used to reconstruct the tail types for four extinct primate species.

RESULTS

Tailless primates had significantly thinner cortical bone than tail-bearing primates. Nonprehensile short-tailed primates had significantly thinner cortical bone than nonprehensile long-tailed primates. Cortical bone cross-sectional thickness did not distinguish between prehensile-tailed and nonprehensile long-tailed taxa. Results are strongly influenced by phylogeny. Corroborating previous studies, Epipliopithecus vindobonensis was reconstructed as tailless, Archaeolemur edwardsi as long-tailed, Megaladapis grandidieri as nonprehensile short-tailed, and Palaeopropithecus kelyus as nonprehensile short-tailed or tailless.

CONCLUSIONS

Results indicate that, in the context of phylogenetic clade, measures of cortical bone cross-sectional thickness can be used to allocate extinct primate species to tail type categories.

Effect of composition and macropore percentage on mechanical and in vitro cell proliferation and differentiation properties of 3D printed HA/β-TCP scaffolds

HA/β-TCP scaffolds were fabricated by 3D printing and exhibited desirable biocompatibilityin vitro.

Bone: A Fertile Soil for Cancer Metastasis

Bone is one of the most common and most dangerous sites for metastatic growth across cancer types, and bone metastasis remains incurable. Unfortunately, the processes by which cancers preferentially metastasize to bone are still not well understood. In this review, we summarize the morphological features, physical properties, and cell signaling events that make bone a unique site for metastasis and bone remodeling. The signaling crosstalk between the tumor cells and bone cells begins a vicious cycle - a self-sustaining feedback loop between the tumor cells and the bone microenvironment composed of osteoclasts, osteoblasts, other bone marrow cells, bone matrix, and vasculature to support both tumor growth and bone destruction. Through this crosstalk, bone provides a fertile microenvironment that can harbor dormant tumor cells, sometimes for long periods, and support their growth by releasing cytokines as the bone matrix is destroyed, similar to providing nutrients for a seed to germinate in soil. However, few models exist to study the late stages of bone colonization by metastatic tumor cells. We describe some of the current methodologies used to study bone metastasis, highlighting the limitations of these methods and alternative future strategies to be used to study bone metastasis. While <i>in vivo</i> animal and patient studies may provide the gold standard for studying metastasis, <i>ex vivo</i> models can be used as an alternative to enable more controlled experiments designed to study the late stages of bone metastasis.

A new multiscale micromechanical model of vertebral trabecular bones

A new three-dimensional (3D) multiscale micromechanical model has been suggested as adept at predicting the overall linear anisotropic mechanical properties of a vertebral trabecular bone (VTB) highly porous microstructure. A nested 3D modeling analysis framework spanning the multiscale nature of the VTB is presented herein. This hierarchical analysis framework employs the following micromechanical methods: the 3D parametric high-fidelity generalized method of cells (HFGMC) as well as the 3D sublaminate model. At the nanoscale level, the 3D HFGMC method is applied to obtain the effective elastic properties of a representative unit cell (RUC) representing the mineral collagen fibrils composite. Next at the submicron scale level, the 3D sublaminate model is used to generate the effective elastic properties of a repeated stack of multilayered lamellae demonstrating the nature of the trabeculae (bone-wall). Thirdly, at the micron scale level, the 3D HFGMC method is used again on a RUC of the highly porous VTB microstructure. The VTB-RUC geometries are taken from microcomputed tomography scans of VTB samples harvested from different vertebrae of human cadavers [Formula: see text]. The predicted anisotropic overall elastic properties for native VTBs are, then, examined as a function of age and sex. The predicted results of the VTBs longitudinal Young's modulus are compared to reported values found in the literature. The proposed 3D nested modeling analysis framework provides a good agreement with reported values of Young's modulus of single trabeculae as well as for VTB-RUC in the literature.

A Computational Efficient Method to Assess the Sensitivity of Finite-Element Models: An Illustration With the Hemipelvis

Assessing the sensitivity of a finite-element (FE) model to uncertainties in geometric parameters and material properties is a fundamental step in understanding the reliability of model predictions. However, the computational cost of individual simulations and the large number of required models limits comprehensive quantification of model sensitivity. To quickly assess the sensitivity of an FE model, we built linear and Kriging surrogate models of an FE model of the intact hemipelvis. The percentage of the total sum of squares (%TSS) was used to determine the most influential input parameters and their possible interactions on the median, 95th percentile and maximum equivalent strains. We assessed the surrogate models by comparing their predictions to those of a full factorial design of FE simulations. The Kriging surrogate model accurately predicted all output metrics based on a training set of 30 analyses (R2 = 0.99). There was good agreement between the Kriging surrogate model and the full factorial design in determining the most influential input parameters and interactions. For the median, 95th percentile and maximum equivalent strain, the bone geometry (60%, 52%, and 76%, respectively) was the most influential input parameter. The interactions between bone geometry and cancellous bone modulus (13%) and bone geometry and cortical bone thickness (7%) were also influential terms on the output metrics. This study demonstrates a method with a low time and computational cost to quantify the sensitivity of an FE model. It can be applied to FE models in computational orthopaedic biomechanics in order to understand the reliability of predictions.

Ex vivo loading of trussed implants for spine fusion induces heterogeneous strains consistent with homeostatic bone mechanobiology

A truss structure was recently introduced as an interbody fusion cage. As a truss system, some of the connected elements may be in a state of compression and others in tension. This study aimed to quantify both the mean and variance of strut strains in such an implant when loaded in a simulated fusion condition with vertebral body or contoured plastic loading platens ex vivo. Cages were each instrumented with 78 fiducial spheres, loaded between platens (vertebral body or contoured plastic), imaged using high resolution micro-CT, and analyzed for deformation and strain of each of the 221 struts. With repeated loading of a cage by vertebral platens, the distribution (variance, indicated by SD) of strut strains widened from 50N control (4±114με, mean±SD) to 1000N (-23±273με) and 2000N (-48±414με), and between 1000N and 2000N. With similar loading of multiple cages, the strain distribution at 2000N (23±389με) increased from 50N control. With repeated loading by contoured plastic platens, induced strains at 2000N had a distribution similar to that induced by vertebral platens (84±426με). In all studies, cages exhibited increases in strut strain amplitude when loaded from 50N to 1000N or 2000N. Correspondingly, at 2000N, 59-64% of struts exhibited strain amplitudes consistent with mechanobiologically-regulated bone homeostasis. At 2000N, vertically-oriented struts exhibited deformation of -2.87±2.04μm and strain of -199±133με, indicating overall cage compression. Thus, using an ex vivo 3-D experimental biomechanical analysis method, a truss implant can have strains induced by physiological loading that are heterogeneous and of amplitudes consistent with mechanobiological bone homeostasis.

SEM and TEM for structure and properties characterization of bacterial cellulose/hydroxyapatite composites

Preparation of composites with different properties and gradient of components is aimed at better performance of materials for bone substitution. Bacterial cellulose-hydroxyapatite (BC-HAP) composites with various mass ratio of the components (BC-25HAP, BC-4HAP, and BC-HAP) were prepared by a novel method of growing HAP nanocrystals (the linear size ≤30 nm) in water solutions in the presence of the BC gel-film micro-fragments. Varying the BC-HAP ratios leads to a gradual change of the physical properties of the materials. It was found that an increase in the BC content results in a decrease of the HAP crystal length and specific surface area, porosity, and pore volume while the values of density and Young's modulus values increase. SCANNING 38:757-765, 2016. © 2016 Wiley Periodicals, Inc.

Inter-trabecular angle: A parameter of trabecular bone architecture in the human proximal femur that reveals underlying topological motifs

Trabecular bone is an intricate 3D network of struts and plates. Although the structure-function relations in trabecular bone have been studied since the time of Julius Wolff, controversy still exists regarding the architectural parameters responsible for its stability and resilience. We present a parameter that measures the angle between two connected trabeculae - the Inter-Trabecular Angle (ITA). We studied the ITA values derived from μCT scans of different regions of the proximal femora of 5 individuals of different age and sex. We show that the ITA angle distribution of nodes with 3 connecting trabeculae has a mean close to 120°, nodes with 4 connecting trabeculae has a mean close to 109° and nodes of higher connectivity have mean ITA values around 100°. This tendency to spread the ITAs around geometrically symmetrical motifs is highly conserved. The implication is that the ITAs are optimized such that the smallest amount of material spans the maximal 3D volume, and possibly by so doing trabecular bone might be better adapted to multidirectional loading. We also draw a parallel between trabecular bone and tensegrity structures - where lightweight, resilient and stable tetrahedron-based shapes contribute to strain redistribution amongst all the elements and to collective impact dampening.

STATEMENT OF SIGNIFICANCE

The Inter-Trabecular Angle (ITA) is a new topological parameter of trabecular bone. The ITA characterizes the way trabeculae connect with each other at nodes, regardless of their thickness and shape. The mean ITA value of nodes with 3 trabeculae is close to 120°, of nodes with 4 trabeculae is just below 109°, and the mean ITA of nodes with 5 and more trabeculae is around 100°. Thus the connections of trabeculae trend towards adopting symmetrical shapes. This implies that trabeculae can maximally span 3D space using the minimal amount of material. We draw a parallel between this motif and the concept of tensegrity - an engineering premise to which many living creatures conform at multiple levels of organization.

Trabecular architecture in the thumb of Pan and Homo: implications for investigating hand use, loading, and hand preference in the fossil record

OBJECTIVES

Humans display an 85-95% cross-cultural right-hand bias in skilled tasks, which is considered a derived behavior because such a high frequency is not reported in wild non-human primates. Handedness is generally considered to be an evolutionary byproduct of selection for manual dexterity and augmented visuo-cognitive capabilities within the context of complex stone tool manufacture/use. Testing this hypothesis requires an understanding of when appreciable levels of right dominant behavior entered the fossil record. Because bone remodels in vivo, skeletal asymmetries are thought to reflect greater mechanical loading on the dominant side, but incomplete preservation of external morphology and ambiguities about past loading environments complicate interpretations. We test if internal trabecular bone is capable of providing additional information by analyzing the thumb of Homo sapiens and Pan.

MATERIALS AND METHODS

We assess trabecular structure at the distal head and proximal base of paired (left/right) first metacarpals using micro-CT scans of Homo sapiens (n = 14) and Pan (n = 9). Throughout each epiphysis we quantify average and local bone volume fraction (BV/TV), degree of anisotropy (DA), and elastic modulus (E) to address bone volume patterning and directional asymmetry.

RESULTS

We find a right directional asymmetry in H. sapiens consistent with population-level handedness, but also report a left directional asymmetry in Pan that may be the result of postural and/or locomotor loading.

CONCLUSION

We conclude that trabecular bone is capable of detecting right/left directional asymmetry, but suggest coupling studies of internal structure with analyses of other skeletal elements and cortical bone prior to applications in the fossil record.

Nonlinear viscoelastic characterization of bovine trabecular bone

The time-independent elastic properties of trabecular bone have been extensively investigated, and several stiffness–density relations have been proposed. Although it is recognized that trabecular bone exhibits time-dependent mechanical behaviour, a property of viscoelastic materials, the characterization of this behaviour has received limited attention. The objective of the present study was to investigate the time-dependent behaviour of bovine trabecular bone through a series of compressive creep–recovery experiments and to identify its nonlinear constitutive viscoelastic material parameters. Uniaxial compressive creep and recovery experiments at multiple loads were performed on cylindrical bovine trabecular bone samples (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 19$$\end{document}n=19). Creep response was found to be significant and always comprised of recoverable and irrecoverable strains, even at low stress/strain levels. This response was also found to vary nonlinearly with applied stress. A systematic methodology was developed to separate recoverable (nonlinear viscoelastic) and irrecoverable (permanent) strains from the total experimental strain response. We found that Schapery’s nonlinear viscoelastic constitutive model describes the viscoelastic response of the trabecular bone, and parameters associated with this model were estimated from the multiple load creep–recovery (MLCR) experiments. Nonlinear viscoelastic recovery compliance was found to have a decreasing and then increasing trend with increasing stress level, indicating possible stiffening and softening behaviour of trabecular bone due to creep. The obtained parameters from MLCR tests, expressed as second-order polynomial functions of stress, showed a similar trend for all the samples, and also demonstrate stiffening–softening behaviour with increasing stress.

Post-euthanasia micro-computed tomography-based strain analysis is able to represent quasi-static in vivo behavior of whole vertebrae

Three-dimensional image-based strain measurement in whole bones allows representation of physiological, albeit quasi-static, loading conditions. However, such work to date has been limited to specimens postmortem. The main purpose of this study is to verify the efficacy of deformable image registration of post-euthanasia strain to characterize the in vivo mechanical behavior of rat vertebrae. A micro-computed tomography-compatible custom loading device was used to apply 75 N load to a three-level caudal motion segment of a healthy rat. Loaded and unloaded micro-computed tomography scans were acquired in vivo and post-sacrifice. A micro-computed tomography-based deformable image registration algorithm was used to calculate vertebral strains live and post-euthanasia. No significant difference was found in the in vivo strains (-0.011 ± 0.001) and ex vivo strains (-0.012 ± 0.001) obtained from the comparisons of loaded and unloaded images (p = 0.3). Comparisons between unloaded-unloaded and loaded-loaded scans yielded significantly lower axial strains, representing the error of the method. Qualitatively, high strains were observed adjacent to growth plate regions in evaluating the loaded-unloaded images. Strain patterns in the loaded-loaded and unloaded-unloaded scans were inconsistent as would be expected in representing noise. Overall, live and dead loaded to unloaded comparisons yielded similar strain patterns and magnitudes. Point-wise differences in axial strain fields also supported this observation. This study demonstrated a proof of concept, suggesting that post-euthanasia micro-computed tomography-based strain analysis is able to represent the in vivo quasi-static behavior of rat tail vertebrae.

Genetics of aging bone

With aging, the skeleton experiences a number of changes, which include reductions in mass and changes in matrix composition, leading to fragility and ultimately an increase of fracture risk. A number of aspects of bone physiology are controlled by genetic factors, including peak bone mass, bone shape, and composition; however, forward genetic studies in humans have largely concentrated on clinically available measures such as bone mineral density (BMD). Forward genetic studies in rodents have also heavily focused on BMD; however, investigations of direct measures of bone strength, size, and shape have also been conducted. Overwhelmingly, these studies of the genetics of bone strength have identified loci that modulate strength via influencing bone size, and may not impact the matrix material properties of bone. Many of the rodent forward genetic studies lacked sufficient mapping resolution for candidate gene identification; however, newer studies using genetic mapping populations such as Advanced Intercrosses and the Collaborative Cross appear to have overcome this issue and show promise for future studies. The majority of the genetic mapping studies conducted to date have focused on younger animals and thus an understanding of the genetic control of age-related bone loss represents a key gap in knowledge.

A Computational Approximation to Model Variation in Cancellous Bone Screw Pull-Out

Cancellous bone screws are used to achieve good pull-out characteristics when connected to cancellous bone. This study examines some screw characteristics, including pitch and inner diameter, using a model of cancellous bone with a range of bone apparent densities (ADs). This was achieved using bone geometry based on microCT-scanned cancellous bone and converted into a geometric model using mimics® software. The finite element (FE) models were produced in ansys®. The calculated reaction force for pull-out of 0.2 mm shows the influence of design parameters. Change in the proximal half angle increased the stiffness by about 15% in line with the experimental findings of others. An increase in pull-out reaction force with an increase in bone AD was also observed. However, when a particular screw geometry in lower AD bone was modeled and then rotated through 180 deg on plan, a significant reduction in reaction force was noted. Further models with screws of similar geometry in the same location showed similar reductions in reaction force and hence pull-out stiffness. Examination of the geometry of the bone/screw interface indicates that in certain positions there is little cancellous bone to support the implant—leading to low pull-out reaction forces, which is very difficult to predict. The study also examined the effect of increasing the bone stiffness adjacent to the implant, and concluded that, even in bone of low AD, increases in pull-out stiffness might be achieved.

Microstructural and compositional contributions towards the mechanical behavior of aging human bone measured by cyclic and impact reference point indentation

The assessment of fracture risk often relies primarily on measuring bone mineral density, thereby accounting for only a single pathology: the loss of bone mass. However, bone's ability to resist fracture is a result of its biphasic composition and hierarchical structure that imbue it with high strength and toughness. Reference point indentation (RPI) testing is designed to directly probe bone mechanical behavior at the microscale in situ, although it remains unclear which aspects of bone composition and structure influence the results at this scale. Therefore, our goal in this study was to investigate factors that contribute to bone mechanical behavior measured by cyclic reference point indentation, impact reference point indentation, and three-point bending. Twenty-eight female cadavers (ages 57-97) were subjected to cyclic and impact RPI in parallel at the unmodified tibia mid-diaphysis. After RPI, the middiaphyseal tibiae were removed, scanned using micro-CT to obtain cortical porosity (Ct.Po.) and tissue mineral density (TMD), then tested using three-point bending, and lastly assayed for the accumulation of advanced glycation end-products (AGEs). Both the indentation distance increase from cyclic RPI (IDI) and bone material strength index from impact RPI (BMSi) were significantly correlated with TMD (r=-0.390, p=0.006; r=0.430, p=0.002; respectively). Accumulation of AGEs was significantly correlated with IDI (r=0.281, p=0.046), creep indentation distance (CID, r=0.396, p=0.004), and BMSi (r=-0.613, p<0.001). There were no significant relationships between tissue TMD or AGEs accumulation with the quasi-static material properties. Toughness decreased with increasing tissue Ct.Po. (r=-0.621, p<0.001). Other three-point bending measures also correlated with tissue Ct.Po. including the bending modulus (r=-0.50, p<0.001) and ultimate stress (r=-0.56, p<0.001). The effects of Ct.Po. on indentation were less pronounced with IDI (r=0.290, p=0.043) and BMSi (r=-0.299, p=0.037) correlated modestly with tissue Ct.Po. These results suggest that RPI may be sensitive to bone quality changes relating to collagen.

μFE models can represent microdamaged regions of healthy and metastatically involved whole vertebrae identified through histology and contrast enhanced μCT imaging

Micro-damage formation within the skeleton is an important stimulant for bone remodeling, however abnormal build-up of micro-damage can lead to skeletal fragility. In this study, µCT imaging based micro finite element (μFE) models were used to evaluate tissue level damage criteria in whole healthy and metastatically-involved vertebrae. T13-L2 spinal segments were excised from osteolytic (n=3) and healthy (n=3) female athymic rnu/rnu rats. Osteolytic metastasis was generated by intercardiac injection of HeLa cancer cells. Micro-mechanical axial loading was applied to the spinal motion segments under μCT imaging. Vertebral samples underwent BaSO4 staining and sequential calcein/fuchsin staining to identify load induced micro-damage. μCT imaging was used generate specimen specific μFE models of the healthy and osteolytic whole rat vertebrae. Model boundary conditions were generated through deformable image registration of loaded and unloaded scans. Elevated stresses and strains were detected in regions of micro-damage identified through histological and BaSO4 staining within healthy and osteolytic vertebral models, as compared to undamaged regions. Additionally, damaged regions of metastatic vertebrae experienced significantly higher local stresses and strains than those in the damaged regions of healthy specimens. Areas identified by BaSO4 staining, however, yielded lower levels of stress and strain in damaged and undamaged regions of healthy and metastatic vertebrae as compared to fuschin staining. The multimodal (experimental, image-based and computational) techniques used in this study demonstrated the ability of local stresses and strains computed through µFE analysis to identify trabecular micro-damage, that can be applied to biomechanical analyses of healthy and diseased whole bones.

Evaluating the Effect of Minimizing Screws on Stabilization of Symphysis Mandibular Fracture by 3D Finite Element Analysis

PURPOSE

The objective of this work is to integrate structural optimization and reliability concepts into mini-plate fixation strategy used in symphysis mandibular fractures. The structural reliability levels are next estimated when considering a single failure mode and multiple failure modes.

PATIENTS AND METHODS

A 3-dimensional finite element model is developed in order to evaluate the ability of reducing the negative effect due to the stabilization of the fracture. Topology optimization process is considered in the conceptual design stage to predict possible fixation layouts. In the detailed design stage, suitable mini-plates are selected taking into account the resulting topology and different anatomical considerations. Several muscle forces are considered in order to obtain realistic predictions. Since some muscles can be cut or harmed during the surgery and cannot operate at its maximum capacity, there is a strong motivation to introduce the loading uncertainties in order to obtain reliable designs. The structural reliability is carried out for a single failure mode and multiple failure modes.

RESULTS

The different results are validated with a clinical case of a male patient with symphysis fracture. In this case while use of the upper plate fixation with four holes, only two screws were applied to protect adjacent vital structure. This behavior does not affect the stability of the fracture.

CONCLUSION

The proposed strategy to optimize bone plates leads to fewer complications and second surgeries, less patient discomfort, and shorter time of healing.

Matrix dimensionality and stiffness cooperatively regulate osteogenesis of mesenchymal stromal cells

Osteogenic potential of mesenchymal stromal cells (MSCs) is mechanosensitive. It's affected by the mechanical properties of the cellular microenvironment, particularly its mechanical modulus. To explore the effect of mechanical modulus on osteogenesis in the third dimension (3D), this study used a novel polyacrylamide (PA) scaffold whose pores are monodisperse and spherical, the mechanical moduli of which can be tuned across a wide range. It was found that MSCs have similar proliferation rates in PA scaffolds independent of the matrix stiffness. The contractile force exerted by MSCs inside PA scaffolds was strong enough to deform the pores of scaffolds made of more compliant PAs (whose shear modulus, G'scaffold<4 kPa). Only scaffolds of the highest stiffness (G'scaffold=12 kPa) can withhold the contraction from MSCs. After osteogenic induction for 21 days, the expression profiles of marker genes showed that PA scaffolds of G'scaffold=12 kPa promoted osteogenesis of MSCs. Confocal image analysis demonstrated that there are more F-actin cytoskeletons and bundled stress fibers at higher matrix moduli in 2D and 3D. Moreover, the 3D porous structure promotes osteogenesis of MSCs more than 2D flat substrates. Together, the differences of cellular behaviors when cultured in 2D and 3D systems are evident. The PA scaffolds developed in the present study can be used for further investigation into the mechanism of MSC mechanosensing in the 3D context.

STATEMENT OF SIGNIFICANCE

Mechanical properties of the microenvironment affect cellular behaviors, such as matrix stiffness. Traditionally, cell biological investigations have mostly employed cells growing on 2D substrates. The 3D porous PA scaffolds with the same topological conformation and pore sizes but different stiffness generated in this study showed that the differences of cellular behaviors in 2D and 3D systems are evident. Our 3D scaffolds provide insights into tissue engineering when stem cells incorporated with 3D scaffolds and support the future studies of cellular mechanobiology as well as the elucidation the role mechanical factor plays on the physiology and fate determination of MSCs in the 3D context.

Three-dimensional spherical spatial boundary conditions differentially regulate osteogenic differentiation of mesenchymal stromal cells

The spatial boundary condition (SBC) arising from the surrounding microenvironment imposes specific geometry and spatial constraints that affect organogenesis and tissue homeostasis. Mesenchymal stromal cells (MSCs) sensitively respond to alterations of mechanical cues generated from the SBC. However, mechanical cues provided by a three-dimensional (3D) environment are deprived in a reductionist 2D culture system. This study investigates how SBC affects osteogenic differentiation of MSCs using 3D scaffolds with monodispersed pores and homogenous spherical geometries. MSCs cultured under SBCs with diameters of 100 and 150 μm possessed the greatest capability of osteogenic differentiation. This phenomenon was strongly correlated with MSC morphology, organization of actin cytoskeleton, and distribution of focal adhesion involving α2 and α5 integrins. Further silencing either α2 or α5 integrin significantly reduced the above mentioned mechanosensitivity, indicating that the α2 and α5 integrins as mechano-sensitive molecules mediate MSCs' ability to provide enhanced osteogenic differentiation in response to different spherical SBCs. Taken together, the findings provide new insights regarding how MSCs respond to mechanical cues from the surrounding microenvironment in a spherical SBC, and such biophysical stimuli should be taken into consideration in tissue engineering and regenerative medicine in conjunction with biochemical cues.

A review of trabecular bone functional adaptation: what have we learned from trabecular analyses in extant hominoids and what can we apply to fossils?

Many of the unresolved debates in palaeoanthropology regarding evolution of particular locomotor or manipulative behaviours are founded in differing opinions about the functional significance of the preserved external fossil morphology. However, the plasticity of internal bone morphology, and particularly trabecular bone, allowing it to respond to mechanical loading during life means that it can reveal greater insight into how a bone or joint was used during an individual's lifetime. Analyses of trabecular bone have been commonplace for several decades in a human clinical context. In contrast, the study of trabecular bone as a method for reconstructing joint position, joint loading and ultimately behaviour in extant and fossil non-human primates is comparatively new. Since the initial 2D studies in the late 1970s and 3D analyses in the 1990 s, the utility of trabecular bone to reconstruct behaviour in primates has grown to incorporate experimental studies, expanded taxonomic samples and skeletal elements, and improved methodologies. However, this work, in conjunction with research on humans and non-primate mammals, has also revealed the substantial complexity inherent in making functional inferences from variation in trabecular architecture. This review addresses the current understanding of trabecular bone functional adaptation, how it has been applied to hominoids, as well as other primates and, ultimately, how this can be used to better interpret fossil hominoid and hominin morphology. Because the fossil record constrains us to interpreting function largely from bony morphology alone, and typically from isolated bones, analyses of trabecular structure, ideally in conjunction with that of cortical structure and external morphology, can offer the best resource for reconstructing behaviour in the past.

Topological design and additive manufacturing of porous metals for bone scaffolds and orthopaedic implants: A review

One of the critical issues in orthopaedic regenerative medicine is the design of bone scaffolds and implants that replicate the biomechanical properties of the host bones. Porous metals have found themselves to be suitable candidates for repairing or replacing the damaged bones since their stiffness and porosity can be adjusted on demands. Another advantage of porous metals lies in their open space for the in-growth of bone tissue, hence accelerating the osseointegration process. The fabrication of porous metals has been extensively explored over decades, however only limited controls over the internal architecture can be achieved by the conventional processes. Recent advances in additive manufacturing have provided unprecedented opportunities for producing complex structures to meet the increasing demands for implants with customized mechanical performance. At the same time, topology optimization techniques have been developed to enable the internal architecture of porous metals to be designed to achieve specified mechanical properties at will. Thus implants designed via the topology optimization approach and produced by additive manufacturing are of great interest. This paper reviews the state-of-the-art of topological design and manufacturing processes of various types of porous metals, in particular for titanium alloys, biodegradable metals and shape memory alloys. This review also identifies the limitations of current techniques and addresses the directions for future investigations.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p<0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p≤0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Pressure and shear stress in trabecular bone marrow during whole bone loading

Skeletal adaptation to mechanical loading is controlled by mechanobiological signaling. Osteocytes are highly responsive to applied strains, and are the key mechanosensory cells in bone. However, many cells residing in the marrow also respond to mechanical cues such as hydrostatic pressure and shear stress, and hence could play a role in skeletal adaptation. Trabecular bone encapsulates marrow, forming a poroelastic solid. According to the mechanical theory, deformation of the pores induces motion in the fluid-like marrow, resulting in pressure and velocity gradients. The latter results in shear stress acting between the components of the marrow. To characterize the mechanical environment of trabecular bone marrow in situ, pore pressure within the trabecular compartment of whole porcine femurs was measured with miniature pressure transducers during stress-relaxation and cyclic loading. Pressure gradients ranging from 0.013 to 0.46 kPa/mm were measured during loading. This range was consistent with calculated pressure gradients from continuum scale poroelastic models with the same permeability. Micro-scale computational fluid dynamics models created from computed tomography images were used to calculate the micromechanical stress in the marrow using the measured pressure differentials as boundary conditions. The volume averaged shear stress in the marrow ranged from 1.67 to 24.55 Pa during cyclic loading, which exceeds the mechanostimulatory threshold for mesenchymal lineage cells. Thus, the loading of bone through activities of daily living may be an essential component of bone marrow health and mechanobiology. Additional studies of cell-level interactions during loading in healthy and disease conditions will provide further incite into marrow mechanobiology.