Computational models for inferring biochemical networks

Biochemical networks are of great practical importance. The interaction of biological compounds in cells has been enforced to a proper understanding by the numerous bioinformatics projects, which contributed to a vast amount of biological information. The construction of biochemical systems (systems of chemical reactions), which include both topology and kinetic constants of the chemical reactions, is NP-hard and is a well-studied system biology problem. In this paper, we propose a hybrid architecture, which combines genetic programming and simulated annealing in order to generate and optimize both the topology (the network) and the reaction rates of a biochemical system. Simulations and analysis of an artificial model and three real models (two models and the noisy version of one of them) show promising results for the proposed method.The Romanian National Authority for Scientific Research, CNDI–UEFISCDI, Project No. PN-II-PT-PCCA-2011-3.2-0917

A structured approach for the engineering of biochemical network models, illustrated for signalling pathways

http://dx.doi.org/10.1093/bib/bbn026Quantitative models of biochemical networks (signal transduction cascades, metabolic pathways, gene regulatory circuits) are a central component of modern systems biology. Building and managing these complex models is a major challenge that can benefit from the application of formal methods adopted from theoretical computing science. Here we provide a general introduction to the field of formal modelling, which emphasizes the intuitive biochemical basis of the modelling process, but is also accessible for an audience with a background in computing science and/or model engineering. We show how signal transduction cascades can be modelled in a modular fashion, using both a qualitative approach { Qualitative Petri nets, and quantitative approaches { Continuous Petri Nets and Ordinary Differential Equations. We review the major elementary building blocks of a cellular signalling model, discuss which critical design decisions have to be made during model building, and present ..

Petri Net modelling approach for analysing the behaviour of Wnt/[inline-formula removed] -catenin and Wnt/Ca 2+ signalling pathways in arrhythmogenic right ventricular cardiomyopathy.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart muscle disease that may result in arrhythmia, heart failure and sudden death. The hallmark pathological findings are progressive myocyte loss and fibro fatty replacement, with a predilection for the right ventricle. This study focuses on the adipose tissue formation in cardiomyocyte by considering the signal transduction pathways including Wnt/[inline-formula removed]-catenin and Wnt/Ca2+ regulation system. These pathways are modelled and analysed using stochastic petri nets (SPN) in order to increase our comprehension of ARVC and in turn its treatment regimen. The Wnt/[inline-formula removed]-catenin model predicts that the dysregulation or absence of Wnt signalling, inhibition of dishevelled and elevation of glycogen synthase kinase 3 along with casein kinase I are key cytotoxic events resulting in apoptosis. Moreover, the Wnt/Ca2+ SPN model demonstrates that the Bcl2 gene inhibited by c-Jun N-terminal kinase protein in the event of endoplasmic reticulum stress due to action potential and increased amount of intracellular Ca2+ which recovers the Ca2+homeostasis by phospholipase C, this event positively regulates the Bcl2 to suppress the mitochondrial apoptosis which causes ARVC

Modular modelling of signalling pathways and their crosstalk

Signalling pathways are well-known abstractions that explain the mechanisms whereby cells respond to signals. Collections of pathways form networks, and interactions between pathways in a network, known as cross-talk, enables further complex signalling behaviours. While there are several formal modelling approaches for signalling pathways, none make cross-talk explicit; the aim of this paper is to define and categorise cross-talk in a rigorous way. We define a modular approach to pathway and network modelling, based on the module construct in the PRISM modelling language, and a set of generic signalling modules. Five different types of cross-talk are defined according to various biologically meaningful combinations of variable sharing, synchronisation labels and reaction renaming. The approach is illustrated with a case-study analysis of cross-talk between the TGF-β, WNT and MAPK pathways

Modeling formalisms in systems biology

Systems Biology has taken advantage of computational tools and high-throughput experimental data to model several biological processes. These include signaling, gene regulatory, and metabolic networks. However, most of these models are specific to each kind of network. Their interconnection demands a whole-cell modeling framework for a complete understanding of cellular systems. We describe the features required by an integrated framework for modeling, analyzing and simulating biological processes, and review several modeling formalisms that have been used in Systems Biology including Boolean networks, Bayesian networks, Petri nets, process algebras, constraint-based models, differential equations, rule-based models, interacting state machines, cellular automata, and agent-based models. We compare the features provided by different formalisms, and discuss recent approaches in the integration of these formalisms, as well as possible directions for the future.Research supported by grants SFRH/BD/35215/2007 and SFRH/BD/25506/2005 from the Fundacao para a Ciencia e a Tecnologia (FCT) and the MIT-Portugal Program through the project "Bridging Systems and Synthetic Biology for the development of improved microbial cell factories" (MIT-Pt/BS-BB/0082/2008)

