Persistence and Stationary Distributions of Biochemical Reaction Networks
|
|
- Kerrie Sharp
- 5 years ago
- Views:
Transcription
1 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 3rd, 2008
2 Overview 1. Biochemical systems can be large and are typically quite complex. 2. There are potentially millions of such systems that could warrant study at one time or another. 3. The dynamics of these systems can be modeled stochastically or deterministically (or a combination of both). 4. Rate constants (key system parameter values) are oftentimes unknown and even vary from cell to cell. 5. Network structure induces equations governing dynamics. 6. Want to infer qualitative properties of the dynamics from the network structure without the need for details of system parameter values. 7. (if time) Simulation methods for stochastic systems.
3 Examples: Gene Regulation and Enzyme Kinetics 1 2 Gene regulatory network within the plant A. thaliana. Pointed arrows reflect an activation relationship, whereas the blunt arrows reflect a repression relationship between genes. S + E 1 SE 1 P + E 1, E 1 A + E 2, E 2. 1 Shahrezaei and Swain, Curr. Opinion in Biotech., Gambin et al., In Silico, 2006.
4 Example: Metabolism
5 Example: Metabolism
6 Chemical Reaction Networks Standard notation for chemical reactions: A + B C is interpreted as a molecule of A combines with a molecule of B to give a molecule of C. We make the association: 2 A + B = ν = , C = ν = Note: each instance of the reaction A + B C or ν ν changes the state of the system by the vector: ν ν.
7 Chemical Reaction Networks: Species, Complexes, and Reactions A chemical reaction network consists of: 1. N 1 chemical species, {S 1,..., S N } = S, undergoing a series of chemical reactions. 2. Vectors, ν k, ν k Z N 0, representing the numbers of molecules of each species consumed and created in the kth reaction, respectively. Each vector may appear as both a source (ν) and product (ν ). We associate with each ν k (and ν k ) a linear combination of the species in which the coefficient of S i is ν ki. 2 3 For example, if ν k = , then ν k is associated with 2S 1 + S 3. Under this association, {ν k } = C are the complexes of the network. 3. R = {ν k ν k} are the reactions. Definition. The triple {S, C, R} is called a chemical reaction network.
8 Reaction Diagrams To each reaction network, {S, C, R}, there is a unique directed graph, called a reaction diagram, constructed in the following manner: 1. The nodes of the graph are the complexes, C. 2. A directed edge (ν, ν ) exists if and only if ν ν R. Each connected component of the resulting graph is termed a linkage class of the graph. Definition. The chemical reaction network is said to be weakly reversible if each linkage class of the corresponding reaction diagram is strongly connected. A network is said to be reversible if ν ν R = ν ν R.
9 Chemical Reaction Networks Example. Consider the reaction network with the following diagram: S + E 1 SE 1 P + E 1, E 1 A + E 2, E S = {S, E 1, SE 1, P, A, E 2 }. 2. C = {S + E 1, SE 1, P + E 1, E 1, A + E 2,, E 2 }. 3. R = {S + E 1 SE 1, SE 1 S + E 1, SE 1 P + E 1, E 1 A + E 2, A + E 2 E 1, E 2, E 2 }. Also, this system is not weakly reversible.
10 Deterministic Dynamics: Mass Action Kinetics We will assume that when the concentration vector of the system is x R N 0, the rate of the reaction ν k ν k, is R k (x) = κ k x ν 1k 1 x ν 2k 2 x ν Nk N = κ k x ν k Example: If S 1 + 2S 2 anything, then R k (x) = κ k x 1 x 2 2. Equations governing dynamics are the nonlinear ODEs: ẋ(t) = X k R k (x(t))(ν k ν k ) = X k κ k x(t) ν k (ν k ν k ), A system with kinetics given by R k (x) = κ k x ν k is said to have deterministic mass action kinetics. Note: up to parameter values, κ k, equations arise solely from the network structure.
11 ω-limit points, Persistence, and a Conjecture Definition. Let x(t) be a trajectory with x(0) R N >0. The set of ω-limit points for this trajectory is: ω(x(0)) = {x R N 0 x(t n) x, for some sequence t n }. Definition. A system is persistent if ω(x(0)) R N 0 =, for all x(0) R N >0. Conjecture (Feinberg ): Let {S, C, R} be a weakly reversible network with deterministic mass action kinetics. Then this system is persistent.
12 Preliminary Observations and Definitions 1. Any solution x(t) satisfies x(t) x(0) = X Z t k 0 «R k (x(s)) ds (ν k ν k ) span{ν k ν k } k. Definition. S = span{ν k ν k } k is the stoichiometric subspace of the network. Let dim(s) = s. We see that for all t 0: x(t) x(0) + S = {z R N z = x(0) + v, for some v S}. 2. Both R N 0 and R N >0 are forward invariant. Definition. We let: P x(0) = P = (x(0) + S) R N 0, be the positive stoichiometric compatibility class associated with x(0).
13 The Faces of P: Boundary Points The set P = (x(0) + S) R N 0 is a polyhedron of dimension dim(s) = s. For W S, let Z W = {x R N x i = 0 if i W }, be the linear space whose zero-set is W. For each face of P, call it F, there is a W S, such that F = F W = P Z W. The dimension of each face satisfies 0 dim(f W ) s. A facet is a face such that dim(f W ) = dim(p) 1 = s 1. A vertex is a face such that dim(f W ) = 0. Note: F W = P Z W = if P does not intersect the face of the positive orthant associated with W. Clearly, no face of P corresponds with such a W. In this case we say that Z W is stoichiometrically unattainable.
14 Example Consider a system with reaction diagram: 2A A + B, B C, S R 3 is the plane spanned by ( 1, 1, 0) and (0, 1, 1). For T > 0, the positive compatibility classes are two-dimensional triangles n o P = (x a, x b, x c) R 3 0 x a + x b + x c = T The three facets of P are the one dimensional segments F {A}, F {B}, F {C}. The three vertices are the three points F {A,B}, F {A,C}, F {B,C}. The set Z {A,B,C} = { 0} is stoichiometrically unattainable.
15 Semilocking Sets Definition. A nonempty subset, W, of the set of species, S, is called a semilocking set if for each reaction in which there is an element of W in the product complex, there is an element of W in the reactant complex. Example: X 1 + X 2 X 3. If X 3 is in a semilocking set, then so too must X 1 and/or X 2. Example: 2A A + B, 2B A + C, D E. Semilocking sets are W 1 = {A, B}, W 2 = {D, E}, W 3 = {A, B, C}, W 4 = {A, B, D, E}. Intuition: If the concentration of each element of W is zero, then any flux which affects the concentrations of the species of W is zero, and the concentrations of the elements of W are locked at zero for all time.
16 Semilocking Sets Theorem (Feinberg , Angeli, De Leenheer, Sontag , A ). Let x(0) R N >0, and let P = (x(0) + S) R N 0 denote the corresponding positive class. If there exists an ω(x(0))-limit point, z ω(x(0)), and a subset of the species, W, such that z is contained within the interior of the face F W, then W is a semilocking set. Corollary. Suppose that F W = P Z W = for each semilocking set W S (i.e. Z W is stoichiometrically unattainable). Then the system is persistent. This corollary has been widely applied to biological examples: 1. Sontag : Model of T-Cell receptor signal transduction. 2. Angeli, De Leenheer, Sontag : Different enzymatic reactions and cascades of reactions (i.e. MAPK). 3. Gnacadja : Reversible binding reactions.
17 Example Consider again the system with reaction diagram: 2A A + B, B C, S R 3 is the plane spanned by ( 1, 1, 0) and (0, 1, 1). For T > 0, the positive compatibility classes are two-dimensional triangles n o P = (x a, x b, x c) R 3 0 x a + x b + x c = T The only semilocking sets are W 1 = {A} and W 2 = {A, B, C}. Z W2 = Z {A,B,C} = { 0} is stoichiometrically unattainable. ω-limit points can only be found on F {A}.
18 ω-limit Points and Facets Theorem (A., Shiu ). Let {S, C, R} be a weakly-reversible network with deterministic mass-action kinetics. Let x(0) R N >0 P. Suppose that z ω(x(0)) lies in the interior of the face F W of P, for some W S. Then F W is not a facet of P; that is, dim(f W ) < dim(p) 1. In words: No ω-limit points exist within the facets of P. Corollary. Let {S, C, R} be a weakly-reversible network with deterministic mass-action kinetics. Suppose that each semilocking set, W, satisfies dim(f W ) = dim(p) 1 or F W =. Then the system is persistent. This shows example on previous slide is persistent.
19 Idea of Proof Proof: Suppose z ω(x(0)) F W with dim(f W ) = dim(p) 1. Will find contradiction. The proof consists of three parts. 1. Argue that the span of {ν k ν k } restricted to W is one-dimensional. Say spanned by v R W. 2. Argue that no subset of W can be conserved and this, combined with above point, implies that v can be chosen so that v R W >0. 3. Argue that v R W >0 repelling. implies there is a neighborhood of z which is
20 Application: Complex-Balanced Systems Definition. An equilibrium c R N 0 is said to be complex-balanced if the following equality holds for each η C: X κ k c ν k = X κ k c νk. {k ν k =η} {k ν k =η} A complex-balanced system is a system that admits a complex-balancing equilibrium. Some chemists/physicists believe that all chemical systems are detailed-balanced, and hence, complex-balanced.
21 Deficiency and complex-balancing Definition. The deficiency of a network is δ = n l s, where n is the number of complexes, l is the number of linkage classes, and s = dim(s) = dim(p). Theorem (Feinberg ). Let {S, C, R} be weakly reversible and satisfy δ = 0. Then, independent of the choice of rate constants, the system with deterministic mass action kinetics is complex-balanced. A partial converse holds: Theorem (Feinberg , Gunawardena ). A complex-balanced system is weakly reversible.
22 Behavior of Complex-Balanced Systems Theorem (Horn, Jackson, Feinberg ). Let {S, C, R} be a complex-balanced system. Then, within the interior of each positive class, P, there is precisely one equilibrium value, c, and that equilibrium is: 1. Complex-balanced. 2. Locally asymptotically stable relative to its compatibility class. Global Attractor Conjecture (GAC): For a complex-balanced system, each of the equilibria c R N >0 is globally asymptotically stable relative to the interior of its compatibility class P.
23 Example Consider the system 2A A + B B + C There are : 3 complexes, 1 linkage classes, and s = 2 = δ = 0. The stoichiometric subspace is of dimension two, and the quantity A + B + C is conserved. Thus, each stoichiometric compatibility class is a triangle. Within each positive stoichiometric compatibility class, there is a unique equilibrium value that is locally asymptotically stable relative to its compatibility class.
24 Partial Results Towards the GAC Can be shown that persistence will guarantee global asymptotic stability. = Semilocking sets and ω-limit points are of interest. There is more structure in complex-balanced systems, so more is already known: Theorem (A , Craciun, Dickenstein, Shiu, Sturmfels ). No vertex of a positive class P can be an ω-limit point. Theorem (Craciun, Dickenstein, Shiu, Sturmfels ). If a system is detailed-balanced, P is two-dimensional, and the network is conservative, then the GAC holds. Can generalize this: Theorem (A., Shiu ). If a system is complex-balanced and P is two-dimensional, then the GAC holds.
25 Towards the GAC Most general statement about GAC to date: Theorem (A., Shiu ). Let {S, C, R} be a complex-balanced system. Suppose that for each semilocking set W S, one of the following holds: 1. F W = Z W P = (Z W is stoichiometrically unattainable). 2. dim(f W ) = 0 (F W is a vertex). 3. dim(f W ) = dim(p) 1 (F W is a facet). Then the GAC holds for this system.
26 Example Consider either of the systems 2A A + B B + C or 2A A + B, B C. Each has δ = 0 and is weakly reversible = they are complex balanced systems. They have same stoichiometric compatibility classes P: dim(p) = 2, so the GAC holds.
27 Dynamics: A Stochastic Model 1. Typically used when abundances of molecules are low. 2. Treats system as a continuous time Markov chain with state X Z N 0 representing the abundance of each species present and with each reaction modeled as a possible transition for the state. 3. If reaction k occurs at time t, then X(t) = X(t ) + ν k ν k. 4. The kth reaction has an associated intensity function, λ k (x), such that P(X(t + t) = X(t) + ν k ν k X(t)) = λ k (X(t)) t + o( t). 5. Kolmogorov s Forward equation (Chemical Master Equation) d dt P(x, t) = X k λ k (x ν k + ν k )P(x ν k + ν k, t) X k λ k (x)p(x, t), describes the time evolution of the distribution of the state of the system.
28 Stochastic Mass Action Kinetics One standard intensity function is stochastic mass-action kinetics: N! Y x i Y N x i! λ k (x) = κ k ν ki! = κ k (x i ν ki )!. i=1 ν ki i=1 Example: If S 1 + S 2 anything, then λ k (x) = κ k x 1 x 2. Example: If S 1 + 2S 2 anything, then λ k (x) = κ k x 1 x 2 (x 2 1). Idea: rate is proportional to number of distinct subsets of the molecules present that can form the inputs for the reaction. Intuitively, this assumption reflects the idea that the system is well-stirred in the sense that all molecules are equally likely to be at any location at any time.
29 Qualitative Properties: Stochastically Modeled Systems Recall, stochastically modeled systems are continuous time Markov chains that satisfy Kolmogorov s Forward equation d dt P(x, t) = X k λ k (x ν k + ν k )P(x ν k + ν k, t) X k λ k (x)p(x, t). We see a stationary distribution, π(x), must satisfy X π(x ν k + ν k )λ k (x ν k + ν k ) = π(x) X k k λ k (x), for all x in some closed, irreducible subset of the state space.
30 Stochastic Version of Complex-Balanced Theorem Theorem (Anderson, Craciun, Kurtz ). Let {S, C, R} be a stochastically modeled chemical reaction network with rate constants κ k. Suppose that 1. the associated deterministic mass-action system with rate constants κ k is complex-balanced with equilibrium c R N >0. Then for any closed, irreducible communicating equivalence class, Γ, the stochastic system has a product-form stationary distribution N π(x) = M cx x! = M Y c x i i x i!, x Γ, i=1 where M is a normalizing constant. Recall: Feinberg showed that weak reversibility and zero deficiency implies complex-balancing. These are easily checked conditions that depend only upon the network structure.
31 Example: Open, first-order system If the network is an open, first-order reaction network, then S = R N, Γ = Z N 0 and so, π(x) = e P N i=1 c i N Y i=1 c x i i x i! = N Y i=1 e c c x i i i x i!, x ZN 0. where c R N >0 is the equilibrium of the associated (linear) deterministic system. Thus, when in distributional equilibrium, the species numbers are independent and have Poisson distributions. This is well known. However: 1. Neither the independence nor the Poisson distribution resulted from the fact that the system under consideration was a first order system. 2. Instead, both facts followed from Γ = Z N 0.
32 Example: Enzyme kinetics Consider the possible model of enzyme kinetics given by E + S ES E + P, E S Easy to see that Γ = Z 4 0. Thus, in distributional equilibrium, the species numbers are independent and have Poisson distributions.
33 Generalizations λ k (x) = κ k N Y i=1 θ i (x i )θ i (x i 1) θ i (x i (ν ik 1)), where θ i ( ) is the rate of association of the ith species. Examples: 1. θ i (x i ) = x i gives mass-action kinetics. 2. θ i (x i ) = v ix i k i + x i gives a stochastic Michaelis-Menten kinetics. 3. Can have different kinetics in same model. Ex: Species A governed by Michaelis-Menten Species B governed by mass-action. Then, A + B ν k has intensity λ k (x) = κ k v 1 x 1 k 1 + x 1 x If (1) ν k {0, 1} and (2) θ i (x i ) = min{n i, x i }, then the system is an M/M/n queueing network.
34 Generalizations λ k (x) = κ k N Y i=1 θ i (x i )θ i (x i 1)θ i (x i (ν ik 1)), (1) Theorem (Anderson, Craciun, Kurtz ). Let {S, C, R} be a stochastically modeled chemical reaction network with intensity functions (1). Suppose that 1. the associated deterministic mass-action system with rate constants κ k is complex-balanced with equilibrium c R N >0. Then for any closed, irreducible communicating equivalence class, Γ, the stochastic system has a product-form stationary distribution π(x) = M NY i=1 c x i i Q xi j=1 θ, x Γ, (2) i(j) where M is a normalizing constant, provided that (2) is summable.
35 Example: Michaelis-Menten kinetics Suppose we have a system with only Michaelis-Menten kinetics. Then,! x Y i x Y i v i j θ i (j) = = v x i i / k i + x i. k i + j x i So, j=1 j=1 π(x) = M NY i=1 c x i i Q xi j=1 θ i(j) = M NY i=1 k i + x i x i! «xi ci. v i Note: k i + x i x i! = O(x k i i ), x i, and so the above π(x) is summable if for each species S i whose abundances are unbounded, we have c i < v i.
36 More Generalizations Note: θ(x) = θ(x) λ k (x) = κ k θ(x ν k ). (3) x NY Y i θ i (j), gives the previous (natural) case. i=1 j=1 Theorem (Anderson, Craciun, Kurtz ). Let {S, C, R} be a stochastically modeled chemical reaction network with intensity functions (3). Suppose that 1. the associated mass-action deterministic system with rate constants κ k has a complex balanced equilibrium c R N >0. Then for any closed, irreducible communicating equivalence class, Γ, the stochastic system has a stationary distribution π(x) = M 1 θ(x) cx = M 1 θ(x) NY i=1 c x i i, x Γ, (4) where M is a normalizing constant, provided that (4) is summable.
37 Simulation of Stochastically Modeled Systems Oftentimes mathematical analysis is not sufficient and/or timely. In such cases, simulation methods are utilized. This is not problematic in deterministic setting. This is problematic in stochastic setting as we want to incorporate Monte Carlo methods. Make use of the random time change representation: X(t) = X(0) + X k Y k Z t 0 «λ k (X(s)) ds (ν k ν k ), where the Y k are independent, unit-rate Poisson processes.
38 Simulation of Stochastically Modeled Systems X(t) = X(0) + X k Y k Z t 0 «λ k (X(s)) ds (ν k ν k ), This representation can be utilized to: (A ) Develop efficient and statistically exact algorithms that are easily adaptable to such typical situations as: 1 Time varying (even stochastically driven) rate constants. 2 Delays in the initiation and completion of reactions (such as is common with gene transcription and/or translation models). (A ) Develop approximate algorithms, that are both more stable and more accurate than existing algorithms. (A., Ganguly, Kurtz ) Perform the first rigorous error analysis of the approximate methods using an appropriate scaling.
39 Tau-Leaping Error Analysis Tau-leaping is equivalent to an Euler approximation of the integral: X(τ) = X(0) + X Z τ «Y k λ k (X(s)) ds (ν k ν k ) k 0 X(0) + X k d = X(0) + X k Y k (τ λ k (X(0))) (ν k ν k ) Poisson(λ k (X(0))τ)(ν k ν k ). Previous error analysis was performed under the scaling τ 0. However, in this scaling, exact methods are faster.
40 Tau-Leaping Error Analysis Using the random time change representation you can perform error analysis in the following way: Let V be order of magnitude of typical abundance. Let τ = c/v β for some c > 0 and 0 < β < 1. As V, τ 0 and λ k (X)τ, which is appropriate regime for the use of approximate methods. In this scaling regime, can show the following: 1. Standard tau-leaping is an order τ = 1/V β method in both the weak and strong sense. 2. Midpoint tau-leaping has a strong order 1/V g(β), where j 2β, β 1/3 g(β) = β + 1, (1 β), β 1/3 2 and a weak order of τ 2 = 1/V 2β. This result is in contrast with previous error analysis (τ 0), but in agreement with numerical examples.
41 Example: Irreversible Isomerization (X 0) Consider the simple example Solutions satisfy X 1 0, X(0) = 10 6 = V. X(t) = X(0) Y 1 Z t 0 «X(s) ds. Simulating 100, 000 sample paths with τ = 1/10 2 = 1/V 1/3 yields
42 Partial List of References 1. David F. Anderson and Anne Shiu, Persistence of deterministic population processes and the Global Attractor Conjecture, in preparation (to be submitted very soon). 2. David F. Anderson, Global asymptotic stability for a class of nonlinear chemical equations, Siam J. Appl. Math., 68(5), , May David F. Anderson, Gheorghe Craciun, and Thomas G. Kurtz, Product-form stationary distributions for deficiency zero chemical reaction networks, submitted, arxiv: , link on my website. 4. David F. Anderson, Arnab Ganguly, and Thomas G. Kurtz, Error Analysis of the tau-leap simulation method for stochastically modeled chemical reaction systems, in preparation. 5. David F. Anderson, Incorporating postleap checks in tau-leaping, Journal of Chemical Physics, 128(5), , February David F. Anderson, A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays, Journal of Chemical Physics, 127(21), , December Martin Feinberg, Lectures on chemical reaction networks, delivered at Univ. of Wisconsin, Madison, Available for download at: feinberg/lecturesonreactionnetworks
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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.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 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 informationFUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models
FUNDAMENTALS of SYSTEMS BIOLOGY From Synthetic Circuits to Whole-cell Models Markus W. Covert Stanford University 0 CRC Press Taylor & Francis Group Boca Raton London New York Contents /... Preface, xi
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 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 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 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 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 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 informationRegulation of metabolism
Regulation of metabolism So far in this course we have assumed that the metabolic system is in steady state For the rest of the course, we will abandon this assumption, and look at techniques for analyzing
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 informationGlobal asymptotic stability for a class of nonlinear chemical equations
arxiv:0708.0319v2 [math.ds] 15 Nov 2007 Global asymptotic stability for a class of nonlinear chemical equations David F. Anderson 1 October 24, 2018 Abstract We consider a class of nonlinear differential
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 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 informationIDENTIFIABILITY OF CHEMICAL REACTION NETWORKS
IDENTIFIABILITY OF CHEMICAL REACTION NETWORKS Gheorghe Craciun a,1 Casian Pantea b 1 Corresponding Author. Address: University of Wisconsin-Madison, Department of Mathematics, 48 Lincoln Drive, Madison
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 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 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 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 informationProgram for the rest of the course
Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling
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 informationNon-homogeneous random walks on a semi-infinite strip
Non-homogeneous random walks on a semi-infinite strip Chak Hei Lo Joint work with Andrew R. Wade World Congress in Probability and Statistics 11th July, 2016 Outline Motivation: Lamperti s problem Our
More informationJOINING AND DECOMPOSING REACTION NETWORKS. 1. Introduction
JOINING AND DECOMPOSING REACTION NETWORKS ELIZABETH GROSS, HEATHER HARRINGTON, NICOLETTE MESHKAT, AND ANNE SHIU Abstract In systems and synthetic biology, much research has focused on the behavior and
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
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 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 informationModelling chemical kinetics
Modelling chemical kinetics Nicolas Le Novère, Institute, EMBL-EBI n.lenovere@gmail.com Systems Biology models ODE models Reconstruction of state variable evolution from process descriptions: Processes
More informationTCOM 501: Networking Theory & Fundamentals. Lecture 6 February 19, 2003 Prof. Yannis A. Korilis
TCOM 50: Networking Theory & Fundamentals Lecture 6 February 9, 003 Prof. Yannis A. Korilis 6- Topics Time-Reversal of Markov Chains Reversibility Truncating a Reversible Markov Chain Burke s Theorem Queues
More informationMonotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005
Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 005 Foreword The aim of this presentation is to give a (primitive) overview of monotone systems and monotone
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 informationA Petri net approach to the study of persistence in chemical reaction networks
A Petri net approach to the study of persistence in chemical reaction networks David Angeli Dip. di Sistemi e Informatica, University of Firenze Patrick De Leenheer Dep. of Mathematics, University of Florida,
More informationNo oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism
No oscillations in the Michaelis-Menten approximation of the dual futile cycle under a sequential and distributive mechanism Luna Bozeman 1 and Adriana Morales 2 1 Department of Mathematical Sciences Clemson
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 informationAn Introduction to Stochastic Simulation
Stephen Gilmore Laboratory for Foundations of Computer Science School of Informatics University of Edinburgh PASTA workshop, London, 29th June 2006 Background The modelling of chemical reactions using
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 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 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 informationStructure and Stability of Certain Chemical Networks and Applications to the Kinetic Proofreading Model of T-Cell Receptor Signal Transduction
1028 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 46, NO 7, JULY 2001 Structure Stability of Certain Chemical Networks Applications to the Kinetic Proofreading Model of T-Cell Receptor Signal Transduction
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 informationMULTIPLE EQUILIBRIA IN COMPLEX CHEMICAL REACTION NETWORKS: II. THE SPECIES-REACTION GRAPH
SIAM J. APPL. MATH. Vol. 66, No. 4, pp. 1321 1338 c 2006 Society for Industrial and Applied Mathematics MULTIPLE EQUILIBRIA IN COMPLEX CHEMICAL REACTION NETWORKS: II. THE SPECIES-REACTION GRAPH GHEORGHE
More informationA Petri Net Approach to Persistence Analysis in Chemical Reaction Networks
A Petri Net Approach to Persistence Analysis in Chemical Reaction Networks David Angeli, Patrick De Leenheer, and Eduardo Sontag 1 Dep. of Systems and Computer Science, University of Florence, Italy angeli@dsi.unifi.it
More informationGlobal Stability of Full Open Reversible Michaelis-Menten Reactions
Preprints of the 8th IFAC Symposium on Advanced Control of Chemical Processes The International Federation of Automatic Control Global Stability of Full Open Reversible Michaelis-Menten Reactions Ismail
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 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 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 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 informationA&S 320: Mathematical Modeling in Biology
A&S 320: Mathematical Modeling in Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.ms.uky.edu/~dmu228/as320/ These slides were modified from Matthew Macauley s lecture
More informationOn the numerical solution of the chemical master equation with sums of rank one tensors
ANZIAM J. 52 (CTAC21) pp.c628 C643, 211 C628 On the numerical solution of the chemical master equation with sums of rank one tensors Markus Hegland 1 Jochen Garcke 2 (Received 2 January 211; revised 4
More informationOverview of Kinetics
Overview of Kinetics [P] t = ν = k[s] Velocity of reaction Conc. of reactant(s) Rate of reaction M/sec Rate constant sec -1, M -1 sec -1 1 st order reaction-rate depends on concentration of one reactant
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 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 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 informationTime Reversibility and Burke s Theorem
Queuing Analysis: Time Reversibility and Burke s Theorem Hongwei Zhang http://www.cs.wayne.edu/~hzhang Acknowledgement: this lecture is partially based on the slides of Dr. Yannis A. Korilis. Outline Time-Reversal
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 informationNetworks in systems biology
Networks in systems biology Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Networks in systems
More informationA First Course on Kinetics and Reaction Engineering. Class 9 on Unit 9
A First Course on Kinetics and Reaction Engineering Class 9 on Unit 9 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going A. Rate Expressions - 4. Reaction Rates and Temperature
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 informationStochastic Transcription Elongation via Rule Based Modelling
Stochastic Transcription Elongation via Rule Based Modelling Masahiro HAMANO hamano@jaist.ac.jp SASB, Saint Malo 2015 1 Purpose of This Talk 2 mechano-chemical TE as Rule Based Modelling 3 Table of Contents
More informationStat 516, Homework 1
Stat 516, Homework 1 Due date: October 7 1. Consider an urn with n distinct balls numbered 1,..., n. We sample balls from the urn with replacement. Let N be the number of draws until we encounter a ball
More informationCS-E5880 Modeling biological networks Gene regulatory networks
CS-E5880 Modeling biological networks Gene regulatory networks Jukka Intosalmi (based on slides by Harri Lähdesmäki) Department of Computer Science Aalto University January 12, 2018 Outline Modeling gene
More informationI/O monotone dynamical systems. Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011
I/O monotone dynamical systems Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011 BEFORE: Santa Barbara, January 2003 Having handed to me a photocopied
More informationExample: physical systems. If the state space. Example: speech recognition. Context can be. Example: epidemics. Suppose each infected
4. Markov Chains A discrete time process {X n,n = 0,1,2,...} with discrete state space X n {0,1,2,...} is a Markov chain if it has the Markov property: P[X n+1 =j X n =i,x n 1 =i n 1,...,X 0 =i 0 ] = P[X
More informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA February 14, 2017 Introduction When we choose to model
More informationIntroduction to Bioinformatics
Systems biology Introduction to Bioinformatics Systems biology: modeling biological p Study of whole biological systems p Wholeness : Organization of dynamic interactions Different behaviour of the individual
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XI Stochastic Stability - H.J. Kushner
STOCHASTIC STABILITY H.J. Kushner Applied Mathematics, Brown University, Providence, RI, USA. Keywords: stability, stochastic stability, random perturbations, Markov systems, robustness, perturbed systems,
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 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 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 informationIntrinsic Noise in Nonlinear Gene Regulation Inference
Intrinsic Noise in Nonlinear Gene Regulation Inference Chao Du Department of Statistics, University of Virginia Joint Work with Wing H. Wong, Department of Statistics, Stanford University Transcription
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 informationOscillations in Michaelis-Menten Systems
Oscillations in Michaelis-Menten Systems Hwai-Ray Tung Brown University July 18, 2017 Hwai-Ray Tung (Brown University) Oscillations in Michaelis-Menten Systems July 18, 2017 1 / 22 Background A simple
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 informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
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 informationLecture 13: Data Analysis for the V versus [S] Experiment and Interpretation of the Michaelis-Menten Parameters
Biological Chemistry Laboratory Biology 3515/Chemistry 3515 Spring 2018 Lecture 13: Data Analysis for the V versus [S] Experiment and Interpretation of the Michaelis-Menten Parameters 20 February 2018
More information1 Periodic stimulations of the Incoherent Feedforward Loop network
1 Periodic stimulations of the Incoherent Feedforward Loop network In this Additional file, we give more details about the mathematical analysis of the periodic activation of the IFFL network by a train
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 informationCycle structure in SR and DSR graphs: implications for multiple equilibria and stable oscillation in chemical reaction networks
Cycle structure in SR and DSR graphs: implications for multiple equilibria and stable oscillation in chemical reaction networks Murad Banaji Department of Mathematics, University of Portsmouth, Lion Gate
More information