The impact of osteoporosis and diabetes on fracture healing under different loading conditions

BACKGROUND

Osteoporosis and diabetes are two prevalent conditions among the elderly population. Each of these conditions can profoundly influence the fracture healing process by disturbing the associated inflammatory process. However, the combined effects of osteoporosis and diabetes on fracture healing remain unclear. Therefore, the purpose of the present study is to investigate the role of osteoporosis and diabetes in fracture healing and the underlying mechanisms by developing numerical models.

METHOD

This study introduces a numerical model that consists of a three-dimensional model of a tibia fracture stabilized by a Locking Compression Plate (LCP), coupled with a two-dimensional axisymmetric model which illustrates the transport and reactions of cells and cytokines throughout the inflammatory phase in early fracture healing. First, the model parameters were calibrated using available experimental data. The model was then implemented to predict the healing outcomes of fractures under five varied conditions, consisting of both osteoporotic and non-osteoporotic bones, each subjected to different physiological loads.

RESULTS

The instability of the fracture callus can significantly escalate in osteoporotic fractures (e.g., when a 150 N physiological load is applied, the unstable region of the osteoporotic fracture callus can reach 26 %, in contrast to 12 % in non-osteoporotic fractures). Additionally, the mesenchymal stem cells (MSCs) proliferation and differentiation can be disrupted in osteoporotic fracture compared to non-osteoporotic fractures (e.g., on the 10th day post-fracture, the decrease in the concentration of MSCs, osteoblasts, and chondrocytes in osteoporotic fractures is nearly double that in non-osteoporotic fractures under a 150 N). Finally, the healing process of fractures can suffer significant impairment when osteoporosis coexists with diabetes (e.g., the concentration of MSCs can be drastically reduced by nearly 37 % in osteoporotic fractures under diabetic conditions when subjected to a load of 200 N) CONCLUSIONS: Fracture calluses destabilized by osteoporosis can negatively affect the fracture healing process by disrupting the proliferation and differentiation of mesenchymal stem cells (MSCs). Moreover, when osteoporosis coexists with diabetes, the fracture healing process can severely impair the fracture healing outcomes.

Depth camera-based model for studying the effects of muscle loading on distal radius fracture healing

BACKGROUND

Distal radius fractures (DRFs) treated with volar locking plates (VLPs) allows early rehabilitation exercises favourable to fracture recovery. However, the role of rehabilitation exercises induced muscle forces on the biomechanical microenvironment at the fracture site remains to be fully explored. The purpose of this study is to investigate the effects of muscle forces on DRF healing by developing a depth camera-based fracture healing model.

METHOD

First, the rehabilitation-related hand motions were captured by a depth camera system. A macro-musculoskeletal model is then developed to analyse the data captured by the system for estimating hand muscle and joint reaction forces which are used as inputs for our previously developed DRF model to predict the tissue differentiation patterns at the fracture site. Finally, the effect of different wrist motions (e.g., from 60° of extension to 60° of flexion) on the DRF healing outcomes will be studied.

RESULTS

Muscle and joint reaction forces in hands which are highly dependent on hand motions could significantly affect DRF healing through imposed compressive and bending forces at the fracture site. There is an optimal range of wrist motion (i.e., between 40° of extension and 40° of flexion) which could promote mechanical stimuli governed healing while mitigating the risk of bony non-union due to excessive movement at the fracture site.

CONCLUSION

The developed depth camera-based fracture healing model can accurately predict the influence of muscle loading induced by rehabilitation exercises in distal radius fracture healing outcomes. The outcomes from this study could potentially assist osteopathic surgeons in designing effective post-operative rehabilitation strategies for DRF patients.

Osteochondral junction leakage and cartilage joint lubrication

BACKGROUND AND OBJECTIVES

Previous studies have shown that there is potentially interstitial fluid exchange between cartilage tissue and the subarticular spongiosa region in the case of injury or disease (e.g., osteoarthritis and osteoporosis). Interstitial flow is also required for cartilage lubrication under joint load. A key question then is how cartilage lubrication is modified by increased interstitial fluid leakage across the osteochondral junction. Thus, the purpose of this study is to develop a numerical model to investigate changes in cartilage lubrication with changes in osteochondral junction leakage.

METHODS

