Modeling the Insulin-Like Growth Factor System in Articular Cartilage

IGF signaling pathway in bone and cartilage development, homeostasis, and disease

The skeleton plays a fundamental role in the maintenance of organ function and daily activities. The insulin-like growth factor (IGF) family is a group of polypeptide substances with a pronounced role in osteoblast differentiation, bone development, and metabolism. Disturbance of the IGFs and the IGF signaling pathway is inextricably linked with assorted developmental defects, growth irregularities, and jeopardized skeletal structure. Recent findings have illustrated the significance of the action of the IGF signaling pathway via growth factors and receptors and its interactions with dissimilar signaling pathways (Wnt/β-catenin, BMP, TGF-β, and Hh/PTH signaling pathways) in promoting the growth, survival, and differentiation of osteoblasts. IGF signaling also exhibits profound influences on cartilage and bone development and skeletal homeostasis via versatile cell-cell interactions in an autocrine, paracrine, and endocrine manner systemically and locally. Our review summarizes the role and regulatory function as well as a potentially integrated gene network of the IGF signaling pathway with other signaling pathways in bone and cartilage development and skeletal homeostasis, which in turn provides an enlightening insight into visualizing bright molecular targets to be eligible for designing effective drugs to handle bone diseases and maladies, such as osteoporosis, osteoarthritis, and dwarfism.

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.

Effect of uncertain clinical conditions on the early healing and stability of distal radius fractures

BACKGROUND AND OBJECTIVES

The healing outcomes of distal radius fracture (DRF) treated with the volar locking plate (VLP) depend on surgical strategies and postoperative rehabilitation. However, the accurate prediction of healing outcomes is challenging due to a range of certainties related to the clinical conditions of DRF patients, including fracture geometry, fixation configuration, and physiological loading. The purpose of this study is to investigate the influence of uncertainty and variability in fracture/fixation parameters on the mechano-biology and biomechanical stability of DRF, using a probabilistic numerical approach based on the results from a series of experimental tests performed in this study.

METHODS

Six composite radius sawboneses fitted with titanium VLP (VLP 2.0, Austofix) were loaded to failure at a rate of 2 N/s. The testing results of the elastic and plastic behaviour of the VLP were used as inputs for a probabilistic-based computational model of DRF, which simulated mechano-regulated tissue differentiation and fixation elastic capacity at the fracture site. Finally, the probability of success in early indirect healing and fracture stabilisation was predicted.

RESULTS

The titanium VLP is a strong and ductile fixation whose flexibility and elastic capacity are governed by flexion working length and bone-to-plate distance, respectively. A fixation with optimised designs and configurations is critical to mechanically stabilising the early fracture site. Importantly, the uncertainty and variability in fracture/fixation parameters could compromise early DRF healing. The physiological loading uncertainty is the most adverse factor, followed by the negative impact of uncertainty in fracture geometry.

CONCLUSIONS

The VRP 2.0 fixation made of grade II titanium is a desirable fixation that is strong enough to resist irreparable deformation during early recovery and is also ductile to deform plastically without implant failure at late rehabilitation.

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.

Histone H3K27 demethylase UTX compromises articular chondrocyte anabolism and aggravates osteoarthritic degeneration

