Assessment of a mechano-regulation theory of skeletal tissue differentiation in an in vivo model of mechanically induced cartilage formation

Development of a subject-specific finite element analysis workflow to assess local biomechanics during segmental bone defect healing

Longitudinal estimation of local strain distributions within the regenerative niche of segmental femoral fractures is important for understanding mechanobiology principles for bone healing to design more effective rehabilitation regimens and mitigate nonunion complications. Finite element (FE) modeling is the standard for investigating these biomechanical parameters, yet most existing models lack clinical relevance due to their use of generic data and computational inefficiency. This study developed a subject-specific FE workflow aimed at accurate biomechanical predictions based on subject-specific data while addressing the limitations of previous approaches. For the experimental study, near-critical-sized segmental bone defects were created in the femurs of Wistar rats and stabilized with internal fixators before rehabilitation. Subject-specific geometries of the defect were generated from in vivo micro-CT scans, which were also used to assign material coefficients. Generalized geometries of the cortical and trabecular bone and fixator were integrated to increase computational efficiency. In addition, axial strain data from strain gauges on the fixators were used to define subject-specific boundary conditions, enabling a longitudinal study of the healing process. Sensitivity analyses revealed that incorporating subject-specific boundary conditions significantly enhanced model accuracy, a factor often overlooked in conventional approaches. The workflow was used to build six defect models to approximate compressive strains within the defect and the joint contact force. Strain distributions correlated with experimentally observed mineralization and better predicted functional bone bridging (union) compared to bone volume metrics. This efficient workflow facilitates the assessment of local biomechanics during bone healing and highlights their influence on adaptive regeneration. Further, the findings support the potential application of the subject-specific modeling workflow to guide clinical decision-making and improve therapeutic outcomes for treating bone fractures.

What is the stable internal fixation for the unstable and osteoporotic supracondylar femoral fractures: a finite element analysis

BACKGROUND

Osteoporotic supracondylar femoral fractures (OSFF) have historically been managed by the lateral anatomical locking plate with reasonable success. However, for some kinds of unstable and osteoporotic supracondylar femoral fractures (UOSFF), especially with bone defects, unilateral locking plate (ULLP) fixation failed or resulted in implant breakage. This paper is going to explore what is the stable internal fixation for UOSFF by adding the bilateral locking plate (BLLP) fixation.

METHODS

OSFF models were divided into two groups according to the fracture line type, which would be further subdivided according to their angle of fracture line, presence of bone defect, location, and degree of bone defect. Thereafter, kinds of locking plate fixation were constructed. A 2010-N load was applied to the femoral head, and a 1086-N load was applied to the greater trochanter. In this condition, the maximum von Mises stress distribution of models were investigated.

RESULTS

Firstly, it was obviously found that the stress concentration in the BLLP group was more dispersed than that in the ULLP group. Secondly, according to the fracture line analysis, the stress value of fracture line type in "\" model group was higher than that of "/" model group. Moreover, with the increase in fracture line angle, the stress value of the model increased. Thirdly, from the bone defect analysis, the stress value of the medial bone defect (MBD) model group was higher than that of the lateral bone defect (LBD) model group. And as the degree of bone defect increased, the stress value increased gradually in the model group.

CONCLUSION

In the following four cases, lateral unilateral locking plate fixation cannot effectively stabilize the fracture end, and double locking plate internal fixation is a necessary choice. First, when the angle of the fracture line is large (30, 45). Second, when the fracture line type is "/." Third, when the bone defect is large. Fourth, when the bone defect is medial.

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

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

Balance Between Mechanical Stability and Mechano-Biology of Fracture Healing Under Volar Locking Plate

The application of volar locking plate (VLP) is promising in the treatment of dorsally comminuted and displaced fracture. However, the optimal balance between the mechanical stability of VLP and the mechanobiology at the fracture site is still unclear. The purpose of this study is to develop numerical models in conjunction with experimental studies to identify the favourable mechanical microenvironment for indirect healing, by optimizing VLP configuration and post-operative loadings for different fracture geometries. The simulation results show that the mechanical behaviour of VLP is mainly governed by the axial compression. In addition, the model shows that, under relatively large gap size (i.e., 3-5 mm), the increase of FWL could enhance chondrocyte differentiation while a large BPD could compromise the mechanical stability of VLP. Importantly, bending moment produced by wrist flexion/extension and torsion moment produced from forearm rotation could potentially hinder endochondral ossification at early stage of healing. The developed model could potentially assist orthopaedic surgeons in surgical pre-planning and designing post-operation physical therapy for treatment of distal radius fractures.

