Exact distributions of threshold crossing times of proteins under post-transcriptional regulation by small RNAs

The timings of several cellular events like cell lysis, cell division, or pore formation in endosomes are regulated by the time taken for the relevant proteins to cross a threshold in number or concentration. Since protein synthesis is stochastic, the threshold crossing time is a first passage problem. The exact distributions of these first passage processes have been obtained recently for unregulated and autoregulated genes. Many proteins are however regulated by post-transcriptional regulation, controlled by small noncoding RNAs (sRNAs). Certain mathematical models of gene expression with post-transcriptional sRNA regulation have been recently exactly mapped to models without sRNA regulation. Utilizing this mapping and the exact distributions, we calculate exact results on fluctuations (full distribution, all cumulants, and characteristic times) of protein threshold crossing times in the presence of sRNA regulation. We derive two interesting predictions from these exact results. We show that the size of the fluctuation of the threshold crossing times have a nonmonotonic U-shaped behavior as a function of the rates of binding and unbinding of the sRNA-mRNA complex. Thus there are optimal parameters that minimize noise. Furthermore, the fluctuations in models with sRNA regulation may be higher or lower compared to the model without regulation, depending on the mean protein burst size.