Epigenome alteration in chondrocytes correlates with osteoarthritis (OA) development. H3K27me3 demethylase UTX regulates tissue homeostasis and deterioration, while its role was not yet studied in articulating joint tissue in situ. We now uncovered that increased UTX and H3K27me3 expression in articular chondrocytes positively correlated with human knee OA. Forced UTX expression upregulated the H3K27me3 enrichment at transcription factor Sox9 promoter, inhibiting key extracellular matrix molecules collagen II, aggrecan, and glycosaminoglycan in articular chondrocytes. Utx overexpression in knee joints aggravated the signs of OA, including articular cartilage damage, synovitis, osteophyte formation, and subchondral bone loss in mice. Chondrocyte-specific Utx knockout mice developed thicker articular cartilage than wild-type mice and showed few gonarthrotic symptoms during destabilized medial meniscus- and collagenase-induced joint injury. In vitro, Utx loss changed H3K27me3-binding epigenomic landscapes, which contributed to mitochondrial activity, cellular senescence, and cartilage development. Insulin-like growth factor 2 (Igf2) and polycomb repressive complex 2 (PRC2) core components Eed and Suz12 were, among others, functional target genes of Utx. Specifically, Utx deletion promoted Tfam transcription, mitochondrial respiration, ATP production and Igf2 transcription but inhibited Eed and Suz12 expression. Igf2 blockade or forced Eed or Suz12 expression increased H3K27 trimethylation and H3K27me3 enrichment at Sox9 promoter, compromising Utx loss-induced extracellular matrix overproduction. Taken together, UTX repressed articular chondrocytic activity, accelerating cartilage loss during OA. Utx loss promoted cartilage integrity through epigenetic stimulation of mitochondrial biogenesis and Igf2 transcription. This study highlighted a novel noncanonical role of Utx, in concert with PRC2 core components, in controlling H3K27 trimethylation and articular chondrocyte anabolism and OA development.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Computational Modeling Intervertebral Disc Pathophysiology: A Review

Lower back pain is a medical condition of epidemic proportion, and the degeneration of the intervertebral disc has been identified as a major contributor. The etiology of intervertebral disc (IVD) degeneration is multifactorial, depending on age, cell-mediated molecular degradation processes and genetics, which is accelerated by traumatic or gradual mechanical factors. The complexity of such intertwined biochemical and mechanical processes leading to degeneration makes it difficult to quantitatively identify cause–effect relationships through experiments. Computational modeling of the IVD is a powerful investigative tool since it offers the opportunity to vary, observe and isolate the effects of a wide range of phenomena involved in the degenerative process of discs. This review aims at discussing the main findings of finite element models of IVD pathophysiology with a special focus on the different factors contributing to physical changes typical of degenerative phenomena. Models presented are subdivided into those addressing role of nutritional supply, progressive biochemical alterations stemming from an imbalance between anabolic and catabolic processes, aging and those considering mechanical factors as the primary source that induces morphological change within the disc. Limitations of the current models, as well as opportunities for future computational modeling work are also discussed.

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.

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.

Modelling the role of redox-related mechanisms in musculoskeletal ageing

The decline in the musculoskeletal system with age is driven at the cellular level by random molecular damage. Cells possess mechanisms to repair or remove damage and many of the pathways involved in this response are regulated by redox signals. However, with ageing there is an increase in oxidative stress which can lead to chronic inflammation and disruption of redox signalling pathways. The complexity of the processes involved has led to the use of computational modelling to help increase our understanding of the system, test hypotheses and make testable predictions. This paper will give a brief background of the biological systems that have been modelled, an introduction to computational modelling, a review of models that involve redox-related mechanisms that are applicable to musculoskeletal ageing, and finally a discussion of the future potential for modelling in this field.

The role of some ADAM-proteins and activation of the insulin growth factor-related pathway in colorectal cancer

Colorectal cancer (CRC) is the third most common malignant neoplasm worldwide. In Poland, colorectal cancer ranks second in tumor incidence regardless of sex; moreover, there has been a steady increase in the incidence of CRC. CRC results from complex interactions between inherited susceptibility, clinical conditions and environmental/lifestyle-related risk factors such as physical inactivity, smoking, alcohol consumption, high-fat/low-fiber diet, and obesity/overweight. The activation of pathways associated with insulin resistance and insulin-like growth factors (IGF) appears to be the epidemiological link between the metabolic syndrome and the development of CRC, which is of particular importance. What is significantly associated with the pathway of IGF is ADAM12 and 28-protein, which belong to a broad family of the adamalysines. These proteins, by adjusting the bioavailability of growth factors, influence the process of carcinogenesis. The aim of this article is to analyze the role of selected adamalysines and activation of the IGF system associated with the formation of colon cancer.

Systems Based Study of the Therapeutic Potential of Small Charged Molecules for the Inhibition of IL-1 Mediated Cartilage Degradation

