Tuesday, 8 April 2008

biochemistry - What does it mean for a distribution to be "consistent with a two rate-limiting stochastic steps"?

WYSIWYG is almost there, but you need one more piece of information to make this explicit.



The distribution cited in the paper is h(t)tet/τ (we're going to ignore normalizing constants today). We can recognize this as a particular case of the Gamma distribution, with PDF:



fGamma(t|k,θ)tk1et/θ.



In particular, h(t) looks like the PDF of a Gamma(2,τ) random variable,
h(t)fGamma(t|2,τ)tet/τ.



So the authors of the paper are claiming that the burst duration τburst is a random variable with a Gamma(2,70sec) distribution.



How do we get from that to "This distribution was consistent with two rate-limiting stochastic steps?"
If we recall from our stats class (or look up properties of the Gamma distribution on Wikipedia), a Gamma distribution with integer values of k is the distribution of the waiting time for k events to occur in a Poisson process. The setup of a Poisson process is that you are tracking the occurrence of unlikely/infrequent events in time, so the fact that the distribution h(t) of burst times looks like a Gamma(2,τ) suggests that the burst duration τburst is determined by the occurrence of 2 stochastic events (with the same "frequency" 1/τ). In other words, this distribution is consistent with the following scenario: a nuclear localization occurs, and then will stop after two particular events, where the timing of the events is governed by a Poisson process (which, as WYSIWYG pointed out, is a reasonable model for the occurrence of chemical reactions with reasonably slow kinetics, i.e. you can count individual reactions in time). If this is true, the distribution of τburst will look like tet/τ.



The authors generalize that statement somewhat by saying that technically there could be more events/reactions required to terminate a localization burst, but all but 2 of those reactions happen extremely rapidly, i.e. there are 2 rate-limiting steps in the decision process (both of which happen to have the same value for τ).




Edit: I realized that one more mathematical point might make this more clear: The statement about the Gamma distribution and Poisson processes above is equivalent to the statement that the sum of k iid random variables with distribution Exponential(τ) is a random variable with a Gamma(k,τ) distribution. Thus h(t) is literally (up to normalization) the distribution of a sum of 2 independent exponential random variables. If you read up on the connection between the exponential distribution and waiting times that might make the claim in the paper seem more transparent.



Source: https://en.wikipedia.org/wiki/Gamma_distribution#Special_cases



http://en.wikipedia.org/wiki/Exponential_distribution

No comments:

Post a Comment