Modeling of signaling pathways in chondrocytes based on phosphoproteomic and cytokine release data

Runx2 deletion in hypertrophic chondrocytes impairs osteoclast mediated bone resorption

Deletion of Runx2 gene in proliferating chondrocytes results in complete failure of endochondral ossification and perinatal lethality. We reported recently that mice with Runx2 deletion specifically in hypertrophic chondrocytes (HCs) using the Col10a1-Cre transgene survive and exhibit enlarged growth plates due to decreased HC apoptosis and cartilage resorption. Bulk of chondrogenesis occurs postnatally, however, the role of Runx2 in HCs during postnatal chondrogenesis is unknown. Despite limb dwarfism, adult homozygous (Runx2HC/HC) mice showed a significant increase in length of growth plate and articular cartilage. Consistent with doubling of the hypertrophic zone, collagen type X expression was increased in Runx2HC/HC mice. In sharp contrast, expression of metalloproteinases and aggrecanases were markedly decreased. Impaired cartilage degradation was evident by the retention of significant amount of safranin-O positive cartilage. Histomorphometry and μCT uncovered increased trabecular bone mass with a significant increase in BV/TV ratio, trabecular number, thickness, and a decrease in trabecular space in Runx2HC/HC mice. To identify if this is due to increased bone synthesis, expression of osteoblast differentiation markers was evaluated and found to be comparable amongst littermates. Histomorphometry confirmed similar number of osteoblasts in the littermates. Furthermore, dynamic bone synthesis showed no differences in mineral apposition or bone formation rates. Surprisingly, three-point-bending test revealed Runx2HC/HC bones to be structurally less strong. Interestingly, both the number and surface of osteoclasts were markedly reduced in Runx2HC/HC littermates. Rankl and IL-17a ligands that promote osteoclast differentiation were markedly reduced in Runx2HC/HC mice. Bone marrow cultures were performed to independently establish Runx2 and hypertrophic chondrocytes role in osteoclast development. The culture from the Runx2HC/HC mice formed significantly fewer and smaller osteoclasts. The expression of mature osteoclast markers, Ctsk and Mmp9, were significantly reduced in the cultures from Runx2HC/HC mice. Thus, Runx2 functions extend beyond embryonic development and chondrocyte hypertrophy by regulating cartilage degradation, osteoclast differentiation, and bone resorption during postnatal endochondral ossification.

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 set-ups 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.

Regulatory network-based model to simulate the biochemical regulation of chondrocytes in healthy and osteoarthritic environments

In osteoarthritis (OA), chondrocyte metabolism dysregulation increases relative catabolic activity, which leads to cartilage degradation. To enable the semiquantitative interpretation of the intricate mechanisms of OA progression, we propose a network-based model at the chondrocyte level that incorporates the complex ways in which inflammatory factors affect structural protein and protease expression and nociceptive signals. Understanding such interactions will leverage the identification of new potential therapeutic targets that could improve current pharmacological treatments. Our computational model arises from a combination of knowledge-based and data-driven approaches that includes in-depth analyses of evidence reported in the specialized literature and targeted network enrichment. We achieved a mechanistic network of molecular interactions that represent both biosynthetic, inflammatory and degradative chondrocyte activity. The network is calibrated against experimental data through a genetic algorithm, and 81% of the responses tested have a normalized root squared error lower than 0.15. The model captures chondrocyte-reported behaviors with 95% accuracy, and it correctly predicts the main outcomes of OA treatment based on blood-derived biologics. The proposed methodology allows us to model an optimal regulatory network that controls chondrocyte metabolism based on measurable soluble molecules. Further research should target the incorporation of mechanical signals.

Evidence-Based Network Modelling to Simulate Nucleus Pulposus Multicellular Activity in Different Nutritional and Pro-Inflammatory Environments

