Oxygen as a critical determinant of bone fracture healing—A multiscale model

Mini-bones: miniaturized bone in vitro models

In bone tissue engineering (TE) and regeneration, miniaturized, (sub)millimeter-sized bone models have become a popular trend since they bring about physiological biomimicry, precise orchestration of concurrent stimuli, and compatibility with high-throughput setups and high-content imaging. They also allow efficient use of cells, reagents, materials, and energy. In this review, we describe the state of the art of miniaturized in vitro bone models, or 'mini-bones', describing these models based on their characteristics of (multi)cellularity and engineered extracellular matrix (ECM), and elaborating on miniaturization approaches and fabrication techniques. We analyze the performance of 'mini-bone' models according to their applications for studying basic bone biology or as regeneration models, disease models, and screening platforms, and provide an outlook on future trends, challenges, and opportunities.

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.

Taohong Siwu Decoction Promotes Osteo-Angiogenesis in Fractures by Regulating the HIF-1α Signaling Pathway

Background

Vascular damage is a major consequence of bone fracture. Taohong Siwu decoction (TSD) can raise the expression of vascular endothelial growth factor (VEGF) in fracture healing. However, its molecular mechanism in promoting angiogenesis is still unknown. The aim of this study was to investigate the potential mechanisms of TSD in the regulation of osteo-angiogenesis in fracture healing.

Methods

A rat tibial fracture model was established. After low- (4.5 g·kg⁻¹), medium- (9 g·kg⁻¹), and high-dose TSD (18 g·kg⁻¹) and panax notoginsenoside (25 mg kg⁻¹) treatment, hematoxylin-eosin staining was employed to visualize pathological changes in bone tissues. The levels of cytokines (interleukin (IL)-2, tumor necrosis factor-α (TNF-α), IL-6, and IL-1β), thromboxane B2 (TXB2), and 6 ketone prostaglandin F1α (6-Keto-PGF1α) were quantified by enzyme-linked immunosorbent assay (ELISA). Immunofluorescence was used to identify the rat aortic endothelial cells (RAECs). Control serum, 10% TSD-containing serum, and 10% TSD-containing serum combined with hypoxia-inducible factor-1α (HIF-1α) inhibitor were used to treat the RAECs and rat osteoblasts. Transwell migration assay was utilized to examine the migration of the RAECs. The Matrigel tubulogenesis assay was used for the assessment of angiogenesis. The expression of angiogenesis- (von Hippel-Lindau tumor suppressor (VHL), HIF-1α, VEGF, angiopoietin-2 (Ang-2), and pVHL) and osteogenesis-related (alkaline phosphatase (ALP), runt-related transcription factor 2 (Runx2), and osteopontin-1 (OPN-1)) protein and gene was detected by western blot and quantitative real-time PCR (qRT-PCR).

Results

Compared with the model group, TSD increased the trabecular bone areas, numbers, and thicknesses in fractured rats. In the plasma, the levels of cytokines and TXB2 in the middle- and high-dose TSD group were significantly lower than those in the model group (P < 0.01). The 6-keto-PGF1α content was increased by middle- and high-dose TSD intervention (P < 0.01). Compared to the control serum group, the angiogenesis and migration of the RAECs were enhanced in the TSD group (P < 0.001). The expression of HIF-1α, VEGF, and Ang-2 in the TSD group upregulated significantly (P < 0.001). VHL and pVHL were inhibited under TSD-containing serum treatment (P < 0.001). ALP, Runx2, and OPN-1 were increased obviously in the TSD group (P < 0.001). Nevertheless, the HIF-1α inhibitor reversed these changes (P < 0.001).

Conclusion

TSD promotes angiogenesis and osteogenesis by regulating the HIF-1α signaling pathway. Meanwhile, it can effectively reduce the risk of inflammation and improve blood circulation.

The vascularization paradox of non-union formation

Despite major research efforts to elucidate mechanisms of non-union formation, failed fracture healing remains a common complication in orthopedic surgery. Adequate vascularization has been recognized as a crucial factor for successful bone regeneration, as newly formed microvessels guarantee the supply of the callus tissue with vital oxygen, nutrients, and growth factors. Accordingly, a vast number of preclinical studies have focused on the development of vascularization strategies to stimulate fracture repair. However, recent evidence suggests that stimulation of blood vessel formation is an oversimplified approach to support bone regeneration. This review discusses the role of vascularization during bone regeneration and delineates a phenomenon, for which we coin the term “the vascularization paradox of non-union-formation”. This view is based on the results of a variety of experimental studies that suggest that the callus tissue of non-unions is indeed densely vascularized and that pro-angiogenic mediators, such as vascular endothelial growth factor, are sufficiently expressed at the facture site. By gaining further insights into the molecular and cellular basis of non-union vascularization, it may be possible to develop more optimized treatment approaches or even prevent the non-union formation in the future.

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.

Self-Oxygenation of Tissues Orchestrates Full-Thickness Vascularization of Living Implants

