Bifurcation Analysis for Phage Lambda with Binding Energy Uncertainty

Abstract

In a phage λ genetic switch model, bistable dynamical behavior can be destroyed due to the bifurcation caused by inappropriately chosen model parameters. Since the values of many parameters with biological significance often cannot be accurately acquired, it is thus of fundamental importance to analyze how and to which extent the system dynamics is influenced by model parameters, especially those parameters pertaining to binding energies. In this paper, we apply a Jacobian method to investigate the relation between bifurcation and parameter uncertainties on a phage λ OR model. By introducing bistable range as a measure of system robustness, we find that RNA polymerase binding energies have the minimum bistable ranges among all the binding energies, which implies that the uncertainties on these parameters tend to demolish the bistability more easily. Moreover, parameters describing mutual prohibition between proteins Cro and CI have finite bistable ranges, whereas those describing self-prohibition have infinity bistable ranges. Hence, the former are more sensitive to parameter uncertainties than the latter. These results help to understand the influence of parameter uncertainties on the bifurcation and thus bistability.