Computational modeling of human bone fracture healing affected by different conditions of initial healing stage

Addressing the challenges of infectious bone defects: a review of recent advances in bifunctional biomaterials

Infectious bone defects present a substantial clinical challenge due to the complex interplay between infection control and bone regeneration. These defects often result from trauma, autoimmune diseases, infections, or tumors, requiring a nuanced approach that simultaneously addresses infection and promotes tissue repair. Recent advances in tissue engineering and materials science, particularly in nanomaterials and nano-drug formulations, have led to the development of bifunctional biomaterials with combined osteogenic and antibacterial properties. These materials offer an alternative to traditional bone grafts, minimizing complications such as multiple surgeries, high antibiotic dosages, and lengthy recovery periods. This review examines the repair mechanisms in the infectious microenvironment and highlights various bifunctional biomaterials that foster both anti-infective and osteogenic processes. Emerging design strategies are also discussed to provide a forward-looking perspective on treating infectious bone defects with clinically significant outcomes.

The effect of early weight-bearing on bone fusion after triple arthrodesis

Triple arthrodesis is an effective method for treating stiff horseshoe feet and severe osteoarthritis. However, it is still a challenge to improve postoperative bone fusion by changing early weight-bearing. This study improved the classical bone remodeling algorithm, established a mathematical relationship between density change rate and mechanical stimulation, and combined it with finite element theory. The proposed algorithm can not only predict the effect of early weight-bearing on triple arthrodesis but also visually demonstrate the change of bone mineral density with time. The analysis results indicated that 2.5% of the initial load was a potential factor leading to bone nonunion, and 50% of the initial load would result in bone resorption. Meanwhile, it was found that 25% of the external load was more conducive to postoperative rehabilitation. The study results have theoretical significance for enhancing the effect of postoperative bone fusion and formulating a more scientific rehabilitation program, thereby supporting patients' postoperative rehabilitation exercise.

Don't mind the gap: reframing the Perren strain rule for fracture healing using insights from virtual mechanical testing

Aims

The "2 to 10% strain rule" for fracture healing has been widely interpreted to mean that interfragmentary strain greater than 10% predisposes a fracture to nonunion. This interpretation focuses on the gap-closing strain (axial micromotion divided by gap size), ignoring the region around the gap where osteogenesis typically initiates. The aim of this study was to measure gap-closing and 3D interfragmentary strains in plated ovine osteotomies and associate local strain conditions with callus mineralization.

Methods

MicroCT scans of eight female sheep with plated mid-shaft tibial osteotomies were used to create image-based finite element models. Virtual mechanical testing was used to compute postoperative gap-closing and 3D continuum strains representing compression (volumetric strain) and shear deformation (distortional strain). Callus mineralization was measured in zones in and around the osteotomy gap.

Results

Gap-closing strains averaged 51% (mean) at the far cortex. Peak compressive volumetric strain averaged 32% and only a small tissue volume (average 0.3 cm3) within the gap experienced compressive strains > 10%. Distortional strains were much higher and more widespread, peaking at a mean of 115%, with a mean of 3.3 cm3 of tissue in and around the osteotomy experiencing distortional strains > 10%. Callus mineralization initiated outside the high-strain gap and was significantly lower within the fracture gap compared to around it at nine weeks.

Conclusion

Ovine osteotomies can heal with high gap strains (> 10%) dominated by shear conditions. High gap strain appears to be a transient local limiter of osteogenesis, not a global inhibitor of secondary fracture repair.

A multi-domain computational framework investigating the short- and long-term viability of bioabsorbable magnesium fixation for tibial fractures

This study presents a multi-domain computational framework to investigate the long-term performance of permanent and bioabsorbable magnesium fixation devices in orthopaedic fracture applications. The framework integrates a coupled model for bone fracture healing and remodeling, with an enhanced surface-based corrosion model to predict the performance of bioabsorbable magnesium devices. It was found that plated fracture fixation enabled fracture healing outcomes compare to non-plated models by facilitating direct fracture healing. During the fracture healing phase, it was found that the stiff titanium plate provided a better healing response compared to the less stiff bioabsorbable magnesium plates. However, in the longer-term remodeling phase, the titanium plate showed evidence of stress-shielding and inhibited bone remodeling. On the other hand, the magnesium plates showed that there was continued remodeling, which meant that the bone tissue gradually returned to the pre-fracture stress state. While the corrosion rate and pit severity heavily influenced the mechanical support provided by the corroding magnesium fixator, the results showed that fixation was only required to provide mechanical stability to the fracture region for approximately the first 30 days for successful fracture union to occur. This coupled computational framework provides a platform to investigate the role of a wide range of magnesium fixation devices and their design and optimisation in orthopaedic applications.

Impacts of post-operation loading and fixation implant on the healing process of fractured tibia