Bioengineering of tissues and organs has the potential to generate functional replacement organs. However, achieving the full-thickness vascularization that is required for long-term survival of living implants has remained a grand challenge, especially for clinically sized implants. During the pre-vascular phase, implanted engineered tissues are forced to metabolically rely on the diffusion of nutrients from adjacent host-tissue, which for larger living implants results in anoxia, cell death, and ultimately implant failure. Here it is reported that this challenge can be addressed by engineering self-oxygenating tissues, which is achieved via the incorporation of hydrophobic oxygen-generating micromaterials into engineered tissues. Self-oxygenation of tissues transforms anoxic stresses into hypoxic stimulation in a homogenous and tissue size-independent manner. The in situ elevation of oxygen tension enables the sustained production of high quantities of angiogenic factors by implanted cells, which are offered a metabolically protected pro-angiogenic microenvironment. Numerical simulations predict that self-oxygenation of living tissues will effectively orchestrate rapid full-thickness vascularization of implanted tissues, which is empirically confirmed via in vivo experimentation. Self-oxygenation of tissues thus represents a novel, effective, and widely applicable strategy to enable the vascularization living implants, which is expected to advance organ transplantation and regenerative medicine applications.

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.

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%).

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.

Fibrous Demineralized Bone Matrix (DBM) Improves Bone Marrow Mononuclear Cell (BMC)-Supported Bone Healing in Large Femoral Bone Defects in Rats

Regeneration of large bone defects is a major objective in trauma surgery. Bone marrow mononuclear cell (BMC)-supported bone healing was shown to be efficient after immobilization on a scaffold. We hypothesized that fibrous demineralized bone matrix (DBM) in various forms with BMCs is superior to granular DBM. A total of 65 male SD rats were assigned to five treatment groups: syngenic cancellous bone (SCB), fibrous demineralized bone matrix (f-DBM), fibrous demineralized bone matrix densely packed (f-DBM 120%), DBM granules (GDBM) and DBM granules 5% calcium phosphate (GDBM5%Ca2+). BMCs from donor rats were combined with different scaffolds and placed into 5 mm femoral bone defects. After 8 weeks, bone mineral density (BMD), biomechanical stability and histology were assessed. Similar biomechanical properties of f-DBM and SCB defects were observed. Similar bone and cartilage formation was found in all groups, but a significantly bigger residual defect size was found in GDBM. High bone healing scores were found in f-DBM (25) and SCB (25). The application of DBM in fiber form combined with the application of BMCs shows promising results comparable to the gold standard, syngenic cancellous bone. Denser packing of fibers or higher amount of calcium phosphate has no positive effect.

Impaired Callus Formation in Pathological Mandibular Fractures in Medication-Related Osteonecrosis of the Jaw and Osteoradionecrosis

PURPOSE

Nonsurgical treatment of mandibular fractures secondary to medication-related osteonecrosis of the jaw (MRONJ) or osteoradionecrosis (ORN) mostly results in nonunion, whereas nonsurgical fracture treatment of atrophic fractures can achieve favorable results in selected cases. The aim of this study was to compare callus formation in pathological mandibular fractures due to MRONJ, ORN, or extreme mandibular atrophy.

METHODS

A retrospective cohort study reviewing the medical records of all MRONJ-, ORN-, or atrophy-related fractures treated at the departments of maxillofacial surgery in the Leuven or Lille university hospitals between 2010 and 2019 was undertaken. The primary predictor variable in this study was disease state (MRONJ, ORN, or extreme mandibular atrophy). The primary outcome variable was callus formation after 1 month of follow-up (present, absent). Additional study variables measured included patient age and gender. T-tests, Fisher exact tests, and multiple logistic regression were used for statistical analysis. The significance level was set at P < .05.

RESULTS

Seventy patients were analyzed (12 MRONJ cases, 54 ORN fractures, 4 atrophic fractures). The callus formation prevalence in nonsurgically approached fractures secondary to ORN and MRONJ after 1 month of follow-up was 3.03% (2/66 cases). In contrast, callus was detected in all patients in the mandibular atrophy-related fracture group. Osteonecrosis was statistically correlated with impaired callus formation (P = .0121).

CONCLUSION

Whereas one would expect indirect fracture healing and thus callus formation to occur in all non-surgically treated fractures, our data demonstrate its absence in the majority of MRONJ- and ORN-related fractures. Multiple plausible explanations for this phenomenon were identified: periosteal damage with loss of callus-forming cells, compromised vasculature, and bacterial colonization.

Development of a tibial experimental non-union model in rats

BACKGROUND

Many non-union animal models have been developed to explore the problems surrounding fracture healing. However, the existing models are not perfect and cannot satisfy all non-union studies. This study aimed to make a non-union model of the tibia in rats by cauterization of the posterior of 2 mm on both sides of the fracture end after open osteotomy of the tibia and fixing the fractured tibia with a Kirschner wire 0.8 mm in diameter.

METHODS

For this study, 96 female adult Sprague-Dawley (SD) rats were used. The rats underwent surgery to produce a tibial open fracture and were fixed with a 0.8-mm diameter Kirschner wire. In 48 of the rats, the periosteum proximal and distal to the fracture end was cauterized.

RESULTS

