Angiogenesis in bone fracture healing: a bioregulatory model

Impacts of permeability and effective diffusivity of porous scaffolds on bone ingrowth: In silico and in vivo analyses

The permeability and the effective diffusivity of a porous scaffold are critical in the bone-ingrowth process. However, design guidelines for porous structures are still lacking due to inadequate understanding of the complex physiological processes involved. In this study, a model integrating the fundamental biological processes of bone regeneration was constructed to investigate the roles of permeability and effective diffusivity in regulating bone deposition in scaffolds. The in silico analysis results were confirmed in vivo by examining bone depositions in three diamond lattice scaffolds manufactured using selective laser melting. The findings show that the scaffolds with better permeability and effective diffusivity had deeper bone ingrowth and greater bone volume. Compared to permeability, effective diffusivity exhibited greater sensitivity to the orientation of porous structures, and bone ingrowth was deeper in the directions with higher effective diffusivity in spite of identical pore size. A 4.8-fold increase in permeability and a 1.6-fold increase in effective diffusivity by changing the porous structure led to a 1.5-fold increase in newly formed bone. The effective diffusivity of the porous scaffold affects the distribution of osteogenic growth factor, which in turn impacts cell migration and bone deposition through chemotaxis effects. Therefore, effective diffusivity may be a more suitable indicator for porous scaffolds because our study shows changes in this parameter determine changes in bone distribution and bone volume.

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.

Angiogenesis is uncoupled from osteogenesis during calvarial bone regeneration

Bone regeneration requires a well-orchestrated cellular and molecular response including robust vascularization and recruitment of mesenchymal and osteogenic cells. In femoral fractures, angiogenesis and osteogenesis are closely coupled during the complex healing process. Here, we show with advanced longitudinal intravital multiphoton microscopy that early vascular sprouting is not directly coupled to osteoprogenitor invasion during calvarial bone regeneration. Early osteoprogenitors emerging from the periosteum give rise to bone-forming osteoblasts at the injured calvarial bone edge. Microvessels growing inside the lesions are not associated with osteoprogenitors. Subsequently, osteogenic cells collectively invade the vascularized and perfused lesion as a multicellular layer, thereby advancing regenerative ossification. Vascular sprouting and remodeling result in dynamic blood flow alterations to accommodate the growing bone. Single cell profiling of injured calvarial bones demonstrates mesenchymal stromal cell heterogeneity comparable to femoral fractures with increase in cell types promoting bone regeneration. Expression of angiogenesis and hypoxia-related genes are slightly elevated reflecting ossification of a vascularized lesion site. Endothelial Notch and VEGF signaling alter vascular growth in calvarial bone repair without affecting the ossification progress. Our findings may have clinical implications for bone regeneration and bioengineering approaches.

An in silico micro-multiphysics agent-based approach for simulating bone regeneration in a mouse femur defect model

Bone defects represent a challenging clinical problem as they can lead to non-union. In silico models are well suited to study bone regeneration under varying conditions by linking both cellular and systems scales. This paper presents an in silico micro-multiphysics agent-based (micro-MPA) model for bone regeneration following an osteotomy. The model includes vasculature, bone, and immune cells, as well as their interaction with the local environment. The model was calibrated by time-lapsed micro-computed tomography data of femoral osteotomies in C57Bl/6J mice, and the differences between predicted bone volume fractions and the longitudinal in vivo measurements were quantitatively evaluated using root mean square error (RMSE). The model performed well in simulating bone regeneration across the osteotomy gap, with no difference (5.5% RMSE, p = 0.68) between the in silico and in vivo groups for the 5-week healing period - from the inflammatory phase to the remodelling phase - in the volume spanning the osteotomy gap. Overall, the proposed micro-MPA model was able to simulate the influence of the local mechanical environment on bone regeneration, and both this environment and cytokine concentrations were found to be key factors in promoting bone regeneration. Further, the validated model matched clinical observations that larger gap sizes correlate with worse healing outcomes and ultimately simulated non-union. This model could help design and guide future experimental studies in bone repair, by identifying which are the most critical in vivo experiments to perform.

External mechanical loading overrules cell-cell mechanical communication in sprouting angiogenesis during early bone regeneration

Sprouting angiogenesis plays a key role during bone regeneration. For example, insufficient early revascularization of the injured site can lead to delayed or non-healing. During sprouting, endothelial cells are known to be mechano-sensitive and respond to local mechanical stimuli. Endothelial cells interact and communicate mechanically with their surroundings, such as outer-vascular stromal cells, through cell-induced traction forces. In addition, external physiological loads act at the healing site, resulting in tissue deformations and impacting cellular arrangements. How these two distinct mechanical cues (cell-induced and external) impact angiogenesis and sprout patterning in early bone healing remains however largely unknown. Therefore, the aim of this study was to investigate the relative role of externally applied and cell-induced mechanical signals in driving sprout patterning at the onset of bone healing. To investigate cellular self-organisation in early bone healing, an in silico model accounting for the mechano-regulation of sprouting angiogenesis and stromal cell organization was developed. Computer model predictions were compared to in vivo experiments of a mouse osteotomy model stabilized with a rigid or a semirigid fixation system. We found that the magnitude and orientation of principal strains within the healing region can explain experimentally observed sprout patterning, under both fixation conditions. Furthermore, upon simulating the selective inhibition of either cell-induced or externally applied mechanical cues, external mechanical signals appear to overrule the mechanical communication acting on a cell-cell interaction level. Such findings illustrate the relevance of external mechanical signals over the local cell-mediated mechanical cues and could be used in the design of fracture treatment strategies for bone regeneration.

In silico assessment of the bone regeneration potential of complex porous scaffolds

Mechanical environment plays a crucial role in regulating bone regeneration in bone defects. Assessing the mechanobiological behavior of patient-specific orthopedic scaffolds in-silico could help guide optimal scaffold designs, as well as intra- and post-operative strategies to enhance bone regeneration and improve implant longevity. Additively manufactured porous scaffolds, and specifically triply periodic minimal surfaces (TPMS), have shown promising structural properties to act as bone substitutes, yet their ability to induce mechanobiologially-driven bone regeneration has not been elucidated. The aim of this study is to i) explore the bone regeneration potential of TPMS scaffolds made of different stiffness biocompatible materials, to ii) analyze the influence of pre-seeding the scaffolds and increasing the post-operative resting period, and to iii) assess the influence of patient-specific parameters, such as age and mechanosensitivity, on outcomes. To perform this study, an in silico model of a goat tibia is used. The bone ingrowth within the scaffold pores was simulated with a mechano-driven model of bone regeneration. Results showed that the scaffold's architectural properties affect cellular diffusion and strain distribution, resulting in variations in the regenerated bone volume and distribution. The softer material improved the bone ingrowth. An initial resting period improved the bone ingrowth but not enough to reach the scaffold's core. However, this was achieved with the implantation of a pre-seeded scaffold. Physiological parameters like age and health of the patient also influence the bone regeneration outcome, though to a lesser extent than the scaffold design. This analysis demonstrates the importance of the scaffold's geometry and its material, and highlights the potential of using mechanobiological patient-specific models in the design process for bone substitutes.

Flavonoid-Loaded Biomaterials in Bone Defect Repair