Healing of tibia demonstrates a complex mechanobiological process as it is stimulated by the major factor of strains applied by body weight. The effect of screw heads and bodies as well as their pressure distribution is often overlooked. Hence, effective mechanical conditions of the healing process of tibia can be categorized into the material of the plate and screws, post-operation loadings, and screw type and pressure. In this paper, a mathematical biodegradation model was used to simulate the PGF/PLA plate-screw device over 8 weeks. The effect of different post-operation loading patterns was studied for both locking and non-locking screws. The aim was to reach the best configuration for the most achievable healing using FEA by computing the healing pattern, trend, and efficiency with the mechano-regulation theory based on deviatoric strain. The biodegradation process of the plate and screws resulted in 82% molecular weight loss and 1.05 GPa decrease in Young's modulus during 8 weeks. The healing efficiency of the cases ranged from 4.72% to 14.75% in the first week and 18.64% to 63.05% in the eighth week. Finally, an optimal case was achieved by considering the prevention of muscle erosion, bone density reduction, and nonunion, according to the obtained results.

A novel strain-based bone-fracture healing algorithm is able to predict a range of healing outcomes

Fracture healing is a complex process which sometimes results in non-unions, leading to prolonged disability and high morbidity. Traditional methods of optimising fracture treatments, such as in vitro benchtop testing and in vivo randomised controlled trials, face limitations, particularly in evaluating the entire healing process. This study introduces a novel, strain-based fracture-healing algorithm designed to predict a wide range of healing outcomes, including both successful unions and non-unions. The algorithm uses principal strains as mechanical stimuli to simulate fracture healing in response to local mechanical environments within the callus region. The model demonstrates good agreement with experimental data from ovine metatarsal osteotomies across six fracture cases with varying gap widths and inter-fragmentary strains, replicates physiological bony growth patterns, and is independent of the initial callus geometry. This computational approach provides a framework for developing new fracture-fixation devices, aid in pre-surgical planning, and optimise rehabilitation strategies.

The Advantages and Shortcomings of Stem Cell Therapy for Enhanced Bone Healing

This review explores the regenerative potential of key progenitor cell types and therapeutic strategies to improve healing of complex fractures and bone defects. We define, summarize, and discuss the differentiation potential of totipotent, pluripotent, and multipotent stem cells, emphasizing the advantages and shortcomings of cell therapy for bone repair and regeneration. The fundamental role of mesenchymal stem cells is highlighted due to their multipotency to differentiate into the key lineage cells including osteoblasts, osteocytes, and chondrocytes, which are crucial for bone formation and remodeling. Hematopoietic stem cells (HSCs) also play a significant role; immune cells such as macrophages and T-cells modulate inflammation and tissue repair. Osteoclasts are multinucleated cells that are important to bone remodeling. Vascular progenitor (VP) cells are critical to oxygen and nutrient supply. The dynamic interplay among these lineages and their microenvironment is essential for effective bone restoration. Therapies involving cells that are more than "minimally manipulated" are controversial and include embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). ESCs, derived from early-stage embryos, possess pluripotent capabilities and have shown promise in preclinical studies for bone healing. iPSCs, reprogrammed from somatic cells, offer personalized medicine applications and can differentiate into various tissue-specific cell lines. Minimally manipulative cell therapy approaches such as the use of bone marrow aspirate concentrate (BMAC), exosomes, and various biomaterials for local delivery are explored for their effectiveness in bone regeneration. BMAC, which contains mostly immune cells but few mesenchymal and VPs, probably improves bone healing by facilitating paracrine-mediated intercellular communication. Exosome isolation harnesses the biological signals and cellular by-products that are a primary source for cell crosstalk and activation. Safe, efficacious, and cost-effective strategies to enhance bone healing using novel cellular therapies are part of a changing paradigm to modulate the inflammatory, repair, and regenerative pathways to achieve earlier more robust tissue healing and improved physical function. Impact Statement Stem cell therapy holds immense potential for bone healing due to its ability to regenerate damaged tissue. Nonmanipulated bone marrow aspirate contains mesenchymal stem cells that promote bone repair and reduce healing time. Induced pluripotent stem cells offer the advantage of creating patient-specific cells that can differentiate into osteoblasts, aiding in bone regeneration. Other delivery methods, such as scaffold-based techniques, enhance stem cell integration and function. Collectively, these approaches can improve treatment outcomes, reduce recovery periods, and advance our understanding of bone healing mechanisms, making them pivotal in orthopedic research and regenerative medicine.

Gelatin Methacryloyl (GelMA) - 45S5 Bioactive Glass (BG) Composites for Bone Tissue Engineering: 3D Extrusion Printability and Cytocompatibility Assessment Using Human Osteoblasts

3D extrusion printing has been widely investigated for low-volume production of complex-shaped scaffolds for tissue regeneration. Gelatin methacryloyl (GelMA) is used as a baseline material for the synthesis of biomaterial inks, often with organic/inorganic fillers, to obtain a balance between good printability and biophysical properties. The present study demonstrates how 45S5 bioactive glass (BG) addition and GelMA concentrations can be tailored to develop GelMA composite scaffolds with good printability and buildability. The experimental results suggest that 45S5 BG addition consistently decreases the compression stiffness, irrespective of GelMA concentration, albeit within 20% of the baseline scaffold (without 45S5 BG). The optimal addition of 2 wt % 45S5 BG in 7.5 wt % GelMA was demonstrated to provide the best combination of printability and buildability in the 3D extrusion printing route. The degradation decreases and the swelling kinetics increases with 45S5 BG addition, irrespective of GelMA concentration. Importantly, the dissolution in simulated body fluid over 3 weeks clearly promoted the nucleation and growth of crystalline calcium phosphate particles, indicating the potential of GelMA-45S5 BG to promote biomineralization. The cytocompatibility assessment using human osteoblasts could demonstrate uncompromised cell proliferation or osteogenic marker expression over 21 days in culture for 3D printable 7.5 wt % GelMA -2 wt % 45S5 BG scaffolds when compared to 7.5 wt % GelMA. The results thus encourage further investigations of the GelMA/45S5 BG composite system for bone tissue engineering applications.