At 2, 4, 6, and 8 weeks after surgery, radiological and histological analysis showed typical physiological healing in the control group, and the healing rate was 100% at 6 weeks. But the non-union group was characterized by resorption of the fracture ends with few callus formations and no bridging callus formation, and the healing rate was 0% at 8 weeks.

CONCLUSIONS

This method represents a reproducible model to create atrophic non-unions. This model provides a new option for studying the basic healing mechanisms and evaluating new therapies for bone regeneration and treatment of non-unions.

Fibrodysplasia Ossificans Progressiva: What Have We Achieved and Where Are We Now? Follow-up to the 2015 Lorentz Workshop

Fibrodysplasia ossificans progressiva (FOP) is an ultra-rare progressive genetic disease effecting one in a million individuals. During their life, patients with FOP progressively develop bone in the soft tissues resulting in increasing immobility and early death. A mutation in the ACVR1 gene was identified as the causative mutation of FOP in 2006. After this, the pathophysiology of FOP has been further elucidated through the efforts of research groups worldwide. In 2015, a workshop was held to gather these groups and discuss the new challenges in FOP research. Here we present an overview and update on these topics.

Inhibition of the mitochondrial permeability transition improves bone fracture repair

Bone fracture is accompanied by trauma, mechanical stresses, and inflammation - conditions known to induce the mitochondrial permeability transition. This phenomenon occurs due to opening of the mitochondrial permeability transition pore (MPTP) promoted by cyclophilin D (CypD). MPTP opening leads to more inflammation, cell death and potentially to disruption of fracture repair. Here we performed a proof-of-concept study and tested a hypothesis that protecting mitochondria from MPTP opening via inhibition of CypD improves fracture repair. First, our in vitro experiments indicated pro-osteogenic and anti-inflammatory effects in osteoprogenitors upon CypD knock-out or pharmacological inhibition. Using a bone fracture model in mice, we observed that bone formation and biomechanical properties of repaired bones were significantly increased in CypD knock-out mice or wild type mice treated with a CypD inhibitor, NIM811, when compared to controls. These effects were evident in young male but not female mice, however in older (13 month-old) female mice bone formation was also increased during fracture repair. In contrast to global CypD knock-out, mesenchymal lineage-specific (Prx1-Cre driven) CypD deletion did not result in improved fracture repair. Our findings implicate MPTP in bone fracture and suggest systemic CypD inhibition as a modality to promote fracture repair.

Preparation and characteristics of a novel oxygen-releasing coating for improved cell responses in hypoxic environment

Affected by environmental factors such as oxygen deficiency, the secretion of growth factor was abnormal in bone injury sites, resulting in the poor responses of osteoblasts and prolonging the healing process. Herein, in this study, we reported an in situ oxygen-releasing porous titanium coating that combines the dual degradability of poly(lactic-co-glycolic acid) with the self-releasing oxygen capacity of the CaO₂ core. The resulting formulation exhibited stable oxygen-releasing capacity as well as the ability to promote proliferation and differentiation of the MC3T3 cell line under hypoxia conditions. According to these results, oxygen-releasing coatings based on improved cellular microenvironment may be a promising bone repair material that would reduce the incidence of difficult bone healing in the future.

Oxygen-releasing nanofibers for breathable bone tissue engineering application

Oxygen is a vital molecule for cell and tissue processes. Electrospun fibers have been extensively used as drug loading carriers due to possibility of well control over drug release with modulating fiber properties. However, they have not been used as depots for oxygen release. In the present study, an oxygen-releasing nanofibrous scaffold has been developed by electrospinning of polylactic acid/nano-calcium peroxide suspension with different polylactic acid concentrations (6.5 and 13% w/v). The electrospun fibers with calcium peroxide cargo provided oxygen content of 30–94 mmHg in a period of 14 days which lies well within the oxygen level of osseous tissue. The release profile of 13% polylactic acid fibers was different with that of 6.5% fibers with respects to the initial content of released oxygen and the release rate. Not only did 13% fibers supply oxygen with a slower rate, but also they resulted in a lower burst release of oxygen. Cell culture studies in hypoxia corroborated that 13% polylactic acid fibers better preserve cell viability comparing 6.5% counterparts as perceived by MTT assay. Moreover, they endowed more favored milieu for adherence, arrangement and migration of mesenchymal stem cells as confirmed by microscopy images. The oxygen-releasing fibers equally affected alkaline phosphatase, osteocalcin, and calcium deposition by mesenchymal stem cells most likely due to interplay between topographical and metabolic cues offered by 6.5 and 13% formulations.

Lipid availability determines fate of skeletal progenitor cells via SOX9

The avascular nature of cartilage makes it a unique tissue^1–4, but whether and how the absence of nutrient supply regulates chondrogenesis remains unknown. Here, we show that obstruction of vascular invasion during bone healing favours chondrogenic over osteogenic differentiation of skeletal progenitor cells. Unexpectedly, this process is driven by a decreased availability of extracellular lipids. When lipids are scarce, skeletal progenitors activate FoxO transcription factors, which bind to the Sox9 promoter and increase its expression. Besides initiating chondrogenesis, SOX9 acts as a regulator of cellular metabolism by suppressing fatty acid oxidation, and thus adapts the cells to an avascular life. Our results define lipid scarcity as an important determinant of chondrogenic commitment, reveal a role for FoxOs during lipid starvation, and identify SOX9 as a critical metabolic mediator. These data highlight the importance of the nutritional microenvironment in the specification of skeletal cell fate.

