It looks like you're new here. If you want to get involved, click one of these buttons!

- All Categories 2.3K
- Chat 501
- Study Groups 21
- Petri Nets 9
- Epidemiology 4
- Leaf Modeling 2
- Review Sections 9
- MIT 2020: Programming with Categories 51
- MIT 2020: Lectures 20
- MIT 2020: Exercises 25
- Baez ACT 2019: Online Course 339
- Baez ACT 2019: Lectures 79
- Baez ACT 2019: Exercises 149
- Baez ACT 2019: Chat 50
- UCR ACT Seminar 4
- General 70
- Azimuth Code Project 110
- Statistical methods 4
- Drafts 9
- Math Syntax Demos 15
- Wiki - Latest Changes 3
- Strategy 113
- Azimuth Project 1.1K
- - Spam 1
- News and Information 148
- Azimuth Blog 149
- - Conventions and Policies 21
- - Questions 43
- Azimuth Wiki 715

Options

From:

- Darren J. Wilkinson,
*Stochastic Modelling for Systems Biology*, Taylor & Francis, New York, 2006. Good introduction to stochastic Petri nets, with applications to gene expression and the Lotka-Volterra equations for predator-prey interactions.

Chapter on case studies has the following applications, with simulation lab exercises:

- Dimerization kinetics
- Michaelis-Menten enzyme kinetics
- An auto-regulatory genetic network
- The lac operon

In this thread we'll look into the one on dimerization kinetics.

## Comments

We'll look at the homodimer case, which is just the reaction \(X + X = 2X \rightarrow X_2\).

For instance, the potassium dimer \(K_2\) is a molecule consisting of two potassium atoms, and the dimerization of potassium would be the reaction where two \(K\) atoms combine into one \(K_2\) molecule: \(2K \rightarrow K_2\).

The reverse reaction, where the dimer splits into two monomers, is called dissociation: \(X_2 \rightarrow 2X\).

`* [Dimer](https://en.wikipedia.org/wiki/Dimer_(chemistry)), Wikipedia. > A dimer is an oligomer consisting of two monomers joined by bonds that can be either strong or weak, covalent or intermolecular. The term homodimer is used when the two molecules are identical and heterodimer when they are not. The reverse of dimerization is often called dissociation. We'll look at the homodimer case, which is just the reaction \\(X + X = 2X \rightarrow X_2\\). For instance, the potassium dimer \\(K_2\\) is a molecule consisting of two potassium atoms, and the dimerization of potassium would be the reaction where two \\(K\\) atoms combine into one \\(K_2\\) molecule: \\(2K \rightarrow K_2\\). The reverse reaction, where the dimer splits into two monomers, is called dissociation: \\(X_2 \rightarrow 2X\\).`

For this lab, the reaction system consists of two reactions, dimerization and dissociation:

Let:

`For this lab, the reaction system consists of two reactions, dimerization and dissociation: * \\(2X \rightarrow X_2\\), with rate coefficient \\(k_1\\) = 50,000 * \\(X_2 \rightarrow 2X\\), with rate coefficient \\(k_2\\) = 0.2 Let: * \\(X[t]\\) be the count of the number of monomers \\(X\\) at time \\(t\\) * \\(X_2[t]\\) be the count of the number of dimers \\(X_2\\) at time \\(t\\)`

So we have a "tug of war" between the two reactions, with dimerization working to build up the count of dimers \(X_2[t]\) and reduce the count of monomers \(X[t]\) - and dissociation doing just the opposite.

`So we have a "tug of war" between the two reactions, with dimerization working to build up the count of dimers \\(X_2[t]\\) and reduce the count of monomers \\(X[t]\\) - and dissociation doing just the opposite.`

The

dynamicsof the system is expressed by the two time series \(X[t]\) and \(X_2[t]\).Intuitively, we would expect that given a starting state \(S[0] = (X[0], X_2[0])\), it will eventually approach some limiting state \(S[{\infty}] = (X', X_2')\) as \(t\) advances - that is the

