A mathematical framework to study the effects of growth factor influences on fracture healing

During fracture healing, multipotential stem cells differentiate into specialized cells responsible for producing the different tissues involved in the bone regeneration process. This cell differentiation has been shown to be regulated by locally expressed growth factors. The details of their regulatory mechanisms need to be understood. In this work, we present a two-dimensional mathematical model of the bone healing process for moderate fracture gap sizes and fracture stability. The inflammatory and tissue regeneration stages of healing are simulated by modeling mesenchymal cell migration; mesenchymal cell, chondrocyte and osteoblast proliferation and differentiation, and extracellular matrix synthesis and degradation over time. The effects of two generic growth factors on cell differentiation are based on the experimentally studied chondrogenic and osteogenic effects of bone morphogenetic proteins-2 and 4 and transforming growth factor-beta-1, respectively. The model successfully simulates the progression of healing and predicts that the rate of osteogenic growth factor production by osteoblasts and the duration of the initial release of growth factors upon injury are particularly important parameters for complete ossification and successful healing. This temporo-spatial model of fracture healing is the first model to consider the effects of growth factors. It will help us understand the regulatory mechanisms involved in bone regeneration and provides a mathematical framework with which to design experiments and understand pathological conditions.     

Mechanistic,Mathematical Model to Predict the Dynamics of Tissue Genesis in Bone Defects via Mechanical Feedback and Mediation of Biochemical Factors

The link between mechanics and biology in the generation and the adaptation of bone has been well studied in context of skeletal development and fracture healing. Yet, the prediction of tissue genesis within - and the spatiotemporal healing of - postnatal defects, necessitates a quantitative evaluation of mechano-biological interactions using experimental and clinical parameters. To address this current gap in knowledge, this study aims to develop a mechanistic mathematical model of tissue genesis using bone morphogenetic protein (BMP) to represent of a class of factors that may coordinate bone healing. Specifically, we developed a mechanistic, mathematical model to predict the dynamics of tissue genesis by periosteal progenitor cells within a long bone defect surrounded by periosteum and stabilized via an intramedullary nail. The emergent material properties and mechanical environment associated with nascent tissue genesis influence the strain stimulus sensed by progenitor cells within the periosteum. Using a mechanical finite element model, periosteal surface strains are predicted as a function of emergent, nascent tissue properties. Strains are then input to a mechanistic mathematical model, where mechanical regulation of BMP-2 production mediates rates of cellular proliferation, differentiation and tissue production, to predict healing outcomes. A parametric approach enables the spatial and temporal prediction of endochondral tissue regeneration, assessed as areas of cartilage and mineralized bone, as functions of radial distance from the periosteum and time. Comparing model results to histological outcomes from two previous studies of periosteum-mediated bone regeneration in a common ovine model, it was shown that mechanistic models incorporating mechanical feedback successfully predict patterns (spatial) and trends (temporal) of bone tissue regeneration. The novel model framework presented here integrates a mechanistic feedback system based on the mechanosensitivity of periosteal progenitor cells, which allows for modeling and prediction of tissue regeneration on multiple length and time scales. Through combination of computational, physical and engineering science approaches, the model platform provides a means to test new hypotheses in silico and to elucidate conditions conducive to endogenous tissue genesis. Next generation models will serve to unravel intrinsic differences in bone genesis by endochondral and intramembranous mechanisms.     

A computational model to explore the role of angiogenic impairment on endochondral ossification during fracture healing

While it is well established that an adequate blood supply is critical to successful bone regeneration, it remains poorly understood how progenitor cell fate is affected by the altered conditions present in fractures with disrupted vasculature. In this study, computational models were used to explore how angiogenic impairment impacts oxygen availability within a fracture callus and hence regulates mesenchymal stem cell (MSC) differentiation and bone regeneration. Tissue differentiation was predicted using a previously developed algorithm which assumed that MSC fate is governed by oxygen tension and substrate stiffness. This model was updated based on the hypothesis that cell death, chondrocyte hypertrophy and endochondral ossification are regulated by oxygen availability. To test this, the updated model was used to simulate the time course of normal fracture healing, where it successfully predicted the observed quantity and spatial distribution of bone and cartilage at 10 and 20 days post-fracture (dpf). It also predicted the ratio of cartilage which had become hypertrophic at 10 dpf. Following this, three models of fracture healing with increasing levels of angiogenic impairment were developed. Under mild impairment, the model predicted experimentally observed reductions in hypertrophic cartilage at 10 dpf as well as the persistence of cartilage at 20 dpf. Models of more severe impairment predicted apoptosis and the development of fibrous tissue. These results provide insight into how factors specific to an ischemic callus regulate tissue regeneration and provide support for the hypothesis that chondrocyte hypertrophy and endochondral ossification during tissue regeneration are inhibited by low oxygen.     