Bone remodeling analysis for a swine skull at continuous scale based on the smoothed finite element method

In order to study bone response during chewing, bone remodeling analysis at a continuous scale is performed to a swine skull obtained using μCT. The smoothed finite element method (S-FEM) is utilized to replace the finite element method (FEM) in bone remodeling as it is solving the "overly-stiff" problem in FEM by introducing strain smoothing technology to soften the stiffness matrix. Three S-FEM models with different levels of softening effects are developed, including node-based, edge-based, and face-based, which leads to various bone remodeling results for a better understanding of the remodeling process. During the remodeling process, the strain energy density is used as the mechanical stimulus, and the surface elements or smoothing domains are regarded as cortical bone. Under the action of mechanical stimuli, cortical bone and cancellous bone have been remodeled. In remodeling progress, ES-FEM shows close results as compared with the experimental μCT in nodal bone density distribution, FEM and FS-FEM are close to the μCT experimental model in average nodal density. In summary, the combined use of several methods provides more angles for the description of bone remodeling.

Mechanical performance and implications on bone healing of different screw configurations for plate fixation of diaphyseal tibia fractures: a computational study

Diaphyseal tibia fractures may require plate fixation for proper healing to occur. Currently, there is no consensus on the number of screws required for proper fixation or the optimal placement of the screws within the plate. Mechanical stability of the construct is a leading criterion for choosing plate and screws configuration. However, number and location of screws have implications on the mechanical environment at the fracture site and, consequently, on bone healing response: The interfragmentary motion attained with a specific plate and screw construct may elicit mechano-transduction signals influencing cell-type differentiation, which in turn affects how well the fracture heals. This study investigated how different screw configurations affect mechanical performance of a tibia plate fixation construct. Three configurations of an eight-hole plate were considered with the fracture in the center of the plate: eight screws-screws at first, fourth, fifth and eighth hole and screws at first, third, sixth and eighth hole. Constructs' stiffness was compared through biomechanical tests on bone surrogates. A finite element model of tibia diaphyseal fracture was used to conduct a stress analysis on the implanted hardware. Finally, the potential for bone regeneration of each screw configuration was assessed via the computational model through the evaluation of the magnitude of mechano-transduction signals at the bone callus. The results of this study indicate that having screws at fourth and fifth holes represents a preferable configuration since it provides mechanical properties similar to those attained by the stiffest construct (eight screws), and elicits an ideal bone regenerative response.

Age-Related Changes in the Mechanical Regulation of Bone Healing Are Explained by Altered Cellular Mechanoresponse

Increasing age is associated with a reduced bone regeneration potential and increased risk of morbidities and mortality. A reduced bone formation response to mechanical loading has been shown with aging, and it remains unknown if the interplay between aging and mechanical stimuli during regeneration is similar to adaptation. We used a combined in vivo/in silico approach to investigate age-related alterations in the mechanical regulation of bone healing and identified the relative impact of altered cellular function on tissue patterns during the regenerative cascade. To modulate the mechanical environment, femoral osteotomies in adult and elderly mice were stabilized using either a rigid or a semirigid external fixator, and the course of healing was evaluated using histomorphometric and micro-CT analyses at 7, 14, and 21 days post-surgery. Computer models were developed to investigate the influence of the local mechanical environment within the callus on tissue formation patterns. The models aimed to identify the key processes at the cellular level that alter the mechanical regulation of healing with aging. Fifteen age-related biological alterations were investigated on two levels (adult and elderly) with a design of experiments setup. We show a reduced response to changes in fixation stability with age, which could be explained by reduced cellular mechanoresponse, simulated as alteration of the ranges of mechanical stimuli driving mesenchymal stem cell differentiation. Cellular mechanoresponse has been so far widely ignored as a therapeutic target in aged patients. Our data hint to mechanotherapeutics as a potential treatment to enhance bone healing in the elderly. © 2019 American Society for Bone and Mineral Research.

Biomechanical Evaluation of the Effect of Mesenchymal Stem Cells on Cartilage Regeneration in Knee Joint Osteoarthritis