The multi-phase coupled model includes domains corresponding to the contact gap, cartilage tissue and subchondral bone plate region (ScBP). Each of these domains are treated as poroelastic systems, with their coupling implemented through mass and pressure continuity. The effects of osteochondral junction leakage on lubrication were investigated with a parametric study on the relative permeability between the ScBP and cartilage tissue.

RESULTS

Significant effects of ScBP permeability were predicted, especially during the early stage of the junction leakage development (early stage of the disease). There is a significant reduction in mixed-mode lubrication duration under the effect of increased junction leakage (the cartilage tissue mixed-mode lubrication duration is about 33% decrease for a relative permeability ratio of 0.1 between ScBP and cartilage tissue, and about 52% decrease under the osteoarthritis condition). In addition, the time for cartilage to reach steady-state consolidation is significantly reduced when ScBP permeability increases (the consolidation time reduces from roughly 2 h to 1.2 h when the relative permeability ratio increases from 0.001 to 0.1, and it reduces to 0.8 h for an advanced osteoarthritis condition). It is predicted that the initial friction coefficient could increase by over 60% when the ScBP permeability is consistent with an advanced osteoarthritis (OA) condition.

CONCLUSION

Increased osteochondral junction leakage induced by joint injury and disease could result in increased cartilage surface wear rates due to more rapid interstitial fluid depressurization within articular cartilage.

The effects of mechanical instability on PDGF mediated inflammatory response at early stage of fracture healing under diabetic condition

BACKGROUND AND OBJECTIVE

Mechanical stability plays an important role in fracture healing process. Excessive interfragmentary movement will continuously damage the tissue and newly formed capillaries at the fracture site, which leads to overproduction of platelet-derived growth factor (PDGF) that attracts more macrophages into fracture callus, ultimately persistent and enhanced inflammatory response happens. For diabetic condition, the impact of mechanical instability of fracture site on inflammatory response could be further compliciated and the relevant research in this field is relatively limited.

METHODS

Building on previous experimental studies, this study presents a numerical model consisting of a system of reactive-transport equations representing the transport as well as interactions of different cells and cytokines within the fracture callus. The model is initially validated by available experimental data, and then implemented to investigate the role of mechanical stability of fracture site in inflammatory response during early stage of healing. It is assumed that there is an increased release of PDGF due to the rupture of blood vessels resulting from mechanical instability, which leads to increased production of inflammatory cytokines (i.e., TNF-α). The bone healing process under three different conditions were investigated, i.e., mechanically stable condition with normal inflammatory response (Control, Case 1), mechanically unstable condition with normal inflammatory response (Case 2) and mechanically unstable condition with diabetes (Case 3).

RESULTS

Mechanical instability can promote the macrophage infiltration and thus induce an enhanced and prolonged inflammatory response, which could impede the MSCs proliferation during the early fracture healing stage (e.g., compared with the control condition, the MSCs concentration in unstable fracture with normal inflammatory response can be reduced by 3.2% and 5.2% on day 2 and day 10 post-fracture, respectively). Under diabetic condition, the mechanical instability of fracture site could lead to a significant increase of TNF-α concentration in fracture callus (Case 3) in comparison to control (Case 1) (e.g., three-fold increase in TNF-α concentration compared to control). In addition, the results show that the mechanical instability affects the cell differentiation and proliferation in fracture callus in a spatially dependent manner, e.g., for diabetic fracture patients, the mechanical instability could potentially decrease the concentration of MSCs, osteoblasts and chondrocytes by around 39%, 30% and 29% in cortical callus, respectively, in comparison to control.

CONCLUSION

The mechanical instability together with diabetic condition can significantly affect the natural resolution of inflammation during early stage of healing by turning acute inflammation into chronic inflammation which is characterized by a continuously upregulated TNF-α pathway.

Computational Modelling for Managing Pathways to Cartilage Failure