Beneficial effects of moderate,early loading and adverse effects of delayed or excessive loading on bone healing

Fracture healing involves the differentiation and proliferation of cells in the callus and the synthesis and degradation of connective, cartilage and bone tissue. These processes are initiated and tightly regulated by growth factors and by the mechanical environment in the callus. In this work we incorporated the effects of mechanical stimulation on cell differentiation and ossification into a previously developed temporal-spatial model of growth factor mediated fracture healing. In particular, the stimulatory and inhibitory effects of dilatational and deviatoric strains were modeled. This predictive model was then calibrated and validated using well-defined in vivo experiments from the literature. As in the experiments, the results of the model demonstrated the beneficial and adverse effects of moderate and excessive loading, respectively, as well as the negative effects of delaying mechanical stimulation of rigidly fixed calluses. In addition, the model examined loading conditions and time points beyond those used in the experiments, providing a more complete and mechanistic characterization of the effects of loading in the biological tissue response associated with fracture healing.     

The spatio-temporal arrangement of different tissues during bone healing as a result of simple mechanobiological rules

During secondary bone healing, different tissue types are formed within the fracture callus depending on the local mechanical and biological environment. Our aim was to understand the temporal succession of these tissue patterns for a normal bone healing progression by means of a basic mechanobiological model. The experimental data stemmed from an extensive, previously published animal experiment on sheep with a 3?mm tibial osteotomy. Using recent experimental data, the development of the hard callus was modelled as a porous material with increasing stiffness and decreasing porosity. A basic phenomenological model was employed with a small number of simulation parameters, which allowed comprehensive parameter studies. The model distinguished between the formation of new bone via endochondral and intramembranous ossification. To evaluate the outcome of the computer simulations, the tissue images of the simulations were compared with experimentally derived tissue images for a normal healing progression in sheep. Parameter studies of the threshold values for the regulation of tissue formation were performed, and the source of the biological stimulation (comprising e.g. stem cells) was varied. It was found that the formation of the hard callus could be reproduced in silico for a wide range of threshold values. However, the bridging of the fracture gap by cartilage on the periosteal side was observed only (i) for a rather specific choice of the threshold values for tissue differentiation and (ii) when assuming a strong source of biological stimulation at the periosteum.     

Collagen IX is indispensable for timely maturation of cartilage during fracture repair in mice.

Fracture repair recapitulates in adult organisms the sequence of cell biological events of endochondral ossification during skeletal development and growth. After initial inflammation and deposition of granulation tissue, a cartilaginous callus is formed which, subsequently, is remodeled into bone. In part, bone formation is influenced also by the properties of the extracellular matrix of the cartilaginous callus. Deletion of individual macromolecular components can alter extracellular matrix suprastructures, and hence stability and organization of mesenchymal tissues. Here, we took advantage of the collagen IX knockout mouse model to better understand the role of this collagen for organization, differentiation and maturation of a cartilaginous template during formation of new bone. Although a seemingly crucial component of cartilage fibrils is missing, collagen IX-deficient mice develop normally, but are predisposed to premature joint cartilage degeneration. However, we show here that lack of collagen IX alters the time course of callus differentiation during bone fracture healing. The maturation of cartilage matrix was delayed in collagen IX-deficient mice calli as judged by collagen X expression during the repair phase and the total amount of cartilage matrix was reduced. Entering the remodeling phase of fracture healing, Col9a1(-/-) calli retained a larger percentage of cartilage matrix than in wild type indicating also a delayed formation of new bone. We concluded that endochondral bone formation can occur in collagen IX knockout mice but is impaired under conditions of stress, such as the repair of an unfixed fractured long bone.     

Differential distribution of V-type H+-ATPase and Na+/K+-ATPase in the branchial chamber of the palaemonid shrimp Macrobrachium amazonicum