Numerous clinical studies have reported cell-based treatments for cartilage regeneration in knee joint osteoarthritis using mesenchymal stem cells (MSCs). However, the post-surgery rehabilitation and weight-bearing times remain unclear. Phenomenological computational models of cartilage regeneration have been only partially successful in predicting experimental results and this may be due to simplistic modeling assumptions and loading conditions of cellular activity. In the present study, we developed a knee joint model of cell and tissue differentiation based on a more mechanistic approach, which was applied to cartilage regeneration in osteoarthritis. First, a phenomenological biphasic poroelastic finite element model was developed and validated according to a previous study. Second, this method was applied to a real knee joint model with a cartilage defect created to simulate the tissue regeneration process. The knee joint model was able to accurately predict several aspects of cartilage regeneration, such as the cell and tissue distributions in the cartilage defect. Additionally, our results indicated that gait cycle loading with flexion was helpful for cartilage regeneration compared to the use of simple weight-bearing loading.

Optimal mechanical properties of a scaffold for cartilage regeneration using finite element analysis

The development of successful scaffolds for bone tissue engineering requires concurrent engineering that combines different research fields. In previous studies, phenomenological computational models predicted the mechanical properties of a scaffold in a simple loading condition using the mechano-regulation theory. Therefore, the aim of this study is to predict the mechanical properties of an optimum scaffold required for cartilage regeneration using three-dimensional knee joint developed from medical imaging and mechano-regulation theory. It was predicted that the scaffold with optimal mechanical properties would result in greater amounts of cartilage tissue formation than without a scaffold. The results demonstrated the ability of the algorithms to design optimized scaffolds with target properties and confirmed the applicability of set techniques for bone tissue engineering. The scaffolds were optimized to suit the site-specific loading requirements, and the results reveal a new approach for computational simulations in tissue engineering.

Far Cortical Locking Fixation of Distal Femur Fractures is Dominated by Shear at Clinically Relevant Bridge Spans

OBJECTIVES

Far cortical locking (FCL) constructs have been shown to increase axial interfragmentary displacement while limiting shear and have been specifically recommended in the treatment of distal femur fractures. However, there is no available data regarding their mechanical behavior within the range of bridge spans typically used for comminuted distal femur fractures. This biomechanical study of distal femur locked plate fixation assessed 4 methods of diaphyseal fixation for associated axial and shear displacement at bridge spans typically used in clinical practice.

METHODS

Distal femur locking plates were used to bridge simulated fractures in femur surrogates with 4 different methods of diaphyseal fixation (bicortical locking, bicortical nonlocking, near cortical locking, and FCL). Axial and shear displacement were assessed at 5 different bridge spans for each fixation method.

RESULTS

Diaphyseal fixation type was associated with the amount of shear (P = 0.04), but not the amount of axial displacement (P = 0.39). Specifically, FCL constructs demonstrated greater shear than bicortical locking (median 4.57 vs. 2.94 mm, P = 0.02) and bicortical nonlocking (median 4.57 vs. 3.41 mm, P = 0.02) constructs.

CONCLUSIONS

Unexpectedly, FCL constructs demonstrated greater shear than bicortical locking and nonlocking constructs and similar axial displacement for all fixation methods. Bridge span had a dominant effect on displacement that interacted negatively with more flexible FCL diaphyseal fixation. Potentially interactive construct features are best studied in concert. Given the complexity of these relationships, computational modeling will likely play an integral role in future mechanotransduction research.

Motion Predicts Clinical Callus Formation: Construct-Specific Finite Element Analysis of Supracondylar Femoral Fractures

BACKGROUND

Mechanotransduction is theorized to influence fracture-healing, but optimal fracture-site motion is poorly defined. We hypothesized that three-dimensional (3-D) fracture-site motion as estimated by finite element (FE) analysis would influence callus formation for a clinical series of supracondylar femoral fractures treated with locking-plate fixation.

METHODS

Construct-specific FE modeling simulated 3-D fracture-site motion for sixty-six supracondylar femoral fractures (OTA/AO classification of 33A or 33C) treated at a single institution. Construct stiffness and directional motion through the fracture were investigated to assess the validity of construct stiffness as a surrogate measure of 3-D motion at the fracture site. Callus formation was assessed radiographically for all patients at six, twelve, and twenty-four weeks postoperatively. Univariate and multivariate linear regression analyses examined the effects of longitudinal motion, shear (transverse motion), open fracture, smoking, and diabetes on callus formation. Construct types were compared to determine whether their 3-D motion profile was associated with callus formation.

