Gary S. Beaupre

Mechanobiology of skeletal regeneration

Skeletal regeneration is accomplished by a cascade of biologic processes that may include differentiation of pluripotential tissue, endochondral ossification, and bone remodeling. It has been shown that all these processes are influenced strongly by the local tissue mechanical loading history. This article reviews some of the mechanobiologic principles that are thought to guide the differentiation of mesenchymal tissue into bone, cartilage, or fibrous tissue during the initial phase of regeneration. Cyclic motion and the associated shear stresses cause cell proliferation and the production of a large callus in the early phases of fracture healing. For intermittently imposed loading in the regenerating tissue: (1) direct intramembranous bone formation is permitted in areas of low stress and strain; (2) low to moderate magnitudes of tensile strain and hydrostatic tensile stress may stimulate intramembranous ossification; (3) poor vascularity can promote chondrogenesis in an otherwise osteogenic environment; (4) hydrostatic compressive stress is a stimulus for chondrogenesis; (5) high tensile strain is a stimulus for the net production of fibrous tissue; and (6) tensile strain with a superimposed hydrostatic compressive stress will stimulate the development of fibrocartilage. Finite element models are used to show that the patterns of tissue differentiation observed in fracture healing and distraction osteogenesis can be predicted from these fundamental mechanobiologic concepts. In areas of cartilage formation, subsequent endochondral ossification normally will proceed, but it can be inhibited by intermittent hydrostatic compressive stress and accelerated by octahedral shear stress (or strain). Later, bone remodeling at these sites can be expected to follow the same mechanobiologic adaptation rules as normal bone.

Mechanical properties of the human achilles tendon

OBJECTIVE To determine whether the human Achilles tendon has higher material properties than other tendons and to test for strain rate sensitivity of the tendon. DESIGN Mechanical testing of excised tendons. BACKGROUND While the human Achilles tendon appears to experience higher in vivo stresses than other tendons, it is not known how the Achilles tendons material properties compare with the properties of other tendons. METHODS Modulus, failure stress, and failure strain were measured for excised human Achilles tendons loaded at strain rates of 1% s(-1) and 10% s(-1). Paired t-tests were used to examine strain rate effects, and average properties from grouped data were used to compare the Achilles tendons properties with properties reported in the literature for other tendons. RESULTS Failure stress and failure strain were higher at the faster strain rate, but no significant difference in modulus was observed. At the 1% s(-1)rate, the mean modulus and failure stress were 816 MPa (SD, 218) and 71 MPa (SD, 17), respectively. The failure strain was 12.8% (SD, 1.7) for the bone-tendon complex and 7.5% (SD, 1.1) for the tendon substance. At the 10% s(-1) rate, the mean modulus and failure stress were 822 MPa (SD, 211) and 86 MPa (SD, 24), respectively. The mean failure strain was 16.1% (SD, 3.6) for the bone-tendon complex and 9.9% (SD, 1.9) for the tendon substance. These properties fall within the range of properties reported in the literature for other tendons. CONCLUSIONS The material properties of the human Achilles tendon measured in this study are similar to the properties of other tendons reported in the literature despite higher stresses imposed on the Achilles tendon in vivo.

The Treatment Mechanism of an Interspinous Process Implant for Lumbar Neurogenic Intermittent Claudication

Study Design. The spinal canal and neural foramina dimensions of cadaver lumbar spines were quantified during flexion and extension using magnetic resonance imaging before and after placement of an interspinous process implant. Objective. To quantify the effect of the implant on the dimensions of the spinal canal and neural foramina during flexion and extension. Summary of the Background Data. Lumbar neurogenic intermittent claudication symptoms are typically exacerbated during extension and relieved during flexion. It is understood that the dimensions of the spinal canal and neural foramen increase in flexion and decrease in extension. The authors hypothesized that an interspinous process implant would significantly prevent narrowing of the canal and foramina in extension and have no significant effect in flexion. Methods. Eight L2–L5 specimens were positioned to 15° of flexion and 15° of extension using a positioning frame. Each specimen was magnetic resonance imaged with and without an interspinous implant (X STOP) placed between the L3–L4 spinous processes. Canal and foramina dimensions were compared between the intact and implanted specimens using a repeated measures analysis of variance with a level of significance of 0.05. Results. In extension, the implant significantly increased the canal area by 18% (231–273 mm2), the subarticular diameter by 50% (2.5–3.7 mm), the canal diameter by 10% (17.8–19.5 mm), the foraminal area by 25% (106–133 mm2), and the foraminal width by 41% (3.4–4.8 mm). Conclusions. The results of this study show that the X STOP interspinous process implant prevents narrowing of the spinal canal and foramina in extension.

Calcaneal loading during walking and running

PURPOSE This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

Mechanical factors in bone growth and development

Mechanobiologic factors strongly influence skeletal ossification and regulate changes in bone geometry and apparent density during ontogeny. We have developed computer models that implement a simple mathematical rule relating cyclic tissue stresses to bone apposition and resorption. Beginning at the fetal stages of the femoral anlage, these models successfully predict the appositional bone growth and modeling observed in the development of the diaphyseal cross section. The same basic mechanobiologic rule can also predict the architectural construction of the proximal cancellous bone formed in regions of endochondral ossification. Geometry and density changes in adult diaphyseal and cancellous bone as a result of changes in physical activity can be simulated by invoking the same rule used during development. Future clinical and experimental work is needed to provide more quantitative data for mechanobiologic rules and elucidate the interactions between chemical and mechanical factors influencing bone biology.

Adaptive bone remodeling incorporating simultaneous density and anisotropy considerations.

