Dynamic analysis of a simplified bone model during the process of fracture healing

The dynamic analysis of fracture healing is tackled numerically by means of a bone model which uses the finite element method. The model is of non-uniform cross-sectional area and moment of inertia. Shear and rotatory inertia are taken into account. Considerable variation of the upper natural frequencies is observed as the healing process progresses. The practical implications, as well as present limitations, of the technique are examined.     

A biomechanical approach to the analysis of methods and procedures of bariatric surgery

Bariatric surgery includes a variety of procedures that are performed on obese people and aim at decreasing the intake of food and calories. This goal is usually pursued by reducing stomach capacity and/or absorbing capability. Adjustable gastric banding is the most common and successful operation. In general, bariatric surgical procedures are effective, but are often associated with major complications.Surgical procedure and post-surgical conformation of the stomach are usually defined on clinical and surgical basis only. Instead, the optimal configuration should be identified by analyzing the mechanical functionality of the stomach and the surrounding structures, and the relationship between food intake, nutrient adsorption, mechanical stimulation of stomach wall and feeling of satiety.A novel approach to bariatric surgery is required, integrating competences in the areas of biomechanics, physiology and surgery, based on a strong interaction between engineers and clinicians. Preliminary results from coupled experimental and computational investigations are here reported. The analyses aim to develop computational tools for the investigation of stomach mechanical functionality in pre- and post-surgical conformations.     

A hybrid bioregulatory model of angiogenesis during bone fracture healing

Bone fracture healing is a complex process in which angiogenesis or the development of a blood vessel network plays a crucial role. In this paper, a mathematical model is presented that simulates the biological aspects of fracture healing including the formation of individual blood vessels. The model consists of partial differential equations, several of which describe the evolution in density of the most important cell types, growth factors, tissues and nutrients. The other equations determine the growth of blood vessels as a result of the movement of leading endothelial (tip) cells. Branching and anastomoses are accounted for in the model. The model is applied to a normal fracture healing case and subjected to a sensitivity analysis. The spatiotemporal evolution of soft tissues and bone, as well as the development of a blood vessel network are corroborated by comparison with experimental data. Moreover, this study shows that the proposed mathematical framework can be a useful tool in the research of impaired healing and the design of treatment strategies.     

Quantitative measures for fracture healing: an in-vitro biomechanical study

A method is presented to assess the functional capabilities of plated osteotomized canine femora as load bearing structures and to quantify their healing status. In this method the osteotomized bones and their intact contralateral controls are tested in nondestructive bending, in twenty-four planes of loading at 15 degree angular increments. Flexural rigidity (EI), in each of the planes, is determined by using classical beam theory. It has been found that the EI values have the expected elliptical distribution. The 24 EI values, obtained for each bone, are curve fitted by regression analysis and the characteristics of each bone are described by the three parameters defining its ellipse. The parameters of the ellipses of the osteotomized bone and of its contralateral control are used to define four new parameters that may serve as measures for the healing efficiency. One of these serves as an indicator for the fragility of the healed bone and the other three add quantitative information on its mechanical state.     

A numerical model of the fracture healing process that describes tissue development and revascularisation

A dynamic model was developed to simulate complex interactions of mechanical stability, revascularisation and tissue differentiation in secondary fracture healing. Unlike previous models, blood perfusion was included as a spatio-temporal state variable to simulate the revascularisation process. A 2D, axisymmetrical finite element model described fracture callus mechanics. Fuzzy logic rules described the following biological processes: angiogenesis, intramembranous ossification, chondrogenesis, cartilage calcification and endochondral ossification, all of which depended on local strain state and local blood perfusion. In order to evaluate how the predicted revascularisation depended on the mechanical environment, we simulated two different healing cases according to two groups of transverse metatarsal osteotomies in sheep with different axial stability. The model predicted slower revascularisation and delayed bony bridging for the less stable case, which corresponded well to the experimental observations. A revascularisation sensitivity analysis demonstrated the potential of the model to account for different conditions regarding the blood supply.     

A numerical model of the fracture healing process that describes tissue development and revascularisation

A dynamic model was developed to simulate complex interactions of mechanical stability, revascularisation and tissue differentiation in secondary fracture healing. Unlike previous models, blood perfusion was included as a spatio-temporal state variable to simulate the revascularisation process. A 2D, axisymmetrical finite element model described fracture callus mechanics. Fuzzy logic rules described the following biological processes: angiogenesis, intramembranous ossification, chondrogenesis, cartilage calcification and endochondral ossification, all of which depended on local strain state and local blood perfusion. In order to evaluate how the predicted revascularisation depended on the mechanical environment, we simulated two different healing cases according to two groups of transverse metatarsal osteotomies in sheep with different axial stability. The model predicted slower revascularisation and delayed bony bridging for the less stable case, which corresponded well to the experimental observations. A revascularisation sensitivity analysis demonstrated the potential of the model to account for different conditions regarding the blood supply.     