RESULTS

Shear disproportionately increased relative to longitudinal motion with increasing bridge span, which was not predicted by our assessment of construct stiffness alone. Callus formation was not associated with open fracture, smoking, or diabetes at six, twelve, or twenty-four weeks. However, callus formation was associated with 3-D fracture-site motion at twelve and twenty-four weeks. Longitudinal motion promoted callus formation at twelve and twenty-four weeks (p = 0.017 for both). Shear inhibited callus formation at twelve and twenty-four weeks (p = 0.017 and p = 0.022, respectively). Titanium constructs with a short bridge span demonstrated greater longitudinal motion with less shear than did the other constructs, and this was associated with greater callus formation (p < 0.001).

CONCLUSIONS

In this study of supracondylar femoral fractures treated with locking-plate fixation, longitudinal motion promoted callus formation, while shear inhibited callus formation. Construct stiffness was found to be a poor surrogate of fracture-site motion. Future implant design and operative fixation strategies should seek to optimize 3-D fracture-site motion rather than rely on surrogate measures such as axial stiffness.

Mechanical microenvironments and protein expression associated with formation of different skeletal tissues during bone healing

Uncovering the mechanisms of the sensitivity of bone healing to mechanical factors is critical for understanding the basic biology and mechanobiology of the skeleton, as well as for enhancing clinical treatment of bone injuries. This study refined an experimental method of measuring the strain microenvironment at the site of a bone injury during bone healing. This method used a rat model in which a well-controlled bending motion was applied to an osteotomy to induce the formation of pseudarthrosis that is composed of a range of skeletal tissues, including woven bone, cartilage, fibrocartilage, fibrous tissue, and clot tissue. The goal of this study was to identify both the features of the strain microenvironment associated with formation of these different tissues and the expression of proteins frequently implicated in sensing and transducing mechanical cues. By pairing the strain measurements with histological analyses that identified the regions in which each tissue type formed, we found that formation of the different tissue types occurs in distinct strain microenvironments and that the type of tissue formed is correlated most strongly to the local magnitudes of extensional and shear strains. Weaker correlations were found for dilatation. Immunohistochemical analyses of focal adhesion kinase and rho family proteins RhoA and CDC42 revealed differences within the cartilaginous tissues in the calluses from the pseudarthrosis model as compared to fracture calluses undergoing normal endochondral bone repair. These findings suggest the involvement of these proteins in the way by which mechanical stimuli modulate the process of cartilage formation during bone healing.

Microcomputed tomography: approaches and applications in bioengineering

Microcomputed tomography (microCT) has become a standard and essential tool for quantifying structure-function relationships, disease progression, and regeneration in preclinical models and has facilitated numerous scientific and bioengineering advancements over the past 30 years. In this article, we recount the early events that led to the initial development of microCT and review microCT approaches for quantitative evaluation of bone, cartilage, and cardiovascular structures, with applications in fundamental structure-function analysis, disease, tissue engineering, and numerical modeling. Finally, we address several next-generation approaches under active investigation to improve spatial resolution, acquisition time, tissue contrast, radiation dose, and functional and molecular information.

A microfluidic-based multi-shear device for investigating the effects of low fluid-induced stresses on osteoblasts

Interstitial fluid flow (IFF) within the extracellular matrix (ECM) produces low magnitude shear stresses on cells. Fluid flow-induced stress (FSS) plays an important role during tissue morphogenesis. To investigate the effect of low FSS generated by IFF on cells, we developed a microfluidic-based cell culture device that can generate multiple low shear stresses. By changing the length and width of the flow-in channels, different continuous low level shear stresses could be generated in individual cell culture chambers. Numerical calculations demonstrate uniform shear stress distributions of the major cell culture area of each chamber. This calculation is further confirmed by the wall shear stress curves. The effects of low FSS on MC3T3-E1 proliferation and differentiation were studied using this device. It was found that FSS ranging from 1.5 to 52.6 µPa promoted MC3T3-E1 proliferation and differentiation, but FSS over 412 µPa inhibited the proliferation and differentiation of MC3T3-E1 cells. FSS ranging from 1.5 to 52.6 µPa also increased the expression of Runx2, a key transcription factor regulating osteoblast differentiation. It is suggested that Runx2 might be an important regulator in low FSS-induced MC3T3-E1 differentiation. This device allows for detailed study of the effect of low FSS on the behaviors of cells; thus, it would be a useful tool for analysis of the effects of IFF-induced shear stresses on cells.