Initiation of intervertebral disc degeneration is thought to be biologically driven. This reflects a process, where biochemical and mechanical stimuli affect cell activity (CA) that compromise the tissue strength over time. Experimental research enhanced our understanding about the effect of such stimuli on different CA, such as protein synthesis or mRNA expression. However, it is still unclear how cells respond to their native environment that consists of a “cocktail” of different stimuli that might locally vary. This work presents an interdisciplinary approach of experimental and in silico research to approximate Nucleus Pulposus CA within multifactorial biochemical environments. Thereby, the biochemical key stimuli glucose, pH, and the proinflammatory cytokines TNF-α and IL1β were considered that were experimentally shown to critically affect CA. To this end, a Nucleus Pulposus multicellular system was modelled. It integrated experimental findings from in vitro studies of human or bovine Nucleus Pulposus cells, to relate the individual effects of targeted stimuli to alterations in CA. Unknown stimulus-CA relationships were obtained through own experimental 3D cultures of bovine Nucleus Pulposus cells in alginate beads. Translation of experimental findings into suitable parameters for network modelling approaches was achieved thanks to a new numerical approach to estimate the individual sensitivity of a CA to each stimulus type. Hence, the effect of each stimulus type on a specific CA was assessed and integrated to approximate a multifactorial stimulus environment. Tackled CA were the mRNA expressions of Aggrecan, Collagen types I & II, MMP3, and ADAMTS4. CA was assessed for four different proinflammatory cell states; non-inflamed and inflamed for IL1β, TNF-α or both IL1β&TNF-α. Inflamed cell clusters were eventually predicted in a multicellular 3D agent-based model. Experimental results showed that glucose had no significant impact on proinflammatory cytokine or ADAMTS4 mRNA expression, whereas TNF-α caused a significant catabolic shift in most explored CA. In silico results showed that the presented methodology to estimate the sensitivity of a CA to a stimulus type importantly improved qualitative model predictions. However, more stimuli and/or further experimental knowledge need to be integrated, especially regarding predictions about the possible progression of inflammatory environments under adverse nutritional conditions. Tackling the multicellular level is a new and promising approach to estimate manifold responses of intervertebral disc cells. Such a top-down high-level network modelling approach allows to obtain information about relevant stimulus environments for a specific CA and could be shown to be suitable to tackle complex biological systems, including different proinflammatory cell states. The development of this methodology required a close interaction with experimental research. Thereby, specific experimental needs were derived from systematic in silico approaches and obtained results were directly used to enhance model predictions, which reflects a novelty in this research field. Eventually, the presented methodology provides modelling solutions suitable for multiscale approaches to contribute to a better understanding about dynamics over multiple spatial scales. Future work should focus on an amplification of the stimulus environment by integrating more key relevant stimuli, such as mechanical loading parameters, in order to better approximate native physiological environments.

A dysfunctional TRPV4-GSK3β pathway prevents osteoarthritic chondrocytes from sensing changes in extracellular matrix viscoelasticity

Changes in the composition and viscoelasticity of the extracellular matrix in load-bearing cartilage influence the proliferation and phenotypes of chondrocytes, and are associated with osteoarthritis. However, the underlying molecular mechanism is unknown. Here we show that the viscoelasticity of alginate hydrogels regulates cellular volume in healthy human chondrocytes (with faster stress relaxation allowing cell expansion and slower stress relaxation restricting it) but not in osteoarthritic chondrocytes. Cellular volume regulation in healthy chondrocytes was associated with changes in anabolic gene expression, in the secretion of multiple pro-inflammatory cytokines, and in the modulation of intracellular calcium regulated by the ion-channel protein transient receptor potential cation channel subfamily V member 4 (TRPV4), which controls the phosphorylation of glycogen synthase kinase 3β (GSK3β), an enzyme with pleiotropic effects in osteoarthritis. A dysfunctional TRPV4-GSK3β pathway in osteoarthritic chondrocytes rendered the cells unable to respond to environmental changes in viscoelasticity. Our findings suggest strategies for restoring chondrocyte homeostasis in osteoarthritis.

Simulating intervertebral disc cell behaviour within 3D multifactorial environments

MOTIVATION

Low back pain is responsible for more global disability than any other condition. Its incidence is closely related to intervertebral disc (IVD) failure, which is likely caused by an accumulation of microtrauma within the IVD. Crucial factors in microtrauma development are not entirely known yet, probably because their exploration in vivo or in vitro remains tremendously challenging. In silico modelling is, therefore, definitively appealing, and shall include approaches to integrate influences of multiple cell stimuli at the microscale. Accordingly, this study introduces a hybrid Agent-based (AB) model in IVD research and exploits network modelling solutions in systems biology to mimic the cellular behaviour of Nucleus Pulposus cells exposed to a 3D multifactorial biochemical environment, based on mathematical integrations of existing experimental knowledge. Cellular activity reflected by mRNA expression of Aggrecan, Collagen type I, Collagen type II, MMP-3 and ADAMTS were calculated for inflamed and non-inflamed cells. mRNA expression over long periods of time is additionally determined including cell viability estimations. Model predictions were eventually validated with independent experimental data.

RESULTS

As it combines experimental data to simulate cell behaviour exposed to a multifactorial environment, the present methodology was able to reproduce cell death within 3 days under glucose deprivation and a 50% decrease in cell viability after 7 days in an acidic environment. Cellular mRNA expression under non-inflamed conditions simulated a quantifiable catabolic shift under an adverse cell environment, and model predictions of mRNA expression of inflamed cells provide new explanation possibilities for unexpected results achieved in experimental research.

