Protein Adsorption Enhances Energy Dissipation in Networks of Lysozyme Amyloid Fibrils.

Hydrogels of amyloid fibrils are a versatile biomaterial for tissue engineering and other biomedical applications. Their suitability for these applications has been partly ascribed to their excellent and potentially engineerable rheological properties. However, while in biomedical applications the gels have to function in compositionally complex physiological solutions, their rheological behavior is typically only characterized in simple buffers. Here we show that the viscoelastic response of networks of amyloid fibrils of the protein lysozyme in biologically relevant solutions substantially differs from the response in simple buffers. We observe enhanced energy dissipation in both cell culture medium and synovial fluid. We attribute this energy dissipation to interactions of the amyloid fibrils with other molecules in these solutions and especially to the adsorption of the abundantly present protein serum albumin. This finding provides the basis for a better understanding of the performance of amyloid hydrogels in biomedical applications.

The Effects of the WNT-Signaling Modulators BIO and PKF118-310 on the Chondrogenic Differentiation of Human Mesenchymal Stem Cells.

Mesenchymal stem cells (MSCs) are multipotent cells, mainly from bone marrow, and an ideal source of cells in bone and cartilage tissue engineering. A study of the chondrogenic differentiation of MSCs is of particular interest for MSCs-based cartilage regeneration. In this study, we aimed to optimize the conditions for the chrondogenic differentiation of MSCs by regulating WNT signaling using the small molecule WNT inhibitor PKF118-310 and activator BIO. Human mesenchymal stem cells (hMSCs) were isolated from bone marrow aspirates and cultured in hMSCs proliferation medium. Pellet culture was subsequently established for three-dimensional chondrogenic differentiation of 5 weeks. WNT signaling was increased by the small molecule glycogen synthase kinase-3 inhibitor 6-bromoindirubin-3-oxim (BIO) and decreased by the WNT inhibitor PKF118-310 (PKF). The effects of BIO and PKF on the chondrogenesis of hMSCs was examined by real-time PCR, histological methods, and ELISA. We found that activation of canonical WNT-signaling by BIO significantly downregulated the expression of cartilage-specific genes SOX9, COL2A1, and ACAN, and matrix metalloproteinase genes MMP1/3/9/13, but increased ADAMTS 4/5. Inhibition of WNT signaling by PKF increased the expression of SOX9, COL2A1, ACAN, and MMP9, but decreased MMP13 and ADAMTS4/5. In addition, a high level of WNT signaling induced the expression of hypertrophic markers COL10A1, ALPL, and RUNX2, the dedifferentiation marker COL1A1, and glycolysis genes GULT1 and PGK1. Deposition of glycosaminoglycan (GAG) and collagen type II in the pellet matrix was significantly lost in the BIO-treated group and increased in the PKF-treated group. The protein level of COL10A1 was also highly induced in the BIO group. Interestingly, BIO decreased the number of apoptotic cells while PKF significantly induced apoptosis during chondrogenesis. The natural WNT antagonist DKK1 and the protein level of MMP1 in the pellet culture medium were decreased after PKF treatment. All of these chondrogenic effects appeared to be mediated through the canonical WNT signaling pathway, since the target gene Axin2 and other WNT members, such as TCF4 and ß-catenin, were upregulated by BIO and downregulated by PKF, respectively, and BIO induced nuclear translocation of ß-catenin while PKF inhibited ß-catenin translocation into the nucleus. We concluded that addition of BIO to a chondrogenic medium of hMSCs resulted in a loss of cartilage formation, while PKF induced chondrogenic differentiation and cartilage matrix deposition and inhibited hypertrophic differentiation. However, BIO promoted cell survival by inhibiting apoptosis while PKF induced cell apoptosis. This result indicates that either an overexpression or overinhibition of WNT signaling to some extent causes harmful effects on chondrogenic differentiation. Cartilage tissue engineering could benefit from the adjustment of the critical level of WNT signaling during chondrogenesis of hMSC.