Trabecular bone fracture healing simulation with finite element analysis and fuzzy logic

Trabecular bone fractures heal through intramembraneous ossification. This process differs from diaphyseal fracture healing in that the trabecular marrow provides a rich vascular supply to the healing bone, there is very little callus formation, woven bone forms directly without a cartilage intermediary, and the woven bone is remodelled to form trabecular bone. Previous studies have used numerical methods to simulate diaphyseal fracture healing or bone remodelling, however not trabecular fracture healing, which involves both tissue differentiation and trabecular formation. The objective of this study was to determine if intramembraneous bone formation and remodelling during trabecular bone fracture healing could be simulated using the same mechanobiological principles as those proposed for diaphyseal fracture healing. Using finite element analysis and the fuzzy logic for diaphyseal healing, the model simulated formation of woven bone in the fracture gap and subsequent remodelling of the bone to form trabecular bone. We also demonstrated that the trabecular structure is dependent on the applied loading conditions. A single model that can simulate bone healing and remodelling may prove to be a useful tool in predicting musculoskeletal tissue differentiation in different vascular and mechanical environments.     

神经生长因子与骨折愈合

神经生长因子主要由来源于神经嵴的神经元支配的靶组织产生,其被这些神经元轴突摄取后逆行运输至胞体,通过多种途径调节神经细胞的基因转录而发挥生物效应,维持神经元的存活、刺激轴突的生长.并对外周神经的发育、营养起重要的作用.在骨组织和骨折骨痴中均可见神经生长因子及其受体的表达,神经生长因子主要是通过促进骨折部位神经的再生参与骨折修复.骨折愈合的机制十分复杂,神经生长因子对骨组织的作用也是多方面、多层次和相互交叉的,其机制尚未完全明确.虽然神经生长因子促进骨折修复作用机制的研究已经取得一些进展,但仍处于初级阶段,其作用机制仍不明确.     

Analyzing the cellular contribution of bone marrow to fracture healing using bone marrow transplantation in mice

The bone marrow is believed to play important roles during fracture healing such as providing progenitor cells for inflammation, matrix remodeling, and cartilage and bone formation. Given the complex nature of bone repair, it remains difficult to distinguish the contributions of various cell types. Here we describe a mouse model based on bone marrow transplantation and genetic labeling to track cells originating from bone marrow during fracture healing. Following lethal irradiation and engraftment of bone marrow expressing the LacZ transgene constitutively, wild type mice underwent tibial fracture. Donor bone marrow-derived cells, which originated from the hematopoietic compartment, did not participate in the chondrogenic and osteogenic lineages during fracture healing. Instead, the donor bone marrow contributed to inflammatory and bone resorbing cells. This model can be exploited in the future to investigate the role of inflammation and matrix remodeling during bone repair, independent from osteogenesis and chondrogenesis.     

Comparison of microCT and an inverse finite element approach for biomechanical analysis: results in a mesenchymal stem cell therapeutic system for fracture healing

An important concern in the study of fracture healing is the ability to assess mechanical integrity in response to candidate therapeutics in small-animal systems. In recent reports, it has been proposed that microCT image-derived densitometric parameters could be used as a surrogate for mechanical property assessment. Recently, we have proposed an inverse methodology that iteratively reconstructs the modulus of elasticity of the lumped soft callus/hard callus region by integrating both intrinsic mechanical property (from biomechanical testing) and geometrical information (from microCT) within an inverse finite element analysis (FEA) to define a callus quality measure. In this paper, data from a therapeutic system involving mesenchymal stem cells is analyzed within the context of comparing traditional microCT densitometric and mechanical property metrics. In addition, a novel multi-parameter regression microCT parameter is analyzed as well as our inverse FEA metric. The results demonstrate that the inverse FEA approach was the only metric to successfully detect both longitudinal and therapeutic responses. While the most promising microCT-based metrics were adequate at early healing states, they failed to track late-stage mechanical integrity. In addition, our analysis added insight to the role of MSCs by demonstrating accelerated healing and was the only metric to demonstrate therapeutic benefits at late-stage healing. In conclusion, the work presented here indicates that microCT densitometric parameters are an incomplete surrogate for mechanical integrity. Additionally, our inverse FEA approach is shown to be very sensitive and may provide a first-step towards normalizing the often challenging process of assessing mechanical integrity of healing fractures.     