Skeletons play an important role in the human body, and can form gaps of varying sizes once damaged. Bone defect healing involves a series of complex physiological processes and requires ideal bone defect implants to accelerate bone defect healing. Traditional grafts are often accompanied by issues such as insufficient donors and disease transmission, while some bone defect implants are made of natural and synthetic polymers, which have characteristics such as good porosity, mechanical properties, high drug loading efficiency, biocompatibility and biodegradability. However, their antibacterial, antioxidant, anti-inflammatory and bone repair promoting abilities are limited. Flavonoids are natural compounds with various biological activities, such as antitumor, anti-inflammatory and analgesic. Their good anti-inflammatory, antibacterial and antioxidant activities make them beneficial for the treatment of bone defects. Several researchers have designed different types of flavonoid-loaded polymer implants for bone defects. These implants have good biocompatibility, and they can effectively promote the expression of angiogenesis factors such as VEGF and CD31, promote angiogenesis, regulate signaling pathways such as Wnt, p38, AKT, Erk and increase the levels of osteogenesis-related factors such as Runx-2, OCN, OPN significantly to accelerate the process of bone defect healing. This article reviews the effectiveness and mechanism of biomaterials loaded with flavonoids in the treatment of bone defects. Flavonoid-loaded biomaterials can effectively promote bone defect repair, but we still need to improve the overall performance of flavonoid-loaded bone repair biomaterials to improve the bioavailability of flavonoids and provide more possibilities for bone defect repair.

Development of numerical model-based machine learning algorithms for different healing stages of distal radius fracture healing

BACKGROUND AND OBJECTIVES

Early therapeutic exercises are vital for the healing of distal radius fractures (DRFs) treated with the volar locking plate. However, current development of rehabilitation plans using computational simulation is normally time-consuming and requires high computational power. Thus, there is a clear need for developing machine learning (ML) based algorithms that are easy for end-users to implement in daily clinical practice. The purpose of the present study is to develop optimal ML algorithms for designing effective DRF physiotherapy programs at different stages of healing.

METHOD

First, a three-dimensional computational model for the healing of DRF was developed by integrating mechano-regulated cell differentiation, tissue formation and angiogenesis. The model is capable of predicting time-dependant healing outcomes based on different physiologically relevant loading conditions, fracture geometries, gap sizes, and healing time. After being validated using available clinical data, the developed computational model was implemented to generate a total of 3600 clinical data for training the ML models. Finally, the optimal ML algorithm for each healing stage was identified.

RESULTS

The selection of the optimal ML algorithm depends on the healing stage. The results from this study show that cubic support vector machine (SVM) has the best performance in predicting the healing outcomes at the early stage of healing, while trilayered ANN outperforms other ML algorithms in the late stage of healing. The outcomes from the developed optimal ML algorithms indicate that Smith fractures with medium gap sizes could enhance the healing of DRF by inducing larger cartilaginous callus, while Colles fractures with large gap sizes may lead to delayed healing by bringing excessive fibrous tissues.

CONCLUSIONS

ML represents a promising approach for developing efficient and effective patient-specific rehabilitation strategies. However, ML algorithms at different healing stages need to be carefully chosen before being implemented in clinical applications.

Landscape of Well-Coordinated Fracture Healing in a Mouse Model Using Molecular and Cellular Analysis

The success of fracture healing relies on overlapping but coordinated cellular and molecular events. Characterizing an outline of differential gene regulation throughout successful healing is essential for identifying crucial phase-specific markers and may serve as the basis for engineering these in challenging healing situations. This study analyzed the healing progression of a standard closed femoral fracture model in C57BL/6N (age = 8 weeks) wild-type male mice. The fracture callus was assessed across various days post fracture (D = days 0, 3, 7, 10, 14, 21, and 28) by microarray, with D0 serving as a control. Histological analyses were carried out on samples from D7 until D28 to support the molecular findings. Microarray analysis revealed a differential regulation of immune response, angiogenesis, ossification, extracellular matrix regulation, mitochondrial and ribosomal genes during healing. In-depth analysis showed differential regulation of mitochondrial and ribosomal genes during the initial phase of healing. Furthermore, the differential gene expression showed an essential role of Serpin Family F Member 1 over the well-known Vascular Endothelial Growth Factor in angiogenesis, especially during the inflammatory phase. The significant upregulation of matrix metalloproteinase 13 and bone sialoprotein from D3 until D21 asserts their importance in bone mineralization. The study also shows type I collagen around osteocytes located in the ossified region at the periosteal surface during the first week of healing. Histological analysis of matrix extracellular phosphoglycoprotein and extracellular signal-regulated kinase stressed their roles in bone homeostasis and the physiological bone-healing process. This study reveals previously unknown and novel candidates, that could serve as a target for specific time points in healing and to remedy cases of impaired healing.

Influence of muscle loading on early-stage bone fracture healing

Designing weight-bearing exercises for patients with lower-limb bone fractures is challenging and requires a systematic approach that accounts for patient-specific loading conditions. However, 'trial-and-error' approaches are commonplace in clinical settings due to the lack of a fundamental understanding of the effect of weight-bearing exercises on the bone healing process. Whilst computational modelling has the potential to assist clinicians in designing effective patient-specific weight-bearing exercises, current models do not explicitly account for the effects of muscle loading, which could play an important role in mediating the mechanical microenvironment of a fracture site. We combined a fracture healing model involving a tibial fracture stabilised with a locking compression plate (LCP) with a detailed musculoskeletal model of the lower limb to determine interfragmentary strains in the vicinity of the fracture site during both full weight-bearing (100% body weight) and partial weight-bearing (50% body weight) standing. We found that muscle loading significantly altered model predictions of interfragmentary strains. For a fractured bone with a standard LCP configuration (bone-plate distance = 2 mm, working length = 30 mm) subject to full weight-bearing, the predicted strains at the near and far cortices were 23% and 11% higher when muscle loading was included compared to the case when muscle loading was omitted. The knee and ankle muscles accounted for 38% of the contact force exerted at the knee joint during quiet standing and contributed significantly to the strains calculated at the fracture site. Thus, models of bone fracture healing ought to account explicitly for the effects of muscle loading. Furthermore, the study indicated that LCP configuration parameters play a crucial role in influencing the fracture site microenvironment. The results highlighted the dominance of working length over bone-plate distance in controlling the flexibility of fracture sites stabilised with LCP devices.

Predicting the formation of different tissue types during Achilles tendon healing using mechanoregulated and oxygen-regulated frameworks

During Achilles tendon healing in rodents, besides the expected tendon tissue, also cartilage-, bone- and fat-like tissue features have been observed during the first twenty weeks of healing. Several studies have hypothesized that mechanical loading may play a key role in the formation of different tissue types during healing. We recently developed a computational mechanobiological framework to predict tendon tissue production, organization and mechanical properties during tendon healing. In the current study, we aimed to explore possible mechanobiological related mechanisms underlying formation of other tissue types than tendon tissue during tendon healing. To achieve this, we further developed our recent framework to predict formation of different tissue types, based on mechanobiological models established in other fields, which have earlier not been applied to study tendon healing. We explored a wide range of biophysical stimuli, i.e., principal strain, hydrostatic stress, pore pressure, octahedral shear strain, fluid flow, angiogenesis and oxygen concentration, that may promote the formation of different tissue types. The numerical framework predicted spatiotemporal formation of tendon-, cartilage-, bone- and to a lesser degree fat-like tissue throughout the first twenty weeks of healing, similar to recent experimental reports. Specific features of experimental data were captured by different biophysical stimuli. Our modeling approach showed that mechanobiology may play a role in governing the formation of different tissue types that have been experimentally observed during tendon healing. This study provides a numerical tool that can contribute to a better understanding of tendon mechanobiology during healing. Developing these tools can ultimately lead to development of better rehabilitation regimens that stimulate tendon healing and prevent unwanted formation of cartilage-, fat- and bone-like tissues.

