The authors have developed a mathematical model of regulation of expression of the Escherichia coli lac operon, and have investigated bistability in its steady-state induction behaviour in the absence of external glucose. Numerical analysis of equations describing regulation by artificial inducers revealed two natural bistability parameters that can be used to control the range of inducer concentrations over which the model exhibits bistability. By tuning these bistability parameters, the authors found a family of biophysically reasonable systems that are consistent with an experimentally determined bistable region for induction by thio-methylgalactoside (TMG) (in Ozbudak et al. Nature, 2004, 427; p. 737). To model regulation by lactose, the authors developed similar equations in which allolactose, a metabolic intermediate in lactose metabolism and a natural inducer of lac, is the inducer. For biophysically reasonable parameter values, these equations yield no bistability in response to induction by lactose - only systems with an unphysically small permease-dependent export effect can exhibit small amounts of bistability for limited ranges of parameter values. These results cast doubt on the relevance of bistability in the lac operon within the natural context of E. coli, and help shed light on the controversy among existing theoretical studies that address this issue. The results also motivate a deeper experimental characterisation of permease-independent transport of lac inducers, and suggest an experimental approach to address the relevance of bistability in the lac operon within the natural context of E. coli. The sensitivity of lac bistability to the type of inducer emphasises the importance of metabolism in determining the functions of genetic regulatory networks.