Over several decades the perception and therefore description of articular cartilage changed substantially. It has transitioned from being described as a relatively inert tissue with limited repair capacity, to a tissue undergoing continuous maintenance and even adaption, through a range of complex regulatory processes. Even from the narrower lens of biomechanics, the engagement with articular cartilage has changed from it being an interesting, slippery material found in the hostile mechanical environment between opposing long bones, to an intriguing example of mechanobiology in action. The progress revealing this complexity, where physics, chemistry, material science and biology are merging, has been described with increasingly sophisticated computational models. Here we describe how these computational models of cartilage as an integrated system can be combined with the approach of structural reliability analysis. That is, causal, deterministic models placed in the framework of the probabilistic approach of structural reliability analysis could be used to understand, predict, and mitigate the risk of cartilage failure or pathology. At the heart of this approach is seeing cartilage overuse and disease processes as a 'material failure', resulting in failure to perform its function, which is largely mechanical. One can then describe pathways to failure, for example, how homeostatic repair processes can be overwhelmed leading to a compromised tissue. To illustrate this 'pathways to failure' approach, we use the interplay between cartilage consolidation and lubrication to analyse the increase in expected wear rates associated with cartilage defects or meniscectomy.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Linear and Nonlinear Biphasic Mechanical Properties of Goat IVDs Under Different Swelling Conditions in Confined Compression

To define technical specifications for artificial substitutes, it is necessary to model their mechanical behaviour. Here we studied the linear and nonlinear biphasic models for Nucleus Pulposus (NP) and Annulus Fibrosus (AF). The associated material parameters were obtained using confined compression stress relaxation tests on goat intervertebral disc (IVD) samples. The first parameter, aggregate modulus H_A0, which essentially describes load-bearing capacity of the solid phase, was larger for AF (H_A0 = 0.53 ± 0.06 MPa) than for NP (H_A0 = 0.26 ± 0.04 MPa). For hydraulic permeability, which quantifies the ability to transmit interstitial fluid, it was the opposite (k₀ = (0.20 ± 0.07) × 10^-15 m⁴/Ns for AF and k₀ = (0.67 ± 0.08)×10^-15 m⁴/Ns for NP). The values of nonlinearity coefficients, nonlinear stiffening coefficient β and non-dimensional nonlinear permeability coefficient M, reflected that these tissues had nonlinear elastic behaviour and permeability. Also, investigating the effect of swelling conditions in sample preparation showed that for both AF and NP, confined-swollen samples had higher aggregate modulus and lower permeability values compared to the free-swollen ones. The quantitative description of the nonlinear properties of AF and NP provided a better understanding of IVD behaviour as well as technical specifications for their artificial substitutes.

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

A Probabilistic Failure Risk Approach to The Problem of Articular Cartilage Lubrication

BACKGROUND AND OBJECTIVE

The geometrical and mechanical properties that characterise the cartilage contact gap are uncertain and spatially varied. To date the effects of such uncertainties on cartilage lubrication have not been explored. Using a probabilistic approach, the purpose of this study is to numerically investigate the influence of surficial cartilage glycoaminoglycan (GAG) content on joint lubrication behaviour. Gap asperity stiffness and polymer brush border (PBB) thickness are affected by the uncertainty of surficial GAG concentration, and so their correlated effects in maintaining hydrodynamic joint lubrication are investigated.

METHODS

Correlated sampling data are first generated by Monte Carlo simulation. These data are used as inputs for the cartilage contact model, which includes three distinctive features of cartilage tissue (tension-compression nonlinearity, aggrecan dependent permeability and compressive modulus) and fluid flow resistance effects of PBB on cartilage surface. The degree of hydrodynamic lubrication after thirty minutes of constant loading is used as an indicator for assessing the lubrication performance at the contact interface.

RESULTS

The increase of PBB thickness with GAG concentration enhances the hydrodynamic lubrication component in the cartilage contact gap, whereas increasing the asperity stiffness with GAG concentration impairs hydrodynamic lubrication. GAG loss rate increases with the rise of GAG concentration. More aggrecan shedding through the surface could result in a thicker and denser PBB, and therefore enhance the lubrication performance in mixed-mode regime. On the other hand, higher GAG content makes the asperities stiffer, which may impede contact gap closure, and thus encourage gap fluid loss and impair the lubrication performance of cartilage.

CONCLUSION

The lubrication performance of cartilage varies with the physiological conditions of the joint. Since a range of variables are internally related, the outcomes on joint lubrication are difficult to predict. A probabilistic approach accounting for the uncertainties can potentially result in more accurate evaluations of joint lubrication performance.

Modeling Riverbank Slope Reinforcement Using Anti-Slide Piles with Geocells

