Derivation of time probability distributions for single molecule kinetics

If we have a protein $$\ce{A}$$ that can bind to protein $$\ce{B}$$:

$$\ce{A + B <=>[k_\mathrm{fwd,1}][k_\mathrm{rev,1}] AB}$$

and the bound complex $$\ce{AB}$$ can be converted into $$\ce{C}$$:

$$\ce{AB <=>[k_\mathrm{fwd,2}][k_\mathrm{rev,2}] C}$$

Assume we have a volume of 1 and some initial amount of $$\ce{A}$$ and $$\ce{B}$$ (which are $$A_0,$$ $$B_0$$) and that $$C_0 = 0.$$ We then measure the amount of time elapsed before the first molecule of $$C$$ is formed. The experiment is repeated many times, and the frequency of occurence of first $$C$$ during each time interval 0 to $$t_1$$, $$t_1$$ to $$t_2$$, $$t_2$$ to $$t_3$$, etc. is plotted, so that the result is a probability distribution of the time elapsed before the first molecule of $$C$$ is formed.

I would like to know what is the mathematical function that describes this distribution and how is it derived?

UPDATE: this is not a homework question and the point is not a specific answer. I would like references to literature/textbooks that derives these distributions analytically for simple examples like the above would be helpful. Have not seen a setup where distribution of duration is derived. Commenters say it's standard. Please give references.

• This is a very standard set-up for a kinetics problem. The typical way involves making the steady-state approximation which says that the change in concentration of AB with time is zero (essentially because it is eliminated immediately). If one can't make an assumption like this, there typically won't be a closed-form solution. One can always solve these types of differential equations numerically though. Basically I'm saying you probably need to give more info for this to be answerable. – jheindel Jun 3 at 3:44
• It would help if you could at least give the equilibrium constants for the two steps to define the distribution of species at equilibrium. "The first C" is a strange thing to ask about - the standard kinetics are about bulk properties. – Karsten Theis Jun 3 at 11:24
• @jheindel happy to assume AB doesn't change with time. do you have a link to an example where durations are computed for a problem of this sort? – user9277 Jun 3 at 17:43
• not a homework question. happy with references instead of answer. closing it as off topic is silly – user9277 Jun 4 at 1:14
• Here's a more basic write-up of the derivation for a unimolecular first-order reaction. Extending it to the more complex system is a matter of a bit more math using the same concepts. makarov.cm.utexas.edu/resources/Lecture-notes,-tutorials-etc./… – Andrew Jun 5 at 16:30

To see how this calculation is done, let us consider the simpler case of just the second step shown above, ie $$\ce{AB ->[k_{f2}] C}$$. Since we are interested in the situation where no $$C$$ has yet formed, the reverse reaction can be ignored. I will assume that this reaction is first order in $$AB$$, so $$k_{f2}$$ has units of $$s^{-1}$$.

In this case, the rate constant $$k_{f2}$$ represents the average frequency of the reaction.

Now let's consider that we establish an infinite number of evenly spaced time intervals of length $$\Delta t$$. When we consider the probability distribution, we need to consider two factors:

1) The probability of reaction occurring in an interval $$t_n + \Delta t$$ and

2) The probability that the reaction has not yet occurred at $$t_n$$

The first value is simply the product of the rate constant $$k_{f2}$$ and the length of the time interval: $$k_{f2}\Delta t$$.

The second probability is slightly more complicated, since it is the chance of the reaction $$not$$ occurring in each of the $$n$$ intervals preceding $$t_n$$. That probability is given by $$1-(k_{f2}\Delta t)^n$$. Conveniently, as our time interval $$\Delta t$$ goes to $$0$$, this probability becomes a continuous function $$e^{-k_{f2}t}$$.

To get our final probability, we simply multiply these two probabilities:

$$[\text{total probability of reaction occurring in interval }t+dt] = [\text{probability of occurring in interval }t+dt] \times [\text{probability that reaction has not occurred by time }t]=k_{f2}dt\times e^{-k_{f2}t}$$

The probability as a continuous function of $$t$$ therefore can be written as

$$p_{\ce{AB->C}}(t)=k_{f2}e^{-k_{f2}t}$$.

As expected, integration from $$t=0$$ to $$t=\infty$$ gives an overall probability of $$1$$.

For the two-step process, if $$\ce{A+B <=> AB}$$ establishes equilibrium very fast relative to the $$\ce{AB->C}$$ reaction, then the above is a reasonable approximation (after adjustment of $$k_{f2}$$ to reflect the equilibrium concentration of $$AB$$). If the equilibrium is not established quickly, then you would use the same reasoning to determine the probability first of $$AB$$ having formed by time $$t$$ and then the probability of it reacting based on how long it has been in existence.
• Yes, instead of a y-intercept of $k_{f2}$, the probability would start at 0 and increase with the concentration of AB. – Andrew Jun 12 at 14:41
• @Andrew: quick clarification: shouldn't you take into account not only the rate $k_{f2}$ but also the amount of $AB$ you have formed thus far in time? the formation of $C$ should depend on the rate and on that abundance – user9277 Jun 13 at 16:30