This paper employs the root locus method to conduct a detailed investigation of the parameter regions that ensure bistability in a well-studied gene regulatory network namely, lac operon of Escherichia coli (E. coli). In contrast to previous works, the parametric bistability conditions observed in this study constitute a complete set of necessary and sufficient conditions. These conditions were derived by applying the root locus method to the polynomial equilibrium equation of the lac operon model to determine the parameter values yielding the multiple real roots necessary for bistability. The lac operon model used was defined as an ordinary differential equation system in a state equation form with a rational right hand side, and it was compatible with the Hill and Michaelis–Menten approaches of enzyme kinetics used to describe biochemical reactions that govern lactose metabolism. The developed root locus method can be used to study the steady-state behavior of any type of convergent biological system model based on mass action kinetics. This method provides a solution to the problem of analyzing gene regulatory networks under parameter uncertainties because the root locus method considers the model parameters as variable, rather than fixed. The obtained bistability ranges for the lac operon model parameters have the potential to elucidate the appearance of bistability for E. coli cells in in vivo experiments, and they could also be used to design robust hysteretic switches in synthetic biology.