Geocells are increasingly used in engineering applications, but the design of riverbank slope reinforcements that use only geocells limits reinforcement performance. Moreover, the design and use of anti-slide piles with geocells are mainly based on experiences that are unsupported by theoretical models. In this paper, by combining the confinement effect and vertical action mechanism of geocells, the horizontal friction mechanism of the geocell layer and the vertical support mechanism of piles, a theoretical model of riverbank slope reinforced by anti-slide piles with geocells was constructed. In addition, to describe the mechanical behavior of a riverbank slope reinforced by anti-slide piles with geocells, the slip-resisting mechanism of the anti-slide pile with interaction between geocells and their internal filler is considered in the model. Furthermore, to investigate the influence of changes in water level on riverbank slope stability, the developed model takes into account settlement, lateral displacement, pile bending moment and pile axial force. The model predications were validated by the field measurement data. The results from a series of parametric studies show that the use of anti-slide pile and geocells can effectively reduce the settlement and the lateral displacement of a riverbank slope. The developed model could contribute to an optimal design of anti-slide pile with geocells for enhancing the stability of a riverbank slope.

A numerical study on tumor-on-chip performance and its optimization for nanodrug-based combination therapy

Microfluidic devices, such as the tumor-on-a-chip (ToC), allow for the delivery of multiple drugs as desired for various therapies such as cancer treatment. Due to the complexity involved, visualizing, and gaining knowledge of the performance of such devices through experimentation alone is difficult if not impossible. In this paper, we performed a numerical simulation study on ToC performance, which focuses on the ability to combine multiple nanodrugs and optimized ToC performance. The numerical simulations of the chip performance were performed based on the typical chip design and operating parameters, as well as the established governing equations, boundary conditions, and fluid-structure interaction. The effect of cell injection time and position, inlet flow rate, number of inlets, medium viscosity, and cell concentration on the chip performance in terms of shear stress and cell distribution were examined. The results illustrate the profound effect of operation parameters, thus allowing for rigorously determining operational parameters to prevent spheroids ejection from microwells and to restrict the shear stresses within a physiological range. Also, the results show that triple-inlets can increase the uniformity of cell distribution in comparison with single or double inlets. Based on the simulation results, the architecture of the primary ToC was further optimized, resulting in a novel design that enables applying multiple, yet simultaneous, nanodrugs with optimal drug combination as desired for an individual patient. Furthermore, our simulations on the optimized chip showed a uniform cell distribution required for uniform-sized tumor spheroids generation, and complete medium exchange. Taken together, this study not only illustrates that numerical simulations are effective to visualize the ToCs performance, but also develops a novel ToC design optimized for nanodrug-based combination therapy.

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

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

A computational model for investigating the effects of changes in bioavailability of insulin-like growth factor-1 on the homeostasis of the intervertebral disc

Insulin-like growth factor-1 (IGF-1) is well-known for upregulating cell proliferation and biosynthesis of the extracellular matrix in the intervertebral disc (IVD). Pathological conditions, such as obesity or chronic kidney disease cause IGF-1 deficiency in plasma. How this deficiency impacts disc homeostasis remains unknown. Pro-anabolic approaches for the treatment of disc degeneration based on enhancing IGF-1 bioavailability to tissue-cells are considered, but knowledge of their effectiveness in enhancing cellular anabolism of a degenerated disc is limited. In this study, we developed a computational model for disc homeostasis specifically addressing the role of IGF-1 in modulating both extracellular matrix biosynthesis and cellularity in the IVD. This model was applied to investigate how changes in IGF-1 bioavailability, namely deficiency or enhancement of growth factor, affect disc health. In this study, it was found that IGF-1 deficiency mainly affects the biosynthesis of ECM components, especially in the most external regions of the IVD such as the cartilage endplates and the outer portion of annulus fibrosus. Also, a total of three approaches for increasing IGF-1 bioavailability as a therapy for degenerated IVDs were investigated. It was found that all these strategies are only beneficial to those disc regions receiving sufficient nutritional supply (i.e., the outmost IVD regions), while they exacerbate tissue degradation in malnourished regions (i.e., inner portion of the disc). This suggests that pro-anabolic growth factor-based therapies are limited in that their success strongly depends on an adequate nutritional supply to the IVD tissue, which is not guaranteed in degenerated discs.

MODEL-BASED ANALYSIS OF IGF-1 EFFECT ON OSTEOBLAST AND OSTEOCLAST REGULATION IN BONE TURNOVER