Case Studies of a Simulation Workflow to Improve Bone Healing Assessment in Impending Non-Unions

Background: The healing potential of a fracture is determined by mechanical and biological factors. Simulation-based workflows can help assess these factors to assist in predicting non-unions. The aim of this study was the introduction of two use cases for a novel patient-specific simulation workflow based on clinically available information. Methods: The used software is an extension of the "Ulm Bone Healing model" and was applied in two cases with non-union development after fracture fixation to show its principal feasibility. The clinical and radiographic information, starting from initial treatment, were used to feed the simulation process. Results: The simulation predicted non-union development and axial deviation in a mechanically driven non-union. In the case of a biological non-union, a slow, incomplete healing course was correctly identified. However, the time offset in callus bridging was discordant between the simulation and the distinctly slower healing response in the clinical case. Conclusions: The simulation workflow presented in the two clinical use cases allowed for the identification of fractures at risk for impending non-union immediately after the initial fixation based on available clinical and radiographic information. Further validation in a large non-union cohort is needed to increase the model's precision, especially in biologically challenging cases, and show its validity as a screening instrument.

Computational Modeling of Bone Fracture Healing Under Different Initial Conditions and Mechanical Load

OBJECTIVE

Computational model of bone healing can replace animal experiments to study the parameters affecting the bone healing process, thus reducing the damage to experimental animals and saving a lot of time. We propose a computational model for continuous simulation of four phases of bone healing to study the effects of mechanical environmental and biological factors, including initial conditions at the fracture site, mechanical stimulus loading, and vascular growth rate.

METHODS

A finite element model of mechanobiological fracture healing containing several pre-determined variables was developed for bone healing after fracture in sheep, which included many relevant parameters and biological effects during fracture healing, such as the effects of mechanical environment, blood supply level in the local fracture area, cell migration and diffusion, and resorption effects of fracture healing. The effects of several parameters on indices such as Young's modulus of the callus during bone healing were obtained by simulation.

RESULTS

The initial geometry of the healing tissue and mechanical loading had the greatest effect on fracture healing, and different preset values were likely to cause delayed or non-healing fractures. Changed initial tissue properties of the healing tissue showed a nonlinear effect on fracture healing rather than a linear delay or advancement. Parameters related to angiogenesis had a greater effect on fracture healing compared to those related to cell migration.

CONCLUSION

This paper quantified the effect of fracture healing pre-determined variables on fracture healing to better understand the application of mechanobiology in fracture healing simulation models and optimization of treatment strategies.

SIGNIFICANCE

The importance of initial conditions and loads on fracture healing has been shown to help physicians treat bone nonunion or delayed bone healing after a fracture.

Computational models of bone fracture healing and applications: a review

Fracture healing is a very complex physiological process involving multiple events at different temporal and spatial scales, such as cell migration and tissue differentiation, in which mechanical stimuli and biochemical factors assume key roles. With the continuous improvement of computer technology in recent years, computer models have provided excellent solutions for studying the complex process of bone healing. These models not only provide profound insights into the mechanisms of fracture healing, but also have important implications for clinical treatment strategies. In this review, we first provide an overview of research in the field of computational models of fracture healing based on CiteSpace software, followed by a summary of recent advances, and a discussion of the limitations of these models and future directions for improvement. Finally, we provide a systematic summary of the application of computational models of fracture healing in three areas: bone tissue engineering, fixator optimization and clinical treatment strategies. The application of computational models of bone healing in clinical treatment is immature, but an inevitable trend, and as these models become more refined, their role in guiding clinical treatment will become more prominent.

Finite element method-based simulation on bone fracture fixation configuration factors for biodegradable embossed locking compression plate

As an evolution, biodegradable implants need to maximize mechanical performance thereby may lead to confusion in selection of the biodegradable material and implant design to the fracture site. This requires selecting a unique fixation configuration to fit within the fractured bone, factors of which can be bone-plate clearance, interfragmentary gap, alteration in screw fixation position and variation in the number of screws whose configuration optimization can re-maximize the mechanical performance of the biodegradable implant. Therefore, these factors have been optimized based on the induced minimum stress using the finite element method-based simulation for which biodegradable embossed locking plates (BELCP) via screws made of Mg-alloy have been fitted over two fragments of femur body (as hollow cylindrical cortical bone). An average human weight of 62 kg is applied to one segment of the femur for all different configurations of each factor, where another segment is assumed to be fixed. By this simulation, the most optimal fixation configuration was found at a minimum induced stress value of 41.96 MPa which is approximately 85%, 18%, 6% and 48% respectively less than all maximum stress induced configurations in each of the factor. This optimized configuration was at the minimum clearance between bone and plate with a 3 mm interfragmentary gap using 8 screws where the locking screw begins to apply from the center of the BELCP. Overall, BELCP may be a better biodegradable implant plate for bone fracture fixation with these optimized fixation configurations as the improved mechanical performance after experimental validation.