Inflammatory cytokines are key drivers of cartilage degradation in post-traumatic osteoarthritis. Cartilage degradation mediated by these inflammatory cytokines has been extensively investigated using in vitro experimental systems. Based on one such study, we have developed a computational model to quantitatively assess the impact of charged small molecules intended to inhibit IL-1 mediated cartilage degradation. We primarily focus on the simplest possible computational model of small molecular interaction with the IL-1 system-direct binding of the small molecule to the active site on the IL-1 molecule itself. We first use the model to explore the uptake and release kinetics of the small molecule inhibitor by cartilage tissue. Our results show that negatively charged small molecules are excluded from the negatively charged cartilage tissue and have uptake kinetics in the order of hours. In contrast, the positively charged small molecules are drawn into the cartilage with uptake and release timescales ranging from hours to days. Using our calibrated computational model, we subsequently explore the effect of small molecule charge and binding constant on the rate of cartilage degradation. The results from this analysis indicate that the small molecules are most effective in inhibiting cartilage degradation if they are either positively charged and/or bind strongly to IL-1α, or both. Furthermore, our results showed that the cartilage structural homeostasis can be restored by the small molecule if administered within six days following initial tissue exposure to IL-1α. We finally extended the scope of the computational model by simulating the competitive inhibition of cartilage degradation by the small molecule. Results from this model show that small molecules are more efficient in inhibiting cartilage degradation by binding directly to IL-1α rather than binding to IL-1α receptors. The results from this study can be used as a template for the design and development of more pharmacologically effective osteoarthritis drugs, and to investigate possible therapeutic options.

Mesenchymal Stem/Stromal Cells seeded on cartilaginous endplates promote Intervertebral Disc Regeneration through Extracellular Matrix Remodeling

Intervertebral disc (IVD) degeneration is characterized by significant biochemical and histomorphological alterations, such as loss of extracellular matrix (ECM) integrity, by abnormal synthesis of ECM main components, resultant from altered anabolic/catabolic cell activities and cell death. Mesenchymal Stem/Stromal Cell (MSC) migration towards degenerated IVD may represent a viable strategy to promote tissue repair/regeneration. Here, human MSCs (hMSCs) were seeded on top of cartilaginous endplates (CEP) of nucleotomized IVDs of bovine origin and cultured ex vivo up to 3 weeks. hMSCs migrated from CEP towards the lesion area and significantly increased expression of collagen type II and aggrecan in IVD, namely in the nucleus pulposus. Concomitantly, hMSCs stimulated the production of growth factors, promoters of ECM synthesis, such as fibroblast growth factor 6 (FGF-6) and 7 (FGF-7), platelet-derived growth factor receptor (PDGF-R), granulocyte-macrophage colony-stimulating factor (GM-CSF) and insulin-like growth factor 1 receptor (IGF-1sR). Overall, our results demonstrate that CEP can be an alternative route to MSC-based therapies for IVD regeneration through ECM remodeling, thus opening new perspectives on endogenous repair capacity through MSC recruitment.

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.

Quantitative analysis of insulin-like growth factor 2 receptor and insulin-like growth factor binding proteins to identify control mechanisms for insulin-like growth factor 1 receptor phosphorylation

Background

The insulin-like growth factor (IGF) system impacts cellular development by regulating proliferation, differentiation, and apoptosis, and is an attractive therapeutic target in cancer. The IGF system is complex, with two ligands (IGF1, IGF2), two receptors (IGF1R, IGF2R), and at least six high affinity IGF-binding proteins (IGFBPs) that regulate IGF ligand bioavailability. While the individual components of the IGF system are well studied, the question of how these different components integrate as a system to regulate cell behavior is less clear.

Results