A finite element inverse analysis to assess functional improvement during the fracture healing process

Assessment of the restoration of load-bearing function is the central goal in the study of fracture healing process. During the fracture healing, two critical aspects affect its analysis: (1) material properties of the callus components, and (2) the spatio-temporal architecture of the callus with respect to cartilage and new bone formation. In this study, an inverse problem methodology is used which takes into account both features and yields material property estimates that can analyze the healing changes. Six stabilized fractured mouse tibias are obtained at two time points during the most active phase of the healing process, respectively 10 days (n=3), and 14 days (n=3) after fracture. Under the same displacement conditions, the inverse procedure estimations of the callus material properties are generated and compared to other fracture healing metrics. The FEA estimated property is the only metric shown to be statistically significant (p=0.0194) in detecting the changes in the stiffness that occur during the healing time points. In addition, simulation studies regarding sensitivity to initial guess and noise are presented; as well as the influence of callus architecture on the FEA estimated material property metric. The finite element model inverse analysis developed can be used to determine the effects of genetics or therapeutic manipulations on fracture healing in rodents.     

Lymphoscintigraphic monitoring of the lower limb lymphatic system response to bone fracture and healing

Closed bone fractures, and torn muscles and tendons are "internal wounds". What kind of reaction do they evoke in the local and systemic immune system? Cellular debris of damaged tissue and extravasated blood cells are removed by scavenger cells. They are transported via lymphatics to the lymph nodes. There elimination of self antigens takes place. Clinically, no enlargement of lymph nodes is observed after closed fractures and soft tissue damage. The question arises whether there is really no enlargement of regional lymph nodes, in other words, no reaction to damaged cell antigens. This question was studied by using lymphoscintigraphy to visualize lymphatics and lymph nodes draining the site of closed bone fracture. The lymphoscintigraphic pictures of two groups of patients, those with a rapid noncomplicated healing of leg fractures, and those with protracted healing and undergoing surgical reconstructions, were evaluated. The surface area of lymphatic pathways and inguinal lymph nodes on the injured and contralateral normal limb were measured. Enlarged superficial lymphatics and inguinal lymph nodes were found in limbs with healed bone fractures, and decreased inguinal lymph nodes and visualization of deep lymphatics and popliteal nodes in the majority of patients with nonhealing fractures. There was a lack of correlation between age of patients, duration of healing, and surgical interventions and the lymphoscintigraphic changes. These findings suggest that the fracture gap tissue is a dominant source of signals to the lymph nodes, releasing cellular and humoral regulatory factors. Taken together, there is a strong immune reaction of lymph node to the fracture, although it cannot be recognized clinically.     

The biomechanical analysis of bone strength: A method and its application to platycnemia

Traditional methods of bone analysis (both metric and topographic) are restricted to external characters. Spacial distribution of material is, however, equally critical to an understanding of a bone's function. Dynamic testing to determine whole bone strength can only be performed on fresh specimens. Methods for the calculation of both bending and torsional strength of other specimens (such as preserved or fossil bones) are developed in this paper. In order to illustrate the methods, the functional significance of tibial shaft cross sectional variation is investigated.     

Computational simulation of bone fracture healing under inverse dynamisation

Adaptive finite element models have allowed researchers to test hypothetical relationships between the local mechanical environment and the healing of bone fractures. However, their predictive power has not yet been demonstrated by testing hypotheses ahead of experimental testing. In this study, an established mechano-biological scheme was used in an iterative finite element simulation of sheep tibial osteotomy healing under a hypothetical fixation regime, “inverse dynamisation”. Tissue distributions, interfragmentary movement and stiffness across the fracture site were compared between stiff and flexible fixation conditions and scenarios in which fixation stiffness was increased at a discrete time-point. The modelling work was conducted blind to the experimental study to be published subsequently. The simulations predicted the fastest and most direct healing under constant stiff fixation, and the slowest healing under flexible fixation. Although low fixation stiffness promoted more callus formation prior to bridging, this conferred little additional stiffness to the fracture in the first 5 weeks. Thus, while switching to stiffer fixation facilitated rapid subsequent bridging of the fracture, no advantage of inverse dynamisation could be demonstrated. In vivo data remains necessary to conclusively test this treatment protocol and this will, in turn, provide an evaluation of the model’s performance. The publication of both hypotheses and their computational simulation, prior to experimental testing, offers an appealing means to test the predictive power of mechano-biological models.     