Oxygen-Dependent Lipid Profiles of Three-Dimensional Cultured Human Chondrocytes Revealed by MALDI-MSI.

Articular cartilage is exposed to a gradient of oxygen levels ranging from 5% at the surface to 1% in the deepest layers. While most cartilage research is performed in supraphysiological oxygen levels (19-21%), culturing chondrocytes under hypoxic oxygen levels (≤8%) promotes the chondrogenic phenotype. Exposure of cells to various oxygen levels alters their lipid metabolism, but detailed studies examining how hypoxia affects lipid metabolism in chondrocytes are lacking. To better understand the chondrocyte's behavior in response to oxygen, we cultured 3D pellets of human primary chondrocytes in normoxia (20% oxygen) and hypoxia (2.5% oxygen) and employed matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI) in order to characterize the lipid profiles and their spatial distribution. In this work we show that chondrocytes cultured in hypoxia and normoxia can be differentiated by their lipid profiles. Among other species, phosphatidylglycerol species were increased in normoxic pellets, whereas phosphatidylinositol species were the most prominent lipids in hypoxic pellets. Moreover, spatial mapping revealed that phospahtidylglyycerol species were less prominent in the center of pellets where the oxygen level is lower. Additional analysis revealed a higher abundance of the mitochondrial-specific lipids, cardiolipins, in normoxic conditions. In conclusion MALDI-MSI described specific lipid profiles that could be used as sensors of oxygen level changes and may especially be relevant for retaining the chondrogenic phenotype, which has important implications for the treatment of bone and cartilage diseases.

Metabolic programming of mesenchymal stromal cells by oxygen tension directs chondrogenic cell fate.

Actively steering the chondrogenic differentiation of mesenchymal stromal cells (MSCs) into either permanent cartilage or hypertrophic cartilage destined to be replaced by bone has not yet been possible. During limb development, the developing long bone is exposed to a concentration gradient of oxygen, with lower oxygen tension in the region destined to become articular cartilage and higher oxygen tension in transient hypertrophic cartilage. Here, we prove that metabolic programming of MSCs by oxygen tension directs chondrogenesis into either permanent or transient hyaline cartilage. Human MSCs chondrogenically differentiated in vitro under hypoxia (2.5% O2) produced more hyaline cartilage, which expressed typical articular cartilage biomarkers, including established inhibitors of hypertrophic differentiation. In contrast, normoxia (21% O2) prevented the expression of these inhibitors and was associated with increased hypertrophic differentiation. Interestingly, gene network analysis revealed that oxygen tension resulted in metabolic programming of the MSCs directing chondrogenesis into articular- or epiphyseal cartilage-like tissue. This differentiation program resembled the embryological development of these distinct types of hyaline cartilage. Remarkably, the distinct cartilage phenotypes were preserved upon implantation in mice. Hypoxia-preconditioned implants remained cartilaginous, whereas normoxia-preconditioned implants readily underwent calcification, vascular invasion, and subsequent endochondral ossification. In conclusion, metabolic programming of MSCs by oxygen tension provides a simple yet effective mechanism by which to direct the chondrogenic differentiation program into either permanent articular-like cartilage or hypertrophic cartilage that is destined to become endochondral bone.

Nitric Oxide Mediates Crosstalk between Interleukin 1ß and WNT Signaling in Primary Human Chondrocytes by Reducing DKK1 and FRZB Expression.