To analyze the relative importance of different mechanisms that control IGF network activity, we developed a mass-action kinetic model incorporating cell surface binding, phosphorylation, and intracellular trafficking events. The model was calibrated and validated using experimental data collected from OVCAR5, an immortalized ovarian cancer cell line. We then performed model analysis to examine the ability of IGF2R or IGFBPs to counteract phosphorylation of IGF1R, a critical step for IGF network activation. This analysis suggested that IGF2R levels would need to be 320-fold greater than IGF1R in order to decrease pIGF1R by 25 %, while IGFBP levels would need to be 390-fold greater. Analysis of The Cancer Genome Atlas (TCGA) data set suggested that this level of overexpression is unlikely for IGF2R in ovarian, breast, and colon cancer. In contrast, IGFBPs can likely reach these levels, suggesting that IGFBPs are the more critical regulator of IGF1R network activity. Levels of phosphorylated IGF1R were insensitive to changes in parameters regulating the IGF2R arm of the network.

Conclusions

Using a mass-action kinetic model, we determined that IGF2R plays a minor role in regulating the activity of IGF1R under a variety of conditions and that due to their high expression levels, IGFBPs are the dominant mechanism to regulating IGF network activation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0263-6) contains supplementary material, which is available to authorized users.

Effects of Tobacco Smoking on the Degeneration of the Intervertebral Disc: A Finite Element Study

Tobacco smoking is associated with numerous pathological conditions. Compelling experimental evidence associates smoking to the degeneration of the intervertebral disc (IVD). In particular, it has been shown that nicotine down-regulates both the proliferation rate and glycosaminoglycan (GAG) biosynthesis of disc cells. Moreover, tobacco smoking causes the constriction of the vascular network surrounding the IVD, thus reducing the exchange of nutrients and anabolic agents from the blood vessels to the disc. It has been hypothesized that both nicotine presence in the IVD and the reduced solute exchange are responsible for the degeneration of the disc due to tobacco smoking, but their effects on tissue homeostasis have never been quantified. In this study, a previously presented computational model describing the homeostasis of the IVD was deployed to investigate the effects of impaired solute supply and nicotine-mediated down-regulation of cell proliferation and biosynthetic activity on the health of the disc. We found that the nicotine-mediated down-regulation of cell anabolism mostly affected the GAG concentration at the cartilage endplate, reducing it up to 65% of the value attained in normal physiological conditions. In contrast, the reduction of solutes exchange between blood vessels and disc tissue mostly affected the nucleus pulposus, whose cell density and GAG levels were reduced up to 50% of their normal physiological levels. The effectiveness of quitting smoking on the regeneration of a degenerated IVD was also investigated, and showed to have limited benefit on the health of the disc. A cell-based therapy in conjunction with smoke cessation provided significant improvements in disc health, suggesting that, besides quitting smoking, additional treatments should be implemented in the attempt to recover the health of an IVD degenerated by tobacco smoking.

Predicting Knee Osteoarthritis

Treatment options for osteoarthritis (OA) beyond pain relief or total knee replacement are very limited. Because of this, attention has shifted to identifying which factors increase the risk of OA in vulnerable populations in order to be able to give recommendations to delay disease onset or to slow disease progression. The gold standard is then to use principles of risk management, first to provide subject-specific estimates of risk and then to find ways of reducing that risk. Population studies of OA risk based on statistical associations do not provide such individually tailored information. Here we argue that mechanistic models of cartilage tissue maintenance and damage coupled to statistical models incorporating model uncertainty, united within the framework of structural reliability analysis, provide an avenue for bridging the disciplines of epidemiology, cell biology, genetics and biomechanics. Such models promise subject-specific OA risk assessment and personalized strategies for mitigating or even avoiding OA. We illustrate the proposed approach with a simple model of cartilage extracellular matrix synthesis and loss regulated by daily physical activity.

Simulation predicts IGFBP2-HIF1α interaction drives glioblastoma growth

