# Optimal first–passage time in gene regulatory networks

## Abstract

The inherent probabilistic nature of the biochemical reactions, and low copy number of species can lead to stochasticity in gene expression across identical cells. As a result, after induction of gene expression, the time at which a specific protein count is reached is stochastic as well. Therefore events taking place at a critical protein level will see stochasticity in their timing. First–passage time (FPT), the time at which a stochastic process hits a critical threshold, provides a framework to model such events. Here, we investigate stochasticity in FPT. Particularly, we consider events for which controlling stochasticity is advantageous. As a possible regulatory mechanism, we also investigate effect of auto–regulation, where the transcription rate of gene depends on protein count, on stochasticity of FPT. Specifically, we investigate for an optimal auto-regulation which minimizes stochasticity in FPT, given fixed mean FPT and threshold.

For this purpose, we model the gene expression at a single cell level. We find analytic formulas for statistical moments of the FPT in terms of model parameters. Moreover, we examine the gene expression model with auto–regulation. Interestingly, our results show that the stochasticity in FPT, for a fixed mean, is minimized when the transcription rate is independent of protein count. Further, we discuss the results in context of lysis time of an E. coli cell infected by a phage virus. An optimal lysis time provides evolutionary advantage to the phage, suggesting a possible regulation to minimize its stochasticity. Our results indicate that there is no auto–regulation of the protein responsible for lysis. Moreover, congruent to experimental evidences, our analysis predicts that the expression of the lysis protein should have a small burst size.

## I Introduction

Gene expression is the process of transcription of genetic information to mRNAs, and translation of each mRNA to proteins. As the copy number of species involved in the process is small, the probabilistic nature of biochemical reactions reflects as stochastcity in gene expression [1, 2, 3, 4, 5, 6].

Stochasticity in gene expression has an important role in several cellular functions. For example, it can lead genetically identical cells to different cell–fates [7, 8, 9, 10, 11, 12]. This helps the cells in responding to the ever–changing environment [13, 14, 15, 16]. On the other hand, stochasticity in expression of housekeeping genes can lead to diseased states [17, 18, 19], and needs to be minimized [20, 21]. Accordingly, different regulatory mechanisms are employed to control stochastic fluctuations [22, 23, 24, 25, 26, 27, 28, 29]. Auto–regulation wherein transcription rate is a function of protein count is an example of one such mechanism. Its effect on stochascticity in gene expression has been a subject of several studies [27, 28, 29].

After onset of gene expression, its stochasticity consequently manifests into stochasticity in the time at which a certain protein level is reached. This implies that the timing of a cellular event which triggers at a critical protein level is stochastic in nature [30, 31]. For instance, lysis time for an E. coli cell infected by a phage virus is stochastic. Lysis of the cell takes place when holin, the protein responsible for lysis, reaches a critical threshold [32, 33, 34].

Further, it has been suggested that optimality in lysis time provides evolutionary advantage to phage virus [35, 36, 37, 38, 39]. This indicates that there could be some regulation of gene expression to ensure lysis at the optimal time, with minimum stochastic fluctuations. In this work, we study stochasticity in first–passage time (FPT), the time it takes for the protein count to reach a fixed threshold for the first time [40], at a single–cell level. We investigate the effect of auto-regulation of transcription on stochasticity of FPT. In particular we seek answer to the question: given the mean FPT (corresponding to optimal lysis time, for example), what auto-regulatory feedback will lead to minimum stochasticity in the FPT?

We first formulate an unregulated gene expression model assuming transcription, translation, and mRNA degradation while considering proteins to be stable. Along the lines of [34], we find expressions for statistical moments of FPT for this model, and discuss their implications with respect to minimizing variance in FPT for given mean FPT. Next, we introduce auto-regulation in the above model and derive the moments for FPT. Then, we deduce the expression for optimal feedback function that minimizes the variance in FPT for a given mean. We show that a negative or positive feedback always results into higher variance in first passage time for a given mean than the case when there is no feedback. The results are validated by carrying out simulations. Also, various notations used in this work are tabulated in Table I.