Over 100 years ago, Wolff hypothesized that cancellous bone altered both its apparent density and trabecular orientation in response to mechanical loads. A mathematical counterpart of this principle is derived by adding a remodeling rule for the rate-of-change of the full anisotropic stiffness tensor (all 21 independent terms) to the density rate-of-change rule adapted from an existing isotropic theory. As a result, anisotropy and density patterns develop such that the local stiffness tensor is optimal for the given series of applied loadings. The method does not rely on additional morphological measures of trabecular orientation. Furthermore, assumptions of material symmetry are not required, and any observed regions of orthotropy, transverse isotropy, or isotropy are a result entirely of the functional adaptation of the bone and not the consequence of a modeling assumption. This approach has been implemented with the finite element method and applied to a two-dimensional model of the proximal femur with encouraging results.

The influence of bone volume fraction and ash fraction on bone strength and modulus

Although bone strength and modulus are known to be influenced by both volume fraction and mineral content (ash fraction), the relative influence of these two parameters remains unknown. Single-parameter power law functions are used widely to relate bone volume or ash fraction to bone strength and elastic modulus. In this study we evaluate the potential for predicting bone mechanical properties with two-parameter power law functions of bone volume fraction (BV/TV) and ash fraction (alpha) of the form y = a(BV/TV)(b) alpha(c) (where y is either ultimate strength or elastic modulus). We derived an expression for bone volume fraction as a function of apparent density and ash fraction to perform a new analysis of data presented by Keller in 1994. Exponents b and c for the prediction of bone strength were found to be 1.92 +/- 0.02 and 2.79 +/- 0.09 (mean +/- SE), respectively, with r(2) = 0.97. The value of b was found to be consistent with that found previously, whereas the value of c was lower than values previously reported. For the prediction of elastic modulus we found b and c to be 2.58 +/- 0.02 and 2.74 +/- 0.13, respectively, with r(2) = 0.97. The exponent related to ash fraction was typically larger than that associated with bone volume fraction, suggesting that a change in mineral content will, in general, generate a larger change in bone strength and stiffness than a similar change in bone volume fraction. These findings are important for interpreting the results of antiresorptive drug treatments that can cause changes in both ash and bone volume fraction.

Mechanobiologic influences in long bone cross-sectional growth

We developed a computer model to simulate the interaction of biological and mechanobiological factors in the development of the cross-sectional morphology of long bones. The model incorporated a strong influence of biologically induced bone formation during early development. In addition, an assumed mechanical loading history during growth and development corresponding to age-related changes in body weight and muscle mass was applied. Based on the bone stress stimulus generated by the assumed loads, mechanically induced apposition and resorption rates were calculated at the periosteal and endosteal surfaces using a previously developed bone modeling theory. These methods successfully emulated the growth-related changes seen in long bone diaphyseal structure as well as changes observed in mature bones during aging. The simulations recreated the rapid increase in bone dimensions during development, stabilizing at maturity, and then the gradual, age-related subperiosteal expansion and cortical thinning. Throughout the growth, development, and aging simulations, the values of the bone radii, area, moments of inertia, and apposition rates corresponded well with measurements documented by other researchers.

Effects of Creep and Cyclic Loading on the Mechanical Properties and Failure of Human Achilles Tendons

AbstractThe Achilles tendon is one of the most frequently injured tendons in humans, and yet the mechanisms underlying its injury are not well understood. This study examines the ex vivo mechanical behavior of excised human Achilles tendons to elucidate the relationships between mechanical loading and Achilles tendon injury. Eighteen tendons underwent creep testing at constant stresses from 35 to 75 MPa. Another 25 tendons underwent sinusoidal cyclic loading at 1 Hz between a minimum stress of 10 MPa and maximum stresses of 30–80 MPa. For the creep specimens, there was no significant relationship between applied stress and time to failure, but time to failure decreased exponentially with increasing initial strain (strain when target stress is first reached) and decreasing failure strain. For the cyclically loaded specimens, secant modulus decreased and cyclic energy dissipation increased over time. Time and cycles to failure decreased exponentially with increasing applied stress, increasing initial strain (peak strain from first loading cycle), and decreasing failure strain. For both creep and cyclic loading, initial strain was the best predictor of time or cycles to failure, supporting the hypothesis that strain is the primary mechanical parameter governing tendon damage accumulation and injury. The cyclically loaded specimens failed faster than would be expected if only time-dependent damage occurred, suggesting that repetitive loading also contributes to Achilles tendon injuries.

Knee muscle forces during walking and running in patellofemoral pain patients and pain-free controls.

One proposed mechanism of patellofemoral pain, increased stress in the joint, is dependent on forces generated by the quadriceps muscles. Describing causal relationships between muscle forces, tissue stresses, and pain is difficult due to the inability to directly measure these variables in vivo. The purpose of this study was to estimate quadriceps forces during walking and running in a group of male and female patients with patellofemoral pain (n = 27, 16 female; 11 male) and compare these to pain-free controls (n = 16, 8 female; 8 male). Subjects walked and ran at self-selected speeds in a gait laboratory. Lower limb kinematics and electromyography (EMG) data were input to an EMG-driven musculoskeletal model of the knee, which was scaled and calibrated to each individual to estimate forces in 10 muscles surrounding the joint. Compared to controls, the patellofemoral pain group had greater co-contraction of quadriceps and hamstrings (p = 0.025) and greater normalized muscle forces during walking, even though the net knee moment was similar between groups. Muscle forces during running were similar between groups, but the net knee extension moment was less in the patellofemoral pain group compared to controls. Females displayed 30-50% greater normalized hamstring and gastrocnemius muscle forces during both walking and running compared to males (p<0.05). These results suggest that some patellofemoral pain patients might experience greater joint contact forces and joint stresses than pain-free subjects. The muscle force data are available as supplementary material.