Computational modelling of cancerous mutations in the EGFR/ERK signalling pathway

This article has been made available through the Brunel Open Access Publishing Fund - Copyright @ 2009 Orton et al.BACKGROUND: The Epidermal Growth Factor Receptor (EGFR) activated Extracellular-signal Regulated Kinase (ERK) pathway is a critical cell signalling pathway that relays the signal for a cell to proliferate from the plasma membrane to the nucleus. Deregulation of the EGFR/ERK pathway due to alterations affecting the expression or function of a number of pathway components has long been associated with numerous forms of cancer. Under normal conditions, Epidermal Growth Factor (EGF) stimulates a rapid but transient activation of ERK as the signal is rapidly shutdown. Whereas, under cancerous mutation conditions the ERK signal cannot be shutdown and is sustained resulting in the constitutive activation of ERK and continual cell proliferation. In this study, we have used computational modelling techniques to investigate what effects various cancerous alterations have on the signalling flow through the ERK pathway. RESULTS: We have generated a new model of the EGFR activated ERK pathway, which was verified by our own experimental data. We then altered our model to represent various cancerous situations such as Ras, B-Raf and EGFR mutations, as well as EGFR overexpression. Analysis of the models showed that different cancerous situations resulted in different signalling patterns through the ERK pathway, especially when compared to the normal EGF signal pattern. Our model predicts that cancerous EGFR mutation and overexpression signals almost exclusively via the Rap1 pathway, predicting that this pathway is the best target for drugs. Furthermore, our model also highlights the importance of receptor degradation in normal and cancerous EGFR signalling, and suggests that receptor degradation is a key difference between the signalling from the EGF and Nerve Growth Factor (NGF) receptors. CONCLUSION: Our results suggest that different routes to ERK activation are being utilised in different cancerous situations which therefore has interesting implications for drug selection strategies. We also conducted a comparison of the critical differences between signalling from different growth factor receptors (namely EGFR, mutated EGFR, NGF, and Insulin) with our results suggesting the difference between the systems are large scale and can be attributed to the presence/absence of entire pathways rather than subtle difference in individual rate constants between the systems.This work was funded by the Department of Trade and Industry (DTI), under their Bioscience Beacon project programme. AG was funded by an industrial PhD studentship from Scottish Enterprise and Cyclacel

Heparan Sulfate Induces Necroptosis in Murine Cardiomyocytes: A Medical-In silico Approach Combining In vitro Experiments and Machine Learning.

Life-threatening cardiomyopathy is a severe, but common, complication associated with severe trauma or sepsis. Several signaling pathways involved in apoptosis and necroptosis are linked to trauma- or sepsis-associated cardiomyopathy. However, the underling causative factors are still debatable. Heparan sulfate (HS) fragments belong to the class of danger/damage-associated molecular patterns liberated from endothelial-bound proteoglycans by heparanase during tissue injury associated with trauma or sepsis. We hypothesized that HS induces apoptosis or necroptosis in murine cardiomyocytes. By using a novel Medical-In silico approach that combines conventional cell culture experiments with machine learning algorithms, we aimed to reduce a significant part of the expensive and time-consuming cell culture experiments and data generation by using computational intelligence (refinement and replacement). Cardiomyocytes exposed to HS showed an activation of the intrinsic apoptosis signal pathway via cytochrome C and the activation of caspase 3 (both p < 0.001). Notably, the exposure of HS resulted in the induction of necroptosis by tumor necrosis factor α and receptor interaction protein 3 (p < 0.05; p < 0.01) and, hence, an increased level of necrotic cardiomyocytes. In conclusion, using this novel Medical-In silico approach, our data suggest (i) that HS induces necroptosis in cardiomyocytes by phosphorylation (activation) of receptor-interacting protein 3, (ii) that HS is a therapeutic target in trauma- or sepsis-associated cardiomyopathy, and (iii) indicate that this proof-of-concept is a first step toward simulating the extent of activated components in the pro-apoptotic pathway induced by HS with only a small data set gained from the in vitro experiments by using machine learning algorithms.This work was supported by an intramural grant to LM (START 46/16) and EZ (START 113/17). LM has received a grant by the Deutsche Forschungsgemeinschaft (DFG, MA 7082/1–1). We thank Dr Claycomb and his coworkers for providing the HL-1 cells and a detailed documentation. The Immunohistochemistry and Confocal Microscopy Unit, a core facility of the Interdisciplinary Center for Clinical Research (IZKF) Aachen, within the Faculty of Medicine at the RWTH Aachen University, supported this work