Methods: We generated a computational model of growth plate chondrocytes (GP model) and a model of articular chondrocytes (AC model) In each model it is possible to set the activity of each of the proteins in the network at a level between 0 and 100 in increments of, for example, 1 or 10. It is also possible to randomly initialize the activities of proteins, or use a combination of random and set initializations. In silico experiments can then be carried out by simulating the evolution of the modelled regulatory network over time. We used ECHO to obtain insight into cross talk of 7 signal transduction pathways important for chondrocyte development and cartilage maintenance: IGF, PTHrP, BMP, FGF, TGFbeta, WNT, IHH. ECHO consists of a network of 123 intracellular signaling molecules (kinases) with 354 interactions. The network faithfully resembles current status of literature. RUNX2+ is the readout for hypertrophy and OA, and SOX9+ the readout of healthy articular chondrocytes. qPCR experiments on human mesenchymal stem cells (hMSCs) and human articular chondrocytes (hACs) were used to validate the model.
Results: Using ECHO we analyzed the most influential pathways, and performed in silico experiments to obtain insight into the molecular mechanisms of cartilage development and disease. For example, we investigated the cell’s response to addition or inhibition of single (Figure 1), or combinations of growth factors or cytokines. Both the GP and AC models showed a dose-dependent response to (combinations of) external stimuli. qPCR validated this dose-dependency in MSC and chondrocytes.We then investigated switching cell state through perturbations of extra-cellular ligands. In this, the models were started from a stable SOX9+ state and were switched to a RUNX2+ state, and vice versa. The in silico experiment predicted a transition from SOX9+ to RUNX2+ in the GP model with single addition of BMP or WNT, whereas in the AC model these factors caused no switch. The in silico experiments predicted a switch in the AC model with inhibition of IHH and simultaneous addition of high levels of either WNT or BMP (figure 2).In addition, we performed in silico double knockout and overexpression of all nodes in the network to predict which signals could interfere with OA development and progression. This resulted in hundreds of possible targets that are currently being investigated further. These experiments provide potential therapeutic clues for treatment of OA.
Conclusions: ECHO was used to mimic biological scenarios during chondrocyte development. We obtained insight into healthy chondrocyte development and OA development. Perturbation of single or double extracellular ligands cause a switch in cell fate, with BMP, WNT and IHH playing important roles. In addition, ECHO was used to identify potential therapeutic targets for OA treatment.
|Number of pages||1|
|Publication status||Published - Apr 2016|
|Event||2016 OARSI World Congress on Osteoarthrtis - Amsterdam RAI, Amsterdam, Netherlands|
Duration: 31 Mar 2016 → 3 Apr 2016
|Conference||2016 OARSI World Congress on Osteoarthrtis|
|Period||31/03/16 → 3/04/16|
- Signaling pathways
- dynamic modeling
- Timed Automata