The osteogenetic activities of mesenchymal stem cells in response to Mg2+ ions and inflammatory cytokines: a numerical approach using fuzzy logic controllers

Magnesium (Mg²⁺) ions are frequently reported to regulate osteogenic activities of mesenchymal stem cells (MSCs). In this study, we propose a numerical model to study the regulatory importance of Mg²⁺ ions on MSCs osteoblastic differentiation in the presence of an inflammatory response. A fuzzy logic controller was formulated to receive the concentrations of Mg²⁺ ions and the inflammatory cytokines of TNF-α, IL-10, IL-1β, and IL-8 as cellular inputs and predict the cells’ early and late differentiation rates. Five sets of empirical data obtained from published cell culture experiments were used to calibrate the model. The model successfully reproduced the empirical data regarding the concentration- and phase-dependent effect of Mg²⁺ ions on the differentiation process. In agreement with the experiments, the model showed the stimulatory role of Mg²⁺ ions on the early differentiation phase, once administered at low concentration, and their inhibitory role on the late differentiation phase. The numerical approach used in this study suggested 6–8 mM as the most effective concentration of Mg²⁺ ions in promoting the early differentiation process. Also, the proposed model sheds light on the fundamental differences in the behavioral properties of cells cultured in different experiments, e.g. differentiation rate and the sensitivity of the cultured cells to stimulatory signals such as Mg²⁺ ions. Thus, it can be used to interpret and compare different empirical findings. Moreover, the model successfully reproduced the nonlinearities in the concentration-dependent role of the inflammatory cytokines in early and late differentiation rates. Overall, the proposed model can be employed in studying the osteogenic properties of Mg-based implants in the presence of an inflammatory response.

Magnesium (Mg) is an attractive material for bone implants as it fully degrades after implantation, saving pain and cost of the second surgery for implant removal. To advance its application in the orthopedic industry, it is paramount to fully understand the biological impact of the degradation products, in particular Mg²⁺ ions. Here, we propose a computer model to study the effects of Mg²⁺ ions on bone regeneration. The model focuses on stem cells and includes both the direct stimulation effects of Mg²⁺ ions on cells and the indirect stimulus through the inflammatory system. The proposed model successfully reproduced the experimental data of five different studies. The model additionally highlighted differences amongst different experiments in terms of the cellular response to Mg²⁺ ions. The proposed system therefore provides an important addition to the field of Mg implant research.

A preclinical model of post-surgery secondary bone healing for subtrochanteric femoral fracture based on fuzzy interpretations

Mechanobiology plays an essential role in secondary bone fracture healing. While the introduction of newer type of plates, e.g. locking plate (LP), is becoming increasingly popular for complex femoral fractures, the conventional technique involving dynamic compression plate (DCP) remains the standard choice. The difference between the two techniques lies primarily in their screw fixation mechanisms. The present study applied 3D dynamic fracture healing scheme modelled on a subtrochanteric femur fracture, regulated by both finite element (FE) analysis and Fuzzy logic control in order to understand the spatio-temporal healing phenomena for both LP and DCP. The study further examined the influence of the two screw fixation mechanisms in determining the comparative progression of fracture healing. The problem was solved iteratively in several healing steps running in loop and accordingly, the local tissue concentrations and material properties were updated. The predicted results accorded well with various previous experimental observations. The study found an initial delay in healing associated with DCP. However, as the healing progressed, there was no significant difference in overall callus modulus. The presented preclinical model may further help predict bone healing for different implantation techniques, and thus can serve as a non-invasive tool for evaluating relative merits of extramedullary plating techniques.

Simulating In Vitro the Bone Healing Potential of a Degradable and Tailored Multifunctional Mg-Based Alloy Platform

This work intended to elucidate, in an in vitro approach, the cellular and molecular mechanisms occurring during the bone healing process, upon implantation of a tailored degradable multifunctional Mg-based alloy. This was prepared by a conjoining anodization of the bare alloy (AZ31) followed by the deposition of a polymeric coating functionalized with hydroxyapatite. Human endothelial cells and osteoblastic and osteoclastic differentiating cells were exposed to the extracts from the multifunctional platform (having a low degradation rate), as well as the underlying anodized and original AZ31 alloy (with higher degradation rates). Extracts from the multifunctional coated alloy did not affect cellular behavior, although a small inductive effect was observed in the proliferation and gene expression of endothelial and osteoblastic cells. Extracts from the higher degradable anodized and original alloys induced the expression of some endothelial genes and, also, ALP and TRAP activities, further increasing the expression of some early differentiation osteoblastic and osteoclastic genes. The integration of these results in a translational approach suggests that, following the implantation of a tailored degradable Mg-based material, the absence of initial deleterious effects would favor the early stages of bone repair and, subsequently, the on-going degradation of the coating and the subjacent alloy would increase bone metabolism dynamics favoring a faster bone formation and remodeling process and enhancing bone healing.

Quantitative Assessment of the Restoration of Original Anatomy after 3D Virtual Reduction of Long Bone Fractures

Background: The purpose of this study was to demonstrate the usefulness of 3D image-based virtual reduction by validating the evaluation criteria according to guidelines suggested by the AO Surgery Reference. Methods: For this experiment, 19 intact radial ORTHObones (ORTHObones radius, 3B Scientific, Germany, Hamburg) without any fractures were prepared. All ORTHObones with six cortical marking holes (three points on the distal part and three points on the proximal part) were scanned using a CT scanner twice (before/after intentional fracture of the ORTHObone). After the virtual reduction of all 19 ORTHObones, accuracy evaluations using the four criteria (length variation, apposition variation, alignment variation, Rotation Variation) suggested in the AO Surgery Reference were performed. Results: The mean (M) length variation was 0.42 mm, with 0.01 mm standard deviation (SD). The M apposition variation was 0.48 mm, with 0.40 mm SD. The M AP angulation variation (for alignment variation) was 3.24°, with 2.95° SD. The M lateral angulation variation (for alignment variation) was 0.09°, with 0.13° SD. The M angle of axial rotation was 1.27° with SD: 1.19°. Conclusions: The method of accuracy evaluation used in this study can be helpful in establishing a reliable plan.

Influence of therapeutic grip exercises induced loading rates in distal radius fracture healing with volar locking plate fixation

BACKGROUND AND OBJECTIVE

Therapeutic exercises could potentially enhance the healing of distal radius fractures (DRFs) treated with volar locking plate (VLP). However, the healing outcomes are highly dependant on the patient-specific fracture geometries (e.g., gap size) and the loading conditions at the fracture site (e.g., loading frequency) resulted from different types of therapeutic exercises. The purpose of this study is to investigate the effects of different loading frequencies induced by therapeutic exercises on the biomechanical microenvironment of the fracture site and the transport of cells and growth factors within the fracture callus, ultimately the healing outcomes. This is achieved through numerical modelling and mechanical testing.

METHODS

Five radius sawbones specimens (Pacific Research Laboratories, Vashon, USA) fixed with VLP (VRP2.0+, Austofix) were mechanically tested using dynamic test instrument (INSTRON E3000, Norwood, MA). The loading protocol used in mechanical testing involved a series of cyclic axial compression tests representing hand and finger therapeutic exercises. The relationship between the dynamic loading rate (i.e., loading frequency) and dynamic stiffness of the construct was established and used as inputs to a developed numerical model for studying the dynamic loading induced cells and growth factors in fracture site and biomechanical stimuli required for healing.