Mechanical regulation of bone regeneration: theories, models, and experiments

How mechanical forces influence the regeneration of bone remains an open question. Their effect has been demonstrated experimentally, which has allowed mathematical theories of mechanically driven tissue differentiation to be developed. Many simulations driven by these theories have been presented, however, validation of these models has remained difficult due to the number of independent parameters considered. An overview of these theories and models is presented along with a review of experimental studies and the factors they consider. Finally limitations of current experimental data and how this influences modeling are discussed and potential solutions are proposed.

Simulation of interstitial fluid flow in ligaments: comparison among Stokes, Brinkman and Darcy models

In this paper, we use Stokes, Brinkman and Darcy equations to approximate the porous continuum media of ligament tissues respectively, simulate the flow field with FLUENT software, and study the shear stress on the cell surface due to the interstitial fluid flow. Since the Brinkman equation approaches Stokes equation well in high hydraulic permeability (k p) condition (k p ≥1.0×10(-8) m(2) in our numerical simulation), and it is an approximation to Darcy model in low k p condition (k p ≤5.0×10(-12) m(2) in our numerical simulation), we used the Brinkman model to simulate the interstitial fluid flow in the ligament where k p is approximately 1.0×10(-16) m(2). It shows k p and anisotropic property have a little effect on the flow field, but have a great effect on the shear stress on the membrane of interstitial cells (τ cell). There is a linear relationship between τ cell and , when k p =1.0×10(-16) m(2) and the maximum τ cell (τ cell,max) is approximately 10 Pa. The anisotropic property will affect τ cell's distribution on the cell surface. When k x/k y>1, low τ cell dominates the cell, while when k x/k y<1, high τ cell dominants the cell.

The emergence of mechanoregulated endochondral ossification in evolution

The differentiation of skeletal tissue phenotypes is partly regulated by mechanical forces. This mechanoregulatory aspect of tissue differentiation has been the subject of many experimental and computational investigations. However, little is known about what factors promoted the emergence of mechanoregulated tissue differentiation in evolution, even though mechanoregulated tissue differentiation, for example during development or healing of adult bone, is crucial for vertebrate phylogeny. In this paper, we use a computational framework to test the hypothesis that the emergence of mechanosensitive genes that trigger endochondral ossification in evolution will stabilise in the population and create a variable mechanoregulated response, if the endochondral ossification process enhances fitness for survival. The model combines an evolutionary algorithm that considers genetic change with a mechanoregulated fracture healing model in which the fitness of animals in a population is determined by their ability to heal their bones. The simulations show that, with the emergence of mechanosensitive genes through evolution enabling skeletal cells to modulate their synthetic activities, novel differentiation pathways such as endochondral ossification could have emerged, which when favoured by natural selection is maintained in a population. Furthermore, the model predicts that evolutionary forces do not lead to a single optimal mechanoregulated response but that the capacity of endochondral ossification exists with variability in a population. The simulations correspond with many existing findings about the mechanosensitivity of skeletal tissues in current animal populations, therefore indicating that this kind of multi-level models could be used in future population based simulations of tissue differentiation.

Interstitial fluid flow: the mechanical environment of cells and foundation of meridians

Using information from the deep dissection, microobservation, and measurement of acupoints in the upper and lower limbs of the human body, we developed a three-dimensional porous medium model to simulate the flow field using FLUENT software and to study the shear stress on the surface of interstitial cells (mast cells) caused by interstitial fluid flow. The numerical simulation results show the following: (i) the parallel nature of capillaries will lead to directional interstitial fluid flow, which may explain the long interstitial tissue channels or meridians observed in some experiments; (ii) when the distribution of capillaries is staggered, increases in the velocity alternate, and the velocity tends to be uniform, which is beneficial for substance exchange; (iii) interstitial fluid flow induces a shear stress, with magnitude of several Pa, on interstitial cell membranes, which will activate cells and lead to a biological response; (iv) capillary and interstitial parameters, such as capillary density, blood pressure, capillary permeability, interstitial pressure, and interstitial porosity, affect the shear stress on cell surfaces. The numerical simulation results suggest that in vivo interstitial fluid flow constitutes the mechanical environment of cells and plays a key role in guiding cell activities, which may explain the meridian phenomena and the acupuncture effects observed in experiments.

