Transitions in a genetic transcriptional regulatory system under Lévy motion

Based on a stochastic differential equation model for a single genetic regulatory system, we examine the dynamical effects of noisy fluctuations, arising in the synthesis reaction, on the evolution of the transcription factor activator in terms of its concentration. The fluctuations are modeled by Brownian motion and α-stable Lévy motion. Two deterministic quantities, the mean first exit time (MFET) and the first escape probability (FEP), are used to analyse the transitions from the low to high concentration states. A shorter MFET or higher FEP in the low concentration region facilitates such a transition. We have observed that higher noise intensities and larger jumps of the Lévy motion shortens the MFET and thus benefits transitions. The Lévy motion activates a transition from the low concentration region to the non-adjacent high concentration region, while Brownian motion can not induce this phenomenon. There are optimal proportions of Gaussian and non-Gaussian noises, which maximise the quantities MFET and FEP for each concentration, when the total sum of noise intensities are kept constant. Because a weaker stability indicates a higher transition probability, a new geometric concept is introduced to quantify the basin stability of the low concentration region, characterised by the escaping behaviour.