A Preclinical Trial Protocol Using an Ovine Model to Assess Scaffold Implant Biomaterials for Repair of Critical-Sized Mandibular Defects

The present work describes a preclinical trial (in silico, in vivo and in vitro) protocol to assess the biomechanical performance and osteogenic capability of 3D-printed polymeric scaffolds implants used to repair partial defects in a sheep mandible. The protocol spans multiple steps of the medical device development pipeline, including initial concept design of the scaffold implant, digital twin in silico finite element modeling, manufacturing of the device prototype, in vivo device implantation, and in vitro laboratory mechanical testing. First, a patient-specific one-body scaffold implant used for reconstructing a critical-sized defect along the lower border of the sheep mandible ramus was designed using on computed-tomographic (CT) imagery and computer-aided design software. Next, the biomechanical performance of the implant was predicted numerically by simulating physiological load conditions in a digital twin in silico finite element model of the sheep mandible. This allowed for possible redesigning of the implant prior to commencing in vivo experimentation. Then, two types of polymeric biomaterials were used to manufacture the mandibular scaffold implants: poly ether ether ketone (PEEK) and poly ether ketone (PEK) printed with fused deposition modeling (FDM) and selective laser sintering (SLS), respectively. Then, after being implanted for 13 weeks in vivo, the implant and surrounding bone tissue was harvested and microCT scanned to visualize and quantify neo-tissue formation in the porous space of the scaffold. Finally, the implant and local bone tissue was assessed by in vitro laboratory mechanical testing to quantify the osteointegration. The protocol consists of six component procedures: (i) scaffold design and finite element analysis to predict its biomechanical response, (ii) scaffold fabrication with FDM and SLS 3D printing, (iii) surface treatment of the scaffold with plasma immersion ion implantation (PIII) techniques, (iv) ovine mandibular implantation, (v) postoperative sheep recovery, euthanasia, and harvesting of the scaffold and surrounding host bone, microCT scanning, and (vi) in vitro laboratory mechanical tests of the harvested scaffolds. The results of microCT imagery and 3-point mechanical bend testing demonstrate that PIII-SLS-PEK is a promising biomaterial for the manufacturing of scaffold implants to enhance the bone-scaffold contact and bone ingrowth in porous scaffold implants. MicroCT images of the harvested implant and surrounding bone tissue showed encouraging new bone growth at the scaffold-bone interface and inside the porous network of the lattice structure of the SLS-PEK scaffolds.

Temporal dynamics of immune-stromal cell interactions in fracture healing

Bone fracture repair is a complex, multi-step process that involves communication between immune and stromal cells to coordinate the repair and regeneration of damaged tissue. In the US, 10% of all bone fractures do not heal properly without intervention, resulting in non-union. Complications from non-union fractures are physically and financially debilitating. We now appreciate the important role that immune cells play in tissue repair, and the necessity of the inflammatory response in initiating healing after skeletal trauma. The temporal dynamics of immune and stromal cell populations have been well characterized across the stages of fracture healing. Recent studies have begun to untangle the intricate mechanisms driving the immune response during normal or atypical, delayed healing. Various in vivo models of fracture healing, including genetic knockouts, as well as in vitro models of the fracture callus, have been implemented to enable experimental manipulation of the heterogeneous cellular environment. The goals of this review are to (1): summarize our current understanding of immune cell involvement in fracture healing (2); describe state-of-the art approaches to study inflammatory cells in fracture healing, including computational and in vitro models; and (3) identify gaps in our knowledge concerning immune-stromal crosstalk during bone healing.

A Murine Delayed-Healing Model Associates Immune Response with Functional Bone Regeneration after Trauma

Nonunion and delayed-union fractures pose a significant clinical challenge, often leading to prolonged morbidity and impaired quality of life. Fracture-induced hematoma and acute inflammation are crucial for establishing the healing cascade. However, aberrant inflammatory phenotypes can suppress healing and cause bone resorption. Elucidating these mechanisms is necessary to develop potent immunomodulatory therapies and prevent nonunion. Here, we report a delayed fracture healing model enabling the modulation of interfragmentary strain that mimics the etiology of hypertrophic nonunions to elucidate the role of dysregulated immune response in poor healing outcomes. High interfragmentary strain (>15%) was associated with larger callus volumes but delayed bone healing, increased inflammation, and inferior healing outcomes, while lower strain levels (<5%) corresponded to normal bone healing. In addition, we found distinct differences in the ossification, chondrification, and fibrosis patterns between high and low-strain groups, underscoring the significant impact of strain on the healing process. A comprehensive analysis of the systemic immune response revealed dynamic alterations in immune cell populations and factors, particularly within the early hours and days post-fracture. Several immune factors exhibited significant correlations with various functional healing outcomes, indicating their potential as predictive markers for assessing fracture healing progression. Our results also highlighted the significance of timely resolution of proinflammatory signals and the elevation of pro-regenerative immune cell phenotypes in promoting bone regeneration. Multivariate analysis revealed that CD25+ T-regulatory cells were influential in predicting proper bone healing, followed by CD206+ macrophages, underscoring the pivotal role of immune cell populations in the bone healing process. In conclusion, our study provides valuable insights into the intricate interplay between interfragmentary strain, immune response, and the ultimate outcomes of fracture healing. By shedding light on the underlying mechanisms that drive hypertrophic nonunion pathogenesis, our research lays the foundation for enhanced surgical management of nonunions and offers a promising avenue for developing targeted therapeutic interventions and personalized treatment strategies for individuals suffering from fracture nonunion.

