Top Banner
Stochastic Analysis of Chemical Reaction Networks using Linear Noise Approximation Luca Cardelli1,2 , Marta Kwiatkowsk1 , Luca Laurent 1 Department of Computer Science, University of Oxford 2 Microsoft Research 1
18

Stochastic Analysis of Chemical Reaction Networks using ...

Jan 27, 2022

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Stochastic Analysis of Chemical Reaction Networks using ...

Stochastic Analysis of Chemical Reaction Networks using Linear Noise

Approximation

Luca Cardell i↑1,2     , Marta Kwiatkowsk 𝑎↑1 , Luca Laurent𝒊↑𝟏  ↑1   Department of Computer Science, University of Oxford

↑2   Microsoft Research

1

Page 2: Stochastic Analysis of Chemical Reaction Networks using ...

•  Motivation•  Background•  Linear Noise Approximation (LNA)•  Stochastic Evolution Logic (SEL)•  Experimental Results

Summary

2

Page 3: Stochastic Analysis of Chemical Reaction Networks using ...

Motivation Ø Biochemical systems are generally analysed considering

deterministic models. However, deterministic models are accurate only when the molecular population is large

Ø When the interacting entities are in low numbers there is a

need of considering a stochastic model

Ø Existing methods for analysis of discrete state space stochastic processes are not scalable and highly dependent on the initial number of molecules

•  Question: Can we derive a formal method to analyse the stochastic semantics of biochemical systems that is scalable and independent of the initial number of molecules?

3

Page 4: Stochastic Analysis of Chemical Reaction Networks using ...

Chemical Reaction Networks (CRNs)

•  A CRN 𝑪=(𝚲,𝑹) is a pair of sets

•  Λ is a finite set of species { 𝜆↓1 ,   𝜆↓2 ,…, 𝜆↓|Λ| }  

•  𝑅 is a finite set of reactions {  𝜏↓1 ,   𝜏↓2 ,…, 𝜏↓|𝑅| }  •  𝑒.𝑔.   𝜏↓𝑖   :   𝜆↓1 + 𝜆↓2    →↑𝑘      𝜆↓3 •  𝑘 is the rate constant

Ø A configuration or state of the system, 𝑥∈ 𝑁↓↑|Λ| , is given by the number of molecules of each species in that configuration

4

Page 5: Stochastic Analysis of Chemical Reaction Networks using ...

Ø Set of autonomous polynomial ODEs:𝑑Φ/𝑑𝑡 =F(t)Φ(𝑡)                                                                              Φ(0)=   𝑥↓0 /𝑁  •  𝑁=𝑉𝑜𝑙𝑢𝑚𝑒⋅ 𝑁↓𝐴  is the Volumetric factor of the

system

•  Φ(𝑡)∈ 𝑅↑|Λ|    represents the species concentration at time t

•  𝐹(𝑡) is determined by mass action kinetics•  Number of differential equations equals number of species•  Valid only for high number of molecules•  Does not take into account the stochastic nature of

molecular interactions 5

Deterministic Semantics of CRNs

Page 6: Stochastic Analysis of Chemical Reaction Networks using ...

•  It is a continous time Markov process ( 𝑋↑𝑁 (𝑡),𝑡≥0) with discrete state space 𝑆 and infinitesimal generator matrix 𝑄 determined by the reactions

•  The transient evoloution of 𝑋↑𝑁   is described by the Chemical Master Equation (CME)o Assuming 𝑃(𝑥,𝑡)=𝑃𝑟𝑜𝑏{𝑋↑𝑁 (𝑡)=𝑥  |   𝑋↑𝑁 (0)= 𝑥↓0 }

then the CME can be written as

𝑑𝑃(𝑥,𝑡)/𝑑𝑡 =𝑃(𝑥,𝑡)𝑄

6

Stochastic Semantics of CRNs

Page 7: Stochastic Analysis of Chemical Reaction Networks using ...

•  One differential equation in the CME for any reachable state•  S highly dependent on the initial number of molecules•  Set of reachable states can be huge or even infinite

•  Solution: USE LINEAR NOISE APPROXIMATION !

7