Substrate stiffness and oxygen as regulators of stem cell differentiation during skeletal tissue regeneration: a mechanobiological model

Extrinsic mechanical signals have been implicated as key regulators of mesenchymal stem cell (MSC) differentiation. It has been possible to test different hypotheses for mechano-regulated MSC differentiation by attempting to simulate regenerative events such as bone fracture repair, where repeatable spatial and temporal patterns of tissue differentiation occur. More recently, in vitro studies have identified other environmental cues such as substrate stiffness and oxygen tension as key regulators of MSC differentiation; however it remains unclear if and how such cues determine stem cell fate in vivo. As part of this study, a computational model was developed to test the hypothesis that substrate stiffness and oxygen tension regulate stem cell differentiation during fracture healing. Rather than assuming mechanical signals act directly on stem cells to determine their differentiation pathway, it is postulated that they act indirectly to regulate angiogenesis and hence partially determine the local oxygen environment within a regenerating tissue. Chondrogenesis of MSCs was hypothesized to occur in low oxygen regions, while in well vascularised regions of the regenerating tissue a soft local substrate was hypothesised to facilitate adipogenesis while a stiff substrate facilitated osteogenesis. Predictions from the model were compared to both experimental data and to predictions of a well established computational mechanobiological model where tissue differentiation is assumed to be regulated directly by the local mechanical environment. The model predicted all the major events of fracture repair, including cartilaginous bridging, endosteal and periosteal bony bridging and bone remodelling. It therefore provides support for the hypothesis that substrate stiffness and oxygen play a key role in regulating MSC fate during regenerative events such as fracture healing.

Simulation of fracture healing in the tibia: mechanoregulation of cell activity using a lattice modeling approach

In this study, a three-dimensional (3D) computational simulation of bone regeneration was performed in a human tibia under realistic muscle loading. The simulation was achieved using a discrete lattice modeling approach combined with a mechanoregulation algorithm to describe the cellular processes involved in the healing process-namely proliferation, migration, apoptosis, and differentiation of cells. The main phases of fracture healing were predicted by the simulation, including the bone resorption phase, and there was a qualitative agreement between the temporal changes in interfragmentary strain and bending stiffness by comparison to experimental data and clinical results. Bone healing was simulated beyond the reparative phase by modeling the transition of woven bone into lamellar bone. Because the simulation has been shown to work with realistic anatomical 3D geometry and muscle loading, it demonstrates the potential of simulation tools for patient-specific pre-operative treatment planning.

Inter-species investigation of the mechano-regulation of bone healing: comparison of secondary bone healing in sheep and rat

Inter-species differences in regeneration exist in various levels. One aspect is the dynamics of bone regeneration and healing, e.g. small animals show a faster healing response when compared to large animals. Mechanical as well as biological factors are known to play a key role in the process. However, it remains so far unknown whether different animals follow at all comparable mechano-biological rules during tissue regeneration, and in particular during bone healing. In this study, we investigated whether differences observed in vivo in the dynamics of bone healing between rat and sheep are only due to differences in the animal size or whether these animals have a different mechano-biological response during the healing process. Histological sections from in vivo experiments were compared to in silico predictions of a mechano-biological computer model for the simulation of bone healing. Investigations showed that the healing processes in both animal models occur under significantly different levels of mechanical stimuli within the callus region, which could explain histological observations of early intramembranous ossification at the endosteal side. A species-specific adaptation of a mechano-biological model allowed a qualitative match of model predictions with histological observations. Specifically, when keeping cell activity processes at the same rate, the amount of tissue straining defining favorable mechanical conditions for the formation of bone had to be increased in the large animal model, with respect to the small animal, to achieve a qualitative agreement of model predictions with histological data. These findings illustrate that geometrical (size) differences alone cannot explain the distinctions seen in the histological appearance of secondary bone healing in sheep and rat. It can be stated that significant differences in the mechano-biological regulation of the healing process exist between these species. Future investigations should aim towards understanding whether these differences are due to differences in cell behavior, material properties of the newly formed tissues within the callus and/or differences in response to the mechanical environment.

Variability observed in mechano-regulated in vivo tissue differentiation can be explained by variation in cell mechano-sensitivity