Interleukin 1 beta (IL1ß) and Wingless-Type MMTV Integration Site Family (WNT) signaling are major players in Osteoarthritis (OA) pathogenesis. Despite having a large functional overlap in OA onset and development, the mechanism of IL1ß and WNT crosstalk has remained largely unknown. In this study, we have used a combination of computational modeling and molecular biology to reveal direct or indirect crosstalk between these pathways. Specifically, we revealed a mechanism by which IL1ß upregulates WNT signaling via downregulating WNT antagonists, DKK1 and FRZB. In human chondrocytes, IL1ß decreased the expression of Dickkopf-1 (DKK1) and Frizzled related protein (FRZB) through upregulation of nitric oxide synthase (iNOS), thereby activating the transcription of WNT target genes. This effect could be reversed by iNOS inhibitor 1400W, which restored DKK1 and FRZB expression and their inhibitory effect on WNT signaling. In addition, 1400W also inhibited both the matrix metalloproteinase (MMP) expression and cytokine-induced apoptosis. We concluded that iNOS/NO play a pivotal role in the inflammatory response of human OA through indirect upregulation of WNT signaling. Blocking NO production may inhibit the loss of the articular phenotype in OA by preventing downregulation of the expression of DKK1 and FRZB.

Correlation between Gene Expression and Osteoarthritis Progression in Human.

Osteoarthritis (OA) is a multifactorial disease characterized by gradual degradation of joint cartilage. This study aimed to quantify major pathogenetic factors during OA progression in human cartilage. Cartilage specimens were isolated from OA patients and scored 0-5 according to the Osteoarthritis Research Society International (OARSI) guidelines. Protein and gene expressions were measured by immunohistochemistry and qPCR, respectively. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assays were used to detect apoptotic cells. Cartilage degeneration in OA is a gradual progress accompanied with gradual loss of collagen type II and a gradual decrease in mRNA expression of SOX9, ACAN and COL2A1. Expression of WNT antagonists DKK1 and FRZB was lost, while hypertrophic markers (RUNX2, COL10A1 and IHH) increased during OA progression. Moreover, DKK1 and FRZB negatively correlated with OA grading, while RUNX2 and IHH showed a significantly positive correlation with OA grading. The number of apoptotic cells was increased with the severity of OA. Taken together, our results suggested that genetic profiling of the gene expression could be used as markers for staging OA at the molecular level. This helps to understand the molecular pathology of OA and may lead to the development of therapies based on OA stage.

Differentiation of mesenchymal stem cells under hypoxia and normoxia: lipid profiles revealed by time-of-flight secondary ion mass spectrometry and multivariate analysis.

Mesenchymal stem cells (MSC) have the ability to self-renew and differentiate into multiple cell types valuable for clinical treatment of rheumatic pathologies. To study the chondrogenic potential of MSC and identify the conditions that recreate the native cartilage environment, we used time-of-flight secondary ion mass spectrometry (TOF-SIMS) for label-free detection of cell-type- and environmental-condition-specific molecular profiles. We observed that coculture of human MSC and chondrocytes under standard culture conditions leads to improved extracellular matrix (ECM) deposition. In marked contrast, this effect was lost under low oxygen tension. This improved extracellular matrix deposition was associated with a significant decrease in lipids and in particular cholesterol under low oxygen tension as revealed by TOF-SIMS coupled to principal component analysis and discriminant analysis. We furthermore demonstrate that the higher cholesterol levels under normoxia might regulate fibroblast growth factor 1 (FGF-1) gene expression which was previously implemented in increased ECM production in the cocultures. In conclusion, our study shows an unexpected role of lipids as orchestrators of chondrogenesis in response to oxygen tension which is, at least in part, mediated through FGF-1.

The Regulatory Role of Signaling Crosstalk in Hypertrophy of MSCs and Human Articular Chondrocytes.

Hypertrophic differentiation of chondrocytes is a main barrier in application of mesenchymal stem cells (MSCs) for cartilage repair. In addition, hypertrophy occurs occasionally in osteoarthritis (OA). Here we provide a comprehensive review on recent literature describing signal pathways in the hypertrophy of MSCs-derived in vitro differentiated chondrocytes and chondrocytes, with an emphasis on the crosstalk between these pathways. Insight into the exact regulation of hypertrophy by the signaling network is necessary for the efficient application of MSCs for articular cartilage repair and for developing novel strategies for curing OA. We focus on articles describing the role of the main signaling pathways in regulating chondrocyte hypertrophy-like changes. Most studies report hypertrophic differentiation in chondrogenesis of MSCs, in both human OA and experimental OA. Chondrocyte hypertrophy is not under the strict control of a single pathway but appears to be regulated by an intricately regulated network of multiple signaling pathways, such as WNT, Bone morphogenetic protein (BMP)/Transforming growth factor-ß (TGFß), Parathyroid hormone-related peptide (PTHrP), Indian hedgehog (IHH), Fibroblast growth factor (FGF), Insulin like growth factor (IGF) and Hypoxia-inducible factor (HIF). This comprehensive review describes how this intricate signaling network influences tissue-engineering applications of MSCs in articular cartilage (AC) repair, and improves understanding of the disease stages and cellular responses within an OA articular joint.

