Stochastic analysis of biochemical reaction networks with absolute concentration robustness
|
|
- Miranda Phoebe Chapman
- 5 years ago
- Views:
Transcription
1 Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson Department of Mathematics University of Wisconsin - Madison Boston University Probability/Statistics seminar Nov. 21st, 2013
2 Big picture Biochemical/population networks can range from simple to very complex. Example 1: A + B C. Example 2: A + B 2B B A Example 3: Gene transcription & translation: G κ 1 G + M M κ 2 M + P M κ 3 P κ 4 G + P κ 5 κ 5 B transcription translation degradation degradation Binding/unbinding of Gene Cartoon representation: 1 1 J. Paulsson, Physics of Life Reviews, 2,
3 Big picture Example 4: EnvZ/OmpR signaling system 2 2 Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, 2010
4 Big picture Hanahan and Weinberg, The Hallmarks of Cancer, Cell, 2000.!
5 Big picture Metabolic Pathways Roche Applied Science!
6 Big picture For complex models, simulation is often used to explore the possible dynamics (for both deterministic and stochastic models). One problem: key system parameters are oftentimes unknown, or known only up to an order of magnitude. Want an alternative approach: discover what pieces of the network architecture determine overall system behavior. Best if results do not depend upon particular choices of rate constants (which are never really known). how can we cut through the bewildering complexity of biological systems and say something sensible? Mathematics. This story is far from complete really just beginning.
7 Big picture Will tell two stories: 1. ( A., Craciun, Kurtz) Provide network conditions that guarantee both an especially stable deterministic model and an especially stable stochastic model. 2. ( A., Enciso, Johnston) Provide network conditions that guarantee deterministic model has component which is Absolutely robust, stochastic model has extinction event! This research is part of chemical reaction network theory, which is part of systems biology. A network property called deficiency will play a key role in the analysis.
8 Networks: {S, C, R} Example: A B S = {A, B}. C = {A, B}. R = {A B}. Example: A + B 2B B A S = {A, B}. C = {A + B, 2B, B, A}. R = {A + B 2B, B A}.
9 Networks: {S, C, R} Example: E + S ES E + P, E S S = {E, S, ES, P}. C = {E + S, ES, E + P, E,, S}. R = {E + S ES, ES E + S,... }.
10 Dynamics: stochastic Example: A + B α 2B B β A (R1) (R2) X(t) = ([ 0 X(0) + R 1 (t) 2 = [ 1 X(0) + R 1 (t) 1 For Markov models can take R 1 (t) = Y 1 (α R 2 (t) = Y 2 (β ] [ 1 1 ] + R 2 (t) t 0 t 0 ]) ([ 1 + R 2 (t) 0 [ 1 1 ]. ) X A (s)x B (s)ds ) X B (s)ds where Y 1, Y 2 are independent unit-rate Poisson processes. ] [ 0 1 ])
11 Dynamics: stochastic For general system, we have S = {X 1,..., X d }, with R : d y ki X i i=1 κ k d y kix i or y k y k i=1 and complexes found at either end of reaction arrow: C = {y k, y k}. The intensity (or propensity) of kth reaction is λ k : Z d 0 R. By analogy with before: X(t) = X(0) + k R k (t)(y k y k ), with X(t) = X(0) + k Y k ( t Y k are independent, unit-rate Poisson processes. 0 ) λ k (X(s))ds (y k y k ),
12 Mass-action kinetics The standard intensity function chosen is mass-action kinetics: λ k (x) = κ k ( ( ) x x i! y ik!) = κ k y k (x i y ik )!. i i Example: If S 1 anything, then λ k (x) = κ k x 1. Example: If S 1 + S 2 anything, then λ k (x) = κ k x 1 x 2. Example: If 2S 2 anything, then λ k (x) = κ k x 2 (x 2 1).
13 Other ways to understand model We can also describe the model as a continuous time Markov chain with infinitesimal generator Af (x) = k λ k (x)(f (x + ζ k ) f (x)). where ζ k = y k y k. Dynkin s formula (See Ethier and Kurtz, 1986, Ch. 1) yields Ef (X(t)) f (X 0 ) = E t 0 (Af )(X(s))ds, Letting f (y) = 1 x(y) above so that E[f (X(t))] = P{X(t) = x} = p x(t), gives Kolmogorov s forward equation (chemical master equation) p t (x) = k λ(x ζ k )p t(x ζ k ) p t(x) k λ k (x)
14 Morris-Lecar Voltage satisfying dv dt = f (V (t), N(t)) = 1 C ( Iapp g Ca m (V (t))(v (t) V Ca ) V K ) g L (V V L ) g o K N(t)(V (t) ). with number of open potassium channels following Open α(v (t)) Closed β(v (t))
15 Morris-Lecar Voltage satisfying dv dt = 1 C ( Iapp g Ca m (V (t))(v (t) V Ca ) V K ) g L (V V L ) g o K N(t)(V (t) ). with number of open potassium channels following Open α(v (t)) Closed β(v (t)) N(t) = N(0) Y close ( t 0 ) ( t ) β(v (s))n(s)ds + Y open α(v (s))(n tot N(s))ds. 0
16 Dynamics: deterministic model Assuming V is a scaling parameter (volume times Avogadro s number), X i = O(V ), and X V (t) = def X(t)/V, λ k (X(t)) V ( κ k X V (t) y k ), Then, X V (t) 1 V X 0 + k LLN for Y k says 1 V Y k(vu) u so a good approximation is solution to ( 1 t ) V Y k V κ k X V (s) y k ds (y k y k ) 0 x(t) = x(0) + k t 0 κ k x(s) y k ds (y k y k ), where is standard mass-action kinetics. u v = u v 1 1 uv d d,
17 Dynamics. Example A + B α 2B B β A (R1) (R2) Stochastic equations X(t) = X(0) + Y 1 (α t 0 ) [ 1 X A (s)x B (s)ds 1 ] ( t ) [ 1 + Y 2 β X B (s)ds 0 1 ]. Deterministic equations or or x(t) = x(0) + α t 0 [ 1 x A (s)x B (s)ds 1 ẋ(t) = αx ax b [ 1 1 ] t [ 1 + β x B (s)ds 0 1 ] + βx B [ 1 1 ẋ A (t) = αx A x B + βx B ẋ B (t) = αx ax B βx B. ] ],
18 Story 1: special stability for both models Network reversibility conditions Definition Let S = {S i }, C = {y k }, and R = {y k y k} denote the sets of species, complexes, and reactions, respectively. The triple {S, C, R} is called the chemical reaction network. Definition Let G be the directed graph with nodes given by the complexes C and directed edges given by the reactions R = {y k y k}, and let G 1,..., G l denote the connected components of G. {G j } are the linkage classes of the reaction network. Example A + B α 2B B β A (R1) (R2) Has two linkage classes.
19 Network reversibility conditions Definition A chemical reaction network, {S, C, R}, is called weakly reversible if each linkage class is strongly connected. A network is called reversible if y k y k R whenever y k y k R. Intuition for probabilists: If the network is weakly reversible, then, thinking of the complexes as states of a Markov chain, the linkage classes are the irreducible communicating equivalence classes of classical Markov chain theory. BUT, these equivalence classes do not correspond to the communicating equivalence classes of the Markov chain model of the reaction network.
20 Network properties Definition S = span {yk y k R}{y k y k } is the stoichiometric subspace of the network. For c R d we say c + S and (c + S) R d >0 are the stoichiometric compatibility classes and positive stoichiometric compatibility classes of the network, respectively. Denote dim(s) = s. Implication: Solutions bound to stoichiometric compatibility class. 1. Deterministic: x(t) x(0) c + S for all t Stochastic: communication classes contained within c + S.
21 Network properties Definition The deficiency of a chemical reaction network, {S, C, R}, is δ = C l s, where C is the number of complexes, l is the number of linkage classes of the network graph, and s is the dimension of the stoichiometric subspace of the network. Lemma (Feinberg) The deficiency of a network is nonnegative. Example A + B 2B B A (R1) (R2) n = 4, l = 2, s = 1 = δ = 1. But, A + B C B A (R1) (R2) n = 4, l = 2, s = 2 = δ = 0.
22 Deficiency zero theorem of Feinberg Theorem (The Deficiency Zero Theorem - Deterministic) Let {S, C, R} be a weakly reversible, deficiency zero chemical reaction network with mass action kinetics. Then, for any choice of rate constants κ k, within each positive stoichiometric compatibility class there is precisely one equilibrium value c, κ k c y k (y k y k ) = 0, k and that equilibrium value is locally asymptotically stable relative to its compatibility class. Actually have stronger result: for each η C, κ k c y k = k:y k =η c is said to be a complex balanced equilibrium. k:y k =η κ k c yk. (1)
23 Deficiency Zero Theorem - stochastic Theorem (A., Craciun, Kurtz, 2010) Let {S, C, R} be a chemical reaction network with rate constants κ k. Suppose that the deterministic system is complex balanced with equilibrium c R m >0. Then, for any irreducible communicating equivalence class, Γ, the stochastic system has a product form stationary measure π(x) = M where M is a normalizing constant. m i=1 c x i i, x Γ, (2) x i! Useful for stochastic averaging.
24 Enzyme kinetics Consider the possible enzyme kinetics given by E + S ES E + P, E S Not too difficult to see that Γ = Z 4 0. Therefore, in distributional equilibrium the specie numbers are independent and have Poisson distributions.
25 Non-independence Consider the simple reversible system S 1 κ 1 κ 2 S 2. Suppose that X 1 (0) + X 2 (0) = N = X 1 (t) + X 2 (t) = N for all t. δ = n l s = = 0. An equilibrium to the system that satisfies the complex balance equation is ( ) κ2 κ 1 c =,, κ 1 + κ 2 κ 1 + κ 2 and the product form stationary distribution for the system is π(x) = M cx 1 1 c x 2 2 x 1! x 2!, x Γ.
26 Non-independence Using that X 1 (t) + X 2 (t) = N for all t yields π 1 (x 1 ) = M cx 1 1 c N x 1 2 x 1! (N x 1 )! = M x 1!(N x 1 )! cx 1 1 (1 c 1) N x 1. After setting M = N!, we see that X 1 is binomially distributed. Similarly, π 2 (x 2 ) = ( N x 2 ) c x 2 2 (1 c 2) N x 2.
27 Enzyme kinetics Consider the slightly different enzyme kinetics given by E + S ES E + P, E 0 We see S + ES + P = N. In distributional equilibrium E has Poisson distribution, S, ES, P have a multinomial distribution, and E is independent from S, ES, and P.
28 Story 2: Absolute Concentration Robustness Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, AIM:
29 Story 2: Absolute Concentration Robustness A + B α 2B B β A (R1) (R2) ẋ A (t) = αx A (t)x B (t) + βx B (t) ẋ B (t) = αx A (t)x B (t) βx B (t) M = def x A (0) + x B (0), Solving for equilibria: x A = β/α, x B = M β/α, Network has absolute concentration robustness in species A.
30 X tot := X + XD + XT + X p + X py + XDY p Y tot := Y + X py + XDY p + Y p.
31 Differing in one species We say that two complexes, y 1, y 2, differ in species S i if for some α > 0. y 1 = y 2 + αe i Examples: 1. differ in species B. A, A + B 2. differ in species Y p. XT, XT + Y p
32 Terminal and non-terminal complexes Viewing complexes as states of a Markov chain, we say two complexes are strongly linked if they are part of the same communication class (of this Markov chain). the communication classes are called strong linkage classes. A strong linkage class is called terminal if the communication class is recurrent. Transient complexes are called non-terminal. k 1 XD X XT X p k 2 [ D] k 6 k 3 k4 X p +Y X p Y X+Y p k 7 k 8 k 9 [ T] XD+Y p XDY p XD+Y k 10 k 11 k 5
33 Conservative system We say there is a global conservation relation if there is a w R d >0 for which w (y k y k ) = 0. Example A + B 2B B A. has a conservation relation of w = [ 1 1 ]. since [ 1 1 [ 1 1 ] [ 1 1 ] [ 1 1 ] = 0 ] = 0.
34 Theorems Theorem (Marty Feinberg and Guy Shinar, Science, 2010 deterministic) Consider a deterministic mass-action system that has a deficiency of one. admits a positive steady state and has two non-terminal complexes that differ only in species S, then the system has absolute concentration robustness in S. Theorem (Anderson, Enciso, Johnston stochastic) Consider a reaction network satisfying the following: has a deficiency of one, the deterministic model admits a positive steady state, has two non-terminal complexes that differ only in species S, (new) is conservative, then for every positive recurrent state of the Markov model the intensity of each non-terminal complex is zero.
35 Do we need a conservation relation? - yes Example A + B α B β A + 2B, have n = 4, l = 2, s = 1 = δ = 1. No conservation relation (but x B x A is conserved), stable fixed point for ODE (if x B (0) > x A (0)) can never run out of B. x A = β α, x B = M + β α,
36 Flavor of proof Consider ẋ = κ k x y k (y k y k ) = Y A κ Ψ(x) k:y k y k Consider the vector space on R C with basis ω y = 1 in yth component. Define Ψ : R d 0 R C by Ψ(x) y = x y = d i=1 x y i i. Consider the mapping A κ : R C R C given by A κ(ψ) = κ k (ω y k ω yk ) Ψ y, y C k:y k =y so that A κ Ψ(x) = y C κ k x y k (ω y k ω yk ) k:y k =y Now define Y : R C R d to be linear mapping with Y (ω y) = y = Y (A κ(ψ(x))) = ẋ
37 Flavor of proof ẋ = f (x) = Y A κ Ψ(x). Note that x is an equilibrium value if Ψ(x) ker(a κ) or Ψ(x) ker(y A κ). So what is A κ? A quick (and easy) calculation shows { A ij κ j i for j i κ = l j κ j l for i = j A T κ is the infinitesimal generator for the complex space CTMC! We know all about the kernel (they correspond to stationary distributions).
38 Flavor of proof ẋ = f (x) = Y A κ Ψ(x). Lemma Let (S, C, R) denote a chemical reaction network with terminal strong linkage classes Λ 1,..., Λ t. Then, for any κ R r >0, ker(a κ) has a basis {b 1,..., b t} where supp(b θ ) = Λ θ for all θ = 1,..., t. Lemma Let (S, C, R) denote a chemical reaction network with deficiency δ and t terminal strong linkage classes. Then dim(ker(ya κ)) δ + t. (3) assumption equilibrium for the ODE model tells ker(y A κ). Not weakly reversible, so can t be linear combination of b i. Deficiency is one, so second lemma tells us that the kernel looks like where c charges all complexes. {c, b 1,..., b t},
39 Flavor of proof Recap: ker(y A κ) has basis {c, b 1,..., b t}, where c charges all complexes supp(b θ ) = Λ θ for all θ = 1,..., t. Rest of proof: consider a modified network with stochastic kinetics, prove that that system satisfies conclusion of theorem by using kernel condition, finally argue the original system does too.
40 An example A + B α 2B B β A (R1) (R2)
41 An example A + B α 2B B β A 2A B + C. (R1) (R2) (R3) Note C = 6, l = 3, s = 2 = δ = 1. w = (1, 1, 1) is conservation relation.
42 flavor of proof Only work with A + B reduced network. Call it {S, C, R }. Call process X. Suppose there is a positive recurrent state X 0 which charges a non-terminal complex. We want to reach a contradiction (with kernel condition). We denote the Nth return time to X 0 as t N. It follows that we have X (t N ) = X (0) for all N so that X (t N ) = X (0) + k ( tn ) Y k κ i (X (s)) y k ds (y k y k ) = X (0), (4) 0 and so ( tn ) Y k κ i (X (s)) y k ds (y k y k ) = 0, (5) k 0 Relatively easy to show that for each y k tn 0 (X (s)) y k ds, as N,
43 flavor of proof and ( tn ) Y k κ i (X (s)) y k ds (y k y k ) = 0, k 0 tn 0 (X (s)) y k ds, as N, multiplying and dividing by appropriate terms, 0 = [ ( 1 tn ) κ k tn Y k κ k (X (s)) y k ds 1 ] tn (X (s)) y k ds (y k κ k 0 (X (s)) y k k ds 0 t y k ). N 0 We have 1 lim N t N tn 0 (X (s)) y k ds = X I π I(X)X y k, (6) Define the vector G π R C via [G π] y := X I π I(X)X y, (7) for y C, and note that [G π] y > 0 for all remaining y. Can then argue that [G π] is in the kernel of Y A κ, which it can not be (contradiction).
44 Generalization Theorem Let (S, C, R) be a conservative chemical reaction network for which the following assumptions hold: 1. There are non-empty sets of non-terminal complexes C and C such that, if y C then there exists a y C such that y y. 2. For some choice of rate constants {k i } r i=1, the following property holds: if Ψ ker(ya k ) and Ψ y = 0 for all y C, then Ψ ȳ = 0 for all non-terminal ȳ. Then, for any choice of stoichiometrically admissible kinetics, all non-terminal complexes of the network are off at each positive recurrent state of the stochastically modeled system. Remark: The second condition in assumption 2 above may also be formulated in the following way, which may be more intuitive to some readers. 2. For some choice of rate constants {k i } r i=1, the following property holds: if Ψ ker(ya k ) has support on a non-terminal complex, then Ψ has support on some y C.
45 Generalization, example XD k 1 X k 3[T ] XT k 5 X p k 2 [D] k 4 X p + Y k 6 X py k 8 X + Y p k 7 k 9 k XD + Y p XDY p 11 XD + Y k 10 (8) Deficiency is two. k 12 k XT + Y p XTY p 14 XT + Yp. k 13 XD + Y p XD, and XT + Y p XT. Ordering non-terminal complexes as {XD, X, XT, X p + Y, X py, XD + Y p, XDY p, XT + Y p, XTY p}. (9) there are two basis vectors of YA k with non-zero support on the complexes as ordered in (9); {[2, 2, 1, 2, 1, 2, 1, 0, 0], [2, 2, 1, 2, 1, 0, 0, 2, 1]}.
46 Quasi-stationary distribution A + B α 2B B β A X A (0) + X B (0) = M, Find π Q M so that for τ absorption time and x transient states, Satisfies lim Pν(X(t) = x τ > t) = t πq M(x). πm(x) Q = P π Q (X(t) = x τ > t). M Can show that quasi-stationary distribution for A converges to Poisson πm(x) Q e (β/α) (β/α) x, as M. x!
47 Quasi-stationary distribution: EnvZ-OmpR signaling system 3 3 Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, 2010
48 Quasi-stationary distribution Quasi stationary probabilities X tot = 100 Y tot = 3500 X tot = 1000 Y tot = X tot = Y tot = Poisson Molecules of Y p
49 That is the story. Thanks! Collaborators: Matthew Johnston, Germán Enciso, Gheorghe Craciun, and Tom Kurtz References: 1. Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, David F. Anderson, Germán Enciso, and Matthew Johnston, Stochastic analysis of biochemical reaction networks with absolute concentration robustness, submitted. 3. David F. Anderson, Gheorghe Craciun, Thomas G. Kurtz, Product-form stationary distributions for deficiency zero chemical reaction networks, Bulletin of Mathematical Biology, Vol. 72, No. 8, , Funding: NSF-DMS , and NSF-DMS
Stochastic and deterministic models of biochemical reaction networks
Stochastic and deterministic models of biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Case Western Reserve November 9th, 215 Goals of the talk
More informationChemical reaction network theory for stochastic and deterministic models of biochemical reaction systems
Chemical reaction network theory for stochastic and deterministic models of biochemical reaction systems University of Wisconsin at Madison anderson@math.wisc.edu MBI Workshop for Young Researchers in
More informationLongtime behavior of stochastically modeled biochemical reaction networks
Longtime behavior of stochastically modeled biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison ASU Math Biology Seminar February 17th, 217 Overview
More informationPersistence and Stationary Distributions of Biochemical Reaction Networks
Persistence and Stationary Distributions of Biochemical Reaction Networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Discrete Models in Systems Biology SAMSI December
More informationStochastic models of biochemical systems
Stochastic models of biochemical systems David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison University of Amsterdam November 14th, 212 Stochastic models
More informationSimulation methods for stochastic models in chemistry
Simulation methods for stochastic models in chemistry David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison SIAM: Barcelona June 4th, 21 Overview 1. Notation
More informationCorrespondence of regular and generalized mass action systems
Correspondence of regular and generalized mass action systems Van Vleck Visiting Assistant Professor University of Wisconsin-Madison Joint Mathematics Meetings (San Antonio, TX) Saturday, January 10, 2014
More informationDiscrepancies between extinction events and boundary equilibria in reaction networks
Discrepancies between extinction events and boundary equilibria in reaction networks David F. nderson Daniele Cappelletti September 2, 208 bstract Reaction networks are mathematical models of interacting
More informationNetwork Analysis of Biochemical Reactions in Complex Environments
1 Introduction 1 Network Analysis of Biochemical Reactions in Complex Environments Elias August 1 and Mauricio Barahona, Department of Bioengineering, Imperial College London, South Kensington Campus,
More informationComputational methods for continuous time Markov chains with applications to biological processes
Computational methods for continuous time Markov chains with applications to biological processes David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Penn.
More informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson, Germán A. Enciso, Matthew D. Johnston October 13, 2013 Abstract It has recently been shown
More informationarxiv: v2 [math.ds] 12 Jul 2011
1 CHEMICAL REACTION SYSTEMS WITH TORIC STEADY STATES MERCEDES PÉREZ MILLÁN, ALICIA DICKENSTEIN, ANNE SHIU, AND CARSTEN CONRADI arxiv:11021590v2 [mathds] 12 Jul 2011 Abstract Mass-action chemical reaction
More information2008 Hotelling Lectures
First Prev Next Go To Go Back Full Screen Close Quit 1 28 Hotelling Lectures 1. Stochastic models for chemical reactions 2. Identifying separated time scales in stochastic models of reaction networks 3.
More informationMathematical and computational methods for understanding the dynamics of biochemical networks
Mathematical and computational methods for understanding the dynamics of biochemical networks Gheorghe Craciun Department of Mathematics and Department of Biomolecular Chemistry University of Wisconsin
More informationA Deficiency-Based Approach to Parametrizing Positive Equilibria of Biochemical Reaction Systems
Bulletin of Mathematical Biology https://doi.org/10.1007/s11538-018-00562-0 A Deficiency-Based Approach to Parametrizing Positive Equilibria of Biochemical Reaction Systems Matthew D. Johnston 1 Stefan
More informationarxiv: v2 [math.na] 21 Jul 2014
A Computational Approach to Steady State Correspondence of Regular and Generalized Mass Action Systems Matthew D. Johnston arxiv:1407.4796v2 [math.na] 21 Jul 2014 Contents Department of Mathematics University
More informationLyapunov functions, stationary distributions, and non-equilibrium potential for chemical reaction networks
Lyapunov functions, stationary distributions, and non-equilibrium potential for chemical reaction networks David F. Anderson, Gheorghe Craciun, Manoj Gopalkrishnan, Carsten Wiuf October 7, 04 Abstract
More informationNotes on Chemical Reaction Networks
Notes on Chemical Reaction Networks JWR May 19, 2009, last update: Nov 27, 2009 These notes are shamelessly cribbed from the extremely helpful online notes [2]. Even most of the notation is the same. I
More informationStatistics 150: Spring 2007
Statistics 150: Spring 2007 April 23, 2008 0-1 1 Limiting Probabilities If the discrete-time Markov chain with transition probabilities p ij is irreducible and positive recurrent; then the limiting probabilities
More informationStochastic Processes (Week 6)
Stochastic Processes (Week 6) October 30th, 2014 1 Discrete-time Finite Markov Chains 2 Countable Markov Chains 3 Continuous-Time Markov Chains 3.1 Poisson Process 3.2 Finite State Space 3.2.1 Kolmogrov
More information8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains
8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains 8.1 Review 8.2 Statistical Equilibrium 8.3 Two-State Markov Chain 8.4 Existence of P ( ) 8.5 Classification of States
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More informationTHE DYNAMICS OF WEAKLY REVERSIBLE POPULATION PROCESSES NEAR FACETS
SIAM J. APPL. MATH. Vol. 70, No. 6, pp. 1840 1858 c 2010 Society for Industrial and Applied Mathematics THE DYNAMICS OF WEAKLY REVERSIBLE POPULATION PROCESSES NEAR FACETS DAVID F. ANDERSON AND ANNE SHIU
More informationThe Transition Probability Function P ij (t)
The Transition Probability Function P ij (t) Consider a continuous time Markov chain {X(t), t 0}. We are interested in the probability that in t time units the process will be in state j, given that it
More informationMARKOV PROCESSES. Valerio Di Valerio
MARKOV PROCESSES Valerio Di Valerio Stochastic Process Definition: a stochastic process is a collection of random variables {X(t)} indexed by time t T Each X(t) X is a random variable that satisfy some
More informationContinuous-Time Markov Chain
Continuous-Time Markov Chain Consider the process {X(t),t 0} with state space {0, 1, 2,...}. The process {X(t),t 0} is a continuous-time Markov chain if for all s, t 0 and nonnegative integers i, j, x(u),
More informationChapter 5. Continuous-Time Markov Chains. Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan
Chapter 5. Continuous-Time Markov Chains Prof. Shun-Ren Yang Department of Computer Science, National Tsing Hua University, Taiwan Continuous-Time Markov Chains Consider a continuous-time stochastic process
More informationCS145: Probability & Computing Lecture 18: Discrete Markov Chains, Equilibrium Distributions
CS145: Probability & Computing Lecture 18: Discrete Markov Chains, Equilibrium Distributions Instructor: Erik Sudderth Brown University Computer Science April 14, 215 Review: Discrete Markov Chains Some
More informationNetwork Theory II: Stochastic Petri Nets, Chemical Reaction Networks and Feynman Diagrams. John Baez, Jacob Biamonte, Brendan Fong
Network Theory II: Stochastic Petri Nets, Chemical Reaction Networks and Feynman Diagrams John Baez, Jacob Biamonte, Brendan Fong A Petri net is a way of drawing a finite set S of species, a finite set
More informationLectures on Markov Chains
Lectures on Markov Chains David M. McClendon Department of Mathematics Ferris State University 2016 edition 1 Contents Contents 2 1 Markov chains 4 1.1 The definition of a Markov chain.....................
More informationStochastic Processes
MTHE/STAT455, STAT855 Fall 202 Stochastic Processes Final Exam, Solutions (5 marks) (a) (0 marks) Condition on the first pair to bond; each of the n adjacent pairs is equally likely to bond Given that
More informationPart I Stochastic variables and Markov chains
Part I Stochastic variables and Markov chains Random variables describe the behaviour of a phenomenon independent of any specific sample space Distribution function (cdf, cumulative distribution function)
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationFock Space Techniques for Stochastic Physics. John Baez, Jacob Biamonte, Brendan Fong
Fock Space Techniques for Stochastic Physics John Baez, Jacob Biamonte, Brendan Fong A Petri net is a way of drawing a finite set S of species, a finite set T of transitions, and maps s, t : T N S saying
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
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 informationKey words. persistence, permanence, global attractor conjecture, mass-action kinetics, powerlaw systems, biochemical networks, interaction networks
PERSISTENCE AND PERMANENCE OF MASS-ACTION AND POWER-LAW DYNAMICAL SYSTEMS GHEORGHE CRACIUN, FEDOR NAZAROV, AND CASIAN PANTEA Abstract. Persistence and permanence are properties of dynamical systems that
More informationKey words. persistence, permanence, global attractor conjecture, mass-action kinetics, powerlaw systems, biochemical networks, interaction networks
PERSISTENCE AND PERMANENCE OF MASS-ACTION AND POWER-LAW DYNAMICAL SYSTEMS GHEORGHE CRACIUN, FEDOR NAZAROV, AND CASIAN PANTEA Abstract Persistence and permanence are properties of dynamical systems that
More informationc 2018 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 78, No. 5, pp. 2692 2713 c 2018 Society for Industrial and Applied Mathematics SOME NETWORK CONDITIONS FOR POSITIVE RECURRENCE OF STOCHASTICALLY MODELED REACTION NETWORKS DAVID
More informationIrreducibility. Irreducible. every state can be reached from every other state For any i,j, exist an m 0, such that. Absorbing state: p jj =1
Irreducibility Irreducible every state can be reached from every other state For any i,j, exist an m 0, such that i,j are communicate, if the above condition is valid Irreducible: all states are communicate
More informationfor all f satisfying E[ f(x) ] <.
. Let (Ω, F, P ) be a probability space and D be a sub-σ-algebra of F. An (H, H)-valued random variable X is independent of D if and only if P ({X Γ} D) = P {X Γ}P (D) for all Γ H and D D. Prove that if
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationSMSTC (2007/08) Probability.
SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................
More informationQuantitative Model Checking (QMC) - SS12
Quantitative Model Checking (QMC) - SS12 Lecture 06 David Spieler Saarland University, Germany June 4, 2012 1 / 34 Deciding Bisimulations 2 / 34 Partition Refinement Algorithm Notation: A partition P over
More informationSTAT STOCHASTIC PROCESSES. Contents
STAT 3911 - STOCHASTIC PROCESSES ANDREW TULLOCH Contents 1. Stochastic Processes 2 2. Classification of states 2 3. Limit theorems for Markov chains 4 4. First step analysis 5 5. Branching processes 5
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Department of Mathematics University of Manitoba Winnipeg Julien Arino@umanitoba.ca 19 May 2012 1 Introduction 2 Stochastic processes 3 The SIS model
More informationDynamic System Properties of Biochemical Reaction Systems
Dynamic System Properties of Biochemical Reaction Systems Balázs Boros Thesis for the Master of Science degree in Applied Mathematics at Institute of Mathematics, Faculty of Science, Eötvös Loránd University,
More informationQUALIFYING EXAM IN SYSTEMS ENGINEERING
QUALIFYING EXAM IN SYSTEMS ENGINEERING Written Exam: MAY 23, 2017, 9:00AM to 1:00PM, EMB 105 Oral Exam: May 25 or 26, 2017 Time/Location TBA (~1 hour per student) CLOSED BOOK, NO CHEAT SHEETS BASIC SCIENTIFIC
More informationMetapopulations with infinitely many patches
Metapopulations with infinitely many patches Phil. Pollett The University of Queensland UQ ACEMS Research Group Meeting 10th September 2018 Phil. Pollett (The University of Queensland) Infinite-patch metapopulations
More informationLecture 20: Reversible Processes and Queues
Lecture 20: Reversible Processes and Queues 1 Examples of reversible processes 11 Birth-death processes We define two non-negative sequences birth and death rates denoted by {λ n : n N 0 } and {µ n : n
More informationOutlines. Discrete Time Markov Chain (DTMC) Continuous Time Markov Chain (CTMC)
Markov Chains (2) Outlines Discrete Time Markov Chain (DTMC) Continuous Time Markov Chain (CTMC) 2 pj ( n) denotes the pmf of the random variable p ( n) P( X j) j We will only be concerned with homogenous
More informationConvergence of Feller Processes
Chapter 15 Convergence of Feller Processes This chapter looks at the convergence of sequences of Feller processes to a iting process. Section 15.1 lays some ground work concerning weak convergence of processes
More informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
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 informationMarkov Chain BSDEs and risk averse networks
Markov Chain BSDEs and risk averse networks Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) 2nd Young Researchers in BSDEs
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationStochastic Chemical Reaction Networks for Robustly Approximating Arbitrary Probability Distributions
Stochastic Chemical Reaction Networks for Robustly Approximating Arbitrary Probability Distributions Daniele Cappelletti Andrés Ortiz-Muñoz David F. Anderson Erik Winfree October 5, 08 Abstract We show
More informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
More informationarxiv: v3 [math.pr] 26 Dec 2014
arxiv:1312.4196v3 [math.pr] 26 Dec 2014 A detailed balanced reaction network is sufficient but not necessary for its Markov chain to be detailed balanced Badal Joshi Abstract Certain chemical reaction
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More information3A Discrete Markov processes and stochastic hybrid systems.
3A Discrete Markov processes and stochastic hybrid systems. So far we have encountered two examples of discrete Markov process: the random walk (Sect. 2.1) and the birth-death process for an ensemble of
More informationIntroduction to Stochastic Processes with Applications in the Biosciences
Introduction to Stochastic Processes with Applications in the Biosciences David F. Anderson University of Wisconsin at Madison Copyright c 213 by David F. Anderson. Contents 1 Introduction 4 1.1 Why study
More informationCDA5530: Performance Models of Computers and Networks. Chapter 3: Review of Practical
CDA5530: Performance Models of Computers and Networks Chapter 3: Review of Practical Stochastic Processes Definition Stochastic ti process X = {X(t), t T} is a collection of random variables (rvs); one
More informationTwo viewpoints on measure valued processes
Two viewpoints on measure valued processes Olivier Hénard Université Paris-Est, Cermics Contents 1 The classical framework : from no particle to one particle 2 The lookdown framework : many particles.
More informationStochastic models for chemical reactions
First Prev Next Go To Go Back Full Screen Close Quit 1 Stochastic models for chemical reactions Formulating Markov models Two stochastic equations Simulation schemes Reaction Networks Scaling limit Central
More informationNoise-induced Mixing and Multimodality in Reaction Networks
Accepted for publication in European Journal of Applied Mathematics (EJAM). 1 Noise-induced Mixing and Multimodality in Reaction Networks TOMISLAV PLESA 1, RADEK ERBAN 1 and HANS G. OTHMER 2 1 Mathematical
More informationHOMOTOPY METHODS FOR COUNTING REACTION NETWORK GHEORGHE CRACIUN, J. WILLIAM HELTON, AND RUTH J. WILLIAMS. August 26, 2008
HOMOTOPY METHODS FOR COUNTING REACTION NETWORK EQUILIBRIA GHEORGHE CRACIUN, J. WILLIAM HELTON, AND RUTH J. WILLIAMS August 26, 2008 Abstract. Dynamical system models of complex biochemical reaction networks
More informationStochastic Processes
Stochastic Processes 8.445 MIT, fall 20 Mid Term Exam Solutions October 27, 20 Your Name: Alberto De Sole Exercise Max Grade Grade 5 5 2 5 5 3 5 5 4 5 5 5 5 5 6 5 5 Total 30 30 Problem :. True / False
More informationIntroduction to self-similar growth-fragmentations
Introduction to self-similar growth-fragmentations Quan Shi CIMAT, 11-15 December, 2017 Quan Shi Growth-Fragmentations CIMAT, 11-15 December, 2017 1 / 34 Literature Jean Bertoin, Compensated fragmentation
More informationMarkov Chains. X(t) is a Markov Process if, for arbitrary times t 1 < t 2 <... < t k < t k+1. If X(t) is discrete-valued. If X(t) is continuous-valued
Markov Chains X(t) is a Markov Process if, for arbitrary times t 1 < t 2
More informationTMA4265 Stochastic processes ST2101 Stochastic simulation and modelling
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 7 English Contact during examination: Øyvind Bakke Telephone: 73 9 8 26, 99 4 673 TMA426 Stochastic processes
More informationA note on Monotone Systems with Positive Translation Invariance
A note on Monotone Systems with Positive Translation Invariance David Angeli Dip. di Sistemi e Informatica Università di Firenze email: angeli@dsi.unifi.it Eduardo D. Sontag Department of Mathematics Rutgers
More informationChapter 22 Examples of Computation of Exact Moment Dynamics for Chemical Reaction Networks
Chapter 22 Examples of Computation of Exact Moment Dynamics for Chemical Reaction Networks Eduardo D. Sontag Abstract The study of stochastic biomolecular networks is a key part of systems biology, as
More informationLecture 20 : Markov Chains
CSCI 3560 Probability and Computing Instructor: Bogdan Chlebus Lecture 0 : Markov Chains We consider stochastic processes. A process represents a system that evolves through incremental changes called
More informationEXAM IN COURSE TMA4265 STOCHASTIC PROCESSES Wednesday 7. August, 2013 Time: 9:00 13:00
Norges teknisk naturvitenskapelige universitet Institutt for matematiske fag Page 1 of 7 English Contact: Håkon Tjelmeland 48 22 18 96 EXAM IN COURSE TMA4265 STOCHASTIC PROCESSES Wednesday 7. August, 2013
More informationMarkov Chains (Part 4)
Markov Chains (Part 4) Steady State Probabilities and First Passage Times Markov Chains - 1 Steady-State Probabilities Remember, for the inventory example we had (8) P &.286 =.286.286 %.286 For an irreducible
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationNecessary and sufficient conditions for strong R-positivity
Necessary and sufficient conditions for strong R-positivity Wednesday, November 29th, 2017 The Perron-Frobenius theorem Let A = (A(x, y)) x,y S be a nonnegative matrix indexed by a countable set S. We
More informationA GLOBAL CONVERGENCE RESULT FOR PROCESSIVE MULTISITE PHOSPHORYLATION SYSTEMS
1 A GLOBAL CONVERGENCE RESULT FOR PROCESSIVE MULTISITE PHOSPHORYLATION SYSTEMS CARSTEN CONRADI AND ANNE SHIU Abstract Multisite phosphorylation plays an important role in intracellular signaling There
More informationData analysis and stochastic modeling
Data analysis and stochastic modeling Lecture 7 An introduction to queueing theory Guillaume Gravier guillaume.gravier@irisa.fr with a lot of help from Paul Jensen s course http://www.me.utexas.edu/ jensen/ormm/instruction/powerpoint/or_models_09/14_queuing.ppt
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationA Simple Protein Synthesis Model
A Simple Protein Synthesis Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 3, 213 Outline A Simple Protein Synthesis Model
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationWinter 2019 Math 106 Topics in Applied Mathematics. Lecture 9: Markov Chain Monte Carlo
Winter 2019 Math 106 Topics in Applied Mathematics Data-driven Uncertainty Quantification Yoonsang Lee (yoonsang.lee@dartmouth.edu) Lecture 9: Markov Chain Monte Carlo 9.1 Markov Chain A Markov Chain Monte
More informationApproximate Counting and Markov Chain Monte Carlo
Approximate Counting and Markov Chain Monte Carlo A Randomized Approach Arindam Pal Department of Computer Science and Engineering Indian Institute of Technology Delhi March 18, 2011 April 8, 2011 Arindam
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
More informationLecture 5. If we interpret the index n 0 as time, then a Markov chain simply requires that the future depends only on the present and not on the past.
1 Markov chain: definition Lecture 5 Definition 1.1 Markov chain] A sequence of random variables (X n ) n 0 taking values in a measurable state space (S, S) is called a (discrete time) Markov chain, if
More informationarxiv: v1 [math.pr] 1 Jan 2013
The role of dispersal in interacting patches subject to an Allee effect arxiv:1301.0125v1 [math.pr] 1 Jan 2013 1. Introduction N. Lanchier Abstract This article is concerned with a stochastic multi-patch
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationEngineering Model Reduction and Entropy-based Lyapunov Functions in Chemical Reaction Kinetics
Entropy 2010, 12, 772-797; doi:10.3390/e12040772 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Engineering Model Reduction and Entropy-based Lyapunov Functions in Chemical Reaction
More informationChapter 7. Markov chain background. 7.1 Finite state space
Chapter 7 Markov chain background A stochastic process is a family of random variables {X t } indexed by a varaible t which we will think of as time. Time can be discrete or continuous. We will only consider
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 informationAlgebraic methods for the study of biochemical reaction networks
Algebraic methods for the study of biochemical reaction networks Alicia Dickenstein Departamento de Matemática, FCEN, Universidad de Buenos Aires, and Instituto de Matemática Luis A. Santaló, UBA CONICET
More informationBirth and Death Processes. Birth and Death Processes. Linear Growth with Immigration. Limiting Behaviour for Birth and Death Processes
DTU Informatics 247 Stochastic Processes 6, October 27 Today: Limiting behaviour of birth and death processes Birth and death processes with absorbing states Finite state continuous time Markov chains
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 2
MATH 56A: STOCHASTIC PROCESSES CHAPTER 2 2. Countable Markov Chains I started Chapter 2 which talks about Markov chains with a countably infinite number of states. I did my favorite example which is on
More informationDiscrete Time Markov Chain of a Dynamical System with a Rest Phase
Discrete Time Markov Chain of a Dynamical System with a Rest Phase Abstract A stochastic model, in the form of a discrete time Markov chain, is constructed to describe the dynamics of a population that
More information1 IEOR 4701: Continuous-Time Markov Chains
Copyright c 2006 by Karl Sigman 1 IEOR 4701: Continuous-Time Markov Chains A Markov chain in discrete time, {X n : n 0}, remains in any state for exactly one unit of time before making a transition (change
More informationA NICE PROOF OF FARKAS LEMMA
A NICE PROOF OF FARKAS LEMMA DANIEL VICTOR TAUSK Abstract. The goal of this short note is to present a nice proof of Farkas Lemma which states that if C is the convex cone spanned by a finite set and if
More informationAlignment processes on the sphere
Alignment processes on the sphere Amic Frouvelle CEREMADE Université Paris Dauphine Joint works with : Pierre Degond (Imperial College London) and Gaël Raoul (École Polytechnique) Jian-Guo Liu (Duke University)
More information