Transcription rate for unregulated gene expression model. | |

Translation rate for both unregulated, and regulated gene expression. models | |

mRNA degradation rate for both unregulated, and regulated gene. expression models | |

Burst size after transcriptional event. | |

Parameter of geometric distribution corresponding to. protein bursts | |

Mean of protein burst size. | |

Protein count at time . | |

Protein count after burst. | |

Transcription rate for auto-regulated gene expression model after transcription event. | |

Threshold for protein count. | |

Minimum number of transcription events for protein count to reach the threshold . | |

Waiting time for transcription event. | |

is an Exponential random variable with parameter . The probability density function of is given by . | |

Probability mass function for minimum number of transcription events to reach the threshold . | |

Probability mass function for protein count after transcription events. | |

Expectation operator. | |

Var | Variance. |

Maximum possible transcription rate in model with feedback implemented using Hill function . | |

Fraction of transcription rate that corresponds to minimum transcription rate in model with feedback implemented using Hill function. | |

Hill coefficient. | |

Coefficient proportional to binding efficiency; decides when half rate concentration is reached. |

## Ii First–Passage Time For Gene Expression Model Without Regulation

In this section, we formulate a stochastic gene expression model (as shown in Fig. 1). Then, we define the FPT for this model and derive expressions for its statistical moments. We also discuss the implications of these expressions in context of minimizing variance of FPT, for fixed mean and threshold.

### Ii-a Model Formulation

In the model under consideration transcription of mRNAs from the gene occurs at a rate , translation of proteins from each mRNA occurs at a rate , and each mRNA degrades at a rate . The time interval between two transcription events is exponentially distributed. We assume proteins to be stable as the lysis protein in phage, i.e. holin, is stable [41]. To further simplify the model, we assume each mRNA molecule degrades instantaneously after producing a burst of random number of protein molecules [42, 43, 44, 45]. Consistent with experimental, and theoretical evidences; we assume that protein burst follows a geometric distribution, and the mean burst size is given by [46, 47]. Thus, the simplified model considers gene expression wherein each burst event (equivalent to transcription event) occurs at an exponentially distributed time with parameter , and size of burst follows a geometric distribution with mean .

Let us denote the size of burst by random variable and the parameter of its distribution by . The probability mass function, therefore, can be written as [48]:

(1) |

The mean burst size, , can be expressed as [48]:

(2) |

Further, let protein count after transcription events be denoted as . It can be expressed as a sum of random variables :

(3) |

Being sum of independent and identically distributed geometric random variables, has a negative binomial distribution with parameters and [49]. The probability mass function of , denoted as , can be expressed as [49]:

(4) |

Also, the cumulative distribution function is given by [50]:

(5) |

where is regularized incomplete beta function:

(6) |

and satisfies the following property:

(7) |

We have determined the distribution for protein population. Next, we defined the first–passage time (FPT) for the protein count to reach a certain threshold.

### Ii-B Expression for First Passage Time

For a random process corresponding to protein count, , with , the first passage time (FPT), for a threshold is defined as:

(8) |

Because in our model, the protein count changes only when a burst occurs (or equivalently, a transcription event occurs); we can calculate the minimum number of transcription events, , it takes for the protein count to reach the threshold and define the FPT as sum of inter–burst arrival times. This has been depicted in Fig. 2.

Let the time between and bursts be denoted by random variable , then:

(9) |

where is given by the following equation:

(10) |

Note that in Eq. (9), are independent, and identically distributed exponential random variables with parameter . We denote this by . Also, each of is independent of .

Using standard results from probability theory, one may write [51]:

(11a) | ||||

(11b) |

It can be noted that to determine statistical moments of FPT in Eq. (11a)–(11b), we need to derive expressions for first two moments of , and .

#### First Two Moments of

The cumulative distribution function for defined in Eq. (10) can be written as:

(12a) | ||||

(12b) |

Since is a negative binomial distribution, we have:

(13a) | ||||

(13b) |

Using the property of incomplete beta function mentioned in Eq. (7), we get:

(14) |

