In the previous lecture we discussed a simple model of protein production with an activator.
Now lets consider a more involved model.
X has its own dynamics
We can write down a rate equation describing the dynamics.
In steady state, the concentration of X is constant.
Transcription factors (in this case it is X) often oligomerize before binding to DNA. Therefore the dynamics must be modified in order to take this into account.
So the rate equation needs to take into account the many X oligomerize and then bind to DNA.
Oligomerization is faster than transcription/translation. So we can safely assume
The above function is known as a Hill function. The function sharpens as increases. Note the following limiting case:
- Oligomerization (also referred to as cooperativity) acts to sharpen gene regulation functions.
- However, when comparing with experiments our model does not agree exactly with experimental results.
- Experiments show that biological systems are very noisy.
- Theoretically there are two possible explanations in regards to gene expression to account for the noise seen in biological systems.
- X is noisy only, and that noise propagates to Y.
- X is noisy and Y is also noisy on its own.
- Turns out that statement 2 is correct.
We will describe the noise in Y as having two contributions. One due to the intrinsic noisiness of Y and the other due to the propagation of noise from X (extrinsic noise).
The noise from X is denoted with subscript “ext” for extrinsic noise. The noise due to Y itself is labelled “int” for intrinsic noise.
Lets first consider only intrinsic noise without the presence of X.
We already wrote down the simple rate equation for this model in the previous lecture. Instead of thinking in terms of rate equations lets now look at this in terms of probabilities. Let denote the number of Y proteins in the cell. Lets formulate the probabilities of going to state with one more or one less Y protein.
Transition Probabilities and Master Equation
Then the master equation is the following:
In stationary state all time derivatives are equal to zero. For n = 0 :
For n = 1 :
A pattern emerges (take my word for it)…
This is a Poisson distribution! We can normalize our distribution in order to solve for which yields:
Like any Poisson distribution the mean and variance are equal.
We now have an expression for the intrinsic noise!
Relative noise is the ratio of the noise to the mean number squared.
So the relative intrinsic noise decreases as the average number of proteins increases.
Written with StackEdit.