Modelling and analysis of biochemical signalling pathway cross-talk

Signalling pathways are abstractions that help life scientists structure the coordination of cellular activity. Cross-talk between pathways accounts for many of the complex behaviours exhibited by signalling pathways and is often critical in producing the correct signal-response relationship. Formal models of signalling pathways and cross-talk in particular can aid understanding and drive experimentation. We define an approach to modelling based on the concept that a pathway is the (synchronising) parallel composition of instances of generic modules (with internal and external labels). Pathways are then composed by (synchronising) parallel composition and renaming; different types of cross-talk result from different combinations of synchronisation and renaming. We define a number of generic modules in PRISM and five types of cross-talk: signal flow, substrate availability, receptor function, gene expression and intracellular communication. We show that Continuous Stochastic Logic properties can both detect and distinguish the types of cross-talk. The approach is illustrated with small examples and an analysis of the cross-talk between the TGF-b/BMP, WNT and MAPK pathways

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

Modelling the cAMP pathway using BioNessie, and the use of BVP techniques for solving ODEs (Poster Presentation)

Copyright @ 2007 Gu et al; licensee BioMed Central LtdBiochemists often conduct experiments in-vivo in order to explore observable behaviours and understand the dynamics of many intercellular and intracellular processes. However an intuitive understanding of their dynamics is hard to obtain because most pathways of interest involve components connected via interlocking loops. Formal methods for modelling and analysis of biochemical pathways are therefore indispensable. To this end, ODEs (ordinary differential equations) have been widely adopted as a method to model biochemical pathways because they have an unambiguous mathematical format and are amenable to rigorous quantitative analysis. BioNessie http://www.bionessie.com webcite is a workbench for the composition, simulation and analysis of biochemical networks which is being developed in by the Systems Biology team at the Bioinformatics Research Centre as a part of a large DTI funded project 'BPS: A Software Tool for the Simulation and Analysis of Biochemical Networks' http://www.brc.dcs.gla.ac.uk/projects/dti_beacon webcite. BioNessie is written in Java using NetBeans Platform libraries that makes it platform independent. The software employs specialised differential equations solvers for stiff and non-stiff systems to produce model simulation traces. BioNessie provides a user-friendly interfact that comes up with an intuitive tree-based graphical layout, an edition function to SBML-compatible models and feature of data output

When kinases meet mathematics: the systems biology of MAPK signalling

The mitogen activated protein kinase/extracellular signal regulated kinase pathway regulates fundamental cellular function such as cell proliferation, survival, differentiation and motility, raising the question how these diverse functions are specified and coordinated. They are encoded through the activation kinetics of the pathway, a multitude of feedback loops, scaffold proteins, subcellular compartmentalisation, and crosstalk with other pathways. These regulatory motifs alone or in combination can generate a multitude of complex behaviour. Systems biology tries to decode this complexity through mathematical modelling and prediction in order to gain a deeper insight into the inner works of signalling networks

Mathematical modelling plant signalling networks

During the last two decades, molecular genetic studies and the completion of the sequencing of the Arabidopsis thaliana genome have increased knowledge of hormonal regulation in plants. These signal transduction pathways act in concert through gene regulatory and signalling networks whose main components have begun to be elucidated. Our understanding of the resulting cellular processes is hindered by the complex, and sometimes counter-intuitive, dynamics of the networks, which may be interconnected through feedback controls and cross-regulation. Mathematical modelling provides a valuable tool to investigate such dynamics and to perform in silico experiments that may not be easily carried out in a laboratory. In this article, we firstly review general methods for modelling gene and signalling networks and their application in plants. We then describe specific models of hormonal perception and cross-talk in plants. This sub-cellular analysis paves the way for more comprehensive mathematical studies of hormonal transport and signalling in a multi-scale setting

A rule-based kinetic model of RNA polymerase II C-terminal domain phosphorylation