#### First Two Moments of

Since , its statistical moments are given by:

(17a) | ||||

(17b) |

We now have expressions for first two moments of , and . The expressions for first two moments of FPT in terms of model parameters can, therefore, be written as:

(18) | ||||

(19) |

where the approximations are valid when . It can be observed a smaller mean burst size would result in smaller variance of FPT. The mean FPT can be kept fixed by a commensurate change in the transcription rate, . Therefore, the variance can independently be reduced by a lower mean burst size . This means adopting a high transcription rate , and a low translation rate (and/or having a higher degradation rate for the mRNAs) results in a lower variance in FPT without affecting its mean.

Further, we note that by using from Eq. (17b), we can deduce the following relationship between and from Eq. (11a) and Eq. (11b):

(20) |

We shall use above relationship in the later part of the paper while deriving expression of the auto-regulation function that minimizes variance in FPT, for given mean FPT.

Next, we introduce auto-regulation of transcription rate by the protein count to investigate how the expressions for statistical moments of FPT change.

## Iii Introducing Auto-regulation in Gene Expression Model

To investigate the effect of auto-regulation on statistical moments of FPT, we assume that transcription rate is a function of protein count, i.e., it changes after each transcription event. We denote the transcription rate after arrival of burst as . Similar to previous section, we need to derive expression for moments of inter–burst arrival times , and minimum number of transcription events in order to derive the expression for FPT moments defined in Eq. (9).

We note that the translation burst size is independent of the transcription rate. Therefore, distribution of to reach a certain threshold is same as gene expression model without any regulation discussed in previous section. However, distribution of each is different and depends upon corresponding rate of transcription.

We derive expressions for first two moments of each to find analytical forms of first two moments of FPT.

### Iii-a Inter–burst arrival time for auto-regulatory gene expression model

It may be noted that if protein count after any burst event is known, arrival time for the next burst will be exponentially distributed. Therefore, the distribution of each can be modelled as a conditional exponential distribution. More specifically, we can write:

(21) |

where , and respectively denote the arrival time for burst and protein count after the burst.

The expressions for mean and variance of can be calculated as follows.

#### Mean

Before arrival of the first burst, there are no protein molecules, i.e., for . Therefore, we can write the mean for arrival time for the first burst as:

(22) |

For , the corresponding arrival times would be conditionally exponential, implying:

(23a) | ||||

(23b) | ||||

(23c) |

#### Second Order Moments

Adopting similar approach as above, we derive the expressions for second order moments of . For , we have:

(24a) | ||||

For : | ||||

(24b) | ||||

(24c) | ||||

(24d) |

Therefore the expression for variance of :

(25) |

For , the expression for will be

(26) |

Moreover, we have following relationship first two moments of the random variable :

(27) |

which alongwith Eq. (26), and (25) yields:

(28) |

We note that the equality above holds for . We will use it in later part of the paper while deducing the expression for optimal auto-regulation that leads to minimum variance in the FPT for fixed mean.

Having derived the expressions for moments of inter–bursts arrival times, we see how the introduction of auto-regulation influences the expressions for FPT moments.

### Iii-B FPT for auto-regulatory gene expression model

We present the expressions for statistical moments of FPT in theorem–proof format. In developing the proofs, we make use of the fact that each will be independent of . Also, are independent of each other. However, they are not identically distributed like the unregulated gene expression case discussed in previous section.

###### Theorem 1 (Mean of First Passage Time)

To prove the result, we first find conditional expectation given then we have:

(30a) | ||||

(30b) | ||||

Unconditioning above expression with respect to : | ||||

(30c) | ||||

(30d) |

This completes the proof.

###### Theorem 2 (Variance of First Passage Time)

Since expression for is known and given by Eq. (29), we need to find expression for , in order to find expression for variance of FPT.

Using the definition of first passage time in Eq. (9), we have:

(32a) | ||||

(32b) | ||||

Since are independent of each other, and are independent of for each ; we can write: | ||||

(32c) | ||||

(32d) | ||||

Using , we have: | ||||

(32e) |