Substrate Stiffness of Bone Microenvironment Controls Functions of Pre-Osteoblasts and Fibroblasts In Vitro

The formation of bone in a bone defect is accomplished by osteoblasts, while the over activation of fibroblasts promotes fibrosis. However, it is not clear how the extracellular matrix stiffness of the bone-regeneration microenvironment affects the function of osteoblasts and fibroblasts. This study aim to investigate the effect of bone-regeneration microenvironment stiffness on cell adhesion, cell proliferation, cell differentiation, synthesizing matrix ability and its potential mechanisms in mechanotransduction, in pre-osteoblasts and fibroblasts. Polyacrylamide substrates mimicking the matrix stiffness of different stages of the bone-healing process (15 kPa, mimic granulation tissue; 35 kPa, mimic osteoid; 150 kPa, mimic calcified bone matrix) were prepared. Mouse pre-osteoblasts MC3T3-E1 and mouse fibroblasts NIH3T3 were plated on three types of substrates, respectively. There were significant differences in the adhesion of pre-osteoblasts and fibroblasts on different polyacrylamide substrates. Runx2 expression increased with increasing substrate stiffness in pre-osteoblasts, while no statistical differences were found in the Acta2 expression in fibroblasts on three substrates. OPN expression in pre-osteoblasts, as well as Fn1 and Col1a1 expression in fibroblasts, decreased with increasing stiffness. The difference between the cell traction force generated by pre-osteoblasts and fibroblasts on substrates was also found. Our results indicated that substrate stiffness is a potent regulator of pre-osteoblasts and fibroblasts with the ability of promoting osteogenic differentiation of pre-osteoblasts, while having no effect on myofibroblast differentiation of fibroblasts.

A coupled computational framework for bone fracture healing and long-term remodelling: Investigating the role of internal fixation on bone fractures

In this study, a coupled computational modelling framework for bone fracture repair is presented that enables predictions of both healing and remodelling phases of the fracture region and is used to investigate the role of an internal fixation plate on the long-term healing performance of a fracture tibia under a range of different conditions. It was found that introduction of a titanium plate allowed the tibia to undergo successful healing at higher loading conditions and fracture gaps, compared with the non-plated versions. While these plated cases showed faster rates of repair in the healing phase, their performance was substantially different once they entered the remodelling phase, with substantial regions of stress shielding predicted. This framework is one of the few implementations of both fracture healing and remodelling phases of bone repair and includes several innovative approaches to smoothing, time-averaging and time incrementation in its implementation, thereby avoiding any unwanted abrupt changes between tissue phenotypes. This provides a better representation of tissue development in the fracture site when compared with fracture healing models alone and provides a suitable platform to investigate the long-term performance of orthopaedic fixation devices. This would enable the more effective design of permanent fixation devices and optimisation of the spatial and temporal performance of bioabsorbable implants.

Tissue-Level Regeneration and Remodeling Dynamics are Driven by Mechanical Stimuli in the Microenvironment in a Post-Bridging Loaded Femur Defect Healing Model in Mice

Bone healing and remodeling are mechanically driven processes. While the generalized response to mechanical stimulation in bone is well-understood, much less is known about the mechanobiology-regulating tissue-scale bone formation and resorption during the reparative and remodeling phases of fracture healing. In this study, we combined computational approaches in the form of finite element analysis and experimental approaches by using a loaded femoral defect model in mice to investigate the role of mechanical stimulation in the microenvironment of bone. Specifically, we used longitudinal micro-computed tomography to observe temporal changes in bone at different densities and micro-finite element analysis to map the mechanics of the microenvironment to tissue-scale formation, quiescence (no change in bone presence between time points), and resorption dynamics in the late reparative and remodeling phases (post bridging). Increasing levels of effective strain led to increasing conditional probability of bone formation, while decreasing levels of effective strain led to increasing probability of bone resorption. In addition, the analysis of mineralization dynamics showed both a temporal and effective strain level-dependent behavior. A logarithmic-like response was displayed, where the conditional probability of bone formation or resorption increased rapidly and plateaued or fell rapidly and plateaued as mechanical strain increased.

Effect of Intramedullary Nailing Patterns on Interfragmentary Strain in a Mouse Femur Fracture: A Parametric Finite Element Analysis