T cell factor 4 is a pro-catabolic and apoptotic factor in human articular chondrocytes by potentiating nuclear factor κB signaling.

T cell factor/lymphoid enhancer factor (TCF/LEF) transcription factors are downstream effectors of Wnt/ß-catenin signaling, which has been implicated in the development and progression of osteoarthritis (OA). This study aimed to investigate the role of TCF/LEF transcription factors in human articular chondrocytes. Primary human osteoarthritic cartilage predominantly expressed TCF4 and to a lesser extent, LEF1 and TCF3 mRNA. Overexpression of TCF4, but not of TCF3 or LEF1, induced MMP-1, -3, and -13 expression and generic MMP activity in human chondrocytes. This was due to potentiating NF-κB signaling by a protein-protein interaction between TCF4 and NF-κB p65 activating established NF-κB target genes such as MMPs and IL-6. LEF1 competed with TCF4 for binding to NF-κB p65. IκB-α was able to counteract the effect of TCF4 on NF-κB target gene expression. Finally, we showed that TCF4 mRNA expression was elevated in OA cartilage compared with healthy cartilage and induced chondrocyte apoptosis at least partly through activating caspase 3/7. Our findings suggest that increased TCF4 expression may contribute to cartilage degeneration in OA by augmenting NF-κB signaling.

Osmolarity-Induced Altered Intracellular Molecular Crowding Drives Osteoarthritis Pathology.

Osteoarthritis (OA) is a multifactorial degenerative joint disease of which the underlying mechanisms are yet to be fully understood. At the molecular level, multiple factors including altered signaling pathways, epigenetics, metabolic imbalance, extracellular matrix degradation, production of matrix metalloproteinases, and inflammatory cytokines, are known to play a detrimental role in OA. However, these factors do not initiate OA, but are mediators or consequences of the disease, while many other factors causing the etiology of OA are still unknown. Here, it is revealed that microenvironmental osmolarity can induce and reverse osteoarthritis-related behavior of chondrocytes via altered intracellular molecular crowding, which represents a previously unknown mechanism underlying OA pathophysiology. Decreased intracellular crowding is associated with increased sensitivity to proinflammatory triggers and decreased responsiveness to anabolic stimuli. OA-induced lowered intracellular molecular crowding could be renormalized via exposure to higher extracellular osmolarity such as those found in healthy joints, which reverse OA chondrocyte's sensitivity to catabolic stimuli as well as its glycolytic metabolism.

WNT signaling and cartilage: of mice and men.

In adult articular cartilage, the extracellular matrix is maintained by a balance between the degradation and the synthesis of matrix components. Chondrocytes that sparsely reside in the matrix and rarely proliferate are the key cellular mediators for cartilage homeostasis. There are indications for the involvement of the WNT signaling pathway in maintaining articular cartilage. Various WNTs are involved in the subsequent stages of chondrocyte differentiation during development, and deregulation of WNT signaling was observed in cartilage degeneration. Even though gene expression and protein synthesis can be activated upon injury, articular cartilage has a limited ability of self-repair and efforts to regenerate articular cartilage have so far not been successful. Because WNT signaling was found to be involved in the development and maintenance of cartilage as well as in the degeneration of cartilage, interfering with this pathway might contribute to improving cartilage regeneration. However, most of the studies on elucidating the role of WNT signaling in these processes were conducted using in vitro or in vivo animal models. Discrepancies have been found in the role of WNT signaling between chondrocytes of mouse and human origin, and extrapolation of results from mouse models to the human situation remains a challenge. Elucidation of detailed WNT signaling functions will provide knowledge to improve cartilage regeneration.