Computational simulations of tissue differentiation have been able to capture the main aspects of tissue formation/regeneration observed in animal experiments-except for the considerable degree of variability reported. Understanding and modelling the source of this variability is crucial if computational tools are to be developed for clinical applications. The objective of this study was to test the hypothesis that differences in cell mechano-sensitivity between individuals can explain the variability of tissue differentiation patterns observed experimentally. Simulations of an experiment of tissue differentiation in a mechanically loaded bone chamber were performed. Finite element analysis was used to determine the biophysical environment, and a lattice-modelling approach was used to simulate cell activity. Differences in cell mechano-sensitivity among individuals were modelled as differences in cell activity rates, with the activation of cell activities regulated by the mechanical environment. Predictions of the tissue distribution in the chambers produced the two different classes of results found experimentally: (i) chambers with a layer of bone across the chamber covered by a layer of cartilage on top and (ii) chambers with almost no bone, mainly fibrous tissue and small islands of cartilage. This indicates that the differing cellular response to the mechanical environment (i.e., subject-specific mechano-sensitivity) could be a reason for the different outcomes found when implants (or tissue engineered constructs) are used in a population.

Correlations between local strains and tissue phenotypes in an experimental model of skeletal healing

Defining how mechanical cues regulate tissue differentiation during skeletal healing can benefit treatment of orthopaedic injuries and may also provide insight into the influence of the mechanical environment on skeletal development. Different global (i.e., organ-level) mechanical loads applied to bone fractures or osteotomies are known to result in different healing outcomes. However, the local stimuli that promote formation of different skeletal tissues have yet to be established. Finite element analyses can estimate local stresses and strains but require many assumptions regarding tissue material properties and boundary conditions. This study used an experimental approach to investigate relationships between the strains experienced by tissues in a mechanically stimulated osteotomy gap and the patterns of tissue differentiation that occur during healing. Strains induced by the applied, global mechanical loads were quantified on the mid-sagittal plane of the callus using digital image correlation. Strain fields were then compared to the distribution of tissue phenotypes, as quantified by histomorphometry, using logistic regression. Significant and consistent associations were found between the strains experienced by a region of the callus and the tissue type present in that region. Specifically, the probability of encountering cartilage increased, and that of encountering woven bone decreased, with increasing octahedral shear strain and, to a lesser extent, maximum principal strain. Volumetric strain was the least consistent predictor of tissue type, although towards the end of the four-week stimulation timecourse, cartilage was associated with increasingly negative volumetric strains. These results indicate that shear strain may be an important regulator of tissue fate during skeletal healing.

Tissue differentiation in an in vivo bioreactor: in silico investigations of scaffold stiffness

Scaffold design remains a main challenge in tissue engineering due to the large number of requirements that need to be met in order to create functional tissues in vivo. Computer simulations of tissue differentiation within scaffolds could serve as a powerful tool in elucidating the design requirements for scaffolds in tissue engineering. In this study, a lattice-based model of a 3D porous scaffold construct derived from micro CT and a mechano-biological simulation of a bone chamber experiment were combined to investigate the effect of scaffold stiffness on tissue differentiation inside the chamber. The results indicate that higher scaffold stiffness, holding pore structure constant, enhances bone formation. This study demonstrates that a lattice approach is very suitable for modelling scaffolds in mechano-biological simulations, since it can accurately represent the micro-porous geometries of scaffolds in a 3D environment and reduce computational costs at the same time.

Schlafen 3 induction by cyclic strain regulates intestinal epithelial differentiation

The intestinal epithelium is subjected to repetitive deformation during normal gut function by peristalsis and villous motility. In vitro, cyclic strain promotes intestinal epithelial proliferation and induces an absorptive phenotype characterized by increased dipeptidyl dipeptidase (DPPIV) expression. Schlafen 3 is a novel gene recently associated with cellular differentiation. We sought to evaluate whether Schlafen 3 mediates the effects of strain on the differentiation of intestinal epithelial cell (IEC)-6 in the absence or presence of cyclic strain. Strain increased Schlafen 3 mRNA and protein. In cells transfected with a control-nontargeting siRNA, strain increased DPPIV-specific activity. However, Schlafen 3 reduction by siRNA decreased basal DPPIV and prevented any stimulation of DPPIV activity by strain. Schlafen 3 reduction also prevented DPPIV induction by sodium butyrate (1 mM) or transforming growth factor (TGF)-beta (0.1 ng/ml), two unrelated differentiating stimuli. However, Schlafen-3 reduction by siRNA did not prevent the mitogenic effect of strain or that of EGF. Blocking Src and phosphatidyl inositol (PI3)-kinase prevented strain induction of Schlafen 3, but Schlafen 3 induction required activation of p38 but not ERK. These results suggest that cyclic strain induces an absorptive phenotype characterized by increased DPPIV activity via Src-, p38-, and PI3-kinase-dependent induction of Schlafen 3 in rat IEC-6 cells on collagen, whereas Schlafen 3 may also be a key factor in the induction of intestinal epithelial differentiation by other stimuli such as sodium butyrate or TGF-beta. The induction of Schlafen 3 or its human homologs may modulate intestinal epithelial differentiation and preserve the gut mucosa during normal gut function.

