Optimized bifurcation diagram for 4 gene model: β1 = 22.186, β2 = 4, β3 = 22.145, β4 = 0.1.