Biochemical simulation by stochastic concurrent constraint programming and hybrid systems
|
|
- Malcolm Campbell
- 5 years ago
- Views:
Transcription
1 Biochemical simulation by stochastic concurrent constraint programming and hybrid systems Luca Bortolussi 1 Alberto Policriti 2,3 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste luca@dmi.units.it 2 Dipartimento di Matematica ed Informatica Università degli studi di Udine 3 Istituto di Genomica Applicata Parco Scientifico Tecnologico, Udine. Siena, 20 th April 2007
2 A biological circuit: MAPKinase
3 A biological circuit: MAPKinase
4 Circadian Clock (J. M. G. Vilar, H. Yuan Kueh, N. Barkai, and S. Leibler. PNAS, 2002.)
5 Circadian Clock p_gate(α A, α A, γ A, θ A, M A, A) p_gate(α R, α R, γ R, θ R, M R, A) reaction(β A, [M A ], [A]) reaction(δ MA, [M A ], []) reaction(β R, [M R ], [R]) reaction(δ MR, [M R ], []) reaction(γ C, [A, R], [AR]) reaction(δ A, [AR], [R]) reaction(δ A, [A], []) reaction(δ R, [R], [])
6 Modeling Biological Systems with Stochastic Process Algebras Pros Cons Simple Language Compositionality Hard to encode general information Lacking computational extensibility
7 How to Increase Expressiveness... Add a computational twist to π-calculus. Add a form of local storage. Keep separated the description of interactions and the description of memory and computations. Introduce more general format for rates. Constraints... why not?
8 Outline 1 Introduction 2 Stochastic Concurrent Constraint Programming 3 Modeling Biological Systems in sccp 4 sccp and Hybrid Automata sccp and ODEs sccp and Hybrid Automata
9 Concurrent Constraint Programming Constraint Store In this process algebra, the main object are constraints, which are formulae over an interpreted first order language (i.e. X = 10, Y > X 3). Constraints can be added to a "container", the constraint store, but can never be removed. Agents Agents can perform two basic operations on this store (asynchronously): Add a constraint (tell ask) Ask if a certain relation is entailed by the current configuration (ask instruction) Syntax of CCP Program = Decl.A D = ε Decl.Decl p(x) : A A = 0 tell(c).a ask(c 1 ).A 1 + ask(c 2 ).A 2 A 1 A 2 x A p(x) V. Saraswat, Concurrent Constraint Programming, MIT press, 1993
10 Syntax of sccp Syntax of Stochastic CCP Program = D.A D = ε D.D p( x ) : A A = 0 tell (c).a M xa A A M = π.g M + M π = tell λ (c) ask λ (c) G = 0 tell (c).g p( y ) M xg G G L. Bortolussi, Stochastic Concurrent Constraint Programming, QAPL, 2006 Stochastic Rates Rates are functions from the constraint store C to positive reals: λ : C R +. Rates can be thought as speed or duration of communications.
11 sccp technical details Operational Semantics There are two transition relations, one instantaneous (finite and confluent) and one stochastic. Show Details Traces are sequences of events with variable time delays among them. Discrete vs. Continuous Semantics Show Details The operational semantics is abstract w.r.t. the notion of time: we can map the labeled transition system into a discrete or a continuous time Markov Chain. Implementation We have an interpreter written in Prolog, using the CLP engine of SICStus to manage the constraint store. Efficiency issues. Stream Variables Quantities varying over time can be represented in sccp as unbounded lists. Hereafter: special meaning of X = X + 1.
12 Continuous Time Markov Chains A Continuous Time Markov Chain (CTMC) is a direct graph with edges labeled by a real number, called the rate of the transition (representing the speed or the frequency at which the transition occurs). In each state, we select the next state according to a probability distribution obtained normalizing rates (from S to S 1 r with prob. 1 ). r 1 +r 2 The time spent in a state is given by an exponentially distributed random variable, with rate given by the sum of outgoing transitions from the actual node (r 1 + r 2 ).
13 General Principles Measurable Entities Stream Variables Logical Entities Processes (Control Variables) Interactions Processes L. Bortolussi, A. Policriti Modeling Biological Systems in Stochastic Concurrent Constraint Programming, 2006
14 Biochemical Arrows to sccp processes R R n k P P m R R n k 1 k2 P P m S E K,V 0 P S E K,V 0,h P reaction(k, [R 1,..., R V n], [P 1,..., P m]) : ask ni=1 rma (k,r 1,...,Rn) (R i > 0). n i=i tell (R i = R i 1) m j=1 tell (P j = P j + 1). reaction(k, [R 1,..., R n], [P 1,..., P m]) reaction(k 1, [R 1,..., R n], [P 1,..., P m]) reaction(k 2, [P 1,..., P m], [R 1,..., R n]) mm_reaction(k, V 0, S, P) : ask rmm (K,V 0,S)(S > 0). (tell (S = S 1) tell (P = P + 1)). mm_reaction(k, V 0, S, P) hill_reaction(k, V 0, h, S, P) : ask rhill (K,V 0,h,S)(S > 0). (tell (S = S h) tell (P = P + h)). Hill_reaction(K, V 0, h, S, P) where r MA (k, X 1,..., X n) = k X 1 X n; r MM (K, V 0, S) = V 0S S + K ; r Hill (k, V 0, h, S) = V h 0S S h + K h
15 Enzymatic reaction S + E k 1 k 1 ES k2 P + E Mass Action Kinetics enz_reaction(k 1, k 1, k 2, S, E, ES, P) :- reaction(k 1, [S, E], [ES]) reaction(k 1, [ES], [E, S]) reaction(k 2, [ES], [E, P]) Michaelis-Menten Equations d[p] = V 0 S dt S+K V 0 = k 2 [E 0 ] K = k 2 +k 1 k 1 Mass Action Equations d[es] = k dt 1 [S][E] k 2 [ES] k 1 [ES] d[e] = k dt 1 [S][E] + k 2 [ES] + k 1 [ES] d[s] = k dt 1 [S][E] d[p] = k dt 2 [ES] Michaelis-Menten Kinetics mm_reaction! k 2 + k 1, k 2 E, S, P k 1
16 Enzymatic reaction S + E k 1 k 1 ES k2 P + E Mass Action Kinetics enz_reaction(k 1, k 1, k 2, S, E, ES, P) :- reaction(k 1, [S, E], [ES]) reaction(k 1, [ES], [E, S]) reaction(k 2, [ES], [E, P]) Michaelis-Menten Kinetics! k 2 + k 1 mm_reaction, k 2 E, S, P k 1
17 MAP-Kinase cascade enz_reaction(k a, k d, k r, KKK, E1, KKKE1, KKKS) enz_reaction(k a, k d, k r, KKKS, E2, KKKSE2, KKK ) enz_reaction(k a, k d, k r, KK, KKKS, KKKKKS, KKP) enz_reaction(k a, k d, k r, KKP, KKP1, KKPKKP1, KK ) enz_reaction(k a, k d, k r, KKP, KKKS, KKPKKKS, KKPP) enz_reaction(k a, k d, k r, KP, KP1, KPKP1, K ) enz_reaction(k a, k d, k r, K, KKPP, KKKPP, KP) enz_reaction(k a, k d, k r, KKPP, KKP1, KKPPKKP1, KKP) enz_reaction(k a, k d, k r, KP, KKPP, KPKKPP, KPP) enz_reaction(k a, k d, k r, KPP, KP1, KPPKP1, KP)
18 The gene machine
19 The instruction set null_gate(k p, X) : tell kp (X = X + 1).null_gate(k p, X) pos_gate(k p, k e, k f, X, Y ) : tell kp (X = X + 1).pos_gate(k p, k e, k f, X, Y ) +ask r(ke,y ) (true).tell ke (X = X + 1).pos_gate(k p, k e, k f, X, Y ) neg_gate(k p, k i, k d, X, Y ) : tell kp (X = X + 1).neg_gate(k p, k i, k d, X, Y ) +ask r(ki,y ) (true).ask k d (true).neg_gate(k p, k i, k d, X, Y ) where r(k, Y ) = k Y. L. Cardelli, A. Phillips, 2005.
20 Repressilator neg_gate(0.1, 1, , A, C) reaction(0.0001, [A], []) neg_gate(0.1, 1, , B, A) reaction(0.0001, [B], []) neg_gate(0.1, 1, , C, B) reaction(0.0001, [C], [])
21 Miscellanea Model Checking Show Details We defined a translation procedure that maps (a restricted version of) sccp into PRISM, a probabilistic symbolic model checker. Hence, we can check if sccp programs satisfy properties specified as Continuous Stochastic Logic formulae.
22 Connecting SPA and ODE models
23 From sccp to ODE What? We want to associate a set of ODE to an sccp program (written with a restricted syntax). Why? ODE can be numerically simulated faster than stochastic processes. On the market... Show Details How things can go wrong There are (syntactic) methods to write set of ODEs for PEPA and stochastic π-calculus, looking at the speed of creation and destruction of terms (We did the same for sccp). However, the ODE can show a behavior different from that of SPA models. J. Hillston, Fluid Flow Approximation of PEPA models, QEST, L. Cardelli, From Processes to ODEs by Chemistry, L. Bortolussi, A. Policriti. Connecting Process Algebras and Differential Equations for systems biology, 2006.
24 The strange beast of repressilator Repressilator with gene gates Repressilator in sccp Repressilator: ODE from sccp Repressilator: average of sccp model Return
25 Hybrid Automata - Intuitively Intuitively, a hybrid automaton is a finite state automaton H with continuous variables Z such that: H evolves from Z to Z in time T according to differential equations Ż = f (v, Z ) Z must always satisfy invariant conditions in each state H can cross e only when guards on Z are true when H crosses e, some variables may be reset.
26 From sccp to Hybrid Automata What? We want to associate an hybrid automata to a sccp network. Why? The mixed discrete/continuous dynamics of HA is more natural, as it can preserve the logical structure of sccp models. Hybrid automata are equipped with well developed analysis methods. How? The separation between constraint store and logical description of agents makes easy to identify (discrete) modes of the automata Activation conditions need to look at the temporal semantics of stochastic actions. L. Bortolussi, A. Policriti. Hybrid approximation of Stochastic Concurrent Constraint Programming, 2007.
27 From sccp to Hybrid Automata: example A :- tell k+ (X = X + 1).A + ask k X (X > 0). tell (X = X 1).A + ask k0 (true).b B :- tell k (X = X + 1).B + ask k+ X (X > 0). tell (X = X 1).B + ask k0 (true).a
28 Hybrid Repressilator Repressilator with gene gates Repressilator in sccp Repressilator: ODE from sccp Hybrid Repressilator from sccp Analysis of Hybrid Repressilator
29 Conclusions We introduced a stochastic version of Concurrent Constraint Programming. We applied it in the field of biological modeling for a wide spectrum of problems: biochemical reactions, genetic networks, but also protein folding. The analysis techniques comprehend: simulation, model checking, solution of related ODEs. A new promising analysis method consists in associating hybrid automata to sccp models. We are currently investigating the use of sccp to simulate implicitly complex protein networks, in order to tame combinatorial explosion
30 Biology, Computation and Information 4 th International School on Biology, Computation, and Information (BCI 2007). 2-7 July, Area Science Park, Trieste.
31 sccp technical details Further biological applications sccp and ODE Operational Semantic: Instantaneous Transition Instantaneous Transition (IR1) tell (c), d, V 0, d c, V (IR2) D p( E D x ), d, V A[ x / E y ], d, V if p( y ) : A (IR3) x A, d, V A[y/x], d, V {y} with y V 2 \ V (IR4) (IR5) A 1, d, V A 1, d, V A 1.A 2, d, V A 1.A 2, d, V A 1, d, V A 1, d, V A 1 A 2, d, V A 1 A 2, d, V Theorem The instantaneous transition is confluent and can be applied only a finite number of steps to each configuration C. Return
32 sccp technical details Further biological applications sccp and ODE Operational Semantic: Stochastic Transition Stochastic Transition (SR1) tell λ (c), d, V = (1,λ(d)) 0, d c, V if d c false (SR2) ask λ (c), d, V = (1,λ(d)) 0, d, V if d c (SR3) π, d, V = (p,λ) 0, d, V π.a, d, V = (p,λ) A, d, V with π = ask or π = tell (SR4) (SR5) (SR6) A 1, d, V = (p,λ) A 1, d, V A 1.A 2, d, V = (p,λ) A 1.A 2, d, V M 1, d, V = (p,λ) A 1, d, V M 1 + M 2, d, V = (p,λ ) A 1, d, V with p pλ = λ+rate( M 2,d,V ) and λ = λ + rate( M 2, d, V ) A 1, d, V = (p,λ) A 1, d, V A 1 A 2, d, V = (p,λ ) A 1 A 2, d, V with p pλ = λ+rate( A 2,d,V ) and λ = λ + rate( A 2, d, V ) rate returns the sum of rates of all active agents.
33 sccp technical details Further biological applications sccp and ODE Operational Semantic: Stochastic Transition Theorem Let A, d, V C be the current configuration. Then the next stochastic transition executes one of the agents prefixed by a guard belonging to the set exec( A, d, V ), call it A. Moreover, the probability of the transition (i.e. the first label in = ) is rate( A, d, V ) rate (exec( A, d, V )), and the rate associated to the transition (the second label in = ) is rate (exec( A, d, V )). Return
34 sccp technical details Further biological applications sccp and ODE Rates Rates can be interpreted as priorities or as frequencies. Rates as Priorities A rate can represent the priority of execution of a process. There is a global scheduler choosing probabilistically between active processes, according to their priority. Discrete time evolution. Rates as Frequencies A rate can represent the frequency or speed of a process. The higher the speed, the higher the probability of seeing a certain process executed. Continuous time evolution. Return
35 sccp technical details Further biological applications sccp and ODE Discrete and Continuous Time Discrete time Discrete time transition can be recovered from stochastic transition = (p,λ) by dropping the second label. Hence we leave only the probability associated to transitions, obtaining a Discrete time Markov Chain. Continuous Time Continuous time transition can be recovered from stochastic transition = (p,λ) by multiplying the two labels. Hence we consider the rate associated to the transition, obtaining a Continuous time Markov Chain. Return
36 sccp technical details Further biological applications sccp and ODE Discrete and Continuous Time Observables Discrete time I/O observables O d ( A, d ) = { (d, p) p = Prob ( A, d 0, d )}. Continuous time I/O observables O c ( A, d )(t) = { (d, p) p = Prob ( A, d 0, d ) (t) }. Theorem lim O c( A, d )(t) = O d ( A, d ). t Return
37 sccp technical details Further biological applications sccp and ODE Continuous Stochastic Logic Return (state formulae) φ ::= true a φ φ φ P p [ψ] S p [φ] E c [φ] (path formulae) ψ ::= Xφ φu I φ φuφ, X: next. U: until. U I : bounded until, second formula must be true within interval I. P p [ψ]: path formula ψ must have probability p, with {<,,, >}. S p [φ]: states satisfying must have steady state probability p, with {<,,, >}. E c [φ]: expected reward of reaching a state satisfying φ is p, with {<,,, >}.
38 sccp technical details Further biological applications sccp and ODE PRISM: a probabilistic model checker Probabilistic Model Checking To model check a formula in CSL, we need essentially to solve linear systems in the size of the state space. Hence model checking is polynomial (but no more linear). PRISM PRISM is a symbolic model checking for CSL (and other probabilistic logics), using both Algebraic Decision Diagrams and sparse matrices to speed up the basic algorithms. Models can be described using a concurrent language having variable update and comparison as basic primitives. Model Checking sccp To model check (restricted) sccp, we defined an encoding of sccp programs into PRISM programs. Return
39 sccp technical details Further biological applications sccp and ODE Example: encoding S-System s Repressilator Ẋ 1 = α 1 X 1 3 β 1 X1 0.5, Ẋ 2 = α 2 X 1 1 β 2 X2 0.5, Ẋ 3 = α 3 X 1 2 β 3 X3 0.5, α i = 0.2, β i = 1. Return subs(x 1 ) subs(x 2 ) subs(x 3 ) ::= (tell(x 1 = X 1 + σ) [α1 X 1 3 ] +tell(x 1 = X 1 σ) [β1 X ] ).subs(x 1) ::= (tell(x 2 = X 2 + σ) [α2 X 1 1 ] +tell(x 2 = X 2 σ) [β2 X ] ).subs(x 2) ::= (tell(x 3 = X 3 + σ) [α3 X 1 2 ] +tell(x 3 = X 3 σ) [β3 X ] ).subs(x 3) σ = 1 Show more on behavioral equivalence
40 sccp technical details Further biological applications sccp and ODE Repressilator gone wild S-System s model of repressilator suffers from an high sensitivity from parameters, differently from the usual PA models. sccp model with variable rates has the same wild behaviour! β i = 0.01 β i = Return
41 sccp technical details Further biological applications sccp and ODE From sccp to ODE: example RW X :- tell 1 (X = X 1).RW X + tell 1 (X = X + 2).RW X ) + ask f (X) (true).( tell 1 (X = X 2).RW X + tell 1 (X = X + 1).RW X ) f (X) = 1 X 2 +1 X P P Ẋ = P 0 P 1 Ṗ 0 = 1 X 2 +1 P 0 + 2P 1 Ṗ 1 = 1 X 2 +1 P 0 2P 1 Return L. Bortolussi, A. Policriti. Stochastic Concurrent Constraint Programming and Differential Equations, 2007.
42 sccp technical details Further biological applications sccp and ODE Analysis of Hybrid Repressilator Determinism vs non-determinism Stability Return
Modeling Biological Systems in Stochastic Concurrent Constraint Programming
Modeling Biological Systems in Stochastic Concurrent Constraint Programming Luca Bortolussi 1 Alberto Policriti 1 1 Department of Mathematics and Computer Science University of Udine, Italy. Workshop on
More informationModeling Biological Systems in Stochastic Concurrent Constraint Programming
Modeling Biological Systems in Stochastic Concurrent Constraint Programming Luca Bortolussi Alberto Policriti Abstract We present an application of stochastic Concurrent Constraint Programming (sccp) for
More informationON THE APPROXIMATION OF STOCHASTIC CONCURRENT CONSTRAINT PROGRAMMING BY MASTER EQUATION
ON THE APPROXIMATION OF STOCHASTIC CONCURRENT CONSTRAINT PROGRAMMING BY MASTER EQUATION Luca Bortolussi 1,2 1 Department of Mathematics and Informatics University of Trieste luca@dmi.units.it 2 Center
More informationLecture 4 The stochastic ingredient
Lecture 4 The stochastic ingredient Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it
More informationStochastic Simulation.
Stochastic Simulation. (and Gillespie s algorithm) Alberto Policriti Dipartimento di Matematica e Informatica Istituto di Genomica Applicata A. Policriti Stochastic Simulation 1/20 Quote of the day D.T.
More informationLecture 1 Modeling in Biology: an introduction
Lecture 1 in Biology: an introduction Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it
More informationStochastic Process Algebra models of a Circadian Clock
Stochastic Process Algebra models of a Circadian Clock Jeremy T. Bradley Thomas Thorne Department of Computing, Imperial College London 180 Queen s Gate, London SW7 2BZ, United Kingdom Email: jb@doc.ic.ac.uk
More informationLogic-based Multi-Objective Design of Chemical Reaction Networks
Logic-based Multi-Objective Design of Chemical Reaction Networks Luca Bortolussi 1 Alberto Policriti 2 Simone Silvetti 2,3 1 DMG, University of Trieste, Trieste, Italy luca@dmi.units.it 2 Dima, University
More informationThe Expressivity of Universal Timed CCP: Undecidability of Monadic FLTL and Closure Operators for Security
The Expressivity of Universal Timed CCP: Undecidability of Monadic FLTL and Closure Operators for Security Carlos Olarte and Frank D. Valencia INRIA /CNRS and LIX, Ecole Polytechnique Motivation Concurrent
More informationSPA for quantitative analysis: Lecture 6 Modelling Biological Processes
1/ 223 SPA for quantitative analysis: Lecture 6 Modelling Biological Processes Jane Hillston LFCS, School of Informatics The University of Edinburgh Scotland 7th March 2013 Outline 2/ 223 1 Introduction
More informationVarieties of Stochastic Calculi
Research is what I'm doing when I don't know what I'm doing. Wernher Von Braun. Artificial Biochemistry Varieties of Stochastic Calculi Microsoft Research Trento, 26-5-22..26 www.luca.demon.co.uk/artificialbiochemistry.htm
More informationIntroduction Probabilistic Programming ProPPA Inference Results Conclusions. Embedding Machine Learning in Stochastic Process Algebra.
Embedding Machine Learning in Stochastic Process Algebra Jane Hillston Joint work with Anastasis Georgoulas and Guido Sanguinetti, School of Informatics, University of Edinburgh 16th August 2017 quan col....
More informationProbabilistic Model Checking for Biochemical Reaction Systems
Probabilistic Model Checing for Biochemical Reaction Systems Ratana Ty Ken-etsu Fujita Ken ichi Kawanishi Graduated School of Science and Technology Gunma University Abstract Probabilistic model checing
More informationSome investigations concerning the CTMC and the ODE model derived from Bio-PEPA
FBTC 2008 Some investigations concerning the CTMC and the ODE model derived from Bio-PEPA Federica Ciocchetta 1,a, Andrea Degasperi 2,b, Jane Hillston 3,a and Muffy Calder 4,b a Laboratory for Foundations
More informationModeling biological systems with delays in Bio-PEPA
Modeling biological systems with delays in Bio-PEPA Giulio Caravagna Dipartimento di Informatica, Università di Pisa, argo Pontecorvo 3, 56127 Pisa, Italy. caravagn@di.unipi.it Jane Hillston aboratory
More informationTimed Automata. Chapter Clocks and clock constraints Clock variables and clock constraints
Chapter 10 Timed Automata In the previous chapter, we have discussed a temporal logic where time was a discrete entities. A time unit was one application of the transition relation of an LTS. We could
More informationModels for Efficient Timed Verification
Models for Efficient Timed Verification François Laroussinie LSV / ENS de Cachan CNRS UMR 8643 Monterey Workshop - Composition of embedded systems Model checking System Properties Formalizing step? ϕ Model
More informationCSL model checking of biochemical networks with Interval Decision Diagrams
CSL model checking of biochemical networks with Interval Decision Diagrams Brandenburg University of Technology Cottbus Computer Science Department http://www-dssz.informatik.tu-cottbus.de/software/mc.html
More informationMASSPA-Modeller: A Spatial Stochastic Process Algebra modelling tool
MASSPA-Modeller: A Spatial Stochastic Process Algebra modelling tool Marcel C. Guenther Jeremy T. Bradley Imperial College London, 180 Queen s Gate, London SW7 2AZ, United Kingdom, Email: {mcg05,jb}@doc.ic.ac.uk
More informationPopulation models from PEPA descriptions
Population models from PEPA descriptions Jane Hillston LFCS, The University of Edinburgh, Edinburgh EH9 3JZ, Scotland. Email: jeh@inf.ed.ac.uk 1 Introduction Stochastic process algebras (e.g. PEPA [10],
More informationModelling Biochemical Pathways with Stochastic Process Algebra
Modelling Biochemical Pathways with Stochastic Process Algebra Jane Hillston. LFCS, University of Edinburgh 13th April 2007 The PEPA project The PEPA project started in Edinburgh in 1991. The PEPA project
More informationFrom CCS to Hybrid π via baby steps. Bill Rounds CSE, U of Michigan
From CCS to Hybrid π via baby steps Bill Rounds CSE, U of Michigan Main idea The hybrid pi-calculus extends pi-calculus by adding a component called the continuous environment, which evolves over time
More informationModelling protein trafficking: progress and challenges
Modelling protein trafficking: progress and challenges Laboratory for Foundations of Computer Science School of Informatics University of Edinburgh 21 May 2012 Outline Protein Trafficking Modelling Biology
More informationModelling in Systems Biology
Modelling in Systems Biology Maria Grazia Vigliotti thanks to my students Anton Stefanek, Ahmed Guecioueur Imperial College Formal representation of chemical reactions precise qualitative and quantitative
More informationRecent results on Timed Systems
Recent results on Timed Systems Time Petri Nets and Timed Automata Béatrice Bérard LAMSADE Université Paris-Dauphine & CNRS berard@lamsade.dauphine.fr Based on joint work with F. Cassez, S. Haddad, D.
More informationThe State Explosion Problem
The State Explosion Problem Martin Kot August 16, 2003 1 Introduction One from main approaches to checking correctness of a concurrent system are state space methods. They are suitable for automatic analysis
More information2 Conceptual Framework Before introducing the probabilistic concurrent constraint (PCCP) language we have to discuss a basic question: What is a proba
On Probabilistic CCP Alessandra Di Pierro and Herbert Wiklicky fadp,herbertg@cs.city.ac.uk City University London, Northampton Square, London EC1V OHB Abstract This paper investigates a probabilistic version
More informationProbabilistic Model Checking Michaelmas Term Dr. Dave Parker. Department of Computer Science University of Oxford
Probabilistic Model Checking Michaelmas Term 2011 Dr. Dave Parker Department of Computer Science University of Oxford Probabilistic model checking System Probabilistic model e.g. Markov chain Result 0.5
More informationCprE 281: Digital Logic
CprE 281: Digital Logic Instructor: Alexander Stoytchev http://www.ece.iastate.edu/~alexs/classes/ Synchronous Sequential Circuits Basic Design Steps CprE 281: Digital Logic Iowa State University, Ames,
More informationQualitative and Quantitative Analysis of a Bio-PEPA Model of the Gp130/JAK/STAT Signalling Pathway
Qualitative and Quantitative Analysis of a Bio-PEPA Model of the Gp13/JAK/STAT Signalling Pathway Maria Luisa Guerriero Laboratory for Foundations of Computer Science, The University of Edinburgh, UK Abstract.
More informationAutomata-based Verification - III
COMP30172: Advanced Algorithms Automata-based Verification - III Howard Barringer Room KB2.20: email: howard.barringer@manchester.ac.uk March 2009 Third Topic Infinite Word Automata Motivation Büchi Automata
More informationAutomata-based Verification - III
CS3172: Advanced Algorithms Automata-based Verification - III Howard Barringer Room KB2.20/22: email: howard.barringer@manchester.ac.uk March 2005 Third Topic Infinite Word Automata Motivation Büchi Automata
More informationA Behavioral Congruence for Concurrent Constraint Programming with Nondeterministic Choice
A Behavioral Congruence for Concurrent Constraint Programming with Nondeterministic Choice Luis Pino*, Filippo Bonchi** and Frank Valencia* (Presented by: Jorge A. Pe rez) *E quipe Come te, LIX, Laboratoire
More informationCoBiC: Context-dependent Bioambient Calculus
QAPL 2009 CoBiC: Context-dependent Bioambient Calculus Luca Bortolussi,3 Dipartimento di Matematica ed Informatica University of Trieste, Trieste, Italy Maria G. Vigliotti 4,2 Department of Computing,
More informationFormal Verification of Mobile Network Protocols
Dipartimento di Informatica, Università di Pisa, Italy milazzo@di.unipi.it Pisa April 26, 2005 Introduction Modelling Systems Specifications Examples Algorithms Introduction Design validation ensuring
More informationHarvard CS 121 and CSCI E-207 Lecture 9: Regular Languages Wrap-Up, Context-Free Grammars
Harvard CS 121 and CSCI E-207 Lecture 9: Regular Languages Wrap-Up, Context-Free Grammars Salil Vadhan October 2, 2012 Reading: Sipser, 2.1 (except Chomsky Normal Form). Algorithmic questions about regular
More informationFormal Models of Timed Musical Processes Doctoral Defense
Formal Models of Timed Musical Processes Doctoral Defense Gerardo M. Sarria M. Advisor: Camilo Rueda Co-Advisor: Juan Francisco Diaz Universidad del Valle AVISPA Research Group September 22, 2008 Motivation
More informationEVOLUTIONARY DISTANCES
EVOLUTIONARY DISTANCES FROM STRINGS TO TREES Luca Bortolussi 1 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste luca@dmi.units.it Trieste, 14 th November 2007 OUTLINE 1 STRINGS:
More informationTowards a Mechanised Denotational Semantics for Modelica
Towards a Mechanised Denotational Semantics for Modelica Simon Foster Bernhard Thiele Jim Woodcock Peter Fritzson Department of Computer Science, University of York PELAB, Linköping University 3rd February
More informationHybrid systems and computer science a short tutorial
Hybrid systems and computer science a short tutorial Eugene Asarin Université Paris 7 - LIAFA SFM 04 - RT, Bertinoro p. 1/4 Introductory equations Hybrid Systems = Discrete+Continuous SFM 04 - RT, Bertinoro
More informationProbabilistic verification and approximation schemes
Probabilistic verification and approximation schemes Richard Lassaigne Equipe de Logique mathématique, CNRS-Université Paris 7 Joint work with Sylvain Peyronnet (LRDE/EPITA & Equipe de Logique) Plan 1
More informationCASPA - A Tool for Symbolic Performance Evaluation and Stochastic Model Checking
CASPA - A Tool for Symbolic Performance Evaluation and Stochastic Model Checking Boudewijn R. Haverkort 1, Matthias Kuntz 1, Martin Riedl 2, Johann Schuster 2, Markus Siegle 2 1 : Universiteit Twente 2
More informationModelling a Circadian Clock with HYPE
Modelling a Circadian Clock with HYPE Vashti Galpin Laboratory for Foundations of Computer Scince School of Informatics, University of Edinburgh Vashti.Galpin@ed.ac.uk HYPE is a process algebra for hybrid
More informationExperiments with Measuring Time in PRISM 4.0
Experiments with Measuring Time in PRISM 4.0 Wolfgang Schreiner Research Institute for Symbolic Computation (RISC) Johannes Kepler University, Linz, Austria Wolfgang.Schreiner@risc.jku.at March 1, 2013
More informationOperational Semantics
Operational Semantics Semantics and applications to verification Xavier Rival École Normale Supérieure Xavier Rival Operational Semantics 1 / 50 Program of this first lecture Operational semantics Mathematical
More informationA Probabilistic Model for Molecular Systems
Fundamenta Informaticae XX (2005) 1 14 1 IOS Press A Probabilistic Model for Molecular Systems Roberto Barbuti Stefano Cataudella Andrea Maggiolo Schettini Paolo Milazzo C Angelo Troina Dipartimento di
More informationLecture 14: State Tables, Diagrams, Latches, and Flip Flop
EE210: Switching Systems Lecture 14: State Tables, Diagrams, Latches, and Flip Flop Prof. YingLi Tian Nov. 6, 2017 Department of Electrical Engineering The City College of New York The City University
More informationEmbedded Systems 5. Synchronous Composition. Lee/Seshia Section 6.2
Embedded Systems 5-1 - Synchronous Composition Lee/Seshia Section 6.2 Important semantic model for concurrent composition Here: composition of actors Foundation of Statecharts, Simulink, synchronous programming
More informationA Tableau-Based Decision Procedure for Right Propositional Neighborhood Logic (RPNL )
A Tableau-Based Decision Procedure for Right Propositional Neighborhood Logic (RPNL ) Davide Bresolin Angelo Montanari Dipartimento di Matematica e Informatica Università degli Studi di Udine {bresolin,
More informationEmbedded Systems Development
Embedded Systems Development Lecture 2 Finite Automata & SyncCharts Daniel Kästner AbsInt Angewandte Informatik GmbH kaestner@absint.com Some things I forgot to mention 2 Remember the HISPOS registration
More informationDES. 4. Petri Nets. Introduction. Different Classes of Petri Net. Petri net properties. Analysis of Petri net models
4. Petri Nets Introduction Different Classes of Petri Net Petri net properties Analysis of Petri net models 1 Petri Nets C.A Petri, TU Darmstadt, 1962 A mathematical and graphical modeling method. Describe
More informationHybrid Automata and ɛ-analysis on a Neural Oscillator
Hybrid Automata and ɛ-analysis on a Neural Oscillator A. Casagrande 1 T. Dreossi 2 C. Piazza 2 1 DMG, University of Trieste, Italy 2 DIMI, University of Udine, Italy Intuitively... Motivations: Reachability
More informationWhat You Must Remember When Processing Data Words
What You Must Remember When Processing Data Words Michael Benedikt, Clemens Ley, and Gabriele Puppis Oxford University Computing Laboratory, Park Rd, Oxford OX13QD UK Abstract. We provide a Myhill-Nerode-like
More informationNew Computational Methods for Systems Biology
New Computational Methods for Systems Biology François Fages, Sylvain Soliman The French National Institute for Research in Computer Science and Control INRIA Paris-Rocquencourt Constraint Programming
More informationUsing Probabilistic Model Checking in Systems Biology
Using Probabilistic Model Checking in Systems Biology Marta Kwiatkowska, Gethin Norman and David Parker Oxford University Computing Laboratory, Wolfson Building, Parks Road, Oxford, OX1 3QD ABSTRACT Probabilistic
More informationAbstracting real-valued parameters in parameterised boolean equation systems
Department of Mathematics and Computer Science Formal System Analysis Research Group Abstracting real-valued parameters in parameterised boolean equation systems Master Thesis M. Laveaux Supervisor: dr.
More informationUndecidable Problems. Z. Sawa (TU Ostrava) Introd. to Theoretical Computer Science May 12, / 65
Undecidable Problems Z. Sawa (TU Ostrava) Introd. to Theoretical Computer Science May 12, 2018 1/ 65 Algorithmically Solvable Problems Let us assume we have a problem P. If there is an algorithm solving
More informationAlternating Time Temporal Logics*
Alternating Time Temporal Logics* Sophie Pinchinat Visiting Research Fellow at RSISE Marie Curie Outgoing International Fellowship * @article{alur2002, title={alternating-time Temporal Logic}, author={alur,
More informationSFM-11:CONNECT Summer School, Bertinoro, June 2011
SFM-:CONNECT Summer School, Bertinoro, June 20 EU-FP7: CONNECT LSCITS/PSS VERIWARE Part 3 Markov decision processes Overview Lectures and 2: Introduction 2 Discrete-time Markov chains 3 Markov decision
More informationEvolutionary Computation
Evolutionary Computation - Computational procedures patterned after biological evolution. - Search procedure that probabilistically applies search operators to set of points in the search space. - Lamarck
More informationSymmetry Reductions. A. Prasad Sistla University Of Illinois at Chicago
Symmetry Reductions. A. Prasad Sistla University Of Illinois at Chicago Model-Checking Concurrent PGM Temporal SPEC Model Checker Yes/No Counter Example Approach Build the global state graph Algorithm
More informationPolynomial-Time Verification of PCTL Properties of MDPs with Convex Uncertainties and its Application to Cyber-Physical Systems
Polynomial-Time Verification of PCTL Properties of MDPs with Convex Uncertainties and its Application to Cyber-Physical Systems Alberto Puggelli DREAM Seminar - November 26, 2013 Collaborators and PIs:
More informationEfficient Computation of Program Equivalence for Confluent Concurrent Constraint Programming
Efficient Computation of Program Equivalence for Confluent Concurrent Constraint Programming Luis F. Pino INRIA/DGA and LIX École Polytechnique 91128 Palaiseau, France luis.pino@lix.polytechnique.fr Filippo
More informationEECS 144/244: Fundamental Algorithms for System Modeling, Analysis, and Optimization
EECS 144/244: Fundamental Algorithms for System Modeling, Analysis, and Optimization Discrete Systems Lecture: Automata, State machines, Circuits Stavros Tripakis University of California, Berkeley Stavros
More informationState Explosion in Almost-Sure Probabilistic Reachability
State Explosion in Almost-Sure Probabilistic Reachability François Laroussinie Lab. Spécification & Vérification, ENS de Cachan & CNRS UMR 8643, 61, av. Pdt. Wilson, 94235 Cachan Cedex France Jeremy Sproston
More informationReachability Results for Timed Automata with Unbounded Data Structures
Acta Informatica manuscript No. (will be inserted by the editor) Reachability Results for Timed Automata with Unbounded Data Structures Ruggero Lanotte Andrea Maggiolo-Schettini Angelo Troina Received:
More informationTimed Automata VINO 2011
Timed Automata VINO 2011 VeriDis Group - LORIA July 18, 2011 Content 1 Introduction 2 Timed Automata 3 Networks of timed automata Motivation Formalism for modeling and verification of real-time systems.
More informationCoBiC: Context-dependent Bioambient Calculus
Electronic Notes in Theoretical Computer Science 253 (2009) 87 20 www.elsevier.com/locate/entcs CoBiC: Context-dependent Bioambient Calculus Luca Bortolussi,3 Dipartimento di Matematica ed Informatica
More informationLet us first give some intuitive idea about a state of a system and state transitions before describing finite automata.
Finite Automata Automata (singular: automation) are a particularly simple, but useful, model of computation. They were initially proposed as a simple model for the behavior of neurons. The concept of a
More informationSymbolic Semantics and Verification of Stochastic Process Algebras. Symbolische Semantik und Verifikation stochastischer Prozessalgebren
Symbolic Semantics and Verification of Stochastic Process Algebras Symbolische Semantik und Verifikation stochastischer Prozessalgebren Der Technischen Fakultät der Universität Erlangen-Nürnberg zur Erlangung
More informationProofs, Strings, and Finite Automata. CS154 Chris Pollett Feb 5, 2007.
Proofs, Strings, and Finite Automata CS154 Chris Pollett Feb 5, 2007. Outline Proofs and Proof Strategies Strings Finding proofs Example: For every graph G, the sum of the degrees of all the nodes in G
More informationLecture 2: From Classical to Quantum Model of Computation
CS 880: Quantum Information Processing 9/7/10 Lecture : From Classical to Quantum Model of Computation Instructor: Dieter van Melkebeek Scribe: Tyson Williams Last class we introduced two models for deterministic
More informationProbabilistic Model Checking and Strategy Synthesis for Robot Navigation
Probabilistic Model Checking and Strategy Synthesis for Robot Navigation Dave Parker University of Birmingham (joint work with Bruno Lacerda, Nick Hawes) AIMS CDT, Oxford, May 2015 Overview Probabilistic
More informationRealization theory for systems biology
Realization theory for systems biology Mihály Petreczky CNRS Ecole Central Lille, France February 3, 2016 Outline Control theory and its role in biology. Realization problem: an abstract formulation. Realization
More informationAn Introduction to Hybrid Systems Modeling
CS620, IIT BOMBAY An Introduction to Hybrid Systems Modeling Ashutosh Trivedi Department of Computer Science and Engineering, IIT Bombay CS620: New Trends in IT: Modeling and Verification of Cyber-Physical
More informationA Graphical Representation for the Stochastic Pi-Calculus
A Graphical Representation for the Stochastic Pi-Calculus Andrew Phillips Luca Cardelli 27th August 2005 Microsoft Research 7 J J Thomson Avenue Cambridge, UK Andrew Phillips - Bioconcur 2005 Introduction
More informationProbabilistic Model Checking Michaelmas Term Dr. Dave Parker. Department of Computer Science University of Oxford
Probabilistic Model Checking Michaelmas Term 2011 Dr. Dave Parker Department of Computer Science University of Oxford Overview CSL model checking basic algorithm untimed properties time-bounded until the
More informationPRISM An overview. automatic verification of systems with stochastic behaviour e.g. due to unreliability, uncertainty, randomisation,
PRISM An overview PRISM is a probabilistic model checker automatic verification of systems with stochastic behaviour e.g. due to unreliability, uncertainty, randomisation, Construction/analysis of probabilistic
More informationLinear Temporal Logic and Büchi Automata
Linear Temporal Logic and Büchi Automata Yih-Kuen Tsay Department of Information Management National Taiwan University FLOLAC 2009 Yih-Kuen Tsay (SVVRL @ IM.NTU) Linear Temporal Logic and Büchi Automata
More informationAnalysis and Optimization of Discrete Event Systems using Petri Nets
Volume 113 No. 11 2017, 1 10 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Analysis and Optimization of Discrete Event Systems using Petri Nets
More informationVerification of Polynomial Interrupt Timed Automata
Verification of Polynomial Interrupt Timed Automata Béatrice Bérard 1, Serge Haddad 2, Claudine Picaronny 2, Mohab Safey El Din 1, Mathieu Sassolas 3 1 Université P. & M. Curie, LIP6 2 ENS Cachan, LSV
More informationThe PITA System for Logical-Probabilistic Inference
The System for Logical-Probabilistic Inference Fabrizio Riguzzi 1 and Terrance Swift 2 1 EDIF University of Ferrara, Via Saragat 1, I-44122, Ferrara, Italy fabrizio.riguzzi@unife.it 2 CETRIA Universidade
More informationA framework for simulation and symbolic state space analysis of non-markovian models
A framework for simulation and symbolic state space analysis of non-markovian models Laura Carnevali, Lorenzo Ridi, Enrico Vicario SW Technologies Lab (STLab) - Dip. Sistemi e Informatica (DSI) - Univ.
More information86 Part 4 SUMMARY INTRODUCTION
86 Part 4 Chapter # AN INTEGRATION OF THE DESCRIPTIONS OF GENE NETWORKS AND THEIR MODELS PRESENTED IN SIGMOID (CELLERATOR) AND GENENET Podkolodny N.L. *1, 2, Podkolodnaya N.N. 1, Miginsky D.S. 1, Poplavsky
More informationBiological Pathways Representation by Petri Nets and extension
Biological Pathways Representation by and extensions December 6, 2006 Biological Pathways Representation by and extension 1 The cell Pathways 2 Definitions 3 4 Biological Pathways Representation by and
More informationDistributed Optimization. Song Chong EE, KAIST
Distributed Optimization Song Chong EE, KAIST songchong@kaist.edu Dynamic Programming for Path Planning A path-planning problem consists of a weighted directed graph with a set of n nodes N, directed links
More informationBounded Model Checking with SAT/SMT. Edmund M. Clarke School of Computer Science Carnegie Mellon University 1/39
Bounded Model Checking with SAT/SMT Edmund M. Clarke School of Computer Science Carnegie Mellon University 1/39 Recap: Symbolic Model Checking with BDDs Method used by most industrial strength model checkers:
More informationComputational Methods in Systems and Synthetic Biology
Computational Methods in Systems and Synthetic Biology François Fages EPI Contraintes, Inria Paris-Rocquencourt, France 1 / 62 Need for Abstractions in Systems Biology Models are built in Systems Biology
More informationConcept of Statistical Model Checking. Presented By: Diana EL RABIH
Concept of Statistical Model Checking Presented By: Diana EL RABIH Introduction Model checking of stochastic systems Continuous-time Markov chains Continuous Stochastic Logic (CSL) Probabilistic time-bounded
More informationPlanning Under Uncertainty II
Planning Under Uncertainty II Intelligent Robotics 2014/15 Bruno Lacerda Announcement No class next Monday - 17/11/2014 2 Previous Lecture Approach to cope with uncertainty on outcome of actions Markov
More informationComputation Tree Logic (CTL) & Basic Model Checking Algorithms
Computation Tree Logic (CTL) & Basic Model Checking Algorithms Martin Fränzle Carl von Ossietzky Universität Dpt. of Computing Science Res. Grp. Hybride Systeme Oldenburg, Germany 02917: CTL & Model Checking
More informationOn simulations and bisimulations of general flow systems
On simulations and bisimulations of general flow systems Jen Davoren Department of Electrical & Electronic Engineering The University of Melbourne, AUSTRALIA and Paulo Tabuada Department of Electrical
More informationProbabilistic model checking with PRISM
Probabilistic model checking with PRISM Marta Kwiatkowska Department of Computer Science, University of Oxford 4th SSFT, Menlo College, May 204 Part 2 Markov decision processes Overview (Part 2) Introduction
More informationHybrid Automata. Lecturer: Tiziano Villa 1. Università di Verona
Hybrid Automata Lecturer: Tiziano Villa 1 1 Dipartimento d Informatica Università di Verona tiziano.villa@univr.it Thanks to Carla Piazza, Dipartimento di Matematica ed Informatica, Università di Udine
More informationDiscrete abstractions of hybrid systems for verification
Discrete abstractions of hybrid systems for verification George J. Pappas Departments of ESE and CIS University of Pennsylvania pappasg@ee.upenn.edu http://www.seas.upenn.edu/~pappasg DISC Summer School
More informationProbabilistic model checking A comparison of tools
Probabilistic model checking A comparison of tools Masters Thesis in Computer Science of H.A. Oldenkamp born on the 23th of September 1979 May, 2007 Chairman: Prof. Dr. Ir. Joost-Pieter Katoen Supervisors:
More informationThe algorithmic analysis of hybrid system
The algorithmic analysis of hybrid system Authors: R.Alur, C. Courcoubetis etc. Course teacher: Prof. Ugo Buy Xin Li, Huiyong Xiao Nov. 13, 2002 Summary What s a hybrid system? Definition of Hybrid Automaton
More informationFormal Verification. Lecture 1: Introduction to Model Checking and Temporal Logic¹
Formal Verification Lecture 1: Introduction to Model Checking and Temporal Logic¹ Jacques Fleuriot jdf@inf.ed.ac.uk ¹Acknowledgement: Adapted from original material by Paul Jackson, including some additions
More informationUmans Complexity Theory Lectures
Complexity Theory Umans Complexity Theory Lectures Lecture 1a: Problems and Languages Classify problems according to the computational resources required running time storage space parallelism randomness
More informationAnalog Electronics Mimic Genetic Biochemical Reactions in Living Cells
Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells Dr. Ramez Daniel Laboratory of Synthetic Biology & Bioelectronics (LSB 2 ) Biomedical Engineering, Technion May 9, 2016 Cytomorphic
More information