Tremendous strides have been made in improving patients’ survival from cancer with one glaring exception: brain cancer. Glioblastoma is the most common, aggressive and highly malignant type of primary brain tumor. The average overall survival remains less than 1 year. Notably, cancer patients with obesity and diabetes have worse outcomes and accelerated progression of glioblastoma. The root cause of this accelerated progression has been hypothesized to involve the insulin signaling pathway. However, while the process of invasive glioblastoma progression has been extensively studied macroscopically, it has not yet been well characterized with regards to intracellular insulin signaling. In this study we connect for the first time microscale insulin signaling activity with macroscale glioblastoma growth through the use of computational modeling. Results of the model suggest a novel observation: feedback from IGFBP2 to HIF1α is integral to the sustained growth of glioblastoma. Our study suggests that downstream signaling from IGFI to HIF1α, which has been the target of many insulin signaling drugs in clinical trials, plays a smaller role in overall tumor growth. These predictions strongly suggest redirecting the focus of glioma drug candidates on controlling the feedback between IGFBP2 and HIF1α.

Current treatment for glioblastoma patients is limited to nonspecific methods: surgery followed by a combination of radio- and chemotherapy. With these methods, glioma patient survival is less than one year post-diagnosis. Targeting specific protein signaling pathways offers potentially more potent therapies. One promising potential target is the insulin signaling pathway, which is known to contribute to glioblastoma progression. However, drugs targeting this pathway have shown mixed results in clinical trials, and the detailed mechanisms of how the insulin signaling pathway promotes glioblastoma growth remain to be elucidated. Here, we developed a computational model of insulin signaling in glioblastoma in order to study this pathway’s role in tumor progression. Using the model, we systematically test contributions of different insulin signaling protein interactions on glioblastoma growth. Our model highlights a key driver for the growth of glioblastoma: IGFBP2-HIF1α feedback. This interaction provides a target that could open the door for new therapies in glioma and other solid tumors.

A computational analysis on the implications of age-related changes in the expression of cellular signals on the role of IGF-1 in intervertebral disc homeostasis

Insulin-like growth factor-1 (IGF-1) is a well-known anabolic agent in intervertebral discs (IVD), promoting both proteoglycan (PG) biosynthesis and cell proliferation. Accordingly, it is believed that IGF-1 plays a central role in IVD homeostasis. The IGF-mediated anabolic activity in IVD occurs when the growth factor, free from binding proteins (IGFBP), binds to IGF cell surface receptors (IGF-1R). Previous studies reported that, with aging, cellular expression of IGFBP increases, while that of IGF-1R decreases. Both changes in cellular signals are thought to be among the factors that are responsible for the age-related decline in IGF-mediated PG biosynthesis, which ultimately leads to disc degeneration. In this study, a computational model describing the role of IGF-1 in the homeostasis of IVD was deployed in a parametric analysis to investigate the effects of age-related changes in expression of IGF-1R and IGFBP on the IGF-mediated upregulation of PG biosynthesis and cellular proliferation. It was found that changes in the expression of IGF-1R and IGFBP mostly affected the nucleus pulposus, while in the most external disc regions (annulus fibrosus and cartilage endplates) the IVD homeostatic balance was unaltered. It was shown that a decrease of IGF-1R expression caused reduction of both PG levels and cell density in the tissue. In contrast, increase in IGFBP expression increased both PG and cell concentration, suggesting that such change in cellular signaling may be a plausible defense mechanism from age-related IVD degeneration.

Analysis of the quantitative balance between insulin-like growth factor (IGF)-1 ligand, receptor, and binding protein levels to predict cell sensitivity and therapeutic efficacy

Background

The insulin-like growth factor (IGF) system impacts cell proliferation and is highly activated in ovarian cancer. While an attractive therapeutic target, the IGF system is complex with two receptors (IGF1R, IGF2R), two ligands (IGF1, IGF2), and at least six high affinity IGF-binding proteins (IGFBPs) that regulate the bioavailability of IGF ligands. We hypothesized that a quantitative balance between these different network components regulated cell response.

Results

OVCAR5, an immortalized ovarian cancer cell line, were found to be sensitive to IGF1, with the dose of IGF1 (i.e., the total mass of IGF1 available) a more reliable predictor of cell response than ligand concentration. The applied dose of IGF1 was depleted by both cell-secreted IGFBPs and endocytic trafficking, with IGFBPs sequestering up to 90% of the available ligand. To explore how different variables (i.e., IGF1, IGFBPs, and IGF1R levels) impacted cell response, a mass-action steady-state model was developed. Examination of the model revealed that the level of IGF1-IGF1R complexes per cell was directly proportional to the extent of proliferation induced by IGF1. Model analysis suggested, and experimental results confirmed, that IGFBPs present during IGF1 treatment significantly decreased IGF1-mediated proliferation. We utilized this model to assess the efficacy of IGF1 and IGF1R antibodies against different network compositions and determined that IGF1R antibodies were more globally effective due to the receptor-limited state of the network.