The main determinant of bone Ca accretion is a bimolecular regulatory network on osteoblast (OB) and osteoclast (OC). Even though IGF-1 is known as an important regulator in bone cell cycle, little has been done to model IGF-1 action in bone cell regulation. Thus, the objective is to develop a mathematical model that depicts the regulatory action of IGF-1 onto the OB and OC interaction, and to evaluate adolescent and adult bone Ca accretion in response to differences in IGF-1 levels. As a result, a dynamic model of OB and OC with two main regulatory systems, i.e., Receptor Activator for Nuclear Factor [Formula: see text]B (RANK)-RANK Ligand (RANKL)-osteoprogerin (OPG) system, and TGF-[Formula: see text], was augmented with the IGF-1, and incorporated into Ca kinetic data to predict exchangeable bone Ca. The developed model could predict a change in OB and OC levels in response to perturbations in regulators, producing results consistent with bone physiology and published experimental data. The model also estimated parametric difference in regulators between adults and adolescents, suggesting that RANKL/OPG in adolescents was about 4 times higher than in adults, while adolescent serum PTH and IGF-1 concentrations were 60% and 220% of those of adults, respectively. This study highlighted the influence of IGF-1 on the regulation of bone cells in positively modulating bone Ca, suggesting that IGF-1 may be an effective target for reducing bone loss by promoting mature OB.

Computational Simulation of Mechanical Microenvironment of Early Stage of Bone Healing under Locking Compression Plate with Dynamic Locking Screws

It is well known that bone healing outcomes highly depend on the mechanical microenvironment of the fracture site, and a certain degree of interfragmentary movement (IFM) is essential for indirect (i.e. natural) bone healing. The application of locking compression plate (LCP) internal fixation in the treatment of bone fracture is a common practice which leads to early mobility and full function of the fractured extremity. However should the fixation configuration be too stiff, it might result in delayed healing or asymmetric tissue development across the fracture site due to the fact that IFM in near cortex area is too small to promote healing. Dynamic locking screw (DLS) has been recently designed to tackle this problem by reducing the stiffness of LCP fixation. However, the actual mechano-regulation mechanisms in which DLS uses to regulate the healing process are still not fully understood. The objective of this paper is to develop a computational model to understand the change of mechanical microenvironment of fracture site under LCP with dynamic locking screw in comparison to standard locking screw, and how this change could potentially regulate tissue development within the fracture callus during the healing process.Keywords: bone healing, locking compression plate, dynamic locking screw, finite element modelling

Computational simulation of the early stage of bone healing under different configurations of locking compression plates

Flexible fixation or the so-called 'biological fixation' has been shown to encourage the formation of fracture callus, leading to better healing outcomes. However, the nature of the relationship between the degree of mechanical stability provided by a flexible fixation and the optimal healing outcomes has not been fully understood. In this study, we have developed a validated quantitative model to predict how cells in fracture callus might respond to change in their mechanical microenvironment due to different configurations of locking compression plate (LCP) in clinical practice, particularly in the early stage of healing. The model predicts that increasing flexibility of the LCP by changing the bone-plate distance (BPD) or the plate working length (WL) could enhance interfragmentary strain in the presence of a relatively large gap size (> 3 mm). Furthermore, conventional LCP normally results in asymmetric tissue development during early stage of callus formation, and the increase of BPD or WL is insufficient to alleviate this problem.

Modeling the Insulin-Like Growth Factor System in Articular Cartilage

IGF signaling is involved in cell proliferation, differentiation and apoptosis in a wide range of tissues, both normal and diseased, and so IGF-IR has been the focus of intense interest as a promising drug target. In this computational study on cartilage, we focus on two questions: (i) what are the key factors influencing IGF-IR complex formation, and (ii) how might cells regulate IGF-IR complex formation? We develop a reaction-diffusion computational model of the IGF system involving twenty three parameters. A series of parametric and sensitivity studies are used to identify the key factors influencing IGF signaling. From the model we predict the free IGF and IGF-IR complex concentrations throughout the tissue. We estimate the degradation half-lives of free IGF-I and IGFBPs in normal cartilage to be 20 and 100 mins respectively, and conclude that regulation of the IGF half-life, either directly or indirectly via extracellular matrix IGF-BP protease concentrations, are two critical factors governing the IGF-IR complex formation in the cartilage. Further we find that cellular regulation of IGF-II production, the IGF-IIR concentration and its clearance rate, all significantly influence IGF signaling. It is likely that negative feedback processes via regulation of these factors tune IGF signaling within a tissue, which may help explain the recent failures of single target drug therapies aimed at modifying IGF signaling.