Increased fragility fracture risk with improper healing is a frequent and severe complication of insulin resistance (IR). The mechanisms impairing bone health in IR are still not fully appreciated, which gives importance to studies on bone pathologies in animal models of diabetes. Mice deficient in leptin signaling are widely used models of IR and its comorbidities. Leptin was first recognized as a hormone, regulating appetite and energy balance; however, recent studies have expanded its role showing that leptin is a link between insulin-dependent metabolism and bone homeostasis. In the light of these findings, it is intriguing to consider the role of leptin resistance in bone regeneration. In this study, we show that obese diabetic mice lacking leptin receptor (db/db) are deficient in postnatal regenerative osteogenesis. We apply an ectopic osteogenesis and a fracture healing model, both showing that db/db mice display compromised bone acquisition and regeneration capacity. The underlying mechanisms include delayed periosteal mesenchymatic osteogenesis, premature apoptosis of the cartilage callus and impaired microvascular invasion of the healing tissue. Our study supports the use of the db/db mouse as a model of IR associated bone-healing deficits and can aid further studies of mesenchymatic cell homing and differentiation, microvascular invasion, cartilage to bone transition and callus remodeling in diabetic fracture healing.     

In silico study of bone tissue regeneration in an idealised porous hydrogel scaffold using a mechano-regulation algorithm

Mechanical stimulation, in the form of fluid perfusion or mechanical strain, enhances osteogenic differentiation and overall bone tissue formation by mesenchymal stems cells cultured in biomaterial scaffolds for tissue engineering applications. In silico techniques can be used to predict the mechanical environment within biomaterial scaffolds, and also the relationship between bone tissue regeneration and mechanical stimulation, and thereby inform conditions for bone tissue engineering experiments. In this study, we investigated bone tissue regeneration in an idealised hydrogel scaffold using a mechano-regulation model capable of predicting tissue differentiation, and specifically compared five loading cases, based on known experimental bioreactor regimes. These models predicted that low levels of mechanical loading, i.e. compression (0.5% strain), pore pressure of 10 kPa and a combination of compression (0.5%) and pore pressure (10 kPa), could induce more osteogenic differentiation and lead to the formation of a higher bone tissue fraction. In contrast greater volumes of cartilage and fibrous tissue fractions were predicted under higher levels of mechanical loading (i.e. compression strain of 5.0% and pore pressure of 100 kPa). The findings in this study may provide important information regarding the appropriate mechanical stimulation for in vitro bone tissue engineering experiments.     

A mechano-regulation model for tissue differentiation during fracture healing: analysis of gap size and loading

Bone has a capability to repair itself when it is fractured. Repair involves the generation of intermediate tissues, such as fibrous connective tissue, cartilage and woven bone, before final bone healing can occur. The intermediate tissues serve to stabilise the mechanical environment and provide a scaffold for differentiation of new tissues. The repair process is fundamentally affected by mechanical loading and by the geometric configuration of the fracture fragments. Biomechanical analyses of fracture healing have previously computed the stress distribution within the callus and identified the components of the stress tensor favouring or inhibiting differentiation of particular tissue phenotypes. In this paper, a biphasic poroelastic finite element model of a fracture callus is used to simulate the time-course of tissue differentiation during fracture healing. The simulation begins with granulation tissue (post-inflammation phase) and finishes with bone resorption. The biomechanical regulatory model assumes that tissue differentiation is controlled by a combination of shear strain and fluid flow acting within the tissue. High shear strain and fluid flows are assumed to deform the precursor cells stimulating formation of fibrous connective tissue, lower levels stimulate formation of cartilage, and lower again allows ossification. This mechano-regulatory scheme was tested by simulating healing in fractures with different gap sizes and loading magnitudes. The appearance and disappearance of the various tissues found in a callus was similar to histological observation. The effect of gap size and loading magnitude on the rate of reduction of the interfragmentary strain was sufficiently close to confirm the hypothesis that tissue differentiation phenomena could be governed by the proposed mechano-regulation model.     

Identification of tissue differentiation rates in a mechanobiological model of fracture healing