RESULTS

There is a strong positive linear relationship between the loading rate and axial stiffness of the construct fixed with VLP. The loading rates induced by the moderate frequencies (i.e., 1-2 Hz) could promote endochondral ossification, whereas relatively high loading frequencies (i.e., over 3 Hz) may hinder the healing outcomes or lead to non-union. In addition, a dynamic loading frequency of 2 Hz in combination of a fracture gap size of 3 mm could produce a better healing outcome by enhancing the transport of cells and growth factors at the fracture site in comparison to free diffusion (i.e. without loading), and thereby produces a biomechanical microenvironment which is favourable for healing.

CONCLUSION

The experimentally validated numerical model presented in this study could potentially contribute to the design of effective patient-specific therapeutic exercises for better healing outcomes. Importantly, the model results demonstrate that therapeutic grip exercises induced dynamic loading could produce a better biomechanical microenvironment for healing without compromising the mechanical stability of the overall volar locking plate fixation construct.

Bone Ingrowth Around an Uncemented Femoral Implant Using Mechanoregulatory Algorithm: A Multiscale Finite Element Analysis

The primary fixation and long-term stability of a cementless femoral implant depend on bone ingrowth within the porous coating. Although attempts were made to quantify the peri-implant bone ingrowth using the finite element (FE) analysis and mechanoregulatory principles, the tissue differentiation patterns on a porous-coated hip stem have scarcely been investigated. The objective of this study is to predict the spatial distribution of evolutionary bone ingrowth around an uncemented hip stem, using a three-dimensional (3D) multiscale mechanobiology-based numerical framework. Multiple load cases representing a variety of daily living activities, including walking, stair climbing, sitting down, and standing up from a chair, were used as applied loading conditions. The study accounted for the local variations in host bone material properties and implant-bone relative displacements of the macroscale implanted FE model, in order to predict bone ingrowth in microscale representative volume elements (RVEs) of 12 interfacial regions. In majority RVEs, 20-70% bone tissue (immature and mature) was predicted after 2 months, contributing toward a progressive increase in average Young's modulus (1200-3000 MPa) of the interbead tissue layer. Higher bone ingrowth (mostly greater than 60%) was predicted in the anterolateral regions of the implant, as compared to the posteromedial side (20-50%). New bone tissue was formed deeper inside the interbead spacing, adhering to the implant surface. The study helps to gain an insight into the degree of osseointegration of a porous-coated femoral implant.

Role of TNF-α in early-stage fracture healing under normal and diabetic conditions

BACKGROUND AND OBJECTIVE

Inflammatory response plays a crucial role in the early stage of fracture healing. Immediately after fracture, the debris and immune cells (e.g., macrophages), recruited into the fracture callus, lead to the secretion of inflammatory cytokines, such as tumor necrosis factor-alpha (TNF-α), which governs the mesenchymal stem cells (MSCs) mediated healing processes. However, it is still unclear how chronic inflammatory diseases (e.g., diabetes) affect the level of TNF-α in fracture callus, ultimately the healing outcomes at the early stage of healing. Therefore, the purpose of this study is to develop a numerical model for investigating TNF-α mediated bone fracture healing.

METHODS

A mathematical model consisting of a system of partial differential equations that represent the reactive transport of cells and cytokines in the fracture callus is developed in this study. The model is first calibrated by using available experimental data and then implemented to study the effect of TNF-α on the early stage of fracture healing under normal and diabetic conditions.

RESULTS

There is a significant elevation of TNF-α level in facture callus during the first 24 h post-fracture in normal condition, and its influence in the concentration of MSCs and cell differentiation becomes significant three days post-fracture (e.g., the absence of TNF-α signaling could reduce the concentration of MSCs more than 20% in cortical callus). In addition, the excessive secretion of TNF-α induced by diabetes could decrease the concentration of MSCs at the initial stage of healing, particularly reduce the concentration of MSCs in cortical callus by around 25%.

CONCLUSION

The model predictions suggested that there should be an optimal concentration of TNF-α in fracture callus, which enhances the early stage of healing, and excessive or insufficient secretion of TNF-α might significantly hinder the healing process.

Mechanobiological model to study the influence of screw design and surface treatment on osseointegration

This study aims at suggesting a new approach to peri-implant healing models, providing a set of taxis-diffusion-reaction equations under the combined influence of mechanical and biochemical factors. Early events of osseointegration were simulated for titanium screw implants inserted into a pre-drilled trabecular bone environment, up to 12 weeks of peri-implant bone healing. Simulations showed the ability of the model to reproduce biological events occurring at the implant interface through osteogenesis. Implants with shallow healing chamber showed higher proportions of lamellar bone, enhanced by the increase of mechanical stimulation. Osteoconduction was observed through the surface treatment model and similar bone healing patterns compared to in vivo studies.

In silico and in vivo studies of the effect of surface curvature on the osteoconduction of porous scaffolds

Recent evidence shows that the curvature of porous scaffold plays a significant role in guiding tissue regeneration. However, the underlying mechanism remains controversial to date. In this study, we developed an in silico model to simulate the effect of surface curvature on the osteoconduction of scaffold implants, which comprises the primary aspects of bone regeneration. Selective laser melting was used to manufacture a titanium scaffold with channels representative of different strut curvatures for in vivo assessment. The titanium scaffold was implanted in the femur condyles of rabbits to validate the mathematical model. Simulation results suggest that the curvature affected the distribution of growth factors and subsequently induced the migration of osteoblast lineage cells and bone deposition to the locations with higher curvature. The predictions of the mathematical model are in good agreement with the in vivo assessment results, in which newly formed bone first appeared adjacent to the vertices of the major axes in elliptical channels. The mechanism of curvature-guided osteoconduction may provide a guide for the design optimization of scaffold implants to achieve enhanced bone ingrowth.

Towards in silico Models of the Inflammatory Response in Bone Fracture Healing

In silico modeling is a powerful strategy to investigate the biological events occurring at tissue, cellular and subcellular level during bone fracture healing. However, most current models do not consider the impact of the inflammatory response on the later stages of bone repair. Indeed, as initiator of the healing process, this early phase can alter the regenerative outcome: if the inflammatory response is too strongly down- or upregulated, the fracture can result in a non-union. This review covers the fundamental information on fracture healing, in silico modeling and experimental validation. It starts with a description of the biology of fracture healing, paying particular attention to the inflammatory phase and its cellular and subcellular components. We then discuss the current state-of-the-art regarding in silico models of the immune response in different tissues as well as the bone regeneration process at the later stages of fracture healing. Combining the aforementioned biological and computational state-of-the-art, continuous, discrete and hybrid modeling technologies are discussed in light of their suitability to capture adequately the multiscale course of the inflammatory phase and its overall role in the healing outcome. Both in the establishment of models as in their validation step, experimental data is required. Hence, this review provides an overview of the different in vitro and in vivo set-ups that can be used to quantify cell- and tissue-scale properties and provide necessary input for model credibility assessment. In conclusion, this review aims to provide hands-on guidance for scientists interested in building in silico models as an additional tool to investigate the critical role of the inflammatory phase in bone regeneration.

A numerical bone regeneration model incorporating angiogenesis, considering oxygen-induced secretion of vascular endothelial growth factor and vascular remodeling

