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.

In Silico Laboratory Experiments Using Statistical Model Checking: A New Model of the Palytoxin-Induced Pump Channel as Case Study

Studying biological systems is a difficult but important task. Traditional methods include laboratory experimentation and computer simulations. However, often researchers need to explore important but potentially rare events that are not easily observed or simulated. We use UPPAAL-SMC, a formal verification tool to support a methodology that allows us to model biological systems, specify events and conditions that we want to analyze, and to explore system executions using controlled simulations. We also describe an efficient way to reproduce laboratory experiments in silico. Unlike traditional simulations, we are able to guide the experiment to explore special events and conditions by expressing these conditions in temporal logic formulas. We have applied this methodology to create a more detailed model of Palytoxin-induced Na ⁺/K ⁺ pump channels than was previously possible. Moreover, we have reproduced experimental protocols and their associated electrophysiological recordings, which has not been done in previous works. As a consequence, we have been able to propose a new diprotomeric model for the PTX-pump complex and study its behaviour. The use of our methodology has enabled us to reduce the effort and time to perform this research. It can be used to model and analyze other complex biological systems, potentially increasing the productivity of such studies.

A methodology for evaluating tooth wear monitoring using timed automata modelling

OBJECTIVES

Tooth wear is a multifactorial condition leading to the loss of dental hard tissues. A counselling/monitoring protocol is of importance in order to keep that loss as limited as possible. Since many factors are involved and a time span of decades is included, research to disentangle all these processes in patients is difficult. Instead, a modelling technique was used that is able to deal with time, costs and probabilistic and stochastic information. The aim was to shed light on the question: does a yearly or a once-in-five-years counselling/monitoring protocol yield better outcome measures?

METHODS

A so-called timed automata model was adopted, analysed with the tool UPPAAL. To our knowledge, this is the first time that formal modelling is applied in dentistry. In this article, a UPPAAL model for the evaluation of tooth wear is described.

RESULTS

Using the UPPAAL model, it was calculated that with a yearly counselling/monitoring protocol the severity of tooth wear at age 74, the total costs per person and the number of restorative treatments were less, and the number of so-called "good years" was higher.

CONCLUSIONS

With the use of the UPPAAL model, it may be concluded that a yearly counselling/monitoring protocol can yield better outcome measures.

CLINICAL SIGNIFICANCE

Regarding dentistry in general and tooth wear in particular, with the use of a timed automata model in UPPAAL, actual research questions can be answered, factors of influence in a multifactorial condition like tooth wear can be clarified, and future research topics can be determined.

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.

Modelling with ANIMO: between fuzzy logic and differential equations

BACKGROUND

Computational support is essential in order to reason on the dynamics of biological systems. We have developed the software tool ANIMO (Analysis of Networks with Interactive MOdeling) to provide such computational support and allow insight into the complex networks of signaling events occurring in living cells. ANIMO makes use of timed automata as an underlying model, thereby enabling analysis techniques from computer science like model checking. Biology experts are able to use ANIMO via a user interface specifically tailored for biological applications. In this paper we compare the use of ANIMO with some established formalisms on two case studies.

RESULTS

ANIMO is a powerful and user-friendly tool that can compete with existing continuous and discrete paradigms. We show this by presenting ANIMO models for two case studies: Drosophila melanogaster circadian clock, and signal transduction events downstream of TNF α and EGF in HT-29 human colon carcinoma cells. The models were originally developed with ODEs and fuzzy logic, respectively.

CONCLUSIONS

Two biological case studies that have been modeled with respectively ODE and fuzzy logic models can be conveniently modeled using ANIMO. The ANIMO models require less parameters than ODEs and are more precise than fuzzy logic. For this reason we position the modelling paradigm of ANIMO between ODEs and fuzzy logic.

BioNSi: A Discrete Biological Network Simulator Tool

Modeling and simulation of biological networks is an effective and widely used research methodology. The Biological Network Simulator (BioNSi) is a tool for modeling biological networks and simulating their discrete-time dynamics, implemented as a Cytoscape App. BioNSi includes a visual representation of the network that enables researchers to construct, set the parameters, and observe network behavior under various conditions. To construct a network instance in BioNSi, only partial, qualitative biological data suffices. The tool is aimed for use by experimental biologists and requires no prior computational or mathematical expertise. BioNSi is freely available at http://bionsi.wix.com/bionsi , where a complete user guide and a step-by-step manual can also be found.

A Semiquantitative Framework for Gene Regulatory Networks: Increasing the Time and Quantitative Resolution of Boolean Networks

Boolean models have been instrumental in predicting general features of gene networks and more recently also as explorative tools in specific biological applications. In this study we introduce a basic quantitative and a limited time resolution to a discrete (Boolean) framework. Quantitative resolution is improved through the employ of normalized variables in unison with an additive approach. Increased time resolution stems from the introduction of two distinct priority classes. Through the implementation of a previously published chondrocyte network and T helper cell network, we show that this addition of quantitative and time resolution broadens the scope of biological behaviour that can be captured by the models. Specifically, the quantitative resolution readily allows models to discern qualitative differences in dosage response to growth factors. The limited time resolution, in turn, can influence the reachability of attractors, delineating the likely long term system behaviour. Importantly, the information required for implementation of these features, such as the nature of an interaction, is typically obtainable from the literature. Nonetheless, a trade-off is always present between additional computational cost of this approach and the likelihood of extending the model’s scope. Indeed, in some cases the inclusion of these features does not yield additional insight. This framework, incorporating increased and readily available time and semi-quantitative resolution, can help in substantiating the litmus test of dynamics for gene networks, firstly by excluding unlikely dynamics and secondly by refining falsifiable predictions on qualitative behaviour.