equilibriumstate towards which it is bound.`The _dynamics_ of the system is expressed by the two time series \\(X[t]\\) and \\(X_2[t]\\). Intuitively, we would expect that given a starting state \\(S[0] = (X[0], X_2[0])\\), it will eventually approach some limiting state \\(S[{\infty}] = (X', X_2')\\) as \\(t\\) advances - that is the _equilibrium_ state towards which it is bound.`

The proportion between the two species in the equilibrium state will depend on the relative strengths of the two reactions - see the reaction coefficients. If dissociation had a coefficient of zero, it would be effectively gone, and the equilibrium state would consist 100% of dimers 0% of monomers. And similarly if it was all dissociation equilibrium would be at 100% momers.

The rate coefficients given above give a might higher strength to dimerization, and so we would expect equilibrium to consist

mostlyof dimers.`The proportion between the two species in the equilibrium state will depend on the relative strengths of the two reactions - see the reaction coefficients. If dissociation had a coefficient of zero, it would be effectively gone, and the equilibrium state would consist 100% of dimers 0% of monomers. And similarly if it was all dissociation equilibrium would be at 100% momers. The rate coefficients given above give a might higher strength to dimerization, and so we would expect equilibrium to consist _mostly_ of dimers.`

Our goal is now to understand the dynamics of the system. This means understanding the evolution of the state function S[t].

`Our goal is now to understand the dynamics of the system. This means understanding the evolution of the state function S[t].`

We have a choice on the menu: deterministic or stochastic modeling of S[t].

In a deterministic model, S[t] is a definite function of the continuous parameter \(t\). Once the initial condition is specified, by giving a value for S[0], then the future values of S[t] are uniquely determined.

This is a "flow." (Semi-flow actually, as S[t] is defined for non-negative t.)

`We have a choice on the menu: deterministic or stochastic modeling of S[t]. In a deterministic model, S[t] is a definite function of the continuous parameter \\(t\\). Once the initial condition is specified, by giving a value for S[0], then the future values of S[t] are uniquely determined. This is a "flow." (Semi-flow actually, as S[t] is defined for non-negative t.)`

In the stochastic model, \(t\) is still a continuous parameter, but the state S[t] makes little "jumps" at random times, when one of the reactions randomly "fires." Every time the dimerization reaction fires, the count of monomers goes down by two, and the count of dimers goes up by one. When the opposite reaction viz. dissociation fires, the dimer count goes down by one and the monomer count goes up by two.

`In the stochastic model, \\(t\\) is still a continuous parameter, but the state S[t] makes little "jumps" at random times, when one of the reactions randomly "fires." Every time the dimerization reaction fires, the count of monomers goes down by two, and the count of dimers goes up by one. When the opposite reaction viz. dissociation fires, the dimer count goes down by one and the monomer count goes up by two.`

Where do the reaction coefficients show up in these models?

In the deterministic case, the semi-flow will be specified by differential equations, and the reaction coefficients will turn up as coefficients in the equations.

In the stochastic model, a reaction coefficient will turn up as a weight affecting the probability that the reaction fires within a small interval of time.

`Where do the reaction coefficients show up in these models? In the deterministic case, the semi-flow will be specified by differential equations, and the reaction coefficients will turn up as coefficients in the equations. In the stochastic model, a reaction coefficient will turn up as a weight affecting the probability that the reaction fires within a small interval of time.`

Either way, we have a second choice on the menu: an analytic approach to the dynamics, or an experimental approach via simulation.

`Either way, we have a second choice on the menu: an analytic approach to the dynamics, or an experimental approach via simulation.`

Let's start with the deterministic case. This applies well to chemical reactions involving large numbers of molecules, where the different species involved in the reaction are "well-mixed." (E.g. no clumps of dimers surrounded by zones of monomers.)

`Let's start with the deterministic case. This applies well to chemical reactions involving large numbers of molecules, where the different species involved in the reaction are "well-mixed." (E.g. no clumps of dimers surrounded by zones of monomers.)`