Stochastic and deterministic models of biochemical reaction networks
|
|
- Bennett Rodgers
- 5 years ago
- Views:
Transcription
1 Stochastic and deterministic models of biochemical reaction networks David F. Anderson Department of Mathematics University of Wisconsin - Madison Case Western Reserve November 9th, 215
2 Goals of the talk Introduce mathematical models from molecular biology reaction networks: Stochastic models continuous time Markov chains with stochastic mass action kinetics. Deterministic models ODEs with deterministic mass action kinetics. Introduce the field of chemical reaction network theory Focus on two results how do models behave in large time-frames? Collaborators: Gheorghe Craciun, Tom Kurtz, German Enciso, Matthew Johnston, Daniele Cappelletti Funding: Army Research Office, grant W911NF
3 Reaction networks Models usually begin with something called a reaction network. A + B 2B B A or E + S ES E + P E The network is a static object. We are not yet choosing a dynamics.
4 Reaction networks Biochemical/population networks can range from simple to very complex. Example 1: A. Example 2: A + B 2B B A Example 3: Gene transcription & translation: G G + M transcription M M + P translation M degradation P degradation G + P B Binding/unbinding Cartoon representation: 1 1 J. Paulsson, Physics of Life Reviews, 2,
5 Reaction networks Example 4: EnvZ/OmpR signaling system 2 2 Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, 21
6 Big picture Hanahan and Weinberg, The Hallmarks of Cancer, Cell, 2.!
7 Big picture Metabolic Pathways Roche Applied Science!
8 Big picture Goal of chemical reaction network theory: discover how the dynamics of the mathematical model depend upon properties of the reaction network. Key point of research: want results that are applicable to whole classes of reaction networks, not studying a single model.
9 Short (incomplete) history of chemical reaction network theory Dates back to at least work of Horn, Jackson, Feinberg in 1972: Detailed balancing is not necessary for many nice results. Developed idea of complex balancing. 197s: development of deficiency theory to predict dynamics of deterministic models based solely on easily checked network conditions (deficiency = and weak reversibility). 2s: many more people started to join this research area 1. Large focus on global attractor conjecture. 2. Focus on possibility of multiple equilibria or oscillations. 3. Algebraists (Sturmfels, Shiu, Miller, Dickenstein) got involved: Toric varieties, Birch s theorem, Beginning to be a focus on stochastic models (green fluorescent proteins, etc.)
10 Reaction Networks: {S, C, R} Example: A S = {A}. C = {, A}. R = { A, A }.
11 Reaction Networks: {S, C, R} Example: A + B 2B B A S = {A, B}. C = {A + B, 2B, B, A}. R = {A + B 2B, B A}.
12 Reaction Networks: {S, C, R} Example: Species: S = {A, B, C, D, E}. Complexes: C = {A, 2B, A + C, D, B + E}. Reactions: R = {A 2B, 2B A, A+C D, D A+C, D B+E, B+E A+C}.
13 Reaction Networks: {S, C, R} Definition A chemical reaction network, {S, C, R}, consists of: 1. Species, S := {S 1,..., S d }: constituent molecules undergoing a series of chemical reactions. 2. Complexes, C: linear combinations of the species representing those used, and produced, in each reaction. 3. A set of reactions, R := {y k y k}. Denote reaction vectors ζ k = y k y k.
14 Dynamics deterministic with mass-action kinetics Example: A + B κ 1 2B B κ 2 A (R1) (R2) Let x(t) R 2 give concentrations of molecules of A and B: [ 1 ẋ = r 1 (x) 1 ] [ 1 + r 2 (x) 1 ]. Deterministic mass-action kinetics says: so r 1 (x) = κ 1 x A x B, and r 2 (x) = κ 2 x B. ẋ A = κ 1 x A x B + κ 2 x B ẋ B = κ 1 x A x B κ 2 x B.
15 Dynamics deterministic Consider a general system with S = {S 1,..., S d }, and kth reaction y k y k The rate of kth reaction is r k : R d R. As before: or ẋ(t) = k x(t) = x() + k r k (x(t))(y k y k ), ( t ) r k (x(s))ds (y k y k ).
16 Deterministic mass-action kinetics Consider reaction Then rate is y k y k r k (x) = κ k x y k = κ k d i=1 x y ki i. Example: If S 1 anything, then r k (x) = κ k x 1. Example: If S 1 + S 2 anything, then r k (x) = κ k x 1 x 2. Example: If 2S 2 anything, then r k (x) = κ k x2 2. Example: If 3S 1 + 2S 2 + S 3 anything, then r k (x) = κ k x1 3 x2 2 x 3. Yields: ẋ = k κ k x y k (y k y k )
17 Dynamics deterministic Example:
18 Dynamics: discrete stochastic Example: α β A (R1/R2) For Markov models can take X(t) = X() + R 1 (t) R 2 (t). R 1 (t) = Y 1 (αt) R 2 (t) = Y 2 (β t ) X A (s)ds where Y 1, Y 2 are independent unit-rate Poisson processes. t ) X(t) = X() + Y 1 (αt) Y 2 (β X A (s)ds.
19 Dynamics: discrete stochastic Example: A + B κ 1 2B B κ 2 A (R1) (R2) For Markov models can take [ 1 X(t) = X() + R 1 (t) 1 R 1 (t) = Y 1 ( κ 1 t R 2 (t) = Y 2 ( κ 2 t ] [ 1 + R 2 (t) 1 ) X A (s)x B (s)ds ) X B (s)ds where Y 1, Y 2 are independent unit-rate Poisson processes. ].
20 Model of viral infection: Srivistava, You, Summers, and Yin, J. Theor. Biol., 218 (22) R 1 ) T 1 T + G, R 2 ) G.25 T, R 3 ) T 1 T + S, R 4 ) T.25, R 5 ) S 2, R 6 ) G + S V, X G (t) = X G () + Y 1 ( t t X S (t) = X S () + Y 3 (1 t X T (t) = X T () + Y 2 (.25 ( t X V (t) = X V () + Y ) ( X T (s)ds Y 2.25 ) ( X T (s)ds Y 5 2 t t ) ( X G (s)ds Y 4.25 ) X G (s)x S (s)ds ) t ) X G (s)ds Y 6 ( X G (s)x S (s)ds ) ) t X S (s)ds Y 6 ( X G (s)x S (s)ds t ) X T (s)ds structure Counts virus genome 1 1 template Time (Days)
21 Dynamics: discrete stochastic For general system, we have S = {S 1,..., S d }, with kth reaction y k y k The rate (or intensity or propensity) of kth reaction is λ k : Z d R. As before: X(t) = X() + k R k (t)(y k y k ), with X(t) = X() + k Y k ( t Y k are independent, unit-rate Poisson processes. ) λ k (X(s))ds (y k y k ),
22 Mass-action kinetics The standard intensity function chosen is stochastic mass-action kinetics: λ k (x) = κ k i x i! (x i y ik )!. 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) κ 2 x 2 2 if x 2 1.
23 Dynamics. Example A + B α 2B B β A (R1) (R2) Stochastic equations X(t) = X() + Y 1 (α t ) [ 1 X A (s)x B (s)ds 1 ] ( t ) [ 1 + Y 2 β X B (s)ds 1 ]. Deterministic equations x(t) = x() + α t [ 1 x A (s)x B (s)ds 1 ] t [ 1 + β x B (s)ds 1 ],
24 Other ways to understand model Could just say that for n Z d, n + y 1 y 1, n + y 2 y 2, n. n + y r y r, with rate λ 1 (n) with rate λ 2 (n) with rate λ r (n) where y k y k Z d. I.e. a continuous time Markov chain with infinitesimal generator Af (x) = k λ k (x)(f (x + y k y k ) f (x)). and forward (chemical master) equation p t (x) = k p t(x y k + y k )λ k (x y k + y k ) k p t(x)λ k (x)
25 First connection (Tom Kurtz 197s) The ODE model is a limit of the stochastic model. Assume V is a scaling parameter: volume times Avogadro s number, or 1. X i = O(V ), and X V (t) = def V 1 X(t), λ k (X(t)) V ( κ k X V (t) y k ),
26 Connection: Law of large numbers (Tom Kurtz 197s) Then, X V (t) 1 V X + k ( 1 t ) V Y k V κ k X V (s) y k ds (y k y k ) LLN for Y k says 1 V Y k(vu) u ( lim sup V 1 Y k (Vu) u =, V u U a.s. ) so a good approximation is solution to x(t) = x() + k t κ k x(s) y k ds (y k y k ), where u v = u v 1 1 uv d d, is standard mass-action kinetics. See Tom Kurtz s works...
27 LLN: Example Stochastic models: A + B 2/V 2B B 1 A (R1) (R2) with X() = [3V, V ] so that [A V, B V ] = X/V satisfies A V () = 3, B V () = 1. ODE model of A + B 2 2B B 1 A, with x() = [3, 1].
28 LLN: Example, A + B 2/V 2B B 1 A V= A B V=1 B A
29 LLN: Example, A + B 2/V 2B B 1 A V= A B V= A B
30 LLN: Example, A + B 2/V 2B B 1 A V=1 A B V= A B
31 Next questions. The limit is: over [, T ] with T finite. X V x So what about infinite time-frames?
32 Story # 1: Complex balanced networks Theorem (Horn, Jackson, Feinberg - 197s) Let {S, C, R} be a chemical reaction network with rate constants κ k. Suppose: 1. the network is weakly reversible, and 2. has a deficiency of zero. Then, for the ODE system (i) within each stoichiometric compatibility class (invariant manifold), there is a unique positive equilibria. (ii) each equilibria is complex-balanced. (iii) each equilibria is globally asymptotically stable relative to the interior of its compatibility class.
33 Connectivity Definition The connected components of the reaction network are called the linkage classes. C 2 C 3 C 5 C 6 C 7 C 1 C 4 C 8 Example Has two linkage classes. A + B α 2B (Linkage Class 1) B β A (Linkage Class 2)
34 Connectivity 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. Weakly Reversible C 1 Reversible C 1 C 3 C 2 C 3 C 2
35 Connectivity The following is not weakly reversible: C 2 C 3 C 5 C 6 C 7 C 1 C 4 C 8 Neither is the following: A + B α 2B (Linkage Class 1) B β A (Linkage Class 2)
36 Network properties Definition S = span {yk y k R}{y k y k } is the stoichiometric subspace of the network. Denote dim(s) = s. Implication: Solutions bound to translations of S. Example: Reaction network A + B 2B B A # B Molecules # A Molecuels
37 Deficiency. deficiency of {S, C, R} = δ = n l s, where 1. n = # of complexes. 2. l = # of linkage classes. 3. s = dimension of span of reaction vectors. So it is easy to check. Example A + B 2B B A (R1) (R2) n = 4, l = 2, s = 1 = δ = 1.
38 Deficiency Example n = 5 l = 2 s = 3 = δ = =.
39 Deficiency Of course, you are thinking I have no idea what deficiency is that was pointless. In reality, ẋ = f (x) = YA κψ(x), where 1. Y and A κ are matrices, 2. and deficiency = δ = dim(ker(y ) imagea κ).
40 Story # 1: Complex balanced networks Theorem (Horn, Jackson, Feinberg - 197s) Let {S, C, R} be a chemical reaction network with rate constants κ k. Suppose: 1. the network is weakly reversible, and 2. has a deficiency of zero. Then, the ODE system is complex-balanced and, (i) within each stoichiometric compatibility class (invariant manifold), there is a unique positive equilibria. (ii) each equilibria is complex-balanced. (iii) each equilibria is globally asymptotically stable relative to the interior of its compatibility class. c is complex-balanced if for each η C, κ k c y k = k:y k =η k:y k =η κ k c yk. (1)
41 Complex-balanced Example Balance for complex A + C means For complex A means κ A+C D c A c C = κ D A+C c D + κ B+E A+C c B c E κ A 2B c A = κ 2B A c 2 B Together gives [ ] f (c) A = κ A 2B c A + κ 2B A cb 2 +[κ D A+C c D + κ B+E A+C c B c E κ A+C D c A c C ] =
42 Story # 1: Deficiency Zero Theorem - stochastic Theorem (A., Craciun, Kurtz, 21) Let {S, C, R} be a chemical reaction network with rate constants κ k. Suppose: 1. the network is weakly reversible, and 2. has a deficiency of zero. Then, for any irreducible communicating equivalence class, Γ, the stochastic system has a product form stationary distribution π(x) = 1 Z V d i=1 e c c x i i i, x Γ, (2) x i! where Z V is a normalizing constant and c is a complexed-balanced equilibrium of the corresponding ODE. David F. Anderson, Gheorghe Craciun, and Thomas G. Kurtz, Product-form stationary distributions for deficiency zero chemical reaction networks, Bulletin of Mathematical Biology, Vol. 72, No. 8, , 21.
43 Example: Enzyme kinetics Consider the possible model of enzyme kinetics given by E + S ES E + P, E S Easy to see that deficiency is and that Γ = Z 4. Thus, in distributional equilibrium, the species numbers are independent and have Poisson distributions.
44 Enzyme kinetics Consider the slightly different enzyme kinetics given by E + S ES E + P, E 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.
45 Differences Story #2: are there times the models behave very differently?
46 Story 2: Absolute Concentration Robustness Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, 21. 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 () + x B (), Solving for equilibria: x A = β α, x B = M β α, Network has absolute concentration robustness in species A.
47
48 Theorems: deterministic and stochastic Theorem (M. Feinberg and G. Shinar, Science, 21 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.
49 Differing in one species Examples: 1. differ in species B. A, A + B 2. differ in species Y p. XT, XT + Y p 3. differ in species G. T, T + G
50 Terminal and non-terminal complexes 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 1 k 11 k 5 The orange complexes are called terminal. The blue complexes are called non-terminal.
51 Theorems: deterministic and stochastic Theorem (M. Feinberg and G. Shinar, Science, 21 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. What about stochastically modeled systems? Theorem (A., 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 with probability one there there is a last time a nonterminal reaction fires. 2 David F. Anderson, Germán Enciso, and Matthew Johnston, Stochastic analysis of biochemical reaction networks with absolute concentration robustness, J. Royal Society Interface, Vol. 11, , February 12, 214.
52 Example Reaction network A + B 2B B A has state space # B Molecules # A Molecuels
53 Example X tot := X + XD + XT + X p + X py + XDY p Y tot := Y + X py + XDY p + Y p.
54 Idea of proof 1. Suppose not true: then there is a positive recurrent state for which a nonterminal complex has positive intensity (rate). 2. This provides a series of stopping times, τ N. 3. Use stochastic equation: = X (τ N ) X () = k = [ ( 1 τn κ k τn κ k k (X (s)) y k ds Y k κ k ( τn ) Y k κ i (X (s)) y k ds (y k y k ) ) (X (s)) y k ds 1 t N ] τn (X (s)) y k ds (y k y k ) 4. In limit N, find something in kernel of YA κ that can not be there.
55 Extinction can be rare event. What s real behavior? A + B α 2B B β A X A () + X B () = M, # B Molecules # A Molecuels Can show that quasi-stationary distribution for A converges to Poisson πm(x) Q e (β/α) (β/α) x, as M. x!
56 Quasi-stationary distribution: EnvZ-OmpR signaling system 3 3 Guy Shinar and Martin Feinberg, Structural Sources of Robustness in Biochemical Reaction Networks, Science, 21
57 Quasi-stationary distribution Quasi stationary probabilities X tot = 1 Y tot = 35 X tot = 1 Y tot = 35 X tot = 1 Y tot = 35 Poisson Molecules of Y p
58 Behavior on compact time intervals Let N be a scaling parameter (total conserved value), and consider N 1. Let q be the value of the ACR equilibrium value and let J Poisson(q). (3) Theorem (A., Cappelletti, Kurtz, 215 being written now) Suppose T > and that some technical assumptions hold on the reaction network and let J be as in (3). Then, for any continuous function ĝ : R χ ACR R with at most polynomial growth rate we have or equivalently, P sup T t 1,t 2 [,T ] t 1 <t 2 t2 t 1 (ĝ(x N ACR (s)) E[ĝ(J)]) ds = N. N (ĝ(xacr (s)) E[ĝ(J)]) ds > ε N
59 That is the story. Thanks! Main References: 1. 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, , 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, J. Royal Society Interface, Vol. 11, , February 12, David F. Anderson, Gheorghe Craciun, Manoj Gopalkrishnan, and Carsten Wiuf, Lyapunov functions, stationary distributions, and non-equilibrium potential for reaction networks, Bulletin of Math. Biol., David F. Anderson and Thomas G. Kurtz, Stochastic analysis of biochemical systems, Springer, Jeremy Gunawardena, Chemical reaction network theory for in-silico biologists, 23,
Longtime 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 informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Boston University
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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationPolynomial Dynamical Systems as Reaction Networks and Toric Differential Inclusions
Polynomial Dynamical Systems as Reaction Networks and Toric Differential Inclusions Gheorghe Craciun Department of Mathematics and Department of Biomolecular Chemistry University of Wisconsin-Madison e-mail:
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 informationALGEBRA: From Linear to Non-Linear. Bernd Sturmfels University of California at Berkeley
ALGEBRA: From Linear to Non-Linear Bernd Sturmfels University of California at Berkeley John von Neumann Lecture, SIAM Annual Meeting, Pittsburgh, July 13, 2010 Undergraduate Linear Algebra All undergraduate
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 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 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 informationarxiv: v1 [q-bio.mn] 25 May 2018
Mathematical Analysis of Chemical Reaction Systems Polly Y. Yu Gheorghe Craciun May 29, 208 arxiv:805.037v [q-bio.mn] 25 May 208 Abstract. The use of mathematical methods for the analysis of chemical reaction
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 informationA geometric approach to the Global Attractor Conjecture
A geometric approach to the Global Attractor Conjecture Manoj Gopalkrishnan, Ezra Miller, and Anne Shiu 16 December 2013 Abstract This paper introduces the class of strongly endotactic networks, a subclass
More informationBasic modeling approaches for biological systems. Mahesh Bule
Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving
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 informationCHEMICAL REACTION NETWORK THEORY
CHEMICAL REACTION NETWORK THEORY Contents 1. Notations 2 2. Network definition 2 3. Deterministic mass-action CRN 2 4. Stochastic mass-action CRN 2 5. Mapping 3 6. Network structure 4 6.1. Linkage classes
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 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 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 informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
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 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 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 informationCybergenetics: Control theory for living cells
Department of Biosystems Science and Engineering, ETH-Zürich Cybergenetics: Control theory for living cells Corentin Briat Joint work with Ankit Gupta and Mustafa Khammash Introduction Overview Cybergenetics:
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 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 informationarxiv: v3 [q-bio.mn] 20 Oct 2011
arxiv:1105.2359v3 [q-bio.mn] 20 Oct 2011 Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class Santiago Akle Onkar Dalal
More informationMatthew Douglas Johnston December 15, 2014 Van Vleck Visiting Assistant Professor
Matthew Douglas Johnston December 15, 2014 Van Vleck Visiting Assistant Professor mjohnston3@wisc.edu Department of Mathematics www.math.wisc.edu/ mjohnston3 Phone: (608) 263-2727 480 Lincoln Dr., Madison,
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationON THE PARAMETRIC UNCERTAINTY OF WEAKLY REVERSIBLE REALIZATIONS OF KINETIC SYSTEMS
HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY VESZPRÉM Vol. 42(2) pp. 3 7 (24) ON THE PARAMETRIC UNCERTAINTY OF WEAKLY REVERSIBLE REALIZATIONS OF KINETIC SYSTEMS GYÖRGY LIPTÁK, GÁBOR SZEDERKÉNYI,2 AND KATALIN
More informationProgress in Chemical Reaction Network Theory
Progress in Chemical Reaction Network Theory Organizers: Carsten Wiuf (University of Copenhagen, DK), Elisenda Feliu (University of Copenhagen, DK) Monday, July 2, 17:15 19:15, Medium Hall B Talks: Casian
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 informationJournal of Symbolic Computation
Journal of Symbolic Computation 44 (2009) 1551 1565 Contents lists available at ScienceDirect Journal of Symbolic Computation journal homepage: wwwelseviercom/locate/jsc Toric dynamical systems Gheorghe
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationMultivariate Risk Processes with Interacting Intensities
Multivariate Risk Processes with Interacting Intensities Nicole Bäuerle (joint work with Rudolf Grübel) Luminy, April 2010 Outline Multivariate pure birth processes Multivariate Risk Processes Fluid Limits
More informationAlgebraic methods for biochemical reaction network theory. Anne Joyce Shiu
Algebraic methods for biochemical reaction network theory by Anne Joyce Shiu A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Mathematics and
More informationprocess on the hierarchical group
Intertwining of Markov processes and the contact process on the hierarchical group April 27, 2010 Outline Intertwining of Markov processes Outline Intertwining of Markov processes First passage times of
More informationExtending the Tools of Chemical Reaction Engineering to the Molecular Scale
Extending the Tools of Chemical Reaction Engineering to the Molecular Scale Multiple-time-scale order reduction for stochastic kinetics James B. Rawlings Department of Chemical and Biological Engineering
More informationNonlinear Dynamical Systems Lecture - 01
Nonlinear Dynamical Systems Lecture - 01 Alexandre Nolasco de Carvalho August 08, 2017 Presentation Course contents Aims and purpose of the course Bibliography Motivation To explain what is a dynamical
More informationComputer Simulation and Data Analysis in Molecular Biology and Biophysics
Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1
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 informationLecture 7: Simple genetic circuits I
Lecture 7: Simple genetic circuits I Paul C Bressloff (Fall 2018) 7.1 Transcription and translation In Fig. 20 we show the two main stages in the expression of a single gene according to the central dogma.
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 informationLIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974
LIMITS FOR QUEUES AS THE WAITING ROOM GROWS by Daniel P. Heyman Ward Whitt Bell Communications Research AT&T Bell Laboratories Red Bank, NJ 07701 Murray Hill, NJ 07974 May 11, 1988 ABSTRACT We study the
More informationAdvanced Computer Networks Lecture 2. Markov Processes
Advanced Computer Networks Lecture 2. Markov Processes Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/28 Outline 2/28 1 Definition
More informationMarkov processes and queueing networks
Inria September 22, 2015 Outline Poisson processes Markov jump processes Some queueing networks The Poisson distribution (Siméon-Denis Poisson, 1781-1840) { } e λ λ n n! As prevalent as Gaussian distribution
More informationLognormal Moment Closures for Biochemical Reactions
Lognormal Moment Closures for Biochemical Reactions Abhyudai Singh and João Pedro Hespanha Abstract In the stochastic formulation of chemical reactions, the dynamics of the the first M -order moments of
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 informationEvolutionary Game Theory: Overview and Recent Results
Overviews: Evolutionary Game Theory: Overview and Recent Results William H. Sandholm University of Wisconsin nontechnical survey: Evolutionary Game Theory (in Encyclopedia of Complexity and System Science,
More informationFrom cell biology to Petri nets. Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE
From cell biology to Petri nets Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE Biology = Concentrations Breitling / 2 The simplest chemical reaction A B irreversible,
More informationInference of genomic network dynamics with non-linear ODEs
Inference of genomic network dynamics with non-linear ODEs University of Groningen e.c.wit@rug.nl http://www.math.rug.nl/ ernst 11 June 2015 Joint work with Mahdi Mahmoudi, Itai Dattner and Ivan Vujačić
More informationMathematical Analysis of Chemical Reaction Systems
DOI:./ijch. Mathematical Analysis of hemical Reaction Systems Polly Y. Yu + * [a] and Gheorghe raciun ++ * [b] Abstract: The use of mathematical methods for the analysis of chemical reaction systems has
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
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 informationGene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back
Gene Expression as a Stochastic Process: From Gene Number Distributions to Protein Statistics and Back June 19, 2007 Motivation & Basics A Stochastic Approach to Gene Expression Application to Experimental
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 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 informationToric Dynamical Systems
Toric Dynamical Systems Gheorghe Craciun Dept. of Mathematics, University of Wisconsin, Madison, WI 53706-1388, USA Alicia Dickenstein Dep. de Matemática, FCEN, Universidad de Buenos Aires, (1428), Argentina
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 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 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 informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
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 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 informationOn Finding Optimal Policies for Markovian Decision Processes Using Simulation
On Finding Optimal Policies for Markovian Decision Processes Using Simulation Apostolos N. Burnetas Case Western Reserve University Michael N. Katehakis Rutgers University February 1995 Abstract A simulation
More informationRandom Locations of Periodic Stationary Processes
Random Locations of Periodic Stationary Processes Yi Shen Department of Statistics and Actuarial Science University of Waterloo Joint work with Jie Shen and Ruodu Wang May 3, 2016 The Fields Institute
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationAn All-Encompassing Global Convergence Result for Processive Multisite Phosphorylation Systems
An All-Encompassing Global Convergence Result for Processive Multisite Phosphorylation Systems Mitchell Eithun Department of Mathematical Sciences, Ripon College, Ripon, WI 54971 Anne Shiu Department of
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
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 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 informationOn oscillating systems of interacting neurons
Susanne Ditlevsen Eva Löcherbach Nice, September 2016 Oscillations Most biological systems are characterized by an intrinsic periodicity : 24h for the day-night rhythm of most animals, etc... In other
More informationarxiv: v2 [q-bio.mn] 31 Aug 2007
A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays David F. Anderson Department of Mathematics, University of Wisconsin-Madison, Madison, Wi 53706
More informationBasic Synthetic Biology circuits
Basic Synthetic Biology circuits Note: these practices were obtained from the Computer Modelling Practicals lecture by Vincent Rouilly and Geoff Baldwin at Imperial College s course of Introduction to
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 informationA Note on the Central Limit Theorem for a Class of Linear Systems 1
A Note on the Central Limit Theorem for a Class of Linear Systems 1 Contents Yukio Nagahata Department of Mathematics, Graduate School of Engineering Science Osaka University, Toyonaka 560-8531, Japan.
More informationCompensating operator of diffusion process with variable diffusion in semi-markov space
Compensating operator of diffusion process with variable diffusion in semi-markov space Rosa W., Chabanyuk Y., Khimka U. T. Zakopane, XLV KZM 2016 September 9, 2016 diffusion process September with variable
More informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More informationDISSIPATIVITY OF NEURAL NETWORKS WITH CONTINUOUSLY DISTRIBUTED DELAYS
Electronic Journal of Differential Equations, Vol. 26(26), No. 119, pp. 1 7. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) DISSIPATIVITY
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
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 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 informationIowa State University. Instructor: Alex Roitershtein Summer Homework #5. Solutions
Math 50 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 205 Homework #5 Solutions. Let α and c be real numbers, c > 0, and f is defined
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 information