No Arabic abstract
We have developed a mathematical model of regulation of expression of the Escherichia coli lac operon, and have investigated bistability in its steady-state induction behavior 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, we found a family of biophysically reasonable systems that are consistent with an experimentally determined bistable region for induction by thio-methylgalactoside (Ozbudak et al. Nature 427:737, 2004). The model predicts that bistability can be abolished when passive transport or permease export becomes sufficiently large; the former case is especially relevant to induction by isopropyl-beta, D-thiogalactopyranoside. To model regulation by lactose, we 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; however, 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 suggest an experimental approach to address the relevance of bistability in the lac operon within the natural context of E. coli.
The lac operon in Escherichia coli has been studied extensively and is one of the earliest gene systems found to undergo both positive and negative control. The lac operon is known to exhibit bistability, in the sense that the operon is either induced or uninduced. Many dynamical models have been proposed to capture this phenomenon. While most are based on complex mathematical formulations, it has been suggested that for other gene systems network topology is sufficient to produce the desired dynamical behavior. We present a Boolean network as a discrete model for the lac operon. We include the two main glucose control mechanisms of catabolite repression and inducer exclusion in the model and show that it exhibits bistability. Further we present a reduced model which shows that lac mRNA and lactose form the core of the lac operon, and that this reduced model also exhibits the same dynamics. This work corroborates the claim that the key to dynamical properties is the topology of the network and signs of interactions.
The lactose operon in Escherichia coli was the first known gene regulatory network, and it is frequently used as a prototype for new modeling paradigms. Historically, many of these modeling frameworks use differential equations. More recently, Stigler and Veliz-Cuba proposed a Boolean network model that captures the bistability of the system and all of the biological steady states. In this paper, we model the well-known arabinose operon in E. coli with a Boolean network. This has several complex features not found in the lac operon, such as a protein that is both an activator and repressor, a DNA looping mechanism for gene repression, and the lack of inducer exclusion by glucose. For 11 out of 12 choices of initial conditions, we use computational algebra and Sage to verify that the state space contains a single fixed point that correctly matches the biology. The final initial condition, medium levels of arabinose and no glucose, successfully predicts the systems bistability. Finally, we compare the state space under synchronous and asynchronous update, and see that the former has several artificial cycles that go away under a general asynchronous update.
The analysis of stress response systems in microorganisms can reveal molecular strategies for regulatory control and adaptation. Here, we focus on the Cad module, a subsystem of E. colis response to acidic stress, which is conditionally activated at low pH only when lysine is available. When expressed, the Cad system counteracts the elevated H+ concentration by converting lysine to cadaverine under the consumption of H+, and exporting cadaverine in exchange for external lysine. Surprisingly, the cad operon displays a transient response, even when the conditions for its induction persist. To quantitatively characterize the regulation of the Cad module, we have experimentally recorded and theoretically modeled the dynamics of important system variables. We establish a quantitative model that adequately describes and predicts the transient expression behavior for various initial conditions. Our quantitative analysis of the Cad system supports a negative feedback by external cadaverine as the origin of the transient response. Furthermore, the analysis puts causal constraints on the precise mechanism of signal transduction via the regulatory protein CadC.
We show how the competition between sensing and adaptation can result in a performance peak in E.coli chemotaxis using extensive numerical simulations in a detailed theoretical model. Receptor clustering amplifies the input signal coming from ligand binding which enhances chemotactic efficiency. But large clusters also induce large fluctuations in total activity since the number of clusters go down. The activity and hence the run-tumble motility now gets controlled by methylation levels which are part of adaptation module, rather than ligand binding. This reduces chemotactic efficiency.
Cellular metabolism, the integrated interconversion of thousands of metabolic substrates through enzyme-catalyzed biochemical reactions, is the most investigated complex intercellular web of molecular interactions. While the topological organization of individual reactions into metabolic networks is increasingly well understood, the principles governing their global functional utilization under different growth conditions pose many open questions. We implement a flux balance analysis of the E. coli MG1655 metabolism, finding that the network utilization is highly uneven: while most metabolic reactions have small fluxes, the metabolisms activity is dominated by several reactions with very high fluxes. E. coli responds to changes in growth conditions by reorganizing the rates of selected fluxes predominantly within this high flux backbone. The identified behavior likely represents a universal feature of metabolic activity in all cells, with potential implications to metabolic engineering.