PREDICTION OF THE TREND OF BONE FRACTURE HEALING BASED ON THE RESULTS OF THE EARLY STAGES SIMULATIONS: A FINITE ELEMENT STUDY

To date, several studies have implied the importance of early stage mechanical stability in the bone fracture healing process. This study aimed at finding a correlation between the predicted different tissue phenotypes in the early stages of healing and the ultimate healing outcome. For this purpose, the process of fracture healing was numerically simulated employing an axisymmetric bi-phasic finite element (FE) model for three initial gap sizes of 1, 3 and 6[Formula: see text]mm and four initial interfragmentary strains (IFS) of 7%, 11%, 15% and 19%. The model was validated with experimental and other numerical studies from the literature. Results of this study showed that the amount of cartilage and fibrous tissue observed in the early stage after fracture can be used to qualitatively assess the outcome of complete bone healing process. Greater amount of cartilage in early stage of healing process yielded faster callus maturation, and delayed maturation of callus was predicted in the case of high fibrous tissue production. Results of this study can be used to provide an estimation of the performance of different fixation systems by considering the amounts of cartilage and fibrous tissues observed in the early stage of healing.

Computational Modeling to Quantify the Contributions of VEGFR1, VEGFR2, and Lateral Inhibition in Sprouting Angiogenesis

Sprouting angiogenesis is a necessary process in regeneration and development as well as in tumorigenesis. VEGF-A is the main pro-angiogenic chemoattractant and it can bind to the decoy receptor VEGFR1 or to VEGFR2 to induce sprouting. Active sprout cells express Dll4, which binds to Notch1 on neighboring cells, in turn inhibiting VEGFR2 expression. It is known that the balance between VEGFR2 and VEGFR1 determines tip selection and network architecture, however the quantitative interrelationship of the receptors and their interrelated balances, also with relation to Dll4-Notch1 signaling, remains yet largely unknown. Here, we present an agent-based computer model of sprouting angiogenesis, integrating VEGFR1 and VEGFR2 in a detailed model of cellular signaling. Our model reproduces experimental data on VEGFR1 knockout. We show that soluble VEGFR1 improves the efficiency of angiogenesis by directing sprouts away from existing cells over a wide range of parameters. Our analysis unravels the relevance of the stability of the active notch intracellular domain as a dominating hub in this regulatory network. Our analysis quantitatively dissects the regulatory interactions in sprouting angiogenesis. Because we use a detailed model of intracellular signaling, the results of our analysis are directly linked to biological entities. We provide our computational model and simulation engine for integration in complementary modeling approaches.

Engineering of gradient osteochondral tissue: From nature to lab

The osteochondral tissue is an interface between two distinct tissues: articular cartilage and bone. These two tissues are significantly diverse with regard to their chemical compositions, mechanical properties, structure, electrical properties, and the amount of nutrient and oxygen consumption. Thus, transition from the surface of the articular cartilage to the subchondral bone needs to face several smooth gradients. These gradients are imperative to study to generate a scaffold suitable for the reconstruction of the cartilaginous and osseous layers of a defected osteochondral tissue, simultaneously. The aim of this review is to peruse the alternation of biochemical, biomechanical, structural, electrical, and metabolic properties of the osteochondral tissue moving from the surface of the articular cartilage to the subchondral bone. Moreover, this review also discusses currently developed approaches and ideal techniques with a focus on gradients present in the interface of the cartilage and bone. STATEMENT OF SIGNIFICANCE: The submitted review paper entitled as "Engineering of the gradient osteochondral tissue: from nature to lab" is a complete review with regard to the osteochondral tissue and transition of different properties between the cartilage and bone tissues. Moreover, previous studies on the osteochondral tissue engineering have been reviewed in this paper. This complete information can be a valuable and useful source for current and future researchers and scientists. Considering the scope of the submitted paper, Acta Biomaterialia would be a suitable journal for publishing this article.

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.

Modeling the Macrophage-Mediated Inflammation Involved in the Bone Fracture Healing Process

A new mathematical model is presented to study the effects of macrophages on the bone fracture healing process. The model consists of a system of nonlinear ordinary differential equations that represents the interactions among classically and alternatively activated macrophages, mesenchymal stem cells, osteoblasts, and pro- and anti-inflammatory cytokines. A qualitative analysis of the model is performed to determine the equilibria and their corresponding stability properties. Numerical simulations are also presented to support the theoretical results, and to monitor the evolution of a broken bone for different types of fractures under various medical interventions. The model can be used to guide clinical experiments and to explore possible medical treatments that accelerate the bone fracture healing process, either by surgical interventions or drug administrations.

Current state of fabrication technologies and materials for bone tissue engineering