As a basis for model-based analysis of the processes in secondary fracture healing, a dynamical model is presented that characterises the physiological status in the fracture area by the location-dependent composition of tissues. Five types of tissue are distinguished: connective tissue, cartilage, bone, haematoma and avascular bone. A rule base is given that describes dynamical tissue differentiation processes. The rules consider not only a mechanical stimulus but also osteogenic and a vasculative factors as biological stimuli. Within this model structure, it is possible, e.g., to distinguish intramembranous from endochondral ossification processes. An objective function is introduced to assess accordance between the model-based simulation results and reference healing stages. By minimising this objective function, relevant tissue differentiation rates can be determined. For a reference process of secondary fracture healing it could be shown that the intramembranous ossification rate of 0.313%/day (from connective tissue to bone) is much smaller than the endochondral ossification rate of 1.136%/day (from cartilage to bone). In order to verify the model approach, it is transferred to simulate long bone distraction. Results show that healing patterns of bone distraction can be predicted. Using this method, it is possible to identify model parameters for individual subjects. This will allow a patient-specific analysis of tissue healing processes in future.     

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

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

Determining the most important cellular characteristics for fracture healing using design of experiments methods

Computational models are employed as tools to investigate possible mechanoregulation pathways for tissue differentiation and bone healing. However, current models do not account for the uncertainty in input parameters, and often include assumptions about parameter values that are not yet established. The objective of this study was to determine the most important cellular characteristics of a mechanoregulatory model describing both cell phenotype-specific and mechanobiological processes that are active during bone healing using a statistical approach. The computational model included an adaptive two-dimensional finite element model of a fractured long bone. Three different outcome criteria were quantified: (1) ability to predict sequential healing events, (2) amount of bone formation at early, mid and late stages of healing and (3) the total time until complete healing. For the statistical analysis, first a resolution IV fractional factorial design (L₆₄) was used to identify the most significant factors. Thereafter, a three-level Taguchi orthogonal array (L₂₇) was employed to study the curvature (non-linearity) of the 10 identified most important parameters. The results show that the ability of the model to predict the sequences of normal fracture healing was predominantly influenced by the rate of matrix production of bone, followed by cartilage degradation (replacement). The amount of bone formation at early stages was solely dependent on matrix production of bone and the proliferation rate of osteoblasts. However, the amount of bone formation at mid and late phases had the rate of matrix production of cartilage as the most influential parameter. The time to complete healing was primarily dependent on the rate of cartilage degradation during endochondral ossification, followed by the rate of cartilage formation. The analyses of the curvature revealed a linear response for parameters related to bone, where higher rates of formation were more beneficial to healing. In contrast, parameters related to fibrous tissue and cartilage showed optimum levels. Some fibrous connective tissue- and cartilage formation was beneficial to bone healing, but too much of either tissue delayed bone formation. The identified significant parameters and processes are further confirmed by in vivo animal experiments in the literature. This study illustrates the potential of design of experiments methods for evaluating computational mechanobiological model parameters and suggests that further experiments should preferably focus at establishing values of parameters related to cartilage formation and degradation.     

The enhancement of bone regeneration by ultrasound

Millions of fractures occur every year worldwide, with nearly 6.2 million fractures reported annually in the United States alone. Even though treatment methods have improved over the last few decades, 5–10% of fractures still show delayed healing. A significant subpopulation of these delayed healings do not heal by nine months and are thus termed non-unions. Experimental studies have shown some evidence that low intensity pulsed ultrasound stimulation (LIPUS) results in enhanced bone regeneration during fracture healing and callus distraction. LIPUS treatment has led to increased callus area and accelerated return of bone strength following fracture. Histological studies suggest that LIPUS influences all major cell types involved in bone healing, including osteoblasts, osteoclasts, chondrocytes and mesenchymal stem cells. The affect of LIPUS seems to be limited to cells in soft tissue, whereas cells in calcified bone seem not to be effected. In vitro cell culture studies as well as tissue culture studies have shown some effects on cell differentiation and protein synthesis. Even though the energy used by LIPUS treatment is extremely low, the effects are evident. The most probable source of the therapeutic benefits observed with LIPUS treatment involves nonthermal mechanisms that influence cell membrane permeability and increase cellular activity. Despite clinical and experimental studies demonstrating the enhancing effect of LIPUS on bone regeneration, the biophysical mechanisms involved in the complex fracture healing process remain unclear and requires further research.     

Mathematical Modeling in Wound Healing,Bone Regeneration and Tissue Engineering

The processes of wound healing and bone regeneration and problems in tissue engineering have been an active area for mathematical modeling in the last decade. Here we review a selection of recent models which aim at deriving strategies for improved healing. In wound healing, the models have particularly focused on the inflammatory response in order to improve the healing of chronic wound. For bone regeneration, the mathematical models have been applied to design optimal and new treatment strategies for normal and specific cases of impaired fracture healing. For the field of tissue engineering, we focus on mathematical models that analyze the interplay between cells and their biochemical cues within the scaffold to ensure optimal nutrient transport and maximal tissue production. Finally, we briefly comment on numerical issues arising from simulations of these mathematical models.     

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

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