State Space Explosion Problem

Ø Not possible to solve the CME for large molecular populations and/or large CRNs

Page 8: Stochastic Analysis of Chemical Reaction Networks using ...

•  Technique pioneered by Van Kampen in his CME espansion

8

•  Φ(𝑡) solution of the deterministic semantics•  𝑍(𝑡) is a Gaussian Process independent of 𝑁

o  𝐸[𝑍(𝑡)]=0    𝑓𝑜𝑟    𝑡≥0o  𝑑𝐶[𝑍(𝑡)]/𝑑𝑡 = 𝐽↓𝐹 (Φ(𝑡))𝐶[𝑍(𝑡)]+𝐶[𝑍(𝑡)]𝐽↓𝐹↑𝑇 (Φ(𝑡))+𝐺(Φ(𝑡))•   𝐽↓𝐹 (Φ(t)) Jacobian of 𝐹(Φ(𝑡)) •  𝐺= 1/𝑁 ∑𝜏∈𝑅↑▒𝜐↓𝜏 𝜐↓𝜏↑𝑇 𝛼↓𝜏  

Linear Noise Approximation (LNA)

𝑋↑𝑁 (𝑡)≈𝑌↑𝑁 (𝑡)=𝑁Φ(𝑡)+√𝑁 𝑍(𝑡)

Page 9: Stochastic Analysis of Chemical Reaction Networks using ...

•  For any CRN, assuming mass action kinetics, the LNA is always accurate at least for a limited time (it is enough to increase 𝑁)

•  Independence of the initial number of moleculeso  The number of differential equations depends only on the

number of species

•  Number of differential equations quadratic in the number of species

•  Still good approximation for a large class of CRNs even for quite small molecular populations o  Not able to handle multinomial distributions

9

Linear Noise Approximation (LNA)

Page 10: Stochastic Analysis of Chemical Reaction Networks using ...

LNA Also Known as Gaussian Approximation

•  𝐵∈ 𝑁↓↑|Λ| , the linear combination of species 𝐵↑𝑇 𝑌↑𝑁 (𝑡) is still Gaussian

o  𝐸[𝐵↑𝑇 𝑌↑𝑁 (𝑡)]= 𝐵↑𝑇 𝐸[𝑌↑𝑁 (𝑡)]=𝑁(𝐵↑T Φ(𝑡))

o  𝐶[𝐵↑𝑇 𝑌↑𝑁 (𝑡)]=𝐵𝐶[𝑌↑𝑁 (𝑡)]𝐵↑𝑇 =𝐵𝐶[𝑍(𝑡)]𝐵↑𝑇 

•  Probability is calculated by solving Gaussian integrals

•  𝑍(𝑡)  is a Gaussian process

Ø  𝑌↑𝑁 (𝑡)=𝑁Φ(𝑡)+√𝑁 𝑍(𝑡) is a Gaussian process

10

Page 11: Stochastic Analysis of Chemical Reaction Networks using ...

•  𝑄={  𝑠𝑢𝑝𝑉,  𝑖𝑛𝑓𝑉,𝑠𝑢𝑝𝐸,𝑖𝑛𝑓𝐸} •  𝐼 set of closed disjoint intervals •  𝐵∈ 𝑍↓≥0  ↑|Λ| 

11

Stochastic Evolution Logic (SEL)

Ø  𝑃↓∼𝑝 [𝐵,𝐼  ]↓[𝑡↓1 , 𝑡↓2 ]    : probabilistic operator

Ø  𝑠𝑢𝑝𝑉/𝑖𝑛𝑓𝑉↓∼𝑣 [𝐵]↓[𝑡↓1 , 𝑡↓2 ]   : supremum/infimum of variance operators

Ø  𝑠𝑢𝑝𝐸/𝑖𝑛𝑓𝐸↓∼𝑣 [𝐵]↓[𝑡↓1 , 𝑡↓2 ]   : supremum/infimum of expected value operators

Page 12: Stochastic Analysis of Chemical Reaction Networks using ...

12

Semantics (SEL)

Page 13: Stochastic Analysis of Chemical Reaction Networks using ...

13