Choice of Protein, Not Its Amyloid-Fold, Determines the Success of Amyloid-Based Scaffolds for Cartilage Tissue Regeneration.

The formation of fibrocartilage during articular cartilage regeneration remains a clinical problem affecting adequate restoration of articular cartilage in joints. To stimulate chondrocytes to form articular cartilage, we investigated the use of amyloid fibril-based scaffolds. The proteins α-synuclein, ß-lactoglobulin, and lysozyme were induced to self-assemble into amyloid fibrils and, during dialysis, formed micrometer scale amyloid networks that resemble the cartilage extracellular matrix. Our results show that lysozyme amyloid micronetworks supported chondrocyte viability and extracellular matrix deposition, while α-synuclein and ß-lactoglobulin maintained cell viability. With this study, we not only confirm the possible use of amyloid materials for tissue regeneration but also demonstrate that the choice of protein, rather than its amyloid-fold per se, affects the cellular response and tissue formation.

Implementing Computational Modeling in Tissue Engineering: Where Disciplines Meet.

In recent years, the mathematical and computational sciences have developed novel methodologies and insights that can aid in designing advanced bioreactors, microfluidic setups or organ-on-chip devices, in optimizing culture conditions, or predicting long-term behavior of engineered tissues in vivo. In this review, we introduce the concept of computational models and how they can be integrated in an interdisciplinary workflow for Tissue Engineering and Regenerative Medicine (TERM). We specifically aim this review of general concepts and examples at experimental scientists with little or no computational modeling experience. We also describe the contribution of computational models in understanding TERM processes and in advancing the TERM field by providing novel insights. Impact Statement Although in recent years the use of mathematical and computational sciences has increased in the Tissue Engineering and Regenerative Medicine (TERM) field, we believe that a further integration of experimental and computational approaches has a huge potential for advancing the field due to the ability of models to explain and predict experimental results and efficiently optimize TERM product and process designs. By providing an overview of existing computational models, how they have contributed to the field, as well as a future perspective, this review represents an important step to help realize TERM's ultimate goal: a cure instead of care.

ERK nuclear translocation is dimerization-independent but controlled by the rate of phosphorylation.

Upon activation, ERKs translocate from the cytoplasm to the nucleus. This process is required for the induction of many cellular responses, yet the molecular mechanisms that regulate ERK nuclear translocation are not fully understood. We have used a mouse embryo fibroblast ERK1-knock-out cell line expressing green fluorescent protein (GFP)-tagged ERK1 to probe the spatio-temporal regulation of ERK1. Real time fluorescence microscopy and fluorescence correlation spectroscopy revealed that ERK1 nuclear accumulation increased upon serum stimulation, but the mobility of the protein in the nucleus and cytoplasm remained unchanged. Dimerization of ERK has been proposed as a requirement for nuclear translocation. However, ERK1-Delta4, the mutant shown consistently to be dimerization-deficient in vitro, accumulated in the nucleus to the same level as wild type (WT), indicating that dimerization of ERK1 is not required for nuclear entry and retention. Consistent with this finding, energy migration Förster resonance energy transfer and fluorescence correlation spectroscopy measurements in living cells did not detect dimerization of GFP-ERK1-WT upon activation. In contrast, the kinetics of nuclear accumulation and phosphorylation of GFP-ERK1-Delta4 were slower than that of GFP-ERK1-WT. These results indicate that the differential shuttling behavior of the mutant is a consequence of delayed phosphorylation of ERK by MEK rather than dimerization. Our data demonstrate for the first time that a delay in cytoplasmic activation of ERK is directly translated into a delay in nuclear translocation.