Conclusions

Changes that affect IGF1R occupancy have predictable effects on IGF1-induced proliferation and our model captured these effects. Analysis of this model suggests that IGF1R antibodies will be more effective than IGF1 antibodies, although the difference was minimal in conditions with low levels of IGF1 and IGFBPs. Examining how different components of the IGF system influence cell response will be critical to improve our understanding of the IGF signaling network in ovarian cancer.

Modeling the role of IGF-1 on extracellular matrix biosynthesis and cellularity in intervertebral disc

The insulin-like growth factor-1 (IGF-1) is a well-known anabolic agent for intervertebral disc (IVD), promoting both proteoglycan (PG) biosynthesis and cell proliferation. Accordingly, it is believed that IGF-1 may play a central role in IVD homeostasis. Furthermore, the exogenous administration of IGF-1 has been proposed as a possible therapeutic strategy for disc degeneration. The objectives of this study were to develop a new computational framework for describing the mechanisms regulating IGF-mediated homeostasis in IVD, and to apply this numerical tool for investigating the effectiveness of exogenous administration of IGF-1 for curing disc degeneration. A diffusive-reactive model was developed for describing competitive binding of IGF-1 to its binding proteins and cell surface receptors, with the latter reaction initiating the intracellular signaling mechanism leading to PG production and cell proliferation. Because PG production increases cell metabolic rate, and cell proliferation increases nutritional demand, nutrients transport and metabolism were also included into the model, and co-regulated, together with IGF-1, IVD cellularity. The sustainability and the effectiveness of IGF-mediated anabolism were investigated for conditions of pathologically insufficient nutrient supply, and for the case of exogenous administration of IGF-1 to degenerated IVD. Results showed that pathological nutrients deprivation, by decreasing cellularity, caused a reduction of PG biosynthesis. Also, exogenous administration of IGF-1 was only beneficial in well-nourished regions of IVD, and exacerbated cell mortality in malnourished regions. These findings remark the central role of nutrition in IVD health, and suggest that adequate nutritional supply is paramount for achieving a successful IGF-based therapy for disc degeneration.

Altered mechano-chemical environment in hip articular cartilage: effect of obesity

The production of extracellular matrix (ECM) components of articular cartilage is regulated, among other factors, by an intercellular signaling mechanism mediated by the interaction of cell surface receptors (CSR) with insulin-like growth factor-1 (IGF-1). In ECM, the presence of binding proteins (IGFBP) hinders IGF-1 delivery to CSR. It has been reported that levels of IGF-1 and IGFBP in obese population are, respectively, lower and higher than those found in normal population. In this study, an experimental-numerical approach was adopted to quantify the effect of this metabolic alteration found in obese population on the homeostasis of femoral hip cartilage. A new computational model, based on the mechano-electrochemical mixture theory, was developed to describe competitive binding kinetics of IGF-1 with IGFBP and CSR, and associated glycosaminoglycan (GAG) biosynthesis. Moreover, a gait analysis was carried out on obese and normal subjects to experimentally characterize mechanical loads on hip cartilage during walking. This information was deployed into the model to account for effects of physiologically relevant tissue deformation on GAG production in ECM. Numerical simulations were performed to compare GAG biosynthesis in femoral hip cartilage of normal and obese subjects. Results indicated that the lower ratio of IGF-1 to IGFBP found in obese population reduces cartilage GAG concentration up to 18 % when compared to normal population. Moreover, moderate physical activity, such as walking, has a modest beneficial effect on GAG production. The findings of this study suggest that IGF-1/IGFBP metabolic unbalance should be accounted for when considering the association of obesity with hip osteoarthritis.