Delayed long bone fracture healing and nonunion continue to be a significant socioeconomic burden. While mechanical stimulation is known to be an important determinant of the bone repair process, understanding how the magnitude, mode, and commencement of interfragmentary strain (IFS) affect fracture healing can guide new therapeutic strategies to prevent delayed healing or nonunion. Mouse models provide a means to investigate the molecular and cellular aspects of fracture repair, yet there is only one commercially available, clinically-relevant, locking intramedullary nail (IMN) currently available for studying long bone fractures in rodents. Having access to alternative IMNs would allow a variety of mechanical environments at the fracture site to be evaluated, and the purpose of this proof-of-concept finite element analysis study is to identify which IMN design parameters have the largest impact on IFS in a murine transverse femoral osteotomy model. Using the dimensions of the clinically relevant IMN as a guide, the nail material, distance between interlocking screws, and clearance between the nail and endosteal surface were varied between simulations. Of these parameters, changing the nail material from stainless steel (SS) to polyetheretherketone (PEEK) had the largest impact on IFS. Reducing the distance between the proximal and distal interlocking screws substantially affected IFS only when nail modulus was low. Therefore, IMNs with low modulus (e.g., PEEK) can be used alongside commercially available SS nails to investigate the effect of initial IFS or stability on fracture healing with respect to different biological conditions of repair in rodents.

Direct electromagnetic coupling to determine diagnostic bone fracture stiffness

Background

Rapid prediction of adverse bone fracture healing outcome (e.g., nonunion and/or delayed union) is essential to advise adjunct therapies to reduce patient suffering and improving healing outcome. Radiographic diagnostic methods remain ineffective during early healing, resulting in average nonunion diagnosis times surpassing six months. To address this clinical deficit, we developed a novel diagnostic device to predict fracture healing outcome by noninvasive telemetric measurements of fracture bending stiffness. This study evaluated the hypothesis that our diagnostic antenna system is capable of accurately measuring temporal fracture healing stiffness, and advises the utility of this data for expedited prediction of healing outcomes during early (≤3 weeks) fracture recovery.

Methods

Fracture repair was simulated, in reverse chronology, by progressively destabilizing cadaveric ovine metatarsals (n=8) stabilized via locking plate fixation. Bending stiffness of each fracture state were predicted using a novel direct electromagnetic coupling diagnostic system, and results were compared to values from material testing (MT) methods. While direct calculation of fracture stiffness in a simplistic cadaver model is possible, comparable analysis of the innumerable permutations of fracture and treatment type is not feasible. Thus, clinical feasibility of direct electromagnetic coupling was explored by parametric finite element (FE) analyses (n=1,632 simulations). Implant mechanics were simulated throughout the course of healing for cases with variations to fracture size, implant type, implant structure, and implant material.

Results

For all fracture states, stiffness values predicted by the direct electromagnetic coupling system were not significantly different than those quantified by in vitro MT methods [P=0.587, P=0.985, P=0.975; for comparing intact, destabilized, and fully fractured (FF) states; respectively]. In comparable models, the total implant deflection reduction (from FF to intact states) was less than 10% different between direct electromagnetic coupling measurements (82.2 µm) and FE predictions (74.7 µm). For all treatment parameters, FE analyses predicted nonlinear reduction in bending induced implant midspan deflections for increasing callus stiffness.

Conclusions

This technology demonstrates potential as a noninvasive clinical tool to accurately quantify healing fracture stiffness to augment and expedite healing outcome predictions made using radiographic imaging.

Evaluation of Bone Consolidation in External Fixation with an Electromechanical System

The monitoring of fracture or osteotomy healing is vital for orthopedists to help advise, if necessary, secondary treatments for improving healing outcomes and minimizing patient suffering. It has been decades since osteotomy stiffness has been identified as one main parameter to quantify and qualify the outcome of a regenerated callus. Still, radiographic imaging remains the current standard diagnostic technique of orthopedists. Hence, with recent technological advancements, engineers need to use the new branches of knowledge and improve or innovate diagnostic technologies. An electromechanical system was developed to help diagnose changes in osteotomy stiffness treated with the external fixator LRS Orthofix®. The concept was evaluated experimentally and numerically during fracture healing simulation using two different models: a simplified model of a human tibia, consisting of a nylon bar with a diameter of 30 mm, and a synthetic tibia with the anatomical model from fourth-generation Sawbones®. Moreover, Sawbones® blocks with different densities simulated the mechanical characteristics of the regenerated bone in many stages of bone callus growth. The experimental measurements using the developed diagnostic were compared to the numerically simulated results. For this external fixator, it was possible to show that the displacement in osteotomy was always lower than the displacement prescribed in the elongator. Nevertheless, a relationship was established between the energy consumption by the electromechanical system used to perform callus stimulus and the degree of osteotomy consolidation. Hence, this technology may lead to methodologies of mechanical stimulation for regenerating bone, which will play a relevant role for bedridden individuals with mobility limitations.

Image-based radiodensity profilometry measures early remodeling at the bone-callus interface in sheep