Unconditioning with respect to , expression for becomes:

(33) |

Therefore, using Eq. (29), and Eq. (33); expression for becomes:

(34) |

This completes the proof. So far we have developed analytical expressions for mean and variance of FPT when there is an auto-regulatory feedback to transcription rate from protein count. In the next section, we make use of these expressions to deduce the optimal auto-regulation function to minimize the variance of FPT assuming fixed mean FPT.

## Iv Minimizing Variance in First Passage Time for Given Mean

Parameter | Unit | Positive feedback | Negative feedback | No feedback |
---|---|---|---|---|

mRNA produced per minute | 19.35 | 84 | 10 | |

protein produced per mRNA per minute | 2.65 | 2.65 | 2.65 | |

per minute | 0.3 | 0.3 | 0.3 | |

molecules | 5000 | 5000 | 5000 | |

- | 0.05 | 0.05 | - | |

per molecule | 0.002 | 0.002 | - | |

- | 2 | 2 | - |

In this section, we find expression for the auto-regulatory feedback function, that gives minimum variance in FPT, given the mean FPT and event threshold are fixed. The result is presented in form of a theorem.

###### Theorem 3 (Optimal feedback for minimum variance)

Let the first passage time be defined as Eq. (9), and its mean and variance, respectively, given by Eq. (29) and Eq. (31). Then, the optimal function to minimize the variance of FPT for a given mean of FPT will be constant, given by following expression:

(35) |

where denotes the minimum number of transcription events required to reach the FPT threshold, and is given by Eq. (16a).

We assume that each burst event adds a perturbation to transcription rate, i.e., can be written as:

(36) |

where is perturbation corresponding to transcription rate after burst. To prove the result, we shall prove that the variance of FPT for given mean will minimize when .

Recalling the expression for from Eq. (29):

(37) |

Using expressions in Eqs. (22), (23c), we can deduce the expressions for as:

(38) |

where is related with by following expression:

(39) |

Substituting expression for from Eq. (38), we have:

(40a) | ||||

(40b) | ||||

(40c) |

Since , we have:

(41) |

Note that for a fixed mean FPT, minimizing the variance of FPT and minimizing the second order moment are equivalent.

Now, we consider the expression for , and use expression in Eq. (41) to deduce the desired optimal function. From Eq. (31), we have:

(42) |

Substituting value of from Eq. (36), we get following expression for :

(43) |

Further simplifying and using relation obtained in Eq. (41) yields:

(44) |

(45a) | ||||

(45b) |

Further, we note that in above expression if (or equivalently ), the expression minimizes and reduces to:

(46) |

Recalling Eq. (20), we observe that equality in above expression holds for unregulated gene expression case, which essentially means . This proves the desired result.

In this section, we proved that having no auto-regulation of transcription rate provides minimum stochasticity in the FPT, if mean FPT and event threshold are kept fixed. However, since our analysis simplified the gene expression model to burst–limit, we are interested in validating whether it is true if we don’t make an approximation. In the next section, we discuss the computer simulations we carried out for this purpose.

## V Simulation Results

In order to verify the result deduced in previous section, we carried out Monte Carlo simulations using Gillespie’s algorithm [52]. We did not specifically assume that production of protein is in geometric bursts with parameter . Instead, we assumed a non–zero half–life for mRNA thereby relaxing the burst approximation.

To simulate, we considered three separate cases: no feedback, negative feedback and positive feedback. The positive feedback is implemented using Hill function as follows:

(47) |

where is maximum transcription rate, represents minimum transcription rate as the fraction of , denotes the Hill coefficient while is coefficient proportional to the binding affinity (when ).

Similarly, the negative feedback is implemented using following function:

(48) |

We carried out the simulations for several sets of parameters assuming a fixed event threshold. Rest of the model parameters were chosen to keep the mean FPT approximately equal. In all of them, we found that no–feedback case has minimum variance in FPT.

In Table II, we present one set of such parameters. We assumed the event threshold . Other parameters are chosen in a way that the mean FPT minutes.