Using FRAP to Quantify Changes in Transcription Factor Dynamics After Cell Stimulation: Cell Culture, FRAP, Data Analysis, and Visualization.

Here we show how to measure the mobility of transcription factors using fluorescence recovery after photobleaching (FRAP). Transcription factors are DNA-binding proteins that, upon binding to specific DNA motifs, regulate transcription of their target genes. FRAP is a simple, fast, and cost-effective method, and is a widely used quantitative method to measure the dynamics of fluorescently labeled molecules in solution, membranes, and inside living cells. Dynamics, specified by the immobile fraction, recovery half-time, diffusion constant, and ratio of molecules contributing to different phases of FRAP recovery, can be quantified by FRAP. This can be useful to understand their function in gene regulation. This tutorial is intended to familiarize the reader with the FRAP procedure to quantify transcription factor dynamics using a standard confocal microscope and analysis using MATLAB (MathWorks®). This article will guide the reader through the preconditions of FRAP, and a detailed and step-by-step procedure of preparing cells, bleaching protocol, data analysis in MATLAB, and visualization of the FRAP data.

Quantitative Molecular Models for Biological Processes: Modeling of Signal Transduction Networks with ANIMO.

Computational modeling of biological networks is increasing in popularity due to the increased demand for understanding biological processes. This understanding requires integration of a variety of experimental data that allows understanding of complex mechanisms regulating cell and tissue function. However, the mathematical complexity of many modeling tools have thusfar prevented broad adaptation and effective use by molecular biologists. In this chapter, we show by example how one can start building a model in ANIMO and how to adapt the model to experimental data. We show how this model can be used for simulating network activities, testing hypotheses, and how to improve the model using wet-lab data.

The Expressions of Dickkopf-Related Protein 1 and Frizzled-Related Protein Are Negatively Correlated to Local Inflammation and Osteoarthritis Severity.

OBJECTIVE: To investigate the presence of WNT antagonists Dickkopf-related protein 1 (DKK1), Frizzled-related protein (FRZB) and BMP antagonist Gremlin 1 (GREM1) in synovial fluid (SF) and serum, respectively, from end-stage knee osteoarthritis (OA) patients, and correlate their expression with other markers of OA. DESIGN: In a cross-sectional study, SF and serum were collected from OA patients (n = 132). The concentrations of DKK1, FRZB and GREM1 in SF and serum were determined using immunoassays. Correlation measurements were performed between groups and previously assessed disease markers, such as synovium nitric oxide (NO), inerleukin-1ß (IL1ß), tumor necrosis factor-α (TNFα), and prostaglandin E2 (PGE2). RESULTS: The OA patients with the celecoxib treatment till surgery have higher median SF FRZB values compared with the control (no treatment); the celecoxib 3-days before surgery stopped treatment group has higher median serum FRZB values than the control and the naproxen treatment group. The combinational analysis of SF DKK1 and SF FRZB negatively correlated with macroscopic cartilage scores and histological synovium scores in OA patients. The expression of DKK1 and FRZB in SF showed the same expression trend as their expression in serum. Furthermore, the SF concentration of DKK1 was positively correlated with FRZB in both SF and serum. In contrast, it was negatively correlated with synovium NO and IL1ß. SF FRZB was negatively correlated with synovium NO, IL1ß, cartilage PGE2, and age. CONCLUSIONS: Our findings suggest DKK1 and FRZB were negatively correlated with OA severity and multiple pro-inflammatory cytokines. Our data indicate that DKK1 and FRZB can be joint disease-specific biomarkers.

Engineering Cartilage Tissue by Co-culturing of Chondrocytes and Mesenchymal Stromal Cells.

Co-culture of chondrocytes and mesenchymal stromal cells (MSCs) has been shown to be beneficial in engineering cartilage tissue in vitro. In these co-cultures, MSCs increase the proliferation and matrix deposition of chondrocytes. The MSCs accomplish this beneficial effect by so-called trophic actions. Thus, large cartilage constructs can be made with a relatively small number of chondrocytes. In this chapter, we describe different methods for making co-cultures of MSCs and chondrocytes. We also provide detailed protocols for analyzing MSC-chondrocyte co-cultures with cell tracking, proliferation assays, species-specific polymerase chain reactions (PCR), rheological analysis, compression analysis, RNA-sequencing analysis, short tandem repeats analysis, and biochemical examination.