AVAILABILITYAND IMPLEMENTATION

The AB model as well as used mathematical functions were built with open source software. Final functions implemented in the AB model and complete AB model parameters are provided as Supplementary Material. Experimental input and validation data were provided through referenced, published papers. The code corresponding to the model can be shared upon request and shall be reused after proper training.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DIPPER, a spatiotemporal proteomics atlas of human intervertebral discs for exploring ageing and degeneration dynamics

The spatiotemporal proteome of the intervertebral disc (IVD) underpins its integrity and function. We present DIPPER, a deep and comprehensive IVD proteomic resource comprising 94 genome-wide profiles from 17 individuals. To begin with, protein modules defining key directional trends spanning the lateral and anteroposterior axes were derived from high-resolution spatial proteomes of intact young cadaveric lumbar IVDs. They revealed novel region-specific profiles of regulatory activities and displayed potential paths of deconstruction in the level- and location-matched aged cadaveric discs. Machine learning methods predicted a ‘hydration matrisome’ that connects extracellular matrix with MRI intensity. Importantly, the static proteome used as point-references can be integrated with dynamic proteome (SILAC/degradome) and transcriptome data from multiple clinical samples, enhancing robustness and clinical relevance. The data, findings, and methodology, available on a web interface (http://www.sbms.hku.hk/dclab/DIPPER/), will be valuable references in the field of IVD biology and proteomic analytics.

The backbone of vertebrate animals consists of a series of bones called vertebrae that are joined together by disc-like structures that allow the back to move and distribute forces to protect it during daily activities. It is common for these intervertebral discs to degenerate with age, resulting in back pain and severely reducing quality of life.

The mechanical features of intervertebral discs are the result of their proteins. These include extracellular matrix proteins, which form the external scaffolding that binds cells together in a tissue, and signaling proteins, which allow cells to communicate. However, how the levels of different proteins in each region of the disc vary with time has not been fully examined.

To establish how protein composition changes with age, Tam, Chen et al. quantified the protein levels and gene activity (which leads to protein production) of intervertebral discs from young and old deceased individuals. They found that the position of different mixtures of proteins in the intervertebral disc changes with age, and that young people have high levels of extracellular matrix proteins and signaling proteins. Levels of these proteins decreased as people got older, as did the amount of proteins produced.

To determine which region of the intervertebral disc different proteins were in, Tam, Chen et al. also performed magnetic resonance imaging (MRI) of the samples to correlate image intensity (which represents water content) with the corresponding protein signature. The data obtained provides a high-quality map of how the location of different proteins changes with age, and is available online under the name DIPPER. This database is an informative resource for research into skeletal biology, and it will likely advance the understanding of intervertebral disc degeneration in humans and animals, potentially leading to the development of new treatment strategies for this condition.

FPGA Accelerated Analysis of Boolean Gene Regulatory Networks

Boolean models are a powerful abstraction for qualitative modeling of gene regulatory networks. With the recent availability of advanced high-throughput technologies, Boolean models have increasingly grown in size and complexity, posing a challenge for existing software simulation tools that have not scaled at the same speed. Field Programmable Gate Arrays (FPGAs) are powerful reconfigurable integrated circuits that can offer massive performance improvements. Due to their highly parallel nature, FPGAs are well suited to simulate complex molecular networks. We present here a new simulation framework for Boolean models, which first converts the model to Verilog, a standardized hardware description language, and then connects it to an execution core that runs on an FPGA coherently attached to a POWER8 processor. We report an order of magnitude speedup over a multi-threaded software simulation tool running on the same processor on a selection of Boolean models. Analysis on a T-cell large granular lymphocyte leukemia (T-LGL) demonstrates that our framework achieves consistent performance improvements resulting in new biological insights. In addition, we show that our solution allows to perform attractor detection at an unprecedented speed, exhibiting a speedup ranging from one to three orders of magnitude compared to alternative software solutions.

A Novel Multiplex Based Platform for Osteoarthritis Drug Candidate Evaluation

Osteoarthritis (OA) is characterized by irreversible cartilage degradation with very limited therapeutic interventions. Drug candidates targeted at prototypic players had limited success until now and systems based approaches might be necessary. Consequently, drug evaluation platforms should consider the biological complexity looking beyond well-known contributors of OA. In this study an ex vivo model of cartilage degradation, combined with measuring releases of 27 proteins, was utilized to study 9 drug candidates. After an initial single drug evaluation step the 3 most promising compounds were selected and employed in an exhaustive combinatorial experiment. The resulting most and least promising treatment candidates were selected and validated in an independent study. This included estimation of mechanical properties via finite element modelling (FEM) and quantification of cartilage degradation as glycosaminoglycan (GAG) release. The most promising candidate showed increase of Young's modulus, decrease of hydraulic permeability and decrease of GAG release. The least promising candidate exhibited the opposite behaviour. The study shows the potential of a novel drug evaluation platform in identifying treatments that might reduce cartilage degradation. It also demonstrates the promise of exhaustive combination experiments and a connection between chondrocyte responses at the molecular level with changes of biomechanical properties at the tissue level.

An ex vivo tissue model of cartilage degradation suggests that cartilage state can be determined from secreted key protein patterns

The pathophysiology of osteoarthritis (OA) involves dysregulation of anabolic and catabolic processes associated with a broad panel of proteins that ultimately lead to cartilage degradation. An increased understanding about these protein interactions with systematic in vitro analyses may give new ideas regarding candidates for treatment of OA related cartilage degradation. Therefore, an ex vivo tissue model of cartilage degradation was established by culturing tissue explants with bacterial collagenase II. Responses of healthy and degrading cartilage were analyzed through protein abundance in tissue supernatant with a 26-multiplex protein profiling assay, after exposing the samples to a panel of 55 protein stimulations present in synovial joints of OA patients. Multivariate data analysis including exhaustive pairwise variable subset selection identified the most outstanding changes in measured protein secretions. MMP9 response to stimulation was outstandingly low in degrading cartilage and there were several protein pairs like IFNG and MMP9 that can be used for successful discrimination between degrading and healthy samples. The discovered changes in protein responses seem promising for accurate detection of degrading cartilage. The ex vivo model seems interesting for drug discovery projects related to cartilage degradation, for example when trying to uncover the unknown interactions between secreted proteins in healthy and degrading tissues.

Primary human chondrocytes respond to compression with phosphoproteomic signatures that include microtubule activation

Chondrocytes are responsible for maintaining the cartilage that helps joints bear load and move smoothly. These cells typically respond to physiological compression with pathways consistent with matrix synthesis, and chondrocyte mechanotransduction is essential for homeostasis. In osteoarthritis (OA), chondrocyte mechanotransduction appears to be dysregulated, yet the mechanisms remain poorly understood. The objective of this study is to document the phosphoproteomic responses of primary osteoarthritic chondrocytes to physiological sinusoidal compression. We show that OA chondrocytes respond to physiological compression by first activating proteins consistent with cytoskeletal remodeling and decreased transcription, and then later activating proteins for transcription. These results show that several microtubule-related proteins respond to compression. Our results demonstrate that compression is a relevant physiological stimulus for osteoarthritic chondrocytes. Future analyses may build on these results to find differences in compression-induced phosphoproteins between normal and OA cells that lead to druggable targets to restore homeostasis to diseased joints.

Computational Modeling and Reverse Engineering to Reveal Dominant Regulatory Interactions Controlling Osteochondral Differentiation: Potential for Regenerative Medicine

The specialization of cartilage cells, or chondrogenic differentiation, is an intricate and meticulously regulated process that plays a vital role in both bone formation and cartilage regeneration. Understanding the molecular regulation of this process might help to identify key regulatory factors that can serve as potential therapeutic targets, or that might improve the development of qualitative and robust skeletal tissue engineering approaches. However, each gene involved in this process is influenced by a myriad of feedback mechanisms that keep its expression in a desirable range, making the prediction of what will happen if one of these genes defaults or is targeted with drugs, challenging. Computer modeling provides a tool to simulate this intricate interplay from a network perspective. This paper aims to give an overview of the current methodologies employed to analyze cell differentiation in the context of skeletal tissue engineering in general and osteochondral differentiation in particular. In network modeling, a network can either be derived from mechanisms and pathways that have been reported in the literature (knowledge-based approach) or it can be inferred directly from the data (data-driven approach). Combinatory approaches allow further optimization of the network. Once a network is established, several modeling technologies are available to interpret dynamically the relationships that have been put forward in the network graph (implication of the activation or inhibition of certain pathways on the evolution of the system over time) and to simulate the possible outcomes of the established network such as a given cell state. This review provides for each of the aforementioned steps (building, optimizing, and modeling the network) a brief theoretical perspective, followed by a concise overview of published works, focusing solely on applications related to cell fate decisions, cartilage differentiation and growth plate biology. Particular attention is paid to an in-house developed example of gene regulatory network modeling of growth plate chondrocyte differentiation as all the aforementioned steps can be illustrated. In summary, this paper discusses and explores a series of tools that form a first step toward a rigorous and systems-level modeling of osteochondral differentiation in the context of regenerative medicine.

MicroRNAs and Osteoarthritis

An imbalance in gene expressional events skewing chondrocyte anabolic and catabolic pathways toward the latter causes an aberrant turnover and loss of extracellular matrix proteins in osteoarthritic (OA) articular cartilage. Thus, catabolism results in the elevated loss of extracellular matrix proteins. There is also evidence of an increase in the frequency of chondrocyte apoptosis that compromises the capacity of articular cartilage to undergo repair. Although much of the fundamental OA studies over the past 20 years identified and characterized many genes relevant to pro-inflammatory cytokines, apoptosis, and matrix metalloproteinases (MMPs)/a disintegrin and metalloproteinase with thrombospondin motif (ADAMTS), more recent studies focused on epigenetic mechanisms and the associated role of microRNAs (miRs) in regulating gene expression in OA cartilage. Thus, several miRs were identified as regulators of chondrocyte signaling pathways, apoptosis, and proteinase gene expression. For example, the reduced expression of miR-146a was found to be coupled to reduced type II collagen (COL2) in OA cartilage, whereas MMP-13 levels were increased, suggesting an association between MMP-13 gene expression and COL2A1 gene expression. Results of these studies imply that microRNAs could become useful in the search for diagnostic biomarkers, as well as providing novel therapeutic targets for intervention in OA.

A device for high-throughput monitoring of degradation in soft tissue samples

This work describes the design and validation of a novel device, the High-Throughput Degradation Monitoring Device (HDD), for monitoring the degradation of 24 soft tissue samples over incubation periods of several days inside a cell culture incubator. The device quantifies sample degradation by monitoring its deformation induced by a static gravity load. Initial instrument design and experimental protocol development focused on quantifying cartilage degeneration. Characterization of measurement errors, caused mainly by thermal transients and by translating the instrument sensor, demonstrated that HDD can quantify sample degradation with <6 μm precision and <10 μm temperature-induced errors. HDD capabilities were evaluated in a pilot study that monitored the degradation of fresh ex vivo human cartilage samples by collagenase solutions over three days. HDD could robustly resolve the effects of collagenase concentration as small as 0.5 mg/ml. Careful sample preparation resulted in measurements that did not suffer from donor-to-donor variation (coefficient of variance <70%). Due to its unique combination of sample throughput, measurement precision, temporal sampling and experimental versality, HDD provides a novel biomechanics-based experimental platform for quantifying the effects of proteins (cytokines, growth factors, enzymes, antibodies) or small molecules on the degradation of soft tissues or tissue engineering constructs. Thereby, HDD can complement established tools and in vitro models in important applications including drug screening and biomaterial development.

Influence of cartilage interstitial fluid on gene expression in cruciate ligament fibroblasts

Loading of articular cartilage during motion squeezes the fluid from the cartilage, termed cartilage interstitial fluid (CIF), which was found to influence gene expression in synovial membrane cells. After crucial ligaments damage, these cells are exposed to synovial fluid containing factors released from articular cartilage; the aim of the present study was to establish the influence of CIF and factors present in CIF (CIF-like cocktails) on crucial ligament fibroblasts. CIF was squeezed from articular-epiphyseal cartilage complexes of newborn rats. Fibroblasts were obtained from crucial ligaments of adult rat knee joints. Cells were cultured in control medium, CIF and CIF-like cocktails, and the expression of selected genes was evaluated using quantitative PCR. CIF stimulated the expression of HAS1, HAS2, aggrecan, lubricin, MMP3, TIMP3 and TGFβ1. Expression of collagen type I, versican, MMP2, TIMP2, TNF and IL1β was inhibited. The CIF-like cocktail stimulated HAS1, HAS2, collagen type I, versican, aggrecan, lubricin, TIMP1, TGFβ1, IL1β, IL6 and inhibited of MMP3 and TNF expression. Both agents exerted similar effects on the expression of HAS2, aggrecan, lubricin, TGFβ1 and TNF. CIF contains inhibitory and stimulatory factors affecting gene expression in crucial ligament fibroblasts and some of them were not included in the CIF-like cocktail. Due to the powerful influence of CIF on crucial ligament fibroblasts and the synovial membrane, further studies on its composition are needed. An improved CIF like-cocktail could be applied in the treatment of various joint or tendon ailments.

Cue-Signal-Response Analysis in 3D Chondrocyte Scaffolds with Anabolic Stimuli

Articular cartilage is an avascular connective tissue responsible for bearing loads. Cell signaling plays a central role in cartilage homeostasis and tissue engineering by directing chondrocytes to synthesize/degrade the extracellular matrix or promote inflammatory responses. The aim of this paper was to investigate anabolic, catabolic and inflammatory pathways of well-known and underreported anabolic stimuli in 3D chondrocyte cultures and connect them to diverse cartilage responses including matrix regeneration and cell communication. A cue-signal-response experiment was performed in chondrocytes embedded in alginate scaffolds subjected to a 9-day treatment with 7 anabolic cues. At the signaling level diverse pathways were measured whereas at the response level glycosaminoglycan (GAG) synthesis and cytokine releases were monitored. A significant increase of GAG was observed for each stimulus and well known anabolic phosphoproteins were activated. In addition, WNK1, an underreported protein of chondrocyte signaling, was uncovered. At the extracellular level, inflammatory and regulating cytokines were measured and DEFB1 and CXCL10 were identified as novel contributors to chondrocyte responses, both closely linked to TLR signaling and inflammation. Finally, two new pro-growth factors with an inflammatory potential, Cadherin-11 and MGP were observed. Interestingly, well-known anabolic stimuli yielded inflammatory responses which pinpoints to the pleiotropic roles of individual stimuli.

Systems approaches in osteoarthritis: Identifying routes to novel diagnostic and therapeutic strategies

Systems orientated research offers the possibility of identifying novel therapeutic targets and relevant diagnostic markers for complex diseases such as osteoarthritis. This review demonstrates that the osteoarthritis research community has been slow to incorporate systems orientated approaches into research studies, although a number of key studies reveal novel insights into the regulatory mechanisms that contribute both to joint tissue homeostasis and its dysfunction. The review introduces both top-down and bottom-up approaches employed in the study of osteoarthritis. A holistic and multiscale approach, where clinical measurements may predict dysregulation and progression of joint degeneration, should be a key objective in future research. The review concludes with suggestions for further research and emerging trends not least of which is the coupled development of diagnostic tests and therapeutics as part of a concerted effort by the osteoarthritis research community to meet clinical needs. © 2017 The Authors. Journal of Orthopaedic Research Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 35:1573-1588, 2017.

A Role for Soluble IL-6 Receptor in Osteoarthritis

Interleukin-6 (IL-6) is one of several pro-inflammatory cytokines present at elevated levels in the synovial fluid of individuals with confirmed clinical diagnosis of rheumatoid arthritis (RA) and osteoarthritis (OA). The mechanism of action of IL-6 was shown to involve its capacity to interact with a membrane-bound IL-6 receptor (mIL-6Rα), also known as the "classical" IL-6 pathway, or through its interaction with a soluble IL-6 receptor (sIL-6R) termed the "trans-signaling" pathway. Activation of downstream signaling is transduced via these IL-6 receptors and principally involves the Janus Kinase/Signal Transduction and Activators of Transcription (JAK/STAT) signaling pathway that is further regulated by glycoprotein-130 (gp130) interacting with the IL-6/mIL-6R complex. Phosphorylation of STAT proteins via JAK activation facilitates STAT proteins to act as transcription factors in inflammation. However, the biological function(s) of the sIL-6R in human chondrocytes requires further elucidation, although we previously showed that exogenous sIL-6R significantly suppressed the synthesis of neutrophil gelatinase-associated lipocalin (NGAL) in the immortalized line of human chondrocytes, C28/I2. NGAL was shown to regulate the activity of matrix metalloproteinase-9 (MMP-9), whose activity is crucial in OA for the destruction of articular cartilage. The "shedding" of sIL-6R from the plasma membrane is carried out by a family of enzymes known as A Distintegrin and Metalloproteinase (ADAM), which are also elevated in OA. In this paper, we have systematically reviewed the role played by IL-6 in OA. We have proposed that sIL-6R may be an important target for future drug development in OA by ameliorating cartilage extracellular protein degradation.

A Qualitative Model of the Differentiation Network in Chondrocyte Maturation: A Holistic View of Chondrocyte Hypertrophy

Differentiation of chondrocytes towards hypertrophy is a natural process whose control is essential in endochondral bone formation. It is additionally thought to play a role in several pathophysiological processes, with osteoarthritis being a prominent example. We perform a dynamic analysis of a qualitative mathematical model of the regulatory network that directs this phenotypic switch to investigate the influence of the individual factors holistically. To estimate the stability of a SOX9 positive state (associated with resting/proliferation chondrocytes) versus a RUNX2 positive one (associated with hypertrophy) we employ two measures. The robustness of the state in canalisation (size of the attractor basin) is assessed by a Monte Carlo analysis and the sensitivity to perturbations is assessed by a perturbational analysis of the attractor. Through qualitative predictions, these measures allow for an in silico screening of the effect of the modelled factors on chondrocyte maintenance and hypertrophy. We show how discrepancies between experimental data and the model’s results can be resolved by evaluating the dynamic plausibility of alternative network topologies. The findings are further supported by a literature study of proposed therapeutic targets in the case of osteoarthritis.

Applications of Chondrocyte-Based Cartilage Engineering: An Overview

Chondrocytes are the exclusive cells residing in cartilage and maintain the functionality of cartilage tissue. Series of biocomponents such as different growth factors, cytokines, and transcriptional factors regulate the mesenchymal stem cells (MSCs) differentiation to chondrocytes. The number of chondrocytes and dedifferentiation are the key limitations in subsequent clinical application of the chondrocytes. Different culture methods are being developed to overcome such issues. Using tissue engineering and cell based approaches, chondrocytes offer prominent therapeutic option specifically in orthopedics for cartilage repair and to treat ailments such as tracheal defects, facial reconstruction, and urinary incontinence. Matrix-assisted autologous chondrocyte transplantation/implantation is an improved version of traditional autologous chondrocyte transplantation (ACT) method. An increasing number of studies show the clinical significance of this technique for the chondral lesions treatment. Literature survey was carried out to address clinical and functional findings by using various ACT procedures. The current study was conducted to study the pharmacological significance and biomedical application of chondrocytes. Furthermore, it is inferred from the present study that long term follow-up studies are required to evaluate the potential of these methods and specific positive outcomes.

Influence of cartilage interstitial fluid on the mRNA levels of matrix proteins, cytokines, metalloproteases and their inhibitors in synovial membrane

Articular cartilage and the synovial membrane both ensure the smooth action of synovial joints; however, the influence of chondrocytes on synovial metabolism remains unclear. The secretory activity of chondrocytes is usually studied in cell cultures and may differ from that in intact cartilage. According to McCutchen's theory of 'weeping' joint lubrication, loading of the articular cartilage during motion squeezes the fluid with lubricating properties from the cartilage. The purpose of the study was to obtain cartilage interstitial fluid (CIF) from intact cartilage and to evaluate its influence on gene expression in the synovial membrane cells. CIF was rinsed out from the cartilage of newborn rats at a pressure of three bar. The chondrocytes survived rinsing and grew in culture. Cytokines in CIF were detected using the enzyme-linked immunosorbent assay (ELISA). The influence of CIF and CIF-like cocktail (all cytokines found in CIF) on gene expression in the synovial membrane cells was studied after a 4 h-incubation, by real-time PCR. Data were analyzed using the Wilcoxon matched-pair test or by the Mann‑Whitney U test. CIF contained basic fibroblast growth factor (bFGF), insulin-like growth factor (IGF)‑1, transforming growth factor β1 (TGFβ1), bone morphogenetic protein 7 (BMP7), macrophage (M)-colony-stimulating factor (CSF), granulocyte (G)-CSF and leukemia inhibitory factor (LIF). CIF stimulated the expression of hyaluronan synthase (HAS)1 and 2, lubricin, collagen I, versican, aggrecan, matrix metalloproteinases (MMPs)2 and 3, tissue inhibitors of metalloproteinases (TIMPs) 1-3, interleukin (IL)-6 and TGFβ1, and decreased the expression of tumor necrosis factor (TNF) and IL-1β. Incubation of the synovial membrane with CIF-like cocktail partially imitated the effects of CIF. Analysis of CIF composition may help to characterize the secretory activity of chondrocytes in their natural environment under various physiological and pathological conditions and to understand the interactions between articular cartilage and the synovial membrane.

Biologic basis of osteoarthritis: state of the evidence

PURPOSE OF REVIEW

Evidence accumulated since 2010 indicates that human osteoarthritis should now be reclassified as a systemic musculoskeletal disease rather than a focal disorder of synovial joints.

RECENT FINDINGS

Inflammation was seen as the key component promoting synovitis as well as progression of cartilage and bone destruction in osteoarthritis. Thus, metabolic-triggered inflammation involving cytokines, adipokines, abnormal metabolites, acute phase reactants and even complement, all appear to play major roles in osteoarthritis pathophysiology. Immune-mediated inflammation involving T cells and B cells as well as macrophages is now considered a common finding in osteoarthritis synovial tissue. Many experimental and clinical analyses showed that the proinflammatory cytokines, which stimulate matrix metalloproteinase and a disintegrin and metalloproteinase with thrombospondin motif gene transcription in normal and osteoarthritis human chondrocyte cultures, are also present at significantly elevated levels in the synovial fluid of osteoarthritis patients compared with nonarthritic synovial fluids.

SUMMARY

Human osteoarthritis is a systemic musculoskeletal disorder involving activation of innate and adaptive immune systems accompanied by inflammation exemplified by the elevated production of proinflammatory cytokines, which play a significant role in the progression of the disease. The future of novel therapies for osteoarthritis should consider developing drug development strategies designed to inhibit proinflammatory cytokine-induced signal transduction. These strategies have been successful in the development of drugs for the treatment of rheumatoid arthritis.

Network reconstruction based on proteomic data and prior knowledge of protein connectivity using graph theory

Modeling of signal transduction pathways is instrumental for understanding cells’ function. People have been tackling modeling of signaling pathways in order to accurately represent the signaling events inside cells’ biochemical microenvironment in a way meaningful for scientists in a biological field. In this article, we propose a method to interrogate such pathways in order to produce cell-specific signaling models. We integrate available prior knowledge of protein connectivity, in a form of a Prior Knowledge Network (PKN) with phosphoproteomic data to construct predictive models of the protein connectivity of the interrogated cell type. Several computational methodologies focusing on pathways’ logic modeling using optimization formulations or machine learning algorithms have been published on this front over the past few years. Here, we introduce a light and fast approach that uses a breadth-first traversal of the graph to identify the shortest pathways and score proteins in the PKN, fitting the dependencies extracted from the experimental design. The pathways are then combined through a heuristic formulation to produce a final topology handling inconsistencies between the PKN and the experimental scenarios. Our results show that the algorithm we developed is efficient and accurate for the construction of medium and large scale signaling networks. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGF/TNFA stimulation against made up experimental data. To avoid the possibility of erroneous predictions, we performed a cross-validation analysis. Finally, we validate that the introduced approach generates predictive topologies, comparable to the ILP formulation. Overall, an efficient approach based on graph theory is presented herein to interrogate protein–protein interaction networks and to provide meaningful biological insights.

Cytokine profiles in the joint depend on pathology, but are different between synovial fluid, cartilage tissue and cultured chondrocytes

Introduction

This study aimed to evaluate whether profiles of several soluble mediators in synovial fluid and cartilage tissue are pathology-dependent and how their production is related to in vitro tissue formation by chondrocytes from diseased and healthy tissue.

Methods

Samples were obtained from donors without joint pathology (n = 39), with focal defects (n = 65) and osteoarthritis (n = 61). A multiplex bead assay (Luminex) was performed measuring up to 21 cytokines: Interleukin (IL)-1α, IL-1β, IL-1RA, IL-4, IL-6, IL-6Rα, IL-7, IL-8, IL-10, IL-13, tumor necrosis factor (TNF)α, Interferon (IFN)γ, oncostatin M (OSM), leukemia inhibitory factor (LIF), adiponectin, leptin, monocyte chemotactic factor (MCP)1, RANTES, basic fibroblast growth factor (bFGF), hepatocyte growth factor (HGF), vascular growth factor (VEGF).

Results

In synovial fluid of patients with cartilage pathology, IL-6, IL-13, IFNγ and OSM levels were higher than in donors without joint pathology (P ≤0.001). IL-13, IFNγ and OSM were also different between donors with cartilage defects and OA (P <0.05). In cartilage tissue from debrided defects, VEGF was higher than in non-pathological or osteoarthritic joints (P ≤0.001). IL-1α, IL-6, TNFα and OSM concentrations (in ng/ml) were markedly higher in cartilage tissue than in synovial fluid (P <0.01). Culture of chondrocytes generally led to a massive induction of most cytokines (P <0.001). Although the release of inflammatory cytokines was also here dependent on the pathological condition (P <0.001) the actual profiles were different from tissue or synovial fluid and between non-expanded and expanded chondrocytes. Cartilage formation was lower by healthy unexpanded chondrocytes than by osteoarthritic or defect chondrocytes.

Conclusions

Several pro-inflammatory, pro-angiogenic and pro-repair cytokines were elevated in joints with symptomatic cartilage defects and/or osteoarthritis, although different cytokines were elevated in synovial fluid compared to tissue or cells. Hence a clear molecular profile was evident dependent on disease status of the joint, which however changed in composition depending on the biological sample analysed. These alterations did not affect in vitro tissue formation with these chondrocytes, as this was at least as effective or even better compared to healthy chondrocytes.