Landscape and global stability of nonadiabatic and adiabatic oscillations in a gene network

We quantify the potential landscape to determine the global stability and coherence of biological oscillations. We explore a gene network motif in our experimental synthetic biology studies of two genes that mutually repress and activate each other with self-activation and self-repression. We find that in addition to intrinsic molecular number fluctuations, there is another type of fluctuation crucial for biological function: the fluctuation due to the slow binding/unbinding of protein regulators to gene promoters. We find that coherent limit cycle oscillations emerge in two regimes: an adiabatic regime with fast binding/unbinding and a nonadiabatic regime with slow binding/unbinding relative to protein synthesis/degradation. This leads to two mechanisms of producing the stable oscillations: the effective interactions from averaging the gene states in the adiabatic regime; and the time delays due to slow binding/unbinding to promoters in the nonadiabatic regime, which can be tested by forthcoming experiments. In both regimes, the landscape has a topological shape of the Mexican hat in protein concentrations that quantitatively determines the global stability of limit cycle dynamics. The oscillation coherence is shown to be correlated with the shape of the Mexican hat characterized by the height from the oscillation ring to the central top. The oscillation period can be tuned in a wide range by changing the binding/unbinding rate without changing the amplitude much, which is important for the functionality of a biological clock. A negative feedback loop with time delays due to slow binding/unbinding can also generate oscillations. Although positive feedback is not necessary for generating oscillations, it can make the oscillations more robust.