An ECHO of Cartilage: In Silico Prediction of Combinatorial Treatments to Switch Between Transient and Permanent Cartilage Phenotypes With Ex Vivo Validation.

A fundamental question in cartilage biology is: what determines the switch between permanent cartilage found in the articular joints and transient hypertrophic cartilage that functions as a template for bone? This switch is observed both in a subset of OA patients that develop osteophytes, as well as in cell-based tissue engineering strategies for joint repair. A thorough understanding of the mechanisms regulating cell fate provides opportunities for treatment of cartilage disease and tissue engineering strategies. The objective of this study was to understand the mechanisms that regulate the switch between permanent and transient cartilage using a computational model of chondrocytes, ECHO. To investigate large signaling networks that regulate cell fate decisions, we developed the software tool ANIMO, Analysis of Networks with interactive Modeling. In ANIMO, we generated an activity network integrating 7 signal transduction pathways resulting in a network containing over 50 proteins with 200 interactions. We called this model ECHO, for executable chondrocyte. Previously, we showed that ECHO could be used to characterize mechanisms of cell fate decisions. ECHO was first developed based on a Boolean model of growth plate. Here, we show how the growth plate Boolean model was translated to ANIMO and how we adapted the topology and parameters to generate an articular cartilage model. In ANIMO, many combinations of overactivation/knockout were tested that result in a switch between permanent cartilage (SOX9+) and transient, hypertrophic cartilage (RUNX2+). We used model checking to prioritize combination treatments for wet-lab validation. Three combinatorial treatments were chosen and tested on metatarsals from 1-day old rat pups that were treated for 6 days. We found that a combination of IGF1 with inhibition of ERK1/2 had a positive effect on cartilage formation and growth, whereas activation of DLX5 combined with inhibition of PKA had a negative effect on cartilage formation and growth and resulted in increased cartilage hypertrophy. We show that our model describes cartilage formation, and that model checking can aid in choosing and prioritizing combinatorial treatments that interfere with normal cartilage development. Here we show that combinatorial treatments induce changes in the zonal distribution of cartilage, indication possible switches in cell fate. This indicates that simulations in ECHO aid in describing pathologies in which switches between cell fates are observed, such as OA.

ECHO, the executable CHOndrocyte: A computational model to study articular chondrocytes in health and disease.

Computational modeling can be used to investigate complex signaling networks in biology. However, most modeling tools are not suitable for molecular cell biologists with little background in mathematics. We have built a visual-based modeling tool for the investigation of dynamic networks. Here, we describe the development of computational models of cartilage development and osteoarthritis, in which a panel of relevant signaling pathways are integrated. In silico experiments give insight in the role of each of the pathway components and reveal which perturbations may deregulate the basal healthy state of cells and tissues. We used a previously developed computational modeling tool Analysis of Networks with Interactive Modeling (ANIMO) to generate an activity network integrating 7 signal transduction pathways resulting in a network containing over 50 nodes and 200 interactions. We performed in silico experiments to characterize molecular mechanisms of cell fate decisions. The model was used to mimic biological scenarios during cell differentiation using RNA-sequencing data of a variety of stem cell sources as input. In a case-study, we wet-lab-tested the model-derived hypothesis that expression of DKK1 (Dickkopf-1) and FRZB (Frizzled related protein, WNT antagonists) and GREM1 (gremlin 1, BMP antagonist) prevents IL1ß (Interleukin 1 beta)-induced MMP (matrix metalloproteinase) expression, thereby preventing cartilage degeneration, at least in the short term. We found that a combination of DKK1, FRZB and GREM1 may play a role in modulating the effects of IL1ß induced inflammation in human primary chondrocytes.