Simulation results for 10000 realizations are shown in Fig. 3. We note that the variance is minimum in no–feedback case, validating our theoretical claims for this set of parameter values.

## Vi Discussion

In this work, we studied stochasticity in event timing at a single cell level. We considered a standard gene expression model without protein degradation. Next, we formulated the FPT problem for this model and derived the formulas for statistical moments of FPT. Further, we introduced auto-regulation in the gene expression wherein the transcription rate is a function of protein count. We derived the formulas for moments of FPT in this case as well, and demonstrated that for a given mean of FPT, the variance in FPT is minimized when there is no auto-regulation of gene expression. The result was verified with simulations as well.

The result can be connected to the phage lysis time. Due to existence of optimal lysis time [35, 36], the phage would possibly like to kill the cell at that time with as much precision as possible. Thus, it should resort to a strategy that would minimize the lysis time variance and hence have no protein–dependent feedback regulation of transcription rate in the expression of holin. In expression from late promoter in phage, which produces holin, has no evidence of a regulation [53, 54].

Recalling that in no auto–regulation case too, the variance of FPT can be independently decreased by lowering the mean burst size . Other studies also reveal that in case of phage, the burst size is indeed small [33, 35]. Also, antiholin, another protein expressed from the same promoter that expresses holin, binds to holin to decrease the effective burst size [55, 34].

In this paper, there is an underlying assumption of protein being stable. In future work, we plan to use a gene expression model with protein degradation, and carry out a similar analysis. This can be further extended to more generalized gene expression models wherein the promoter can also switch between on and off states [43, 12].

## Acknowledgment

AS is supported by the National Science Foundation Grant DMS-1312926, University of Delaware Research Foundation (UDRF) and Oak Ridge Associated Universities (ORAU).

### References