A range of traditional and free-form fabrication technologies have been investigated and, in numerous occasions, commercialized for use in the field of regenerative tissue engineering (TE). The demand for technologies capable of treating bone defects inherently difficult to repair has been on the rise. This quest, accompanied by the advent of functionally tailored, biocompatible, and biodegradable materials, has garnered an enormous research interest in bone TE. As a result, different materials and fabrication methods have been investigated towards this end, leading to a deeper understanding of the geometrical, mechanical and biological requirements associated with bone scaffolds. As our understanding of the scaffold requirements expands, so do the capability requirements of the fabrication processes. The goal of this review is to provide a broad examination of existing scaffold fabrication processes and highlight future trends in their development. To appreciate the clinical requirements of bone scaffolds, a brief review of the biological process by which bone regenerates itself is presented first. This is followed by a summary and comparisons of commonly used implant techniques to highlight the advantages of TE-based approaches over traditional grafting methods. A detailed discussion on the clinical and mechanical requirements of bone scaffolds then follows. The remainder of the manuscript is dedicated to current scaffold fabrication methods, their unique capabilities and perceived shortcomings. The range of biomaterials employed in each fabrication method is summarized. Selected traditional and non-traditional fabrication methods are discussed with a highlight on their future potential from the authors' perspective. This study is motivated by the rapidly growing demand for effective scaffold fabrication processes capable of economically producing constructs with intricate and precisely controlled internal and external architectures. STATEMENT OF SIGNIFICANCE: The manuscript summarizes the current state of fabrication technologies and materials used for creating scaffolds in bone tissue engineering applications. A comprehensive analysis of different fabrication methods (traditional and free-form) were summarized in this review paper, with emphasis on recent developments in the field. The fabrication techniques suitable for creating scaffolds for tissue engineering was particularly targeted and their use in bone tissue engineering were articulated. Along with the fabrication techniques, we emphasized the choice of materials in these processes. Considering the limitations of each process, we highlighted the materials and the material properties critical in that particular process and provided a brief rational for the choice of the materials. The functional performance for bone tissue engineering are summarized for different fabrication processes and the choice of biomaterials. Finally, we provide a perspective on the future of the field, highlighting the knowledge gaps and promising avenues in pursuit of effective scaffolds for bone tissue engineering. This extensive review of the field will provide research community with a reference source for current approaches to scaffold preparation. We hope to encourage the researchers to generate next generation biomaterials to be used in these fabrication processes. By providing both advantages and disadvantage of each fabrication method in detail, new fabrication techniques might be devised that will overcome the limitations of the current approaches. These studies should facilitate the efforts of researchers interested in generating ideal scaffolds, and should have applications beyond the repair of bone tissue.

A Review on Properties of Natural and Synthetic Based Electrospun Fibrous Materials for Bone Tissue Engineering

Bone tissue engineering is an interdisciplinary field where the principles of engineering are applied on bone-related biochemical reactions. Scaffolds, cells, growth factors, and their interrelation in microenvironment are the major concerns in bone tissue engineering. Among many alternatives, electrospinning is a promising and versatile technique that is used to fabricate polymer fibrous scaffolds for bone tissue engineering applications. Copolymerization and polymer blending is a promising strategic way in purpose of getting synergistic and additive effect achieved from either polymer. In this review, we summarize the basic chemistry of bone, principle of electrospinning, and polymers that are used in bone tissue engineering. Particular attention will be given on biomechanical properties and biological activities of these electrospun fibers. This review will cover the fundamental basis of cell adhesion, differentiation, and proliferation of the electrospun fibers in bone tissue scaffolds. In the last section, we offer the current development and future perspectives on the use of electrospun mats in bone tissue engineering.

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.

Biocompatibility and osteoconductivity of PLCL coated and noncoated xenografts: An in vitro and preclinical trial

BACKGROUND

Cells, scaffolds, and growth factors are the key components in bone tissue engineering. Scaffold composition, topography, and architecture influence the amount of regenerated bone in the implantation site. The aims of the study were to compare viability and proliferation of mesenchymal stem cells (MSCs) seeded onto two commercial xenografts: Bio-Oss (BO) and bioactive bone bovine (BB). Next, these materials were compared for histomorphometric bone formation in a socket preservation model in rats.

MATERIALS AND METHODS

MSCs were seeded onto monolayers of BO or BB granules. Cell viability and proliferation were evaluated after incubation of 0, 2, 20, and 48 h. A total of 24 Sprague Dawley rats underwent unilateral extraction of maxillary molars. Rats were randomly divided into three groups: natural healing (nongrafted socket) or socket preservation with either BO or BB. Rats were sacrificed after 8 weeks, and histomorphometric analysis was done to evaluate bone formation and residual scaffold at the extraction site.

RESULTS

Differences in the metabolic activity of MSCs that were seeded onto BO or BB was observed at 2 h after seeding: the metabolic activity was elevated compared to baseline in the BB (P = .046) and not changed in the BO wells (P = .84). After 20 h, the metabolic activity of MSCs seeded onto BO was decreasing (P = .005), while cell viability was not changed in the BB group (P = .356). Intergroup comparison revealed higher metabolic activity of MSCs seeded on BB after 48 h compared with BO (P = .016). The in vivo results demonstrated differences in socket healing between the groups: percentage of new bone was higher in the BB compared to BO group (39.1 ± 14.3 vs. 23.7 ± 10.8%, respectively, P = .096). Connective tissue portion was higher in the BO group compared with BB (73.7 ± 11.1 vs. 49.6 ± 13.7%, respectively, P = .018). Residual grafting martial was higher in the BB (11.34 ± 4.18 vs. 2.62 ± 1.23%, P = .011).

