Most of the biological systems including gene regulatory networks can be described well by ordinary differential equation models with rational non-linearities. These models are derived either based on the reaction kinetics or by curve fitting to experimental data. This study demonstrates the applicability of the root-locus-based bifurcation analysis method for studying the complex dynamics of such models. The effectiveness of the bifurcation analysis in determining the exact parameter regions in each of which the system shows a certain dynamical behaviour, such as bistability, oscillation, and asymptotically equilibrium dynamics is shown by considering two mostly studied gene regulatory networks, namely Gardner's genetic toggle switch and p53 gene network possessing two-phase (mono-stable/oscillation) dynamics.