Ø We use the LNA for a numerical approximate model checking algorithm of SEL

o  𝑋↑𝑁  is approximated by the Gaussian Process 𝑌↑𝑁 

o  The probability that 𝐵↑𝑇 𝑌↑𝑁  is within the interval [𝑙,𝑟] at time t is:

∫𝑙↑𝑟▒𝑔𝑥 𝐸[𝐵↑𝑇 𝑌↑𝑁 (𝑡)],𝐶[𝐵↑𝑇 𝑌↑𝑁 (𝑡)] 𝑑𝑡  where 𝑔(𝑥|𝐸,𝐶) is the Gaussian distribution with

expected value 𝐸  and variance 𝐶.

o  𝐸[𝐵↑𝑇 𝑌↑𝑁 (𝑡)] and C[𝐵↑𝑇 𝑌↑𝑁 (𝑡)] are obtained by solving the LNA for the given initial condition

Approximate Model Checking

Page 14: Stochastic Analysis of Chemical Reaction Networks using ...

Property to check:

Species = {𝐿1,𝐿1𝑝,𝐿2,𝐿2𝑝,𝐿3,𝐿3𝑝,𝐵} 𝜏↓1 :  𝐿1+𝐵   →↓↑𝑘↓1  𝐵+𝐿1𝑝 𝜏↓2 :  𝐿1𝑝+𝐿2   →↑𝑘↓2    𝐿1+𝐿2𝑝 𝜏↓3 :  𝐿2𝑝+𝐿3   →↑𝑘↓2    𝐿2+𝐿3𝑝 𝜏↓4 :  𝐿3𝑝→↑𝑘↓3    𝐿3

Comparison with standard Uniformization

Initial Condition: x↓0 (𝐿1)= 𝑥↓0 (𝐿2)= 𝑥↓0 (𝐿3)=𝐼𝑛𝑖𝑡;   𝑥↓0 (𝐵)=3⋅𝐼𝑛𝑖𝑡;             where  𝐼𝑛𝑖𝑡 is a variable with values in 𝑁

14

Phosphorelay Network

Page 15: Stochastic Analysis of Chemical Reaction Networks using ...

•  CRN composed of more than 50 reactions and species!o  Initial condition such that all species with non zero concentration

have 105 moleculeso  Exploration of state space infeasibleo  Simulations are time consuming for such a biochemical system

𝑠𝑢𝑝𝐸↓=? [#𝑠𝑟𝑐:𝐹𝑅𝑆2]↓[𝑇,𝑇]  𝑠𝑢𝑝𝑉↓=? [#𝑆𝑟𝑐:𝐹𝑅𝑆2]↓[𝑇,𝑇]       For    𝑇∈[0,8000]

Comparison of SEL and single stochastic simulation

15

FGF Pathways

Page 16: Stochastic Analysis of Chemical Reaction Networks using ...

•  We have presented SEL with an approximate model checking algorithm based on the LNA

•  Our method can be useful for a fast stochastic characterization of biochemical systems or for stochastic analysis of systems too large to be checked with standard techniques.

•  Increasing the number of molecules the LNA is always a valid model assuming mass action kinetics, but can be accurate even far from the thermodynamic limit for a large class of CRNs

16

Conclusion

Page 17: Stochastic Analysis of Chemical Reaction Networks using ...

•  Wallace, E. W. J., et al. "Linear noise approximation is valid over limited times for any chemical system that is sufficiently large." IET systems biology 6.4 (2012): 102-115.

•  Van Kampen, Nicolaas Godfried. Stochastic processes in physics and chemistry. Vol. 1. Elsevier, 1992. e

•  Gillespie, Daniel T. Deterministic limit of stochastic chemical kinetics. The Journal of Physical Chemistry B 113.6 (2009): 1640-1644.

•  Luca Cardelli, Marta Kwiatkowska, Luca Laurenti. Stochastic Analysis of Chemical Reaction Networks Using Linear Noise Approximation.  arXiv preprint arXiv:1506.07861 (2015).

17

Some References

Page 18: Stochastic Analysis of Chemical Reaction Networks using ...

THANK YOU!

18