Angiogenesis is considered playing an important role in bone regeneration. Studies have shown that angiogenesis is affected by biological factors, oxygen tension, and blood flow. In this paper, we propose a bone regeneration model with angiogenesis based on the theories of mechanobiology regulation, vascular network modeling, oxygen-induced secretion of vascular endothelial growth factor (VEGF), and vascular remodeling. The results showed that this model can describe the distribution and concentration of vascular endothelial growth factor induced by oxygen tension during bone regeneration, the growth and remodeling of vascular tissue under the influence of vascular endothelial growth factor and mechanical loading, and the correspondence between vascular tissue and bone regeneration.

Magnesium ions regulate mesenchymal stem cells population and osteogenic differentiation: A fuzzy agent-based modeling approach

Mesenchymal stem cells (MSCs) are proliferative and multipotent cells that play a key role in the bone regeneration process. Empirical data have repeatedly shown the bioregulatory importance of magnesium (Mg) ions in MSC growth and osteogenesis. In this study, we propose an agent-based model to predict the spatiotemporal dynamics of the MSC population and osteogenic differentiation in response to Mg2+ ions. A fuzzy-logic controller was designed to govern the decision-making process of cells by predicting four cellular processes of proliferation, differentiation, migration, and mortality in response to several important bioregulatory factors such as Mg2+ ions, pH, BMP2, and TGF-β1. The model was calibrated using the empirical data obtained from three sets of cell culture experiments. The model successfully reproduced the empirical observations regarding live cell count, viability, DNA content, and the differentiation-related markers of alkaline phosphate (ALP) and osteocalcin (OC). The simulation results, in agreement with the empirical data, showed that Mg2+ ions within 3-6 mM concentration have the highest stimulation effect on cell population growth. The model also correctly reproduced the stimulatory effect of Mg2+ ions on ALP and its inhibitory effect on OC as the early and late differentiation markers, respectively. Besides, the numerical simulation shed light on the innate cellular differences of the cells cultured in different experiments in terms of the proliferative capacity as well as sensitivity to Mg2+ ions. The proposed model can be adopted in the study of the osteogenesis around Mg-based implants where ions released due to degradation interact with local cells and regulate bone regeneration.

A probabilistic approach for modelling bone fracture healing under Ilizarov circular fixator

Bone fracture treatments using Ilizarov circular fixator (ICF) involve dealing with uncertainties about a range of critical factors that control the mechanical microenvironment of the fracture site such as ICF configuration, fracture gap size, physiological loading etc. To date, the effects of the uncertainties about these critical factors on the mechanical microenvironment of the fracture site have not been fully understood. The purpose of this study is to tackle this challenge by using computational modelling in conjunction with engineering reliability analysis. Particularly, the effects of uncertainties in fracture gap size (GS), level of weight-bearing (P), ICF wire pretension (T) and wire diameter (WD) on the fracture site mechanical microenvironment at the beginning of the reparative phase of healing was investigated in this study. The results show that the mechanical microenvironment of fracture site stabilised with ICF is very sensitive to the uncertainties in P and GS. For example, an increase in the coefficient of variation of P (COV_P ) from 0.1 to 0.9 (i.e., an increase in the uncertainty in P) could reduce the probability of achieving a favourable mechanical microenvironment within the fracture site (i.e., Probability of Success, PoS) by more than 50%, while an increase in the coefficient of variation of GS (COV_GS ) from 0.1 to 0.9 could decrease PoS by around 30%. In contrast, an increase in the uncertainties in T and WD (COV increase from 0.1 to 0.9) has little influence on the fracture site mechanical microenvironment (PoS changes <5%).

Targeting angiogenesis for fracture nonunion treatment in inflammatory disease

Atrophic fracture nonunion poses a significant clinical problem with limited therapeutic interventions. In this study, we developed a unique nonunion model with high clinical relevance using serum transfer-induced rheumatoid arthritis (RA). Arthritic mice displayed fracture nonunion with the absence of fracture callus, diminished angiogenesis and fibrotic scar tissue formation leading to the failure of biomechanical properties, representing the major manifestations of atrophic nonunion in the clinic. Mechanistically, we demonstrated that the angiogenesis defect observed in RA mice was due to the downregulation of SPP1 and CXCL12 in chondrocytes, as evidenced by the restoration of angiogenesis upon SPP1 and CXCL12 treatment in vitro. In this regard, we developed a biodegradable scaffold loaded with SPP1 and CXCL12, which displayed a beneficial effect on angiogenesis and fracture repair in mice despite the presence of inflammation. Hence, these findings strongly suggest that the sustained release of SPP1 and CXCL12 represents an effective therapeutic approach to treat impaired angiogenesis and fracture nonunion under inflammatory conditions.

Optimal time-dependent levels of weight-bearing for bone fracture healing under Ilizarov circular fixators

It is known that weight-bearing exercises under Ilizarov circular fixators (ICF) could enhance bone fracture healing by mechano-regulation. However, interfragmentary movements at the fracture site induced by weight-bearing may inhibit angiogenesis and ultimately delay the healing process. To tackle this challenge, a computational model is presented in this study which considers the spatial and temporal changes in mechanical properties of fracture callus to predict optimal levels of weight-bearing during fracture healing under ICF. The study takes sheep fractures as example and shows that the developed model has the capability of predicting patient specific, time-dependent optimal levels of weight-bearing which enhances mechano-regulation mediated healing without hindering the angiogenesis process. The results demonstrate that allowable level of weight-bearing and timings depend on fracture gap size. For normal body weights (BW) and moderate fracture gap sizes (e.g. 3 mm), weight-bearing with 30% BW could start by week 4 post-operation and gradually increase to 100% BW by week 11. In contrast, for relatively large fracture gap sizes (i.e. 6 mm), weight-bearing is recommended to commence in later stages of healing (e.g. week 11 post-operation). Furthermore, increasing ICF stiffness (e.g. using half pins instead of pretension wires) can increase the level of weight-bearing significantly in the early stages up to a certain time point (e.g. week 8 post-operation) beyond which no noticeable benefits could be achieved. The findings of this study have potential applications in designing post-operative weight bearing exercises.

Digital Twins for Tissue Culture Techniques—Concepts, Expectations, and State of the Art

Techniques to provide in vitro tissue culture have undergone significant changes during the last decades, and current applications involve interactions of cells and organoids, three-dimensional cell co-cultures, and organ/body-on-chip tools. Efficient computer-aided and mathematical model-based methods are required for efficient and knowledge-driven characterization, optimization, and routine manufacturing of tissue culture systems. As an alternative to purely experimental-driven research, the usage of comprehensive mathematical models as a virtual in silico representation of the tissue culture, namely a digital twin, can be advantageous. Digital twins include the mechanistic of the biological system in the form of diverse mathematical models, which describe the interaction between tissue culture techniques and cell growth, metabolism, and the quality of the tissue. In this review, current concepts, expectations, and the state of the art of digital twins for tissue culture concepts will be highlighted. In general, DT’s can be applied along the full process chain and along the product life cycle. Due to the complexity, the focus of this review will be especially on the design, characterization, and operation of the tissue culture techniques.

Mechano-driven regeneration predicts response variations in large animal model based on scaffold implantation site and individual mechano-sensitivity

