Stochastic Gene Expression: Modeling, Analysis, and Identification
|
|
- Ronald Sullivan
- 5 years ago
- Views:
Transcription
1 Munsky; q-bio Stochastic Gene Expression: Modeling, Analysis, and Identification Mustafa Khammash University of California, Santa Barbara Center for Control, Dynamical Systems and Computation CC DC
2 Outline Randomness at the molecular level Exploiting the randomness Stochastic modeling framework Simulation and Analysis Methods Identification from Noise Conclusions
3 Randomness at the Molecular Level
4 Stochastic Influences on Phenotype Fingerprints of identical twins Cc, the first cloned cat and her genetic mother J. Raser and E. O Shea, Science, J. Raser and E. O Shea, Science, gen gen gen.. gen variability in gene expression Elowitz et al, Science 2002 Piliated Unpiliated
5 Capturing Randomness in Gene Expression Models Deterministic model protein γ p φ k p mrna γ r φ k r DNA
6 Capturing Randomness in Gene Expression Models Stochastic model protein γ p φ Probability a single mrna is transcribed in time dt is k r dt. Probability a single mrna is degraded in time dt is (#mrna) γ r dt k p mrna γ r φ k r DNA...
7 Fluctuations at Small Copy Numbers protein k p γ p φ Protein Molecules Time (s) mrna DNA k r γ r φ mrna Molecules Time (s)
8 Fluctuations at Small Copy Numbers γ p φ (protein) protein k p mrna γ r φ (mrna) DNA k r C v = coefficient of variation = standard deviation mean
9 Exploiting the Randomness
10 Noise Induced Oscillations Circadian rhythm Vilar, Kueh, Barkai, Leibler, PNAS 2002 Oscillations disappear from deterministic model after a small reduction in deg. of repressor (Coherence resonance) Regularity of noise induced oscillations can be manipulated by tuning the level of noise [El-Samad, Khammash]
11 Stochastic Focusing: Fluctuation Enhanced Sensitivity Signaling Circuit stochastic Number of Molecules of P deterministic 100 P Time mean k s reduced by 50% Stochastic mean value different from deterministic steady state Noise enhances signal! after shift S before shift Johan Paulsson, Otto G. Berg, and Måns Ehrenberg, PNAS 2000
12 Bacterial Competence ¼ Competence is a process by which bacteria takes up foreign DNA Only a fraction of cells become competent SSA runs for 40 hours competence ComS ComK Suel et al, 2006
13 Stochastic Modeling Framework
14 A Simple Example mrna k γ φ mrna copy number N(t) is a random variable Transcription: Probability a single mrna is transcribed in time dt is k r dt DNA Degradation: Probability a single mrna is degraded in time dt is nγdt k k k k k k N n 1 n n γ 2γ 3γ (n 1)γ nγ (n + 1)γ
15 Key Question: k k k k k k n 1 n n + 1 γ 2γ 3γ (n 1)γ nγ (n + 1)γ. Find p(n, t), the probability that N(t) = n. P (n, t + dt) = P (n 1, t) kdt + P (n + 1, t) (n + 1)γdt Prob.{N(t) = n 1 and mrna created in [t,t+dt)} Prob.{N(t) = n + 1 and mrna degraded in [t,t+dt)} + P (n, t) (1 kdt)(1 nγdt) Prob.{N(t) = n and mrna not created nor degraded in [t,t+dt)} P (n, t + dt) P (n, t) = P (n 1, t)kdt + P (n + 1, t)(n + 1)γdt P (n, t)(k + nγ)dt +O(dt 2 ) Dividing by dt and taking the limit as dt 0 The Chemical Master Equation d P (n, t) = kp (n 1, t) + (n + 1)γP (n + 1, t) (k + nγ)p (n, t) dt
16 mrna Stationary Distribution We look for the stationary distribution P (n, t) = p(n) t The stationary solution satisfies: d dtp (n, t) = 0 From the Master Equation... (k + nγ)p(n) = kp(n 1) + (n + 1)γp(n + 1) n = 0 n = 1 n = 2. kp(0) = γp(1) kp(1) = 2γp(2) kp(2) = 3γp(3) kp(n 1) = nγ p(n)
17 kp(n 1) = nγ p(n) p(n) = k γ =. = 1 n 2 k 1 γ k γ We can solve for p(0) using the fact n k 1 1 = γ n! p(0) n=0 We can express p(n) as a function of p(0): p(n 1) n 1 n 1 n 1 n! p(0) p(n 2) n=0 p(n) = 1 = e k/γ p(0) p(0) = e k/γ p(n) = e aan n! a = k γ Poisson Distribution
18 Poisson, a =3!#&$"!#&" probability!#%$"!#%"!#!$"!" mrna count (n) Stationary distribution: P (n) =e aan n! a = k γ Poisson Distribution
19 Formulation of Stochastic Chemical Kinetics Reaction volume=ω Key Assumptions (Well-Mixed) The probability of finding any molecule in a region dω is given by dω Ω. (Thermal Equilibrium) The molecules move due to the thermal energy. The reaction volume is at a constant temperature T. The velocity of a molecule is determined according to a Boltzman distribution: f vx (v) = f vy (v) = f vz (v) = m 2πk B T e m 2k B T v2
20 Stochastic Chemical Kinetics population of S (N-species) S 1,...,S N. Population of each is an integer r.v.: X(t) =[X 1 (t),...,x N (t)] T (M-reactions) The system s state can change through any one of M reaction: R k : k {1, 2,...,M}. (State transition) Firing of reaction R k causes a state transition from X(t) =x to X(t + )=x + s k population of S 1 Stoich. matrix: S = s 1 s M (Transition Probability) The probability that reaction R k fires in the next dt time units is: w k (x)dt. Example: w 1 (x) =c 1 ; w 2 (x) =c 2 x 1 x 2 ; w 3 (x) =c 3 x 1 ;
21 The Chemical Master Equation X(t) is Continuous-time discrete-state Markov Chain p(x, t) := prob(x(t) = x) The Chemical Master Equation (Forward Kolmogorov Equation) dp(x, t) dt = p(x, t) k w k (x)+ k p(x s k,t)w k (x s k )
22 From Stochastic to Deterministic Define X Ω (t) = X(t) Ω. Question: How does X Ω (t) relate to Φ(t)? Fact: Let Φ(t) be the deterministic solution to the reaction rate equations dφ dt = Sf(Φ), Φ(0) = Φ 0. Let X Ω (t) be the stochastic representation of the same chemical systems with X Ω (0) = Φ 0. Then for every t 0: lim sup Ω s t X Ω (s) Φ(s) = 0 a.s.
23 Simulation and Analysis Tools
24 1. Sample Paths Computation Gillespie s Stochastic Simulation Algorithm: To each of the reactions {R 1,...,R M } we associate a RV τ i : τ i is the time to the next firing of reaction R i Fact 0: τ i is exponentially distributed with parameter w i We define two new RVs: τ =min{τ i } i (Time to the next reaction) µ =argmin{τ i } i (Index of the next reaction) Fact 1: τ is exponentially distributed with parameter i w i Fact 2: P (µ = k) = w k i w i
25 Stochastic Simulation Algorithm Step 0 Initialize time t and state population x Step 1 Draw a sample τ from the distribution of τ r 1 U([0, 1]) 1 Cumulative distribution of τ : F (t) =1 exp( k w kt) 0 τ = 1 k w k log 1 1 r 1 time (s) Step 2 Draw a sample µ from the distribution of µ r 2 U([0, 1]) 1 0 Cumulative distribution of µ reaction index µ (w 1 + w 2 + w 3 + w 4 )/ k w k (w 1 + w 2 + w 3 )/ k w k (w 1 + w 2 )/ k w k w 1 / k w k Step 3 Update time: t t + τ. Update state: x x + s µ.
26 2. Moment Computations Let w(x) =[w 1 (x),...,w M (x)] T be the vector of propensity functions Moment Dynamics de[x] dt de[xx T ] dt = S E[w(X)] = SE[w(X)X T ] + E[Xw T (X)]S T + S diag(e[w(x)]) S T Affine propensity. Closed moment equations. Quadratic propensity. Not generally closed. Mass Fluctuation Kinetics (Gomez-Uribe, Verghese) Derivative Matching (Singh, Hespanha)
27 3. SDE Approximation Let X Ω (t) := X(t) Ω Write X Ω = Φ 0 (t) + 1 Ω V Ω where Φ 0 (t) solves the deterministic RRE Linear Noise Approximation dφ dt = Sf(Φ) V Ω (t) V (t) as Ω, where dv (t) = A(t)V (t)dt + B(t)dW t A(t) = d[sf(φ)] dφ (Φ 0(t)), B(t) := S diag[f(φ 0 (t))] Linear Noise Approximation: X Ω (t) Φ(t) + 1 Ω V (t)
28 ω (white gaussian noise) Ẏ = AY + ΩB ω Y (t) = ΩV (t) Ω φ (mean) + X(t)
29 Density Computation Goal: Compute p(x, t), the probability that X(t) = x. Enumerate the state space: X = {x 1,x 2,x 3,...} Form the probability density state vector P (X, ) :R + 1 P (X,t):=[p(x 1,t) p(x 2,t) p(x 3,t)... ] T The Chemical Master Equation (CME): dp(x, t) dt = p(x, t) k w k (x)+ k p(x s k,t)w k (x s k ) can now be written in matrix form: Ṗ (X,t)=A P (X,t)
30 The Finite State Projection Approach
31 The Finite State Projection Approach A finite subset is appropriately chosen
32 The Finite State Projection Approach A finite subset is appropriately chosen The remaining (infinite) states are projected onto a single state (red)
33 The Finite State Projection Approach A finite subset is appropriately chosen The remaining (infinite) states are projected onto a single state (red) Only transitions into removed states are retained The projected system can be solved exactly!
34 Finite Projection Bounds Notation: For a matrix A, leta J to be the principle submatrix of A indexed by J, where J =[m 1...m n ]. Projection Error Bounds Consider any Markov process described by the Forward Kolmogorov Equation: Ṗ (X ; t) =A P (X ; t). If for an indexing vector J: 1 T exp(a J T )P (X J ; 0) 1, then P (XJ ; t) exp(aj t)p (X J ; 0) < t [0,T] P (X J ; t) 0 1 Munsky and Khammash, Journal of Chemical Physics, 2006
35 Example: Analysis of A Synthetic Stochastic Switch Two repressors, u and v. s 2 v s 2 s 2 Gene Promoter s 1 Promoter s 1 Gene Gardner, et al., Nature 403, (2000) s 1 u v inhibits the production of u: a 1 (u, v) = α 1 1+v β ν 1 = [ 1 0 ] u inhibits the production of v: a 3 (u, v) = α 2 1+u γ ν 3 = [ 0 1 ] u and v degrade exponentially: a 2 (u, v) =u a 4 (u, v) =v [ 1 ν 2 = 0 [ 0 ν 4 = 1 ] ] α 1 =50 α 2 =16 β =2.5 γ =1 u(0) = v(0) = 0
36 Using Noise to Identify Model Parameters
37 Why use noise? A?? B 0 hr 3 hr 4 hr 5 hr Noise provides an excitation source for the network dynamics Resulting distributions of proteins can be measured Probability Density q =0 LacI laci Promoter laci laci lac GFP Such distributions q =0 provide a lot of information about the dynamics Can they be used to identify model parameters? q =0 q =1 q =2 q =5 q =2 q =2 q =4 q =2 q =2 q =3 Noise has been used to discriminate among competing models Dunlop et. al (2008). Nature Genetics. Regulatory activity revealed by dynamic correlations in gene expression noise. Total Florescence: GFP + Background (Arbitrary Units)???
38 Identification from Moment Information y y protein γ 2 φ x mrna k 2 γ 1 φ k 1 k 12 y DNA Identifiability Can one identify the parameters λ = {k 1, γ 1,k 2, γ 2,k 21 } from measurements of the moments v(t)?
39 Identifying Using Steady-State Moments y protein γ 2 φ Can the stationary distribution be used to identify all the parameters? x mrna k 2 γ 1 φ k 1 k 12 y DNA Full Identifiability with Stationary Moments Full identifiability is impossible using only v. Munsky et. al, MSB, 2009 Identifiability is possible if lim t E[x(t)x(t + s)] is available. Cinquemani et al, lect. notes comp. sci, 2009
40 Identifiability from Transient Time-Measurements y protein γ 2 φ x mrna k 2 γ 1 φ Multiple Measurements Suppose v j := v(t j ) has been measured at equally separated points in time {t 0,t 1,...,t m } k 1 k 12 y DNA Identifiability with Multiple Moment Measurements For m = 6 the model parameters are identifiable. I 0 A = 1 τ log(g)
41 Identification with Two Measurements x mrna γ φ Identifiability of Transcription Parameters Suppose the mean and variance are known at two times t 0 <t 1 <, and define (µ 0, σ 0 ) := (µ(t 0 ), σ(t 0 )) and (µ 1, σ 1 ) := (µ(t 1 ), σ(t 1 )). DNA k Then the transcription parameters are identifiable, and γ = 1 σ 2 2τ log 1 µ 1 σ0 2 µ k = γ µ 1 exp( γτ)µ exp( γτ) (τ := t 1 t 0 ) protein y γ 2 φ Identifiability of Transcription & Translation Parameters x mrna k 2 γ 1 φ Given v(t 0 ) and v(t 1 ), identifiability of all parameters k 1,k 2, γ 1, γ 2 is generically possible. k 1 An expression exists for finding the parameters. DNA
42 Using Densities to Identify Network Parameters Moment equations can be written only in special cases. Densities (distributions) contain much more information than first two moments. Using the Chemical Master Equation, we propose to use density measurements for model identification. Using Density: Suppose we measure P at different times: P(t 0 ), P(t 1 ),...,P(t N 1 ) We can use these to identify unknown network parameters λ: Find λ subject to Ṗ FSP = A(λ)P FSP P FSP (t 0 ) = P(t 0 ) P FSP (t 1 ) = P(t 1 ). P FSP (t N 1 ) = P(t N 1 )
43 Tsien Lab, UCSD Agar Plate of Fluorescent Bacteria Colonies Flow cytometry protein A histogram joint pdf protein B histogram
44 Identification of lac Induction A Model φ k L B LacI LacI δ L φ Probability Density φ w G GFP δ G φ Experiment laci LacI laci Promoter laci lac 0 hr 3 hr 4 hr 5 hr q =0 GFP w G = q =0 q =1 q =2 q =5 δ L = δ (0) L + δ(1) L [IPTG] IN 1+α[LacI] η, q =2 q =2 q =4 IPTG IN = IPTG OUT (1 e rt ) q =0 q =2 q =2 q =3 k G GFP 9 unknown parameters! E. coli strain Total DL5905 Florescence: GFP + Background (Arbitrary Units) 0 hr 3 hr 4 hr 5 hr Induced Cwith different IPTG concentrations: 5,10, 20, 40, 100 um Induction times: 0, 1, 2, 3, 4, 5 hours before flow cytometry
45 Identified Parameters Identified Model vs. Experiment Model Predictions A B Probability Density C laci can specify the rate δ G = N s, corresponding to a half life of thirty minutes. The remaining seven parameters can then be identified as: k L = s 1 k G =1.0 LacI 10 1 s 1 η =2.1 δ (0) L = N 1 s 1 δ (1) L =5.0 laci Promoter 10 2 (µm N) 1 s 1 α = N η r = s 1 µ GFP = 220 AU σ GFP = 390 AU laci lac GFP where N refers to molecule number. Since the 0 hr assumed hr model 3 hr represents hr a simplified 4 hr hr description 5 hr hr of multiple events (folding dynamics, 10elongation, 2 4 spatial 10 2 motion, 10 4 etc...), 10 2 these 10 4 parameters are 4 best viewed as empirical measurements in q the =0context of q the =1assumed model. q =2Still, it is possible q =5 to make some comparisons 0.1 between the identified parameters and previous analyses. First, the production and degradation rates of 0 LacI yield a mean number of k L /γ (0) L 5 molecules per cell at steady state in the absence of IPTG, on the same magnitude of wild-type levels of about ten per cell. Second, the level 0.1 of LacI q required =0 for half q =2 occupancy of q the =2lac operon q is =4[LacI] 1/2 = (1/α) 1/η =0.012 which compares well to values molecules ( M, Oehler et al., 1990). Third, 0a Hill coefficient of 2.1 is reasonable considering that LacI binds to the operon as a tetramer. Finally, the degradation rate LacI, δ (0) L 0.1 q =0 q =2 q =2 is close to q the =3 dilution rate of , reflecting the high stability of that protein. In addition to comparing the parameters to values in the literature, we have used the parameter set identified from 4, 10, and 20 µm IPTG induction to predict 0 the fluorescence under 40µM IPTG. Fig. 3C shows that these predictions match the subsequent experimental measurements very well despite the vastly different shapes observed at the high Total induction Florescence: levels. GFP + Background (Arbitrary Units) 0 hr hr 3 hr hr 4 hr hr 5 hr hr With single cell experimental techniques such as flow cytometry, it has become possible to efficiently 0.1 measure the fluctuations in cellular species. When properly extracted and processed 0.05with rapidly improving computational tools, these measurements contain sufficiently rich information 0 as to enable the unique identification of parameters. In principle, this can be accomplished when 10 2 accurate distributions are 4 measured at 4 only two distinct time points. More time points are needed if the distributions are poorly measured, but the idea remains the 0.1 same. In this study we have used experimental measurements of cell variability to identify 0.05 the parameters of candidate models, and we have shown that noise and transient dynamics are 0 important to 10 this effort It is10 10 easy 44 to10 10 envision the 4 next iterative step in the scientific process, where one will use these identified models and design the next set of experiments to improve Total Florescence: GFP + Background (Arbitrary Units) the identification. Hence, the proposed exploitation of single cell measurements and stochastic B. Munsky, B. Trinh, M. Khammash, Nature Molecular Systems Biology, in press., Model Experiment
46 Conclusions Randomness noise leads to cell-cell variability Stochastic models are necessary Some stochastic analysis tools available (more needed) Kinetic Monte Carlo Moment approximation Linear noise approximation (van Kampen) Density computation (FSP) Noise reveals network parameters Enabling technologies: flow cytometry and FISH/microscopy A small number of transient measurements suffices FSP exploits full pdf measurements Cellular noise (process noise) vs. measurement noise (output noise)
47 Acknowledgement Brian Munsky, UCSB, LANL Brooke Trinh UCSB (lac induction)
Modeling and Systems Analysis of Gene Regulatory Networks
Modeling and Systems Analysis of Gene Regulatory Networks Mustafa Khammash Center for Control Dynamical-Systems and Computations University of California, Santa Barbara Outline Deterministic A case study:
More information2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements
1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic
More informationStochastic Gene Expression: Modeling, Analysis, and Identification
Stochastic Gene Expression: Modeling, Analysis, and Identification Mustafa Khammash Center for Control, Dynamical Systems, and Computations Department of Mechanical Engineering Engineering II Bldg. Rm.
More informationthe noisy gene Biology of the Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB)
Biology of the the noisy gene Universidad Autónoma de Madrid Jan 2008 Juan F. Poyatos Spanish National Biotechnology Centre (CNB) day III: noisy bacteria - Regulation of noise (B. subtilis) - Intrinsic/Extrinsic
More informationLecture 4: Transcription networks basic concepts
Lecture 4: Transcription networks basic concepts - Activators and repressors - Input functions; Logic input functions; Multidimensional input functions - Dynamics and response time 2.1 Introduction The
More informationA synthetic oscillatory network of transcriptional regulators
A synthetic oscillatory network of transcriptional regulators Michael B. Elowitz & Stanislas Leibler, Nature, 403, 2000 igem Team Heidelberg 2008 Journal Club Andreas Kühne Introduction Networks of interacting
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 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 informationStochastic dynamics of small gene regulation networks. Lev Tsimring BioCircuits Institute University of California, San Diego
Stochastic dynamics of small gene regulation networks Lev Tsimring BioCircuits Institute University of California, San Diego Nizhni Novgorod, June, 2011 Central dogma Activator Gene mrna Protein Repressor
More information56:198:582 Biological Networks Lecture 8
56:198:582 Biological Networks Lecture 8 Course organization Two complementary approaches to modeling and understanding biological networks Constraint-based modeling (Palsson) System-wide Metabolism Steady-state
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 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 informationWhy Should I Care About. Stochastic Hybrid Systems?
Research Supported by NSF & ARO Why Should I Care About Stochastic Hybrid Systems? João Hespanha 1 Why Care Should I Care About SHSs? Eamples Feedback over shared networks Estimation using remote sensors
More informationA Synthetic Oscillatory Network of Transcriptional Regulators
A Synthetic Oscillatory Network of Transcriptional Regulators Michael Elowitz & Stanislas Leibler Nature, 2000 Presented by Khaled A. Rahman Background Genetic Networks Gene X Operator Operator Gene Y
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 informationAnalog Electronics Mimic Genetic Biochemical Reactions in Living Cells
Analog Electronics Mimic Genetic Biochemical Reactions in Living Cells Dr. Ramez Daniel Laboratory of Synthetic Biology & Bioelectronics (LSB 2 ) Biomedical Engineering, Technion May 9, 2016 Cytomorphic
More 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 simulations
Stochastic simulations Application to molecular networks Literature overview Noise in genetic networks Origins How to measure and distinguish between the two types of noise (intrinsic vs extrinsic)? What
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech DRAFT v0.4a, January 16, 2011 c California Institute of Technology All rights reserved. This manuscript is for review purposes
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 informationStochastic Simulation of Biochemical Reactions
1 / 75 Stochastic Simulation of Biochemical Reactions Jorge Júlvez University of Zaragoza 2 / 75 Outline 1 Biochemical Kinetics 2 Reaction Rate Equation 3 Chemical Master Equation 4 Stochastic Simulation
More informationBioinformatics: Network Analysis
Bioinformatics: Network Analysis Kinetics of Gene Regulation COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 The simplest model of gene expression involves only two steps: the
More informationSTOCHASTIC MODELING OF BIOCHEMICAL REACTIONS
STOCHASTIC MODELING OF BIOCHEMICAL REACTIONS Abhyudai Singh and João Pedro Hespanha* Department of Electrical and Computer Engineering University of California, Santa Barbara, CA 93101. Abstract The most
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 informationSTOCHASTIC REDUCTION METHOD FOR BIOLOGICAL CHEMICAL KINETICS USING TIME-SCALE SEPARATION
STOCHASTIC REDUCTION METHOD FOR BIOLOGICAL CHEMICAL KINETICS USING TIME-SCALE SEPARATION CHETAN D. PAHLAJANI 12, PAUL J. ATZBERGER 23, AND MUSTAFA KHAMMASH 3 Abstract. Many processes in cell biology encode
More informationWhen do diffusion-limited trajectories become memoryless?
When do diffusion-limited trajectories become memoryless? Maciej Dobrzyński CWI (Center for Mathematics and Computer Science) Kruislaan 413, 1098 SJ Amsterdam, The Netherlands Abstract Stochastic description
More informationBoulder 07: Modelling Stochastic Gene Expression
Boulder 7: Modelling Stochastic Gene Expression Peter Swain, Physiology, McGill University swain@cnd.mcgill.ca The whys and wherefores of stochasticity A system evolves stochastically if its dynamics is
More informationSTOCHASTIC CHEMICAL KINETICS
STOCHASTIC CHEICAL KINETICS Dan Gillespie GillespieDT@mailaps.org Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC))
More informationEngineering of Repressilator
The Utilization of Monte Carlo Techniques to Simulate the Repressilator: A Cyclic Oscillatory System Seetal Erramilli 1,2 and Joel R. Stiles 1,3,4 1 Bioengineering and Bioinformatics Summer Institute,
More informationGLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data
GLOBEX Bioinformatics (Summer 2015) Genetic networks and gene expression data 1 Gene Networks Definition: A gene network is a set of molecular components, such as genes and proteins, and interactions between
More informationBi 1x Spring 2014: LacI Titration
Bi 1x Spring 2014: LacI Titration 1 Overview In this experiment, you will measure the effect of various mutated LacI repressor ribosome binding sites in an E. coli cell by measuring the expression of a
More informationStochastic simulations!
Stochastic simulations! Application to biomolecular networks! Literature overview Noise in genetic networks! Origins! How to measure the noise and distinguish between the two sources of noise (intrinsic
More informationURL: <
Citation: ngelova, Maia and en Halim, sma () Dynamic model of gene regulation for the lac operon. Journal of Physics: Conference Series, 86 (). ISSN 7-696 Published by: IOP Publishing URL: http://dx.doi.org/.88/7-696/86//7
More informationNumerical Solution of the Fokker- Planck Approximation of the Chemical Master Equation
IT Licentiate theses 2005-010 Numerical Solution of the Fokker- Planck Approximation of the Chemical Master Equation PAUL SJÖBERG UPPSALA UNIVERSITY Department of Information Technology Numerical Solution
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 informationStochastically driven genetic circuits
Stochastically driven genetic circuits CHAOS 16, 026103 2006 L. S. Tsimring Institute for Nonlinear Science, University of California, San Diego, La Jolla, California 92093-0402 D. Volfson Institute for
More informationStochastic Processes around Central Dogma
Stochastic Processes around Central Dogma Hao Ge haoge@pu.edu.cn Beijing International Center for Mathematical Research Biodynamic Optical Imaging Center Peing University, China http://www.bicmr.org/personal/gehao/
More informationBiomolecular Feedback Systems
Biomolecular Feedback Systems Domitilla Del Vecchio MIT Richard M. Murray Caltech Version 1.0b, September 14, 2014 c 2014 by Princeton University Press All rights reserved. This is the electronic edition
More informationProblem Set 5. 1 Waiting times for chemical reactions (8 points)
Problem Set 5 1 Waiting times for chemical reactions (8 points) In the previous assignment, we saw that for a chemical reaction occurring at rate r, the distribution of waiting times τ between reaction
More informationTopic 4 - #14 The Lactose Operon
Topic 4 - #14 The Lactose Operon The Lactose Operon The lactose operon is an operon which is responsible for the transport and metabolism of the sugar lactose in E. coli. - Lactose is one of many organic
More informationBi 8 Lecture 11. Quantitative aspects of transcription factor binding and gene regulatory circuit design. Ellen Rothenberg 9 February 2016
Bi 8 Lecture 11 Quantitative aspects of transcription factor binding and gene regulatory circuit design Ellen Rothenberg 9 February 2016 Major take-home messages from λ phage system that apply to many
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process https://www.youtube.com/
More informationA(γ A D A + γ R D R + γ C R + δ A )
Title: Mechanisms of noise-resistance in genetic oscillators Authors: José M. G. Vilar 1,2, Hao Yuan Kueh 1, Naama Barkai 3, and Stanislas Leibler 1,2 1 Howard Hughes Medical Institute, Departments of
More informationBistability and a differential equation model of the lac operon
Bistability and a differential equation model of the lac operon Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Fall 2016 M. Macauley
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 informationSystem Biology - Deterministic & Stochastic Dynamical Systems
System Biology - Deterministic & Stochastic Dynamical Systems System Biology - Deterministic & Stochastic Dynamical Systems 1 The Cell System Biology - Deterministic & Stochastic Dynamical Systems 2 The
More information7.32/7.81J/8.591J: Systems Biology. Fall Exam #1
7.32/7.81J/8.591J: Systems Biology Fall 2013 Exam #1 Instructions 1) Please do not open exam until instructed to do so. 2) This exam is closed- book and closed- notes. 3) Please do all problems. 4) Use
More informationChapter 15 Active Reading Guide Regulation of Gene Expression
Name: AP Biology Mr. Croft Chapter 15 Active Reading Guide Regulation of Gene Expression The overview for Chapter 15 introduces the idea that while all cells of an organism have all genes in the genome,
More informationReaction time distributions in chemical kinetics: Oscillations and other weird behaviors
Introduction The algorithm Results Summary Reaction time distributions in chemical kinetics: Oscillations and other weird behaviors Ramon Xulvi-Brunet Escuela Politécnica Nacional Outline Introduction
More informationMultistability in the lactose utilization network of Escherichia coli
Multistability in the lactose utilization network of Escherichia coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Agenda Motivation Intro to multistability Purpose
More informationCHAPTER 13 PROKARYOTE GENES: E. COLI LAC OPERON
PROKARYOTE GENES: E. COLI LAC OPERON CHAPTER 13 CHAPTER 13 PROKARYOTE GENES: E. COLI LAC OPERON Figure 1. Electron micrograph of growing E. coli. Some show the constriction at the location where daughter
More informationTesting the transition state theory in stochastic dynamics of a. genetic switch
Testing the transition state theory in stochastic dynamics of a genetic switch Tomohiro Ushikubo 1, Wataru Inoue, Mitsumasa Yoda 1 1,, 3, and Masaki Sasai 1 Department of Computational Science and Engineering,
More informationIntroduction. Gene expression is the combined process of :
1 To know and explain: Regulation of Bacterial Gene Expression Constitutive ( house keeping) vs. Controllable genes OPERON structure and its role in gene regulation Regulation of Eukaryotic Gene Expression
More informationBistability in ODE and Boolean network models
Bistability in ODE and Boolean network models Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson)
More informationBiomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems
Biomolecular reaction networks: gene regulation & the Repressilator Phys 682/ CIS 629: Computational Methods for Nonlinear Systems Processing information to coordinate activity regulatory & signaling networks
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 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 informationTopic 4: Equilibrium binding and chemical kinetics
Topic 4: Equilibrium binding and chemical kinetics Outline: Applications, applications, applications use Boltzmann to look at receptor-ligand binding use Boltzmann to look at PolII-DNA binding and gene
More informationModeling Multiple Steady States in Genetic Regulatory Networks. Khang Tran. problem.
Modeling Multiple Steady States in Genetic Regulatory Networks Khang Tran From networks of simple regulatory elements, scientists have shown some simple circuits such as the circadian oscillator 1 or the
More informationarxiv: v2 [q-bio.qm] 12 Jan 2017
Approximation and inference methods for stochastic biochemical kinetics - a tutorial review arxiv:1608.06582v2 [q-bio.qm] 12 Jan 2017 David Schnoerr 1,2,3, Guido Sanguinetti 2,3, and Ramon Grima 1,3,*
More informationStochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology
research supported by NSF Stochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology João P. Hespanha Center for Control Engineering and Computation University of California
More informationGene regulation I Biochemistry 302. Bob Kelm February 25, 2005
Gene regulation I Biochemistry 302 Bob Kelm February 25, 2005 Principles of gene regulation (cellular versus molecular level) Extracellular signals Chemical (e.g. hormones, growth factors) Environmental
More informationMultifidelity Approaches to Approximate Bayesian Computation
Multifidelity Approaches to Approximate Bayesian Computation Thomas P. Prescott Wolfson Centre for Mathematical Biology University of Oxford Banff International Research Station 11th 16th November 2018
More informationLecture 1 Modeling in Biology: an introduction
Lecture 1 in Biology: an introduction Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it
More informationHow to Build a Living Cell in Software or Can we computerize a bacterium?
How to Build a Living Cell in Software or Can we computerize a bacterium? Tom Henzinger IST Austria Turing Test for E. coli Fictional ultra-high resolution video showing molecular processes inside the
More informationEuropean Conference on Mathematical and Theoretical Biology Gothenburg
European Conference on Mathematical and Theoretical Biology 2014 - Gothenburg Minisymposium: Numerical methods for high-dimensional problems in biology Date: Monday, 16 th of June Time: 16:00-19:00 Organizers:
More informationGene Regulation and Expression
THINK ABOUT IT Think of a library filled with how-to books. Would you ever need to use all of those books at the same time? Of course not. Now picture a tiny bacterium that contains more than 4000 genes.
More informationName: SBI 4U. Gene Expression Quiz. Overall Expectation:
Gene Expression Quiz Overall Expectation: - Demonstrate an understanding of concepts related to molecular genetics, and how genetic modification is applied in industry and agriculture Specific Expectation(s):
More informationNoise is often perceived as being undesirable and unpredictable;
Intrinsic noise in gene regulatory networks Mukund Thattai and Alexander van Oudenaarden* Department of Physics, Room 13-2010, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge,
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationGenetic transcription and regulation
Genetic transcription and regulation Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process DNA codes for DNA replication
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 informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and
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 informationSimulation of Chemical Reactions
Simulation of Chemical Reactions Cameron Finucane & Alejandro Ribeiro Dept. of Electrical and Systems Engineering University of Pennsylvania aribeiro@seas.upenn.edu http://www.seas.upenn.edu/users/~aribeiro/
More informationBrief contents. Chapter 1 Virus Dynamics 33. Chapter 2 Physics and Biology 52. Randomness in Biology. Chapter 3 Discrete Randomness 59
Brief contents I First Steps Chapter 1 Virus Dynamics 33 Chapter 2 Physics and Biology 52 II Randomness in Biology Chapter 3 Discrete Randomness 59 Chapter 4 Some Useful Discrete Distributions 96 Chapter
More information3.B.1 Gene Regulation. Gene regulation results in differential gene expression, leading to cell specialization.
3.B.1 Gene Regulation Gene regulation results in differential gene expression, leading to cell specialization. We will focus on gene regulation in prokaryotes first. Gene regulation accounts for some of
More informationStochastic Process Algebra models of a Circadian Clock
Stochastic Process Algebra models of a Circadian Clock Jeremy T. Bradley Thomas Thorne Department of Computing, Imperial College London 180 Queen s Gate, London SW7 2BZ, United Kingdom Email: jb@doc.ic.ac.uk
More informationComplete all warm up questions Focus on operon functioning we will be creating operon models on Monday
Complete all warm up questions Focus on operon functioning we will be creating operon models on Monday 1. What is the Central Dogma? 2. How does prokaryotic DNA compare to eukaryotic DNA? 3. How is DNA
More informationRegulation of Gene Expression
Chapter 18 Regulation of Gene Expression PowerPoint Lecture Presentations for Biology Eighth Edition Neil Campbell and Jane Reece Lectures by Chris Romero, updated by Erin Barley with contributions from
More informationLecture 4 The stochastic ingredient
Lecture 4 The stochastic ingredient Luca Bortolussi 1 Alberto Policriti 2 1 Dipartimento di Matematica ed Informatica Università degli studi di Trieste Via Valerio 12/a, 34100 Trieste. luca@dmi.units.it
More informationStochastic simulations Application to molecular networks
Stochastic simulations Application to molecular networks Didier Gonze May 1, 27 Les Houches - Spring School - April 9-2, 27 Contents 1 Introduction 3 1.1 Noise in biology.................................
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 informationBioControl - Week 6, Lecture 1
BioControl - Week 6, Lecture 1 Goals of this lecture Large metabolic networks organization Design principles for small genetic modules - Rules based on gene demand - Rules based on error minimization Suggested
More informationStochastic Processes around Central Dogma
Stochastic Processes around Central Dogma Hao Ge haoge@pku.edu.cn Beijing International Center for Mathematical Research Biodynamic Optical Imaging Center Peking University, China http://www.bicmr.org/personal/gehao/
More informationModelling Stochastic Gene Expression
Modelling Stochastic Gene Expression Peter Swain Centre for Systems Biology at Edinburgh University of Edinburgh peter.swain@ed.ac.uk The whys and wherefores of stochasticity A system evolves stochastically
More informationarxiv: v2 [stat.me] 11 Nov 2017
arxiv:1610.07213v2 [stat.me] 11 Nov 2017 Stochastic Modeling and Statistical Inference of Intrinsic Noise in Gene Regulation System via the Chemical Master Equation Chao Du, and Wing Hong Wong November
More informationSlide 1 / 7. Free Response
Slide 1 / 7 Free Response Slide 2 / 7 Slide 3 / 7 1 The above diagrams illustrate the experiments carried out by Griffith and Hershey and Chaserespectively. Describe the hypothesis or conclusion that each
More informationProkaryotic Gene Expression (Learning Objectives)
Prokaryotic Gene Expression (Learning Objectives) 1. Learn how bacteria respond to changes of metabolites in their environment: short-term and longer-term. 2. Compare and contrast transcriptional control
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationnatural development from this collection of knowledge: it is more reliable to predict the property
1 Chapter 1 Introduction As the basis of all life phenomena, the interaction of biomolecules has been under the scrutiny of scientists and cataloged meticulously [2]. The recent advent of systems biology
More informationSimulating ribosome biogenesis in replicating whole cells
Simulating ribosome biogenesis in replicating whole cells Tyler Earnest / PI: Zan Luthey-Schulten (UIUC) NCSA Blue Waters Symposium for Petascale Science and Beyond: June 15, 2016 Multi-scale modeling
More information13.4 Gene Regulation and Expression
13.4 Gene Regulation and Expression Lesson Objectives Describe gene regulation in prokaryotes. Explain how most eukaryotic genes are regulated. Relate gene regulation to development in multicellular organisms.
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
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 informationKinetic Monte Carlo. Heiko Rieger. Theoretical Physics Saarland University Saarbrücken, Germany
Kinetic Monte Carlo Heiko Rieger Theoretical Physics Saarland University Saarbrücken, Germany DPG school on Efficient Algorithms in Computational Physics, 10.-14.9.2012, Bad Honnef Intro Kinetic Monte
More informationVilla et al. (2005) Structural dynamics of the lac repressor-dna complex revealed by a multiscale simulation. PNAS 102:
Villa et al. (2005) Structural dynamics of the lac repressor-dna complex revealed by a multiscale simulation. PNAS 102: 6783-6788. Background: The lac operon is a cluster of genes in the E. coli genome
More informationA Simple Protein Synthesis Model
A Simple Protein Synthesis Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 3, 213 Outline A Simple Protein Synthesis Model
More informationUNIT 6 PART 3 *REGULATION USING OPERONS* Hillis Textbook, CH 11
UNIT 6 PART 3 *REGULATION USING OPERONS* Hillis Textbook, CH 11 REVIEW: Signals that Start and Stop Transcription and Translation BUT, HOW DO CELLS CONTROL WHICH GENES ARE EXPRESSED AND WHEN? First of
More informationMechanisms for Precise Positional Information in Bacteria: The Min system in E. coli and B. subtilis
Mechanisms for Precise Positional Information in Bacteria: The Min system in E. coli and B. subtilis Martin Howard Imperial College London Bacterial Organization Many processes where bacterial cell needs
More information