Resilience of the autocatalytic feedback loop for gene regulation

Abstract: Gene expression in response to stimuli is regulated by transcription factors (TFs) through feedback loop motifs, aimed at maintaining the desired TF concentration despite uncertainties and perturbations. In this work, we consider a stochastic model of the autocatalytic feedback loop for gene regulation and we probabilistically quantify its resilience, \textit{i.e.}, its ability to preserve the equilibrium associated with a prescribed concentration of TFs, and the corresponding basin of attraction, in the presence of noise. We show that the formation of larger oligomers, corresponding to larger Hill coefficients of the regulation function, and thus to sharper non-linearities, improves the system resilience, even close to critical concentrations of TFs. Our formal results within a stochastic formulation relying on the Fokker-Planck equation are accompanied by numerical simulations.