It is well founded that the mechanical environment may regulate bone regeneration in orthopedic applications. The purpose of this study is to investigate the mechanical contributions of the scaffold and the host to bone regeneration, in terms of subject specificity, implantation site and sensitivity to the mechanical environment. Using a computational approach to model mechano-driven regeneration, bone ingrowth in porous titanium scaffolds was simulated in the distal femur and proximal tibia of three goats and compared to experimental results. The results showed that bone ingrowth shifted from a homogeneous distribution pattern, when scaffolds were in contact with trabecular bone (max local ingrowth 12.47%), to a localized bone ingrowth when scaffolds were implanted in a diaphyseal location (max local ingrowth 20.64%). The bone formation dynamics revealed an apposition rate of 0.37±0.28%/day in the first three weeks after implantation, followed by limited increase in bone ingrowth until the end of the experiment (12 weeks). According to in vivo data, we identified one animal whose sensitivity to mechanical stimulation was higher than the other two. Moreover, we found that the stimulus initiating bone formation was consistently higher in the femur than in the tibia for all the individuals. Overall, the dependence of the osteogenic response on the host biomechanics means that, from a mechanical perspective, the regenerative potential depends on both the scaffold and the host environment. Therefore, this work provides insights on how the mechanical conditions of both the recipient and the scaffold contribute to meet patient and location-specific characteristics.

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.

Amygdalin Promotes Fracture Healing through TGF-β/Smad Signaling in Mesenchymal Stem Cells

Chondrogenesis and subsequent osteogenesis of mesenchymal stem cells (MSCs) and angiogenesis at injured sites are crucial for bone fracture healing. Amygdalin, a cyanogenic glycoside compound derived from bitter apricot kernel, has been reported to inhibit IL-1β-induced chondrocyte degeneration and to stimulate blood circulation, suggesting a promising role of amygdalin in fracture healing. In this study, tibial fractures in C57BL/6 mice were treated with amygdalin. Fracture calluses were then harvested and subjected to radiographic, histological, and biomechanical testing, as well as angiography and gene expression analyses to evaluate fracture healing. The results showed that amygdalin treatment promoted bone fracture healing. Further experiments using MSC-specific transforming growth factor- (TGF-) β receptor 2 conditional knockout (KO) mice (Tgfbr2^Gli1-Cre) and C3H10 T1/2 murine mesenchymal progenitor cells showed that this effect was mediated through TGF-β/Smad signaling. We conclude that amygdalin could be used as an alternative treatment for bone fractures.

Modeling osteoinduction in titanium bone scaffold with a representative channel structure

Optimizing scaffold architecture for perfect osteointegration depends on good understanding of bone ingrowth in the porous space of implants. This study developed an immunoregulatory agent-based model to discover the osteoinduction mechanism in porous scaffolds. Immunoreaction, macrophage polarization, and the corresponding growth factors were combined in the model, and all played critical roles in recruiting osteogenic cells that migrated into the scaffolds. Angiogenesis was also considered in this model. The bone ingrowth predicted by the model coincides with results from published in vivo experiments. Simulation results suggested that the pore architecture affected the diffusion process of chemotactic factors in the scaffolds, subsequently influencing the complex reactions of diverse cells and the osteoinduction location. In flexural pore spaces, bone formation spread from the periphery into the center of scaffolds due to larger M2 phenotype macrophage populations colonizing boundary regions and the distribution of corresponding growth factors concentration. In straight channels, osteogenic cells migrated further inward and osteoinduction initiated in deeper position as a result of the deeper distribution of osteogenic cytokines concentration field.

Bone marrow-derived mesenchymal stem cells promote healing of rabbit tibial fractures via JAK-STAT signaling pathway

Influence of bone marrow-derived mesenchymal stem cells (BMMSCs) on the healing of rabbit tibial fractures and the role of the Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway in fracture healing were explored. Rabbit BMMSCs were isolated and cultured in vitro, and their purity was determined using flow cytometry. Rabbit fracture models were established, and injected with BMMSCs, BMMSCs + TG101348 or TG101348, with those injected with an equal volume of normal saline as control group, and the repair of fracture ends was evaluated via X-ray examination 3 weeks later. The BMMSCs isolated in vitro grew well, and flow cytometry assay results showed that the positive expression rates of cluster of differentiation (CD)90 and CD105 in cells were 99.21 and 99.56%, respectively, with no CD45 expressed. According to the results of CCK-8 assay, TG101348 lowered the proliferation level of BMMSCs, and the wound healing assay revealed that the migration ability of BMMSCs at 24 and 48 h was substantially weaker than that in control group (P<0.05). After induction of osteogenic differentiation for 3 weeks, alizarin red staining results manifested that osteogenic induction group had notably more calcium nodules than TG101348 group (P<0.05). Compared with those in control group, the protein expression levels of p-JAK2 and p-STAT3 were remarkably raised by osteogenic induction (P<0.05), but the protein expression levels of JAK2, p-JAK2 and p-STAT3 were considerably decreased by TG101348 (P<0.05). It was found through the X-ray examination that the rabbits in control group and those injected with BMMSCs recovered well, and the latter had larger external calluses at fracture ends than the former, while the fracture ends of those injected with TG101348 and BMMSCs + TG101348 were not healed completely. BMMSCs promote healing of rabbit tibial fractures through the JAK-STAT signaling pathway.

FRACTURE HEALING SIMULATION REGULATED BY FORCE AND OXYGEN

According to the mechanical conditions of fracture fixation and the oxygen levels in the tissues, a simulation model of fracture healing process was built to describe the relationship among mechanical stability, oxygen levels in tissues and tissue differentiation during the second fracture healing. Different from the previous simulation model, in this paper, we took the three-dimensional model as the research object, solved the mechanical stimulation by finite element method, established the partial differential equation to solve the spatial and temporal variation of the oxygen in tissues. The process of tissue differentiation was described by fuzzy control method. The initial stage of fracture healing, intramembranous ossification, chondrogenesis, cartilage calcification and endochondral ossification during the fracture healing process were simulated, and the properties of tissue materials were continuously updated to complete the iterative process. The simulation program of fracture healing process was independently developed in Eclipse environment, and the simulation results were compared with experimental data and those of other fracture healing simulation models to verify the simulation program in this paper. Finally, the processes of transverse fracture healing in rats with different axial stability under normoxic, hypoxic and hyperoxic conditions was simulated, and the effects of different tissue oxygen levels and interosseous stabilities on fracture healing were analyzed. It is concluded by simulation that the delayed healing or non-union of bone will occur when in state of tissue hypoxia or interosseous instability, normal healing will occur when in state of tissue normoxia, and the healing will be accelerated when in state of tissue hyperoxia.

The Size of Intramedullary Fixation Affects Endochondral-Mediated Angiogenesis During Fracture Repair

OBJECTIVES

To explore the effect of intramedullary pin size on the biology of a healing fracture, specifically endochondral angiogenesis. We hypothesized that fracture fixation with a smaller pin would permit greater interfragmentary strain resulting in increased total amount of vascular endothelial growth factor within the callus and greater angiogenesis compared to fixation with a larger pin.

METHODS

Transverse mid-shaft femur fractures in 8-week-old mice were fixed with either a 23-gauge (G) or 30-G pin. Differences in interfragmentary strain at the fracture site were estimated between cohorts. A combination of histology, gene expression, serial radiography, and microcomputed tomography with and without vascular contrast agent were used to assess fracture healing and vascularity for each cohort.

RESULTS

