Given a stable, steady solution for a gene system with a set of nominal parameters, a relevant engineering question is how to construct a genetic oscillator. For the current system, the analytical solutions for the bifurcation manifolds are available via the Routh-Hurwitz stability criterion ([26], supplementary data). However, such analysis cannot be easily carried out for more complex models. Instead, computational methods have to be used.