Mechanoregulation modeling of bone healing in realistic fracture geometries

Biomechanics and Modeling in Mechanobiology - In bone fracture healing, new tissue gradually forms, ossifies, and eventually remodels itself to restore mechanical stiffness and strength across...     

Biglycan modulates angiogenesis and bone formation during fracture healing

Matrix proteoglycans such as biglycan (Bgn) dominate skeletal tissue and yet its exact role in regulating bone function is still unclear. In this paper we describe the potential role of (Bgn) in the fracture healing process. We hypothesized that Bgn could regulate fracture healing because of previous work showing that it can affect normal bone formation. To test this hypothesis, we created fractures in femurs of 6-week-old male wild type (WT or Bgn+/0) and Bgn-deficient (Bgn-KO or Bgn-/0) mice using a custom-made standardized fracture device, and analyzed the process of healing over time. The formation of a callus around the fracture site was observed at both 7 and 14 days post-fracture in WT and Bgn-deficient mice and immunohistochemistry revealed that Bgn was highly expressed in the fracture callus of WT mice, localizing within woven bone and cartilage. Micro-computed tomography (μCT) analysis of the region surrounding the fracture line showed that the Bgn-deficient mice had a smaller callus than WT mice. Histology of the same region also showed the presence of less cartilage and woven bone in the Bgn-deficient mice compared to WT mice. Picrosirius red staining of the callus visualized under polarized light showed that there was less fibrillar collagen in the Bgn-deficient mice, a finding confirmed by immunohistochemistry using antibodies to type I collagen. Interestingly, real time RT-PCR of the callus at 7 days post-fracture showed a significant decrease in relative vascular endothelial growth factor A (VEGF) gene expression by Bgn-deficient mice as compared to WT. Moreover, VEGF was shown to bind directly to Bgn through a solid-phase binding assay. The inability of Bgn to directly enhance VEGF-induced signaling suggests that Bgn has a unique role in regulating vessel formation, potentially related to VEGF storage or stabilization in the matrix. Taken together, these results suggest that Bgn has a regulatory role in the process of bone formation during fracture healing, and further, that reduced angiogenesis could be the molecular basis.     

Molecular signaling in bone fracture healing and distraction osteogenesis

The process of fracture healing has been described in detail in many histological studies. Recent work has focused on the mechanisms by which growth and differentiation factors regulate the fracture healing process. Rapid progress in skeletal cellular and molecular biology has led to the identification of many signaling molecules associated with the formation of skeletal tissues, including members of the transforming growth factor-beta (TGF-beta) superfamily and the insulin-like growth factor (IGF) family. Increasing evidence indicates that they are critical regulators of cellular proliferation, differentiation, extracellular matrix biosynthesis and mineralization. Limb lengthening procedure (distraction osteogenesis) is a relevant model to investigate the in vivo correlation between mechanical stimulation and biological responses as the callus is stretched by a proper rate and rhythm of mechanical strain. This model also provides additional insights into the molecular and cellular events during bone fracture repair. TGF-beta 1 was significantly increased in both the distracted callus and the fracture callus. The increased level of TGF-beta 1, together with a low concentration of calcium and an enhanced level of collagen synthesis, was maintained in the distracted callus as long as mechanical strain was applied. Less mineralization is also associated with a low level of osteocalcin production. These observations provide further insights into the molecular basis for the cellular events during distraction osteogenesis.     

A methodological approach for the biomechanical cause analysis of golf-related lumbar spine injuries

A new methodological approach employing mechanical work (MW) determination and relative portion of its elemental analysis was applied to investigate the biomechanical causes of golf-related lumbar spine injuries. Kinematic and kinetic parameters at the lumbar and lower limb joints were measured during downswing in 18 golfers. The MW at the lumbar joint (LJ) was smaller than at the right hip but larger than the MWs at other joints. The contribution of joint angular velocity (JAV) to MW was much greater than that of net muscle moment (NMM) at the LJ, whereas the contribution of NMM to MW was greater rather than or similar to that of JAV at other joints. Thus, the contribution of JAV to MW is likely more critical in terms of the probability of golf-related injury than that of NMM. The MW-based golf-related injury index (MWGII), proposed as the ratio of the contribution of JAV to MW to that of NMM, at the LJ (1.55) was significantly greater than those at other joints ( < 1.05). This generally corresponds to the most frequent occurrence of golf-related injuries around the lumbar spine. Therefore, both MW and MWGII should be considered when investigating the biomechanical causes of lumbar spine injuries.