Larger soft-tissue callus formation increased vascular endothelial growth factor-A expression, and a corresponding increase in vascular volume was observed in the higher strain, 30-G cohort. Radiographic analysis demonstrated earlier hard callus formation with greater initial interfragmentary strain, similar rates of union between pin size cohorts, yet delayed callus remodeling in mice with the larger pin size.

CONCLUSIONS

These findings suggest that the stability conferred by an intramedullary nail influences endochondral angiogenesis at the fracture.

Exosomes in the Repair of Bone Defects: Next-Generation Therapeutic Tools for the Treatment of Nonunion

Nonunion with bone defects, a common complication after long bone fracture, is a major challenge for orthopaedic surgeons worldwide because of the high incidence rate and difficulties in achieving successful treatment. Bone defects are the main complications of nonunion. The conventional biological treatments for nonunion with bone defects involve the use of autologous bone grafts or bone graft substitutes and cell-based therapy. Traditional nonunion treatments have always been associated with safety issues and various other complications. Bone grafts have limited autologous cancellous bone and there is a risk of infection. Additionally, problems with bone graft substitutes, including rejection and stimulation of bone formation, have been noted, and the health of the stem cell niche is a major consideration in cell-based therapy. In recent years, researchers have found that exosomes can be used to deliver functional RNA and mediate cell-to-cell communication, suggesting that exosomes may repair bone defects by regulating cells and cytokines involved in bone metabolism. In this review, we highlight the possible relationships between risk factors for nonunion and exosomes. Additionally, we discuss the roles of exosomes in bone metabolism and bone regeneration.

Limiting angiogenesis to modulate scar formation

Angiogenesis, the process of new blood vessel formation from existing blood vessels, is a key aspect of virtually every repair process. During wound healing an extensive, but immature and leaky vascular plexus forms which is subsequently reduced by regression of non-functional vessels. More recent studies indicate that uncontrolled vessel growth or impaired vessel regression as a consequence of an excessive inflammatory response can impair wound healing, resulting in scarring and dysfunction. However, in order to elucidate targetable factors to promote functional tissue regeneration we need to understand the molecular and cellular underpinnings of physiological angiogenesis, ranging from induction to resolution of blood vessels. Especially for avascular tissues (e.g. cornea, tendon, ligament, cartilage, etc.), limiting rather than boosting vessel growth during wound repair potentially is beneficial to restore full tissue function and may result in favourable long-term healing outcomes.

The Effects of Dynamic Loading on Bone Fracture Healing Under Ilizarov Circular Fixators

Early weight bearing appears to enhance bone fracture healing under Ilizarov circular fixators (ICFs). However, the role of early weight bearing in the healing process remains unclear. This study aims to provide insights into the effects of early weight bearing on healing of bone fractures stabilized with ICFs, with the aid of mathematical modeling. A computational model of fracture site was developed using poro-elastic formulation to simulate the transport of mesenchymal stem cells (MSCs), fibroblasts, chondrocytes, osteoblasts, osteogenic growth factor (OGF), and chondrogenic growth factor (CGF) and MSC differentiation during the early stage of healing, under various combinations of fracture gap sizes (GS), ICF wire pretension forces, and axial loads. 1 h of physiologically relevant cyclic axial loading followed by 23 h of rest in the post-inflammation phase (i.e., callus with granulation tissue) was simulated. The results show that physiologically relevant dynamic loading could significantly enhance cell and growth factor concentrations in the fracture site in a time and spatially dependent manner. 1 h cyclic loading (axial load with amplitude, PA, of 200 N at 1 Hz) increased the content of chondrocytes up to 37% (in all zones of callus), CGF up to 28% (in endosteal and periosteal callus) and OGF up to 50% (in endosteal and cortical callus) by the end of the 24 h period simulated. This suggests that the synergistic effect of dynamic loading-induced advective transport and mechanical stimuli due to early weight bearing is likely to enhance secondary healing. Furthermore, the study suggests that relatively higher PA values or lower ICF wire pretension forces or smaller GS could result in increased chondrocyte and GF content within the callus.

Effect of ultrasound on bone fracture healing: A computational mechanobioregulatory model

Bone healing process is a complicated phenomenon regulated by biochemical and mechanical signals. Experimental studies have shown that ultrasound (US) accelerates bone ossification and has a multiple influence on cell differentiation and angiogenesis. In a recent work of the authors, a bioregulatory model for providing bone-healing predictions was addressed, taking into account for the first time the salutary effect of US on the involved angiogenesis. In the present work, a mechanobioregulatory model of bone solidification under the US presence incorporating also the mechanical environment on the regeneration process, which is known to affect cellular processes, is presented. An iterative procedure is adopted, where the finite element method is employed to compute the mechanical stimuli at the linear elastic phases of the poroelastic callus region and a coupled system of partial differential equations to simulate the enhancement by the US cell angiogenesis process and thus the oxygen concentration in the fractured area. Numerical simulations with and without the presence of US that illustrate the influence of progenitor cells' origin in the healing pattern and the healing rate and simultaneously demonstrate the salutary effect of US on bone repair are presented and discussed.

Differential osteopontin expression in human osteoblasts derived from iliac crest and alveolar bone and its role in early stages of angiogenesis

In our previous study, we revealed significant differences of osteopontin (OPN) gene expression in primary human osteoblasts (HOBs) derived from iliac crest bone (iHOBs) and alveolar bone (aHOBs). The present study aims at assigning this discriminative expression to a possible biologic function. OPN is known to be involved in several pathologic and physiologic processes, among others angiogenesis. Therefore, we studied the reaction of human umbilical vein endothelial cells (HUVECs) to HOB-derived OPN regarding angiogenesis. To this end, human primary explant cultures of both bone entities from ten donors were established. Subsequent transcription analysis detected higher gene expression of OPN in iHOBs compared to aHOBs, thereby confirming the results of our previous study. This difference was particularly apparent when cultures were derived from female donors. Hence, OPN protein expression as well as the angiogenic potential of OPN was analyzed, originating from HOBs of one female donor. In accordance to the gene expression level, secreted OPN was more abundant in the supernatant of iHOBs than in aHOBs. Moreover, secreted OPN was found to stimulate migration of HUVECs, but not proliferation or tube formation. These results indicate an involvement in very early stages of angiogenesis and a functional distinction of OPN from HOBs derived from different bone entities.

Bone Fracture Acute Phase Response-A Unifying Theory of Fracture Repair: Clinical and Scientific Implications

Bone fractures create five problems that must be resolved: bleeding, risk of infection, hypoxia, disproportionate strain, and inability to bear weight. There have been enormous advancements in our understanding of the molecular mechanisms that resolve these problems after fractures, and in best clinical practices of repairing fractures. We put forth a modern, comprehensive model of fracture repair that synthesizes the literature on the biology and biomechanics of fracture repair to address the primary problems of fractures. This updated model is a framework for both fracture management and future studies aimed at understanding and treating this complex process. This model is based upon the fracture acute phase response (APR), which encompasses the molecular mechanisms that respond to injury. The APR is divided into sequential stages of "survival" and "repair." Early in convalescence, during "survival," bleeding and infection are resolved by collaborative efforts of the hemostatic and inflammatory pathways. Later, in "repair," avascular and biomechanically insufficient bone is replaced by a variable combination of intramembranous and endochondral ossification. Progression to repair cannot occur until survival has been ensured. A disproportionate APR-either insufficient or exuberant-leads to complications of survival (hemorrhage, thrombosis, systemic inflammatory response syndrome, infection, death) and/or repair (delayed- or non-union). The type of ossification utilized for fracture repair is dependent on the relative amounts of strain and vascularity in the fracture microenvironment, but any failure along this process can disrupt or delay fracture healing and result in a similar non-union. Therefore, incomplete understanding of the principles herein can result in mismanagement of fracture care or application of hardware that interferes with fracture repair. This unifying model of fracture repair not only informs clinicians how their interventions fit within the framework of normal biological healing but also instructs investigators about the critical variables and outputs to assess during a study of fracture repair.