Bone healing has been traditionally described as a four-phase process: inflammatory response, soft callus formation, hard callus development, and remodeling. The remodeling phase has been largely neglected in most numerical mechanoregulation models of fracture repair in favor of capturing early healing using a pre-defined callus domain. However, in vivo evidence suggests that remodeling occurs concurrently with repair and causes changes in cortical bone adjacent to callus that are typically neglected in numerical models of bone healing. The objective of this study was to use image processing techniques to quantify this early-stage remodeling in ovine osteotomies. To accomplish this, we developed a numerical method for radiodensity profilometry with optimization-based curve fitting to mathematically model the bone density gradients in the radial direction across the cortical wall and callus. After assessing data from 26 sheep, we defined a dimensionless density fitting function that revealed significant remodeling occurring in the cortical wall adjacent to callus during early healing, a 23% average reduction in density compared to intact. This fitting function is robust for modeling radial density gradients in both intact bone and fracture repair scenarios and can capture a wide variety of the healing responses. The fitting function can also be scaled easily for comparison to numerical model predictions and may be useful for validating future mechanoregulatory models of coupled fracture repair and remodeling.

Computational model of endochondral ossification: Simulating growth of a long bone

Mechanical loading is a crucial factor in joint and bone development. Using a computational model, we investigated the role of mechanics on cartilage growth rate, ossification of the secondary center, formation of the growth plate, and overall bone shape. A computational algorithm was developed and implemented into finite element models to simulate the endochondral ossification for symmetric and asymmetric motion in a generic diarthrodial joint. Under asymmetric loading condition the secondary center ossifies asymmetrically leaning toward the external load and results in tilted growth plate. Also the mechanics seems to have greater influence in the early onset of the ossification of the secondary center rather than later progression of the center. While previous models have simulated select stages of skeletal development, our model can simulate growth and ossification during the entirety of post-natal development. Such computational models of skeletal development may provide insight into specific loading conditions that cause bone and joint deformities, and the required timing for rehabilitative repair.

Bone fracture healing within a continuum bone remodelling framework

Bone fracture healing is a complex process which is still under research. Computer-aided patient-specific prediction of bone development, fracture risk, prevention and treatment approaches promises a significant milestone in clinical practice. With this long-term goal in mind, a novel model is presented and examined in this work in the context of continuum bone remodelling. Therein, a clear distinction is made between external mechanical stimulation and the biological healing process of an injured bone tissue. The model is implemented within a finite element framework and investigated for the example of a fractured proximal femur head. The results show promising perspectives for further application. Besides, the model offers the possibility of easily integrating other factors like age-dependency and the availability of nutrition. For the future, further studies with large clinical datasets are essential for validation.

Multiscale modeling of bone tissue mechanobiology

Mechanical environment has a crucial role in our organism at the different levels, ranging from cells to tissues and our own organs. This regulatory role is especially relevant for bones, given their importance as load-transmitting elements that allow the movement of our body as well as the protection of vital organs from load impacts. Therefore bone, as living tissue, is continuously adapting its properties, shape and repairing itself, being the mechanical loads one of the main regulatory stimuli that modulate this adaptive behavior. Here we review some key results of bone mechanobiology from computational models, describing the effect that changes associated to the mechanical environment induce in bone response, implant design and scaffold-driven bone regeneration.

Diagnostic prediction of ovine fracture healing outcomes via a novel multi-location direct electromagnetic coupling antenna

Background

Expedient prediction of adverse bone fracture healing (delayed- or non-union) is necessary to advise secondary treatments for improving healing outcome to minimize patient suffering. Radiographic imaging, the current standard diagnostic, remains largely ineffective at predicting nonunions during the early stages of fracture healing resulting in mean nonunion diagnosis times exceeding six months. Thus, there remains a clinical deficit necessitating improved diagnostic techniques. It was hypothesized that adverse fracture healing expresses impaired biological progression at the fracture site, thus resulting in reduced temporal progression of fracture site stiffness which may be quantified prior to the appearance of radiographic indicators of fracture healing (i.e., calcified tissue).

Methods

A novel multi-location direct electromagnetic coupling antenna was developed to diagnose relative changes in the stiffness of fractures treated by metallic orthopaedic hardware. The efficacy of this diagnostic was evaluated during fracture healing simulated by progressive destabilization of cadaveric ovine metatarsals treated by locking plate fixation (n=8). An ovine in vivo comparative fracture study (n=8) was then utilized to better characterize the performance of the developed diagnostic in a clinically translatable setting. In vivo measurements using the developed diagnostic were compared to weekly radiographic images and postmortem biomechanical, histological, and micro computed tomography analyses.

Results

For all cadaveric samples, the novel direct electromagnetic coupling antenna displayed significant differences at the fracture site (P<0.05) when measuring a fully fractured sample versus partially intact and fully intact fracture states. In subsequent in vivo fracture models, this technology detected significant differences (P<0.001) in fractures trending towards delayed healing during the first 30 days post-fracture.

Conclusions

This technology, relative to traditional X-ray imaging, exhibits potential to greatly expedite clinical diagnosis of fracture nonunion, thus warranting additional technological development.