- W. J. Blake, M. Kaern, C. R. Cantor, and J. J. Collins, “Noise in eukaryotic gene expression,” Nature, vol. 422, pp. 633–637, 2003.
- J. M. Raser and E. K. O’Shea, “Noise in gene expression: origins, consequences, and control,” Science, vol. 309, pp. 2010–2013, 2005.
- A. Raj and A. van Oudenaarden, “Nature, nurture, or chance: Stochastic gene expression and its consequences,” Cell, vol. 135, pp. 216 – 226, 2008.
- B. Munsky, G. Neuert, and A. van Oudenaarden, “Using gene expression noise to understand gene regulation,” Science, vol. 336, pp. 183–187, 2012.
- M. Kaern, T. C. Elston, W. J. Blake, and J. J. Collins, “Stochasticity in gene expression: from theories to phenotypes,” Nature Review Genetics, vol. 6, pp. 451–64, 2005.
- A. Singh and M. Soltani, “Quantifying intrinsic and extrinsic variability in stochastic gene expression models,” PLoS ONE, vol. 8, p. e84301, 12 2013.
- R. Losick and C. Desplan, “Stochasticity and cell fate,” Science, vol. 320, pp. 65–68, 2008.
- A. Arkin, J. Ross, and H. McAdams, “Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda–infected escherichia coli cells,” Genetics, vol. 149, pp. 1633–1648, 1998.
- L. S. Weinberger, J. C. Burnett, J. E. Toettcher, A. P. Arkin, and D. V. Schaffer, “Stochastic gene expression in a lentiviral positive-feedback loop: Hiv-1 tat fluctuations drive phenotypic diversity,” Cell, vol. 122, pp. 169–182, 2005.
- J.-W. Veening, W. K. Smits, and O. P. Kuipers, “Bistability, epigenetics, and bet-hedging in bacteria,” Annu. Rev. Microbiol., vol. 62, pp. 193–210, 2008.
- J. Hasty, J. Pradines, M. Dolnik, and J. J. Collins, “Noise-based switches and amplifiers for gene expression,” Proceedings of the National Academy of Sciences, vol. 97, pp. 2075–2080, 2000.
- A. Singh, B. Razooky, C. D. Cox, M. L. Simpson, and L. S. Weinberger, “Transcriptional bursting from the hiv-1 promoter is a significant source of stochastic noise in hiv-1 gene expression,” Biophysical Journal, vol. 98, pp. L32–L34, 2010.
- A. Eldar and M. B. Elowitz, “Functional roles for noise in genetic circuits,” Nature, vol. 467, pp. 167–173, Sept. 2010.
- E. Kussell and S. Leibler, “Phenotypic diversity, population growth, and information in fluctuating environments,” Science, vol. 309, pp. 2075–2078, 2005.
- N. Balaban, J. Merrin, R. Chait, L. Kowalik, and S. Leibler, “Bacterial persistence as a phenotypic switch,” Science, vol. 305, pp. 1622–1625, 2004.
- M. Acar, J. T. Mettetal, and A. van Oudenaarden, “Stochastic switching as a survival strategy in fluctuating environments,” Nature Genetics, vol. 40, pp. 471–475, 2008.
- R. Kemkemer, S. Schrank, W. Vogel, H. Gruler, and D. Kaufmann, “Increased noise as an effect of haploinsufficiency of the tumor-suppressor gene neurofibromatosis type 1 in vitro,” Proceedings of the National Academy of Sciences, vol. 99, pp. 13 783–13 788, 2002.
- D. L. Cook, A. N. Gerber, and S. J. Tapscott, “Modeling stochastic gene expression: implications for haploinsufficiency,” Proceedings of the National Academy of Sciences, vol. 95, pp. 15 641–15 646, 1998.
- R. Bahar, C. H. Hartmann, K. A. Rodriguez, A. D. Denny, R. A. Busuttil, M. E. Dollé, R. B. Calder, G. B. Chisholm, B. H. Pollock, C. A. Klein, et al., “Increased cell-to-cell variation in gene expression in ageing mouse heart,” Nature, vol. 441, pp. 1011–1014, 2006.
- B. Lehner, “Selection to minimise noise in living systems and its implications for the evolution of gene expression,” Molecular systems biology, vol. 4, 2008.
- H. B. Fraser, A. E. Hirsh, G. Giaever, J. Kumm, and M. B. Eisen, “Noise minimization in eukaryotic gene expression,” PLoS biology, vol. 2, p. e137, 2004.
- U. Alon, “Network motifs: theory and experimental approaches,” Nature Reviews Genetics, vol. 8, pp. 450–461, 2007.
- A. Becskei and L. Serrano, “Engineering stability in gene networks by autoregulation,” Nature, vol. 405, pp. 590–593, 2000.
- H. El-Samad and M. Khammash, “Regulated degradation is a mechanism for suppressing stochastic fluctuations in gene regulatory networks,” Biophysical journal, vol. 90, pp. 3749–3761, 2006.
- P. S. Swain, “Efficient attenuation of stochasticity in gene expression through post-transcriptional control,” Journal of Molecular Biology, vol. 344, pp. 965 – 976, 2004.
- D. Orrell and H. Bolouri, “Control of internal and external noise in genetic regulatory networks,” Journal of theoretical biology, vol. 230, pp. 301–312, 2004.
- A. Singh and J. P. Hespanha, “Optimal feedback strength for noise suppression in autoregulatory gene networks,” Biophysical journal, vol. 96, pp. 4013 – 4023, 2009.
- Y. Tao, X. Zheng, and Y. Sun, “Effect of feedback regulation on stochastic gene expression,” Journal of Theoretical Biology, vol. 247, pp. 827 – 836, 2007.
- A. Singh, “Negative feedback through mrna provides the best control of gene-expression noise,” NanoBioscience, IEEE Transactions on, vol. 10, pp. 194–200, 2011.
- A. Amir, O. Kobiler, A. Rokney, A. B. Oppenheim, and J. Stavans, “Noise in timing and precision of gene activities in a genetic cascade,” Molecular Systems Biology, vol. 3, 2007.
- R. Murugan and G. Kreiman, “On the minimization of fluctuations in the response times of autoregulatory gene networks,” Biophysical Journal, vol. 101, pp. 1297–1306, 2011.
- R. White, S. Chiba, T. Pang, J. S. Dewey, C. G. Savva, A. Holzenburg, K. Pogliano, and R. Young, “Holin triggering in real time,” Proceedings of the National Academy of Sciences, vol. 108, pp. 798–803, 2011.
- J. Dennehy and I.-N. Wang, “Factors influencing lysis time stochasticity in bacteriophage lambda,” BMC Microbiology, vol. 11, no. 1, p. 174, 2011.
- A. Singh and J. Dennehy, “Stochastic holin expression can account for lysis time variation in the bacteriophage ,” Journal of the Royal Society Interface (to appear), 2014.
- I.-N. Wang, “Lysis timing and bacteriophage fitness,” BMC Microbiology, vol. 172, pp. 17–26, January 2006.
- I.-N. Wang, D. E. Dykhuizen, and L. B. Slobodkin, “The evolution of phage lysis timing,” Evolutionary Ecology, vol. 10, pp. 545–558, 1996.
- R. Heineman and J. Bull, “Testing optimality with experimental evolution: lysis time in a bacteriophage,” Evolution, vol. 61, pp. 169;5–1709, 2007.
- Y. Shao and I.-N. Wang, “Bacteriophage adsorption rate and optimal lysis time,” Genetics, vol. 180, pp. 471–482, 2008.
- J. A. Bonachela and S. A. Levin, “Evolutionary comparison between viral lysis rate and latent period,” Journal of Theoretical Biology, vol. 345, pp. 32 – 42, 2014.
- S. Redner, A guide to first-passage processes. Cambridge University Press, 2001.
- Y. Shao and N. Wang, “Effect of late promoter activity on bacteriophage fitness,” Genetics, vol. 181, pp. 1467–1475, 2009.
- N. Friedman, L. Cai, and X. S. Xie, “Linking stochastic dynamics to population distribution: An analytical framework of gene expression,” Phys. Rev. Lett., vol. 97, p. 168302, 2006.
- V. Shahrezaei and P. S. Swain, “Analytical distributions for stochastic gene expression,” Proceedings of the National Academy of Sciences, vol. 105, pp. 17 256–17 261, 2008.
- J. Paulsson, “Models of stochastic gene expression,” Physics of Life Reviews, vol. 2, pp. 157 – 175, 2005.
- O. G. Berg, “A model for the statistical fluctuations of protein numbers in a microbial population,” Journal of Theoretical Biology, vol. 71, pp. 587 – 603, 1978.
- J. Yu, J. Xiao, X. Ren, K. Lao, and X. S. Xie, “Probing gene expression in live cells, one protein molecule at a time,” Science, vol. 311, pp. 1600–1603, 2006.
- V. Elgart, T. Jia, A. T. Fenley, and R. Kulkarni, “Connecting protein and mrna burst distributions for stochastic models of gene expression,” Physical biology, vol. 8, p. 046001, 2011.
- M. H. DeGroot and M. J. Schervish, Probability and Statistics, 4th ed. Pearson, 2012.
- A. Papoulis and S. U. Pillai, Probability, Random Variables and Stochastic Processes, 4th ed. McGraw Hill, 2002.
- M. R. Spiegel, Theory and Problems of Probability and Statistics. McGraw-Hill, 1992.
- S. K. Ross, Introduction to Probability Models, 10th ed. Academic Press, 2010.
- D. Gillespie, “Exact stochastic simulation of coupled chemical reactions,” J Phys Chem, vol. 81, pp. 2340–2361, 1977.
- M. Ptashne, A Genetic Switch – Phage and Higher Organisms, 2nd ed. Cell Press & Blackwell Scientific Publications, 1991.
- A. B. Oppenheim, O. Kobiler, J. Stavans, D. L. Court, and S. Adhya, “Switches in bacteriophage lambda development,” Annual Review of Genetics, vol. 39, pp. 409–429, 2005.
- D. L. Smith, U. Blasi, and R. Young, “Dimerization between the holin and holin inhibitor of phage ,” Journal of Bacteriology, vol. 182, pp. 6075–6081, 2000.