CONCLUSIONS

The results of this study demonstrating higher vitality and proliferation of MSCs seeded onto BB. Furthermore, following ridge preservation, higher percentage of new bone and lower residual scaffold were found in the BB compared with BO. This enhanced regenerative response might be the result of an enhancement of metabolic activity in cells attached to it. Further research will be needed to understand the precise mechanism.

In silico clinical trials for pediatric orphan diseases

To date poor treatment options are available for patients with congenital pseudarthrosis of the tibia (CPT), a pediatric orphan disease. In this study we have performed an in silico clinical trial on 200 virtual subjects, generated from a previously established model of murine bone regeneration, to tackle the challenges associated with the small, pediatric patient population. Each virtual subject was simulated to receive no treatment and bone morphogenetic protein (BMP) treatment. We have shown that the degree of severity of CPT is significantly reduced with BMP treatment, although the effect is highly subject-specific. Using machine learning techniques we were also able to stratify the virtual subject population in adverse responders, non-responders, responders and asymptomatic. In summary, this study shows the potential of in silico medicine technologies as well as their implications for other orphan diseases.

3D-printed biphasic calcium phosphate scaffolds coated with an oxygen generating system for enhancing engineered tissue survival

Tissue engineering scaffolds with oxygen generating elements have shown to be able to increase the level of oxygen and cell survivability in specific conditions. In this study, biphasic calcium phosphate (BCP) scaffolds with the composition of 60% hydroxyapatite (HA) and 40% beta-tricalcium phosphate (β-TCP), which have shown a great potential for bone tissue engineering applications, were fabricated by a direct-write assembly (robocasting) technique. Then, the three-dimensional (3D)-printed scaffolds were coated with different ratios of an oxygen releasing agent, calcium peroxide (CPO), which encapsulated within a polycaprolactone (PCL) matrix through dip-coating, and used for in situ production of oxygen in the implanted sites. The structure, composition and morphology of the prepared scaffolds were characterized by different techniques. The oxygen release kinetics and biological investigations of the scaffolds were also studied in vitro. The results showed that oxygen release behaviour was sustained and dependant on the concentration of CPO encapsulated in the PCL coating matrix. It was also demonstrated that the coated scaffolds, having 3% CPO in the coating system, could provide a great potential for promoting bone ingrowth with improving osteoblast cells viability and proliferation under hypoxic conditions. The findings indicated that the prepared scaffolds could play a significant role in engineering of large bone tissue implants with limitations in oxygen diffusion.

The Virtual Physiological Human: Ten Years After

Biomedical research and clinical practice are struggling to cope with the growing complexity that the progress of health care involves. The most challenging diseases, those with the largest socioeconomic impact (cardiovascular conditions; musculoskeletal conditions; cancer; metabolic, immunity, and neurodegenerative conditions), are all characterized by a complex genotype-phenotype interaction and by a "systemic" nature that poses a challenge to the traditional reductionist approach. In 2005 a small group of researchers discussed how the vision of computational physiology promoted by the Physiome Project could be translated into clinical practice and formally proposed the term Virtual Physiological Human. Our knowledge about these diseases is fragmentary, as it is associated with molecular and cellular processes on the one hand and with tissue and organ phenotype changes (related to clinical symptoms of disease conditions) on the other. The problem could be solved if we could capture all these fragments of knowledge into predictive models and then compose them into hypermodels that help us tame the complexity that such systemic behavior involves. In 2005 this was simply not possible-the necessary methods and technologies were not available. Now, 10 years later, it seems the right time to reflect on the original vision, the results achieved so far, and what remains to be done.

Bone fracture healing in mechanobiological modeling: A review of principles and methods

Bone fracture is a very common body injury. The healing process is physiologically complex, involving both biological and mechanical aspects. Following a fracture, cell migration, cell/tissue differentiation, tissue synthesis, and cytokine and growth factor release occur, regulated by the mechanical environment. Over the past decade, bone healing simulation and modeling has been employed to understand its details and mechanisms, to investigate specific clinical questions, and to design healing strategies. The goal of this effort is to review the history and the most recent work in bone healing simulations with an emphasis on both biological and mechanical properties. Therefore, we provide a brief review of the biology of bone fracture repair, followed by an outline of the key growth factors and mechanical factors influencing it. We then compare different methodologies of bone healing simulation, including conceptual modeling (qualitative modeling of bone healing to understand the general mechanisms), biological modeling (considering only the biological factors and processes), and mechanobiological modeling (considering both biological aspects and mechanical environment). Finally we evaluate different components and clinical applications of bone healing simulation such as mechanical stimuli, phases of bone healing, and angiogenesis.

The Use of Finite Element Analyses to Design and Fabricate Three-Dimensional Scaffolds for Skeletal Tissue Engineering