Design and analysis of biodegradable buttress threaded screws for fracture fixation in orthopedics: a finite element analysis

Screws made up of non-biodegradable materials (Ti-alloy, etc.) have been used since long for temporary joining/fixation in applications involving skeleton damage or bone fracture. These screws need to be removed after complete healing as their sustained presence results in many complications, such as - micro-fracturing, stress shielding, etc. The removal of these screws is a little difficult too as it may result in the healed bone getting broken/damaged again. These problems can be overcome by employing metallic implants (plate, screws, etc.) made up of biodegradable metallic materials (Mg-alloy, etc.). Such implants exhibit optimal mechanical performance, are biocompatible, have adequate biodegradation rates, and rely on a unique design. Internal fracture fixation makes usage of screws with or without an accompanying plate. Buttress-threaded screws are the most frequently used ones. These screws must have the capacity to bear usually occurring loads and hold fractured segments of bone all through the process of healing. Finite element analysis (FEA) is an effective technique used for testing and validation of desired characteristics for Mg-based biodegradable buttress-threaded screw (BBTS). The characteristics of interest include maximum possible pullout resistance to tightly hold segments of bone, torsional ability for tightening or tapping, bending ability during providing plate support by screw head, and resistance to combined loading (tensile/compressive and bending) during the self-support stage using merely the screw(s). According to test results and subsequent validation through discretization error and convergence plot, BBTS made up of Mg-alloy are found safe for regular applications under usually encountered impact loads. Topological optimization and vibration analysis are also performed wherein it is observed that design of BBTS is good enough for possible usage in fracture fixation in orthopaedics.

Domain-independent simulation of physiologically relevant callus shape in mechanoregulated models of fracture healing

Mechanoregulatory models have been used to predict the progression of bone fracture healing for more than two decades. However, many published studies share the same fundamental limitation: callus development proceeds within a pre-defined domain that both restricts and directs healing and leads to some non-physiologic healing patterns. To address this limitation, we added two spatial proximity functions to an existing mechanoregulatory model of fracture healing to control the localization of callus within the healing domain. We tested the performance of the new model in an idealized ovine tibial osteotomy with medial plate fixation using three sizes of healing domains and multiple variations of the spatial proximity functions. All model variations produced outward callus growth and bridging weighted toward the far cortex, which is consistent with in vivo healing. With and without the proximity functions, there were marked differences in the predicted callus volume and shape. With no proximity functions, the callus produced was strongly domain dependent, with a 15% difference in volume between the smallest and largest initialization domains. With proximity function control, callus growth was restricted to near the fracture line and there was only 2% difference in volume between domain sizes. Superimposing both proximity functions - one to control outward growth and one representing a decay in periosteal activity away from the fracture - produced a predicted callus size that was within the physiologic range for sheep and had a realistic morphology when compared with fluorescent dye co-localization with calcium deposition over time and histology.

Numerical Analysis of the Influence of Porosity and Pore Geometry on Functionality of Scaffolds Designated for Orthopedic Regenerative Medicine

BACKGROUND

Scaffolds are vital for orthopedic regenerative medicine. Therefore, comprehensive studies evaluating their functionality with consideration of variable parameters are needed. The research aim was to evaluate pore geometry and scaffold porosity influence on first, cell culture efficiency in a perfusion bioreactor and second, osteogenic cell diffusion after its implantation.

METHODS

For the studies, five pore geometries were selected (triangular prism with a rounded and a flat profile, cube, octagonal prism, sphere) and seven porosities (up to 80%), on the basis of which 70 models were created for finite element analyses. First, scaffolds were placed inside a flow channel to estimate growth medium velocity and wall shear stress. Secondly, scaffolds were placed in a bone to evaluate osteogenic cell diffusion.

RESULTS

In terms of fluid minimal velocity (0.005 m/s) and maximal wall shear stress (100 mPa), only cubic and octagonal pores with 30% porosity and spherical pores with 20% porosity fulfilled the requirements. Spherical pores had the highest osteogenic cell diffusion efficiency for porosities up to 30%. For higher porosities, the octagonal prism's pores gave the best results up to 80%, where no differences were noted.

CONCLUSIONS

The data obtained allows for the appropriate selection of pore geometry and scaffold porosity for orthopedic regenerative medicine.

[Numerical simulation of fracture healing]

Fracture is a common physical injury. Its healing process involves complex biological activities at tissue, cellular and molecular levels and is affected by mechanical and biological factors. Over recent years, numerical simulation methods have been widely used to explore the mechanisms of fracture healing, design fixators and develop novel treatment strategies, etc. This paper mainly recommend the numerical methods used for simulating fracture healing and their latest research progress, which helps people better understand the mechanism of fracture healing, and also provides direction and guidance for the numerical simulation research of fracture healing in the future. First, the fracture healing process and its relationship with mechanical stimulation and biological factors are described. Then, the numerical models used for simulating fracture healing (including mechano-regulatory model, biological regulatory model and mechano-biological regulatory model) and corresponding modeling techniques (mainly including agent-based techniques and fuzzy logic controlling method) were summarized in particular. Finally, the future research directions in numerical simulation of fracture healing were preliminarily prospected.