A multichamber fluidic device for 3D cultures under interstitial flow with live imaging: development, characterization, and applications

Interstitial flow is an important biophysical cue that can affect capillary morphogenesis, tumor cell migration, and fibroblast remodeling of the extracellular matrix, among others. Current models that incorporate interstitial flow and that are suitable for live imaging lack the ability to perform multiple simultaneous experiments, for example, to compare effects of growth factors, extracellular matrix composition, etc. We present a nine-chamber radial flow device that allows simultaneous 3D fluidic experiments for relatively long-term culture with live imaging capabilities. Flow velocity profiles were characterized by fluorescence recovery after photobleaching (FRAP) for flow uniformity and estimating the hydraulic conductivity. We demonstrate lymphatic and blood capillary morphogenesis in fibrin gels over 10 days, comparing flow with static conditions as well as the effects of an engineered variant of VEGF that binds fibrin via Factor XIII. We also demonstrate the culture of contractile fibroblasts and co-cultures with tumor cells for modeling the tumor microenvironment. Therefore, this device is useful for studies of capillary morphogenesis, cell migration, contractile cells like fibroblasts, and multicellular cultures, all under interstitial flow.

Mechano-regulation of mesenchymal stem cell differentiation and collagen organisation during skeletal tissue repair

A number of mechano-regulation theories have been proposed that relate the differentiation pathway of mesenchymal stem cells (MSCs) to their local biomechanical environment. During spontaneous repair processes in skeletal tissues, the organisation of the extracellular matrix is a key determinant of its mechanical fitness. In this paper, we extend the mechano-regulation theory proposed by Prendergast et al. (J Biomech 30(6):539-548, 1997) to include the role of the mechanical environment on the collagen architecture in regenerating soft tissues. A large strain anisotropic poroelastic material model is used in a simulation of tissue differentiation in a fracture subject to cyclic bending (Cullinane et al. in J Orthop Res 20(3):579-586, 2002). The model predicts non-union with cartilage and fibrous tissue formation in the defect. Predicted collagen fibre angles, as determined by the principal decomposition of strain- and stress-type tensors, are similar to the architecture seen in native articular cartilage and neoarthroses induced by bending of mid-femoral defects in rats. Both stress and strain-based remodelling stimuli successfully predicted the general patterns of collagen fibre organisation observed in vivo. This provides further evidence that collagen organisation during tissue differentiation is determined by the mechanical environment. It is envisioned that such predictive models can play a key role in optimising MSC-based skeletal repair therapies where recapitulation of the normal tissue architecture is critical to successful repair.

Corroboration of mechanobiological simulations of tissue differentiation in an in vivo bone chamber using a lattice-modeling approach

It is well established that the mechanical environment modulates tissue differentiation, and a number of mechanoregulatory theories for describing the process have been proposed. In this study, simulations of an in vivo bone chamber experiment were performed that allowed direct comparison with experimental data. A mechanoregulation theory for mesenchymal stem cell differentiation based on a combination of fluid flow and shear strain (computed using finite element analysis) was implemented to predict tissue differentiation inside mechanically controlled bone chambers inserted into rat tibae. To simulate cell activity, a lattice approach with stochastic cell migration, proliferation, and selected differentiation was adopted; because of its stochastic nature, each run of the simulation gave a somewhat different result. Simulations predicted the load-dependency of the tissue differentiation inside the chamber and a qualitative agreement with histological data; however, the full variability found between specimens in the experiment could not be predicted by the mechanoregulation algorithm. This result raises the question whether tissue differentiation predictions can be linked to genetic variability in animal populations.