Computational modeling has been increasingly applied to the field of tissue engineering and regenerative medicine. Where in early days computational models were used to better understand the biomechanical requirements of targeted tissues to be regenerated, recently, more and more models are formulated to combine such biomechanical requirements with cell fate predictions to aid in the design of functional three-dimensional scaffolds. In this review, we highlight how computational modeling has been used to understand the mechanisms behind tissue formation and can be used for more rational and biomimetic scaffold-based tissue regeneration strategies. With a particular focus on musculoskeletal tissues, we discuss recent models attempting to predict cell activity in relation to specific mechanical and physical stimuli that can be applied to them through porous three-dimensional scaffolds. In doing so, we review the most common scaffold fabrication methods, with a critical view on those technologies that offer better properties to be more easily combined with computational modeling. Finally, we discuss how modeling, and in particular finite element analysis, can be used to optimize the design of scaffolds for skeletal tissue regeneration.

Multiscale Modeling of Bone Healing: Toward a Systems Biology Approach

Bone is a living part of the body that can, in most situations, heal itself after fracture. However, in some situations, healing may fail. Compromised conditions, such as large bone defects, aging, immuno-deficiency, or genetic disorders, might lead to delayed or non-unions. Treatment strategies for those conditions remain a clinical challenge, emphasizing the need to better understand the mechanisms behind endogenous bone regeneration. Bone healing is a complex process that involves the coordination of multiple events at different length and time scales. Computer models have been able to provide great insights into the interactions occurring within and across the different scales (organ, tissue, cellular, intracellular) using different modeling approaches [partial differential equations (PDEs), agent-based models, and finite element techniques]. In this review, we summarize the latest advances in computer models of bone healing with a focus on multiscale approaches and how they have contributed to understand the emergence of tissue formation patterns as a result of processes taking place at the lower length scales.

Stimulating Fracture Healing in Ischemic Environments: Does Oxygen Direct Stem Cell Fate during Fracture Healing?

Bone fractures represent an enormous societal and economic burden as one of the most prevalent causes of disability worldwide. Each year, nearly 15 million people are affected by fractures in the United States alone. Data indicate that the blood supply is critical for fracture healing; as data indicate that concomitant bone and vascular injury are major risk factors for non-union. However, the various role(s) that the vasculature plays remains speculative. Fracture stabilization dictates stem cell fate choices during repair. In stabilized fractures stem cells differentiate directly into osteoblasts and heal the injury by intramembranous ossification. In contrast, in non-stable fractures stem cells differentiate into chondrocytes and the bone heals through endochondral ossification, where a cartilage template transforms into bone as the chondrocytes transform into osteoblasts. One suggested role of the vasculature has been to participate in the stem cell fate decisions due to delivery of oxygen. In stable fractures, the blood vessels are thought to remain intact and promote osteogenesis, while in non-stable fractures, continual disruption of the vasculature creates hypoxia that favors formation of cartilage, which is avascular. However, recent data suggests that non-stable fractures are more vascularized than stable fractures, that oxygen does not appear associated with differentiation of stem cells into chondrocytes and osteoblasts, that cartilage is not hypoxic, and that oxygen, not sustained hypoxia, is required for angiogenesis. These unexpected results, which contrast other published studies, are indicative of the need to better understand the complex, spatio-temporal regulation of vascularization and oxygenation in fracture healing. This work has also revealed that oxygen, along with the promotion of angiogenesis, may be novel adjuvants that can stimulate healing in select patient populations.

Multiscale modeling methods in biomechanics

More and more frequently, computational biomechanics deals with problems where the portion of physical reality to be modeled spans over such a large range of spatial and temporal dimensions, that it is impossible to represent it as a single space-time continuum. We are forced to consider multiple space-time continua, each representing the phenomenon of interest at a characteristic space-time scale. Multiscale models describe a complex process across multiple scales, and account for how quantities transform as we move from one scale to another. This review offers a set of definitions for this emerging field, and provides a brief summary of the most recent developments on multiscale modeling in biomechanics. Of all possible perspectives, we chose that of the modeling intent, which vastly affect the nature and the structure of each research activity. To the purpose we organized all papers reviewed in three categories: 'causal confirmation,' where multiscale models are used as materializations of the causation theories; 'predictive accuracy,' where multiscale modeling is aimed to improve the predictive accuracy; and 'determination of effect,' where multiscale modeling is used to model how a change at one scale manifests in an effect at another radically different space-time scale. Consistent with how the volume of computational biomechanics research is distributed across application targets, we extensively reviewed papers targeting the musculoskeletal and the cardiovascular systems, and covered only a few exemplary papers targeting other organ systems. The review shows a research subdomain still in its infancy, where causal confirmation papers remain the most common. WIREs Syst Biol Med 2017, 9:e1375. doi: 10.1002/wsbm.1375 For further resources related to this article, please visit the WIREs website.

In silico bone mechanobiology: modeling a multifaceted biological system

