View PDF HTML (experimental)

Abstract:The logistic two-gene negative-feedback oscillator is locally asymptotically stable for all biological parameter values, since the trace of the Jacobian is uniformly negative. Real biological oscillators (circadian rhythms, the segmentation clock, Hes1, p53) nevertheless rely on delays. We extend the logistic two-gene model to a delay-differential system with transcriptional delays $\tau_1$ and $\tau_2$, and prove that the equilibrium loses stability through a Hopf bifurcation as the total delay $\tau=\tau_1+\tau_2$ crosses an explicit critical value $\tau_c$. The Hopf frequency $\omega_c$ and $\tau_c$ are computed in closed form from the logistic derivatives; the loop-gain condition $AB>\gamma_1\gamma_2$ is necessary and sufficient; the transversality $\mathrm{Re}(d\mu/d\tau)|_{\tau_c}>0$ admits a parameter-uniform positive lower bound; and the bifurcation persists globally. A sum-of-delays symmetry reduces the analysis to the scalar parameter $\tau$. Numerical simulations confirm three regimes (damped, small limit cycle, relaxation), the supercritical amplitude scaling $A\sim c\sqrt{\tau-\tau_c}$, and the deep-relaxation period asymptote $T\sim 2\tau+C_\infty$ with closed-form offset $C_\infty$. For the symmetric-threshold loop, supercriticality is proved by a Lindstedt--Poincaré reduction yielding closed-form amplitude and frequency laws; for the general asymmetric loop it delivers a closed-form first Lyapunov coefficient and an explicit criticality criterion. Calibrated to p53--Mdm2 data, the closed-form Hopf period matches the observed oscillation within $3\%$, and the standard Hill-function model within a few percent. The analysis extends to cyclic $N$-gene loops, with a closed-form transversality rate valid for every $N$ and -- in the symmetric case -- an explicit delay-induced-Hopf window $\gamma^N<\Lambda<\gamma^N\sec^N(\pi/N)$.

Submission history

From: Ismail Belgacem [view email]
[v1] Fri, 22 May 2026 15:02:52 UTC (244 KB)