The complexity ofmany RNA processing pathways is such that a conventional systemsmodelling approach is inadequate to represent all themolecular species involved. We demonstrate that rule-based modelling permits a detailed model of a complex RNA signalling pathway to be defined. Phosphorylation of the RNApolymerase II (RNAPII)C-terminal domain (CTD; a flexible tail-like extension of the largest subunit) couples pre-messenger RNA capping, splicing and 30 end maturation to transcriptional elongation and termination, and plays a central role in integrating these processes. The phosphorylation states of the serine residues of many heptapeptide repeats of the CTD alter along the coding region of genes as a function of distance from the promoter. From a mechanistic perspective, both the changes in phosphorylation and the location atwhich they take place on the genes are a function of the time spent byRNAPII in elongation as this interval provides the opportunity for the kinases and phosphatases to interactwith theCTD.On this basis,we synthesize the available data to create a kinetic model of the action of the known kinases and phosphatases to resolve the phosphorylation pathways and their kinetics.</p

A predictive computational model reveals that GIV/girdin serves as a tunable valve for EGFR-stimulated cyclic AMP signals.

Cellular levels of the versatile second messenger cyclic (c)AMP are regulated by the antagonistic actions of the canonical G protein → adenylyl cyclase pathway that is initiated by G-protein-coupled receptors (GPCRs) and attenuated by phosphodiesterases (PDEs). Dysregulated cAMP signaling drives many diseases; for example, its low levels facilitate numerous sinister properties of cancer cells. Recently, an alternative paradigm for cAMP signaling has emerged in which growth factor-receptor tyrosine kinases (RTKs; e.g., EGFR) access and modulate G proteins via a cytosolic guanine-nucleotide exchange modulator (GEM), GIV/girdin; dysregulation of this pathway is frequently encountered in cancers. In this study, we present a network-based compartmental model for the paradigm of GEM-facilitated cross-talk between RTKs and G proteins and how that impacts cellular cAMP. Our model predicts that cross-talk between GIV, Gαs, and Gαi proteins dampens ligand-stimulated cAMP dynamics. This prediction was experimentally verified by measuring cAMP levels in cells under different conditions. We further predict that the direct proportionality of cAMP concentration as a function of receptor number and the inverse proportionality of cAMP concentration as a function of PDE concentration are both altered by GIV levels. Taking these results together, our model reveals that GIV acts as a tunable control valve that regulates cAMP flux after growth factor stimulation. For a given stimulus, when GIV levels are high, cAMP levels are low, and vice versa. In doing so, GIV modulates cAMP via mechanisms distinct from the two most often targeted classes of cAMP modulators, GPCRs and PDEs

Compositional modelling of signalling pathways in timed concurrent constraint programming

International audienceThe biological data regarding the signalling pathways often consider single pathways or a small number of them. We propose a methodology for composing this kind of data in a coherent framework, in order to be able to investigate a bigger number of signalling pathways. We specify a biological system by means of a set of stoichiometric-like equations resembling the essential features of molecular interactions. We represent these equations by a timed concurrent constraint (ntcc) language, which can deal with partial information and the time for a reaction to occur. We describe a freely available prototypical implementation of our framework

Reverse engineering of drug induced DNA damage response signalling pathway reveals dual outcomes of ATM kinase inhibition

The DNA Damage Response (DDR) pathway represents a signalling mechanism that is activated in eukaryotic cells following DNA damage and comprises of proteins involved in DNA damage detection, DNA repair, cell cycle arrest and apoptosis. This pathway consists of an intricate network of signalling interactions driving the cellular ability to recognise DNA damage and recruit specialised proteins to take decisions between DNA repair or apoptosis. ATM and ATR are central components of the DDR pathway. The activities of these kinases are vital in DNA damage induced phosphorylational induction of DDR substrates. Here, firstly we have experimentally determined DDR signalling network surrounding the ATM/ATR pathway induced following double stranded DNA damage by monitoring and quantifying time dependent inductions of their phosphorylated forms and their key substrates. We next involved an automated inference of unsupervised predictive models of time series data to generate in silico (molecular) interaction maps. We characterized the complex signalling network through system analysis and gradual utilisation of small time series measurements of key substrates through a novel network inference algorithm. Furthermore, we demonstrate an application of an assumption-free reverse engineering of the intricate signalling network of the activated ATM/ATR pathway. We next studied the consequences of such drug induced inductions as well as of time dependent ATM kinase inhibition on cell survival through further biological experiments. Intermediate and temporal modelling outcomes revealed the distinct signaling profile associated with ATM kinase activity and inhibition and explained the underlying signalling mechanism for dual ATM functionality in cytotoxic and cytoprotective pathways