Mineral and matrix contributions to rigidity in fracture healing

The purpose of this study was to investigate the relationships among selected properties of fracture callus: bending rigidity, tissue density, mineral density, matrix density and mineral-to-matrix ratio. The experimental model was an osteotomized canine radius in which the development of the fracture callus was modified by electrical stimulation with various levels of direct current. This resulted in a range of values for the selected properties of the callus, determined post mortem at 7 weeks after osteotomy. We found that the rigidity (R) of the bone-callus combination obeyed relationships of the form R = axb, where x is the tissue density, mineral density, matrix density or the mineral-to-matrix ratio of the repair tissue. These are analogous to power-law relationships found in studies of compact and cancellous bone. The results suggest that fracture callus at 7 weeks after osteotomy in canine radius behaves more like immature compact bone than cancellous bone in its mechanical and physicochemical properties. The present study demonstrates the feasibility of developing non-invasive in vivo densitometric methods to monitor fracture healing, since models may be developed that can predict mechanical properties from densitometric data. Further studies are needed to develop a refined model based on experimental data on the mechanical and physicochemical properties and microstructure of fracture callus at different stages of healing.     

Expansion on a matrix deposited by nonchondrogenic urine stem cells strengthens the chondrogenic capacity of repeated-passage bone marrow stromal cells

Human urine-derived stem cells (hUSCs) are a newly found type of stem cell with a potential for therapeutic application in urology. The aim of this study is to investigate whether hUSCs contribute to cartilage regeneration. Despite their characterization with multi-lineage differentiation capacities, in terms of osteogenesis, adipogenesis and myogenesis, hUSCs do not show the ability to differentiate into chondrocytes. Human bone marrow stromal cells (hBMSCs) are a tissue-specific stem cell for endochondral bone formation; however, repeated-passage hBMSCs have a lower capacity for chondrogenic differentiation. We found that the extracellular matrix (ECM) deposited by hUSCs (UECM) can greatly recharge repeated-passage hBMSCs toward chondrogenic differentiation, a result that might be explained by trophic factors released from hUSCs being immobilized in UECM. We also found that ECM from repeated-passage hBMSCs (BECM) have a limited rejuvenation effect. The Wnt11-mediated noncanonical signaling pathway might be responsible for UECM-mediated hBMSC rejuvenation and subsequent chondrogenic differentiation. Our data indicate that commercially available UECM from young healthy donors might represent a simple and promising approach for autologous hBMSC rejuvenation. This study also provides an excellent model for investigating the effect of trophic factors released by stem cells on tissue regeneration without interference by stem cell differentiation.     

Regulation of fracture repair by growth factors.

Fractured bones heal by a cascade of cellular events in which mesenchymal cells respond to unknown regulators by proliferating, differentiating, and synthesizing extracellular matrix. Current concepts suggest that growth factors may regulate different steps in this cascade (10). Recent studies suggest regulatory roles for PDGF, aFGF, bFGF, and TGF-beta in the initiation and the development of the fracture callus. Fracture healing begins immediately following injury, when growth factors, including TGF-beta 1 and PDGF, are released into the fracture hematoma by platelets and inflammatory cells. TGF-beta 1 and FGF are synthesized by osteoblasts and chondrocytes throughout the healing process. TGF-beta 1 and PDGF appear to have an influence on the initiation of fracture repair and the formation of cartilage and intramembranous bone in the initiation of callus formation. Acidic FGF is synthesized by chondrocytes, chondrocyte precursors, and macrophages. It appears to stimulate the proliferation of immature chondrocytes or precursors, and indirectly regulates chondrocyte maturation and the expression of the cartilage matrix. Presumably, growth factors in the callus at later times regulate additional steps in repair of the bone after fracture. These studies suggest that growth factors are central regulators of cellular proliferation, differentiation, and extracellular matrix synthesis during fracture repair. Abnormal growth factor expression has been implicated as causing impaired or abnormal healing in other tissues, suggesting that altered growth factor expression also may be responsible for abnormal or delayed fracture repair. As a complete understanding of fracture-healing regulation evolves, we expect new insights into the etiology of abnormal or delayed fracture healing, and possibly new therapies for these difficult clinical problems.