Mechanobiology, the study of the influence of mechanical loads on biological processes through signaling to cells, is fundamental to the inherent ability of bone tissue to adapt its structure in response to mechanical stimulation. The immense contribution of computational modeling to the nascent field of bone mechanobiology is indisputable, having aided in the interpretation of experimental findings and identified new avenues of inquiry. Indeed, advances in computational modeling have spurred the development of this field, shedding new light on problems ranging from the mechanical response to loading by individual cells to tissue differentiation during events such as fracture healing. To date, in silico bone mechanobiology has generally taken a reductive approach in attempting to answer discrete biological research questions, with research in the field broadly separated into two streams: (1) mechanoregulation algorithms for predicting mechanobiological changes to bone tissue and (2) models investigating cell mechanobiology. Future models will likely take advantage of advances in computational power and techniques, allowing multiscale and multiphysics modeling to tie the many separate but related biological responses to loading together as part of a larger systems biology approach to shed further light on bone mechanobiology. Finally, although the ever‐increasing complexity of computational mechanobiology models will inevitably move the field toward patient‐specific models in the clinic, the determination of the context in which they can be used safely for clinical purpose will still require an extensive combination of computational and experimental techniques applied to in vitro and in vivo applications. WIREs Syst Biol Med 2016, 8:485–505. doi: 10.1002/wsbm.1356

For further resources related to this article, please visit the WIREs website.

An intramembranous ossification model for the in silico analysis of bone tissue formation in tooth extraction sites

The accurate modeling of biological processes allows us to predict the spatiotemporal behavior of living tissues by computer-aided (in silico) testing, a useful tool for the development of medical strategies, avoiding the expenses and potential ethical implications of in vivo experimentation. A model for bone healing in mouth would be useful for selecting proper surgical techniques in dental procedures. In this paper, the formulation and implementation of a model for Intramembranous Ossification is presented aiming to describe the complex process of bone tissue formation in tooth extraction sites. The model consists in a mathematical description of the mechanisms in which different types of cells interact, synthesize and degrade extracellular matrices under the influence of biochemical factors. Special attention is given to angiogenesis, oxygen-dependent effects and growth factor-induced apoptosis of fibroblasts. Furthermore, considering the depth-dependent vascularization of mandibular bone and its influence on bone healing, a functional description of the cell distribution on the severed periodontal ligament (PDL) is proposed. The developed model was implemented using the finite element method (FEM) and successfully validated by simulating an animal in vivo experiment on dogs reported in the literature. A good fit between model outcome and experimental data was obtained with a mean absolute error of 3.04%. The mathematical framework presented here may represent an important tool for the design of future in vitro and in vivo tests, as well as a precedent for future in silico studies on osseointegration and mechanobiology.

Role of oxygen as a regulator of stem cell fate during the spontaneous repair of osteochondral defects

The complexity of the in vivo environment makes it is difficult to isolate the effects of specific cues on regulating cell fate during regenerative events such as osteochondral defect repair. The objective of this study was to develop a computational model to explore how joint specific environmental factors regulate mesenchymal stem cell (MSC) fate during osteochondral defect repair. To this end, the spontaneous repair process within an osteochondral defect was simulated using a tissue differentiation algorithm which assumed that MSC fate was regulated by local oxygen levels and substrate stiffness. The developed model was able to predict the main stages of tissue formation observed by a number of in vivo studies. Following this, a parametric study was conducted to better understand why interventions that modulate angiogenesis dramatically impact the outcome of osteochondral defect healing. In the simulations where angiogenesis was reduced, by week 12, the subchondral plate was predicted to remain below the native tidemark, although the chondral region was composed entirely of cartilage and fibrous tissue. In the simulations where angiogenesis was increased, more robust cell proliferation and cartilage formation were observed during the first 4 weeks, however, by week 12 the subchondral plate had advanced above the native tidemark although any remaining tissue was either hypertrophic cartilage or fibrous tissue. These results suggest that osteochondral defect repair could be enhanced by interventions where angiogenesis is promoted but confined to within the subchondral region of the defect. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1026-1036, 2016.

Capturing the wide variety of impaired fracture healing phenotypes in Neurofibromatosis Type 1 with eight key factors: a computational study

Congenital pseudarthrosis of the tibia (CPT) is a rare disease which normally presents itself during early childhood by anterolateral bowing of the tibia and spontaneous tibial fractures. Although the exact etiology of CPT is highly debated, 40–80% of CPT patients are carriers of a mutation in the Neurofibromatosis Type 1 (NF1) gene, which can potentially result in an altered phenotype of the skeletal cells and impaired bone healing. In this study we use a computational model of bone regeneration to examine the effect of the Nf1 mutation on bone fracture healing by altering the parameter values of eight key factors which describe the aberrant cellular behaviour of Nf1 haploinsufficient and Nf1 bi-allelically inactivated cells. We show that the computational model is able to predict the formation of a hamartoma as well as a wide variety of CPT phenotypes through different combinations of altered parameter values. A sensitivity analysis by “Design of Experiments” identified the impaired endochondral ossification process and increased infiltration of fibroblastic cells as key contributors to the degree of severity of CPT. Hence, the computational model results have added credibility to the experimental hypothesis of a genetic cause (i.e. Nf1 mutation) for CPT.