Comparative evaluation of the osteogenic capacity of human mesenchymal stem cells from bone marrow and umbilical cord tissue

Mesenchymal stem cells (MSCs) have been extensively investigated in the field of regenerative medicine. Human bone MSCs (BMSCs) have become a common type of seed cell for bone tissue engineering. However, the viability and cell number of BMSCs are negatively correlated with donor age, and as the extraction process is painful, this method has not been widely used. As human umbilical cord MSCs (UCMSCs) may be harvested inexpensively and inexhaustibly, the present study evaluated and compared the regenerative potential of UCMSCs and BMSCs to determine whether UCMSCs may be used as a novel cell type for bone regeneration. In the present study, the proliferation and osteogenic capacity of BMSCs and UCMSCs was compared in vitro. BMSCs and UCMSCs were respectively combined with biofunctionalized macroporous calcium phosphate cement, and their bone regenerative potentials were determined by investigating their capacity for ectopic bone formation in a nude mouse model as well as their efficacy in a rat model of tibia bone defect. The extent of bone regeneration was examined by X-ray, histological and immunohistochemical analyses. The results revealed that UCMSCs exhibited a good osteogenic differentiation potential, similarly to that of BMSCs, and that UCMSCs were able to contribute to the regeneration of bone and blood vessels. Furthermore, no significant differences were identified between BMSCs and UCMSCs in terms of their bone regenerative effect.

Inflammatory-Driven Angiogenesis in Bone Augmentation with Bovine Hydroxyapatite, B-Tricalcium Phosphate, and Bioglasses: A Comparative Study

Introduction

The clinical use of bioactive materials for bone augmentation has remained a challenge because of predictability and effectiveness concerns, as well as increased costs. The purpose of this study was to analyse the ability to integrate bone substitutes by evaluating the immunohistochemical expression of the platelet endothelial cell adhesion molecules, vascular endothelial growth factor, collagen IV, laminin, and osteonectin, in the vicinity of bone grafts, enabling tissue revascularization and appearance of bone lamellae. There is a lack of in vivo studies of inflammatory-driven angiogenesis in bone engineering using various grafts.

Methods

The study was performed in animal experimental model on the standardized monocortical defects in the tibia of 20 New Zealand rabbits. The defects were augmented with three types of bone substituents. The used bone substituents were beta-tricalcium phosphate, bovine hydroxyapatite, and bioactive glasses. After a period of 6 months, bone fragments were harvested for histopathologic examination. Endothelial cell analysis was done by analysing vascularization with PECAM/CD31 and VEGF and fibrosis with collagen IV, laminin, and osteonectin stains. Statistical analysis was realized by descriptive analysis which was completed with the kurtosis and skewness as well as the Kruskal-Wallis and Mann-Whitney statistical tests.

Results

The discoveries show that the amount of bone that is formed around beta-tricalcium phosphate and bovine hydroxyapatite is clearly superior to the bioactive glasses. Both the lumen diameter and the number of vessels were slightly increased in favor of beta-tricalcium phosphate.

Conclusion

We can conclude that bone substitutes as bovine bone and beta-tricalcium phosphate have significant increased angiogenesis (and subsequent improved osteogenesis) compared to the bioactive glass. In our study, significant angiogenesis is linked with a greater tissue formation, indicating that in bone engineering with the allografts we used, inflammation has more benefic effects, the catabolic action being exceeded by the tissue formation.

Titanium-based implant comprising a porous microstructure assembled with nanoleaves and controllable silicon-ion release for enhanced osseointegration

Osseointegration is crucial for early fixation as well as for long-term implantation success, hence numerous efforts have been made to tune the surface topography or chemical composition of biomedical implants to improve osseointegration. In this work, various nanostructures, including nanoflocs, nanobundles, nanorods, and nanoleaves, were introduced to the surface of silicon (Si)-incorporated microporous structure to form Si-incorporated micro-nano hierarchical structures on titanium (Ti)-based implants. The osseointegration of the implants were systemically assessed in vivo and in vitro. The in vitro evaluations showed that the nanostructures promoted the protein adsorption, thus modulating the early cellular responses, including the attachment and spreading of osteoblasts and human endothelial cells (HUVECs), and subsequent cell proliferation and differentiation. Furthermore, compared with the single microporous structure, the nanostructures located over the microporous structure protected the Si ions from quick release and allowed the long-term sustained Si-ions release, which further contributed to the proliferation and differentiation of osteoblasts and vascular endothelial growth factor (VEGF) secretion as well as the tube formation of HUVECs. Collectively, the favorable nano-surface structures, especially the nanoleaves structure, and the constant Si-ion release together led to robust osteogenic and angiogenic activities. More importantly, in vivo micro-CT evaluation and histological observations further verified that the Si-incorporated micro-nano hierarchical implant with nanoleaves structure could efficiently promote new bone formation, thus indicating it was an attractive candidate as a next-generation bone-implant material.

Engineering vascularized and innervated bone biomaterials for improved skeletal tissue regeneration

Blood vessels and nerve fibers are distributed throughout the entirety of skeletal tissue, and play important roles during bone development and fracture healing by supplying oxygen, nutrients, and cells. However, despite the successful development of bone mimetic materials that can replace damaged bone from a structural point of view, most of the available bone biomaterials often do not induce sufficient formation of blood vessels and nerves. In part, this is due to the difficulty of integrating and regulating multiple tissue types within artificial materials, which causes a gap between native skeletal tissue. Therefore, understanding the anatomy and underlying interaction mechanisms of blood vessels and nerve fibers in skeletal tissue is important to develop biomaterials that can recapitulate its complex microenvironment. In this perspective, we highlight the structure and osteogenic functions of the vascular and nervous system in bone, in a coupled manner. In addition, we discuss important design criteria for engineering vascularized, innervated, and neurovascularized bone implant materials, as well as recent advances in the development of such biomaterials. We expect that bone implant materials with neurovascularized networks can more accurately mimic native skeletal tissue and improve the regeneration of bone tissue.

Three-dimensional computational model simulating the fracture healing process with both biphasic poroelastic finite element analysis and fuzzy logic control

A dynamic model regulated by both biphasic poroelastic finite element analysis and fuzzy logic control was established. Fuzzy logic control was an easy and comprehensive way to simulate the tissue differentiation process, and it is convenient for researchers and medical experts to communicate with one another to change the fuzzy logic rules and improve the simulation of the tissue differentiation process. In this study, a three-dimensional fracture healing model with two different interfragmentary movements (case A: 0.25 mm and case B: 1.25 mm) was analysed with the new set-up computational model. As the healing process proceeded, both simulated interfragmentary movements predicted a decrease and the time that the decrease started for case B was later than that for case A. Compared with experimental results, both cases corresponded with experimental data well. The newly established dynamic model can simulate the healing process under different mechanical environments and has the potential to extend to the multiscale healing model, which is essential for reducing the animal experiments and helping to characterise the complex dynamic interaction between tissue differentiations within the callus region.