Deterministic and stochastic population-level simulations of an artificial lac operon genetic network

Enhanced electrochemical measurement of β-galactosidase activity in whole cells by coexpression of lactose permease, LacY

Whole-cell biosensing links the sensing and computing capabilities of microbes to the generation of a detectable reporter. Whole cells enable dynamic biological computation (filtered noise, amplified signals, logic gating etc.). Enzymatic reporters enable in situ signal amplification. Electrochemical measurements are easily quantified and work in turbid environments. In this work we show how the coexpression of the lactose permease, LacY, dramatically improves electrochemical sensing of β-galactosidase (LacZ) expressed as a reporter in whole cells. The permease facilitates transport of the LacZ substrate, 4-aminophenyl β-d-galactopyranoside, which is converted to redox active p-aminophenol, which, in turn, is detected via cyclic voltammetry or chronocoulometry. We show a greater than fourfold improvement enabled by lacY coexpression in cells engineered to respond to bacterial signal molecules, pyocyanin and quorum-sensing autoinducer-2.

Inherent regulatory asymmetry emanating from network architecture in a prevalent autoregulatory motif

Predicting gene expression from DNA sequence remains a major goal in the field of gene regulation. A challenge to this goal is the connectivity of the network, whose role in altering gene expression remains unclear. Here, we study a common autoregulatory network motif, the negative single-input module, to explore the regulatory properties inherited from the motif. Using stochastic simulations and a synthetic biology approach in E. coli, we find that the TF gene and its target genes have inherent asymmetry in regulation, even when their promoters are identical; the TF gene being more repressed than its targets. The magnitude of asymmetry depends on network features such as network size and TF-binding affinities. Intriguingly, asymmetry disappears when the growth rate is too fast or too slow and is most significant for typical growth conditions. These results highlight the importance of accounting for network architecture in quantitative models of gene expression.

Phenotypic instability in fungi

Fungi are prone to phenotypic instability, that is, the vegetative phase of these organisms, be they yeasts or molds, undergoes frequent switching between two or more behaviors, often with different morphologies, but also sometime having different physiologies without any obvious morphological outcome. In the context of industrial utilization of fungi, this can have a negative impact on the maintenance of strains and/or on their productivity. Instabilities have been shown to result from various mechanisms, either genetic or epigenetic. This chapter will review different types of instabilities and discuss some lesser-known ones, mostly in filamentous fungi, while it will direct readers to additional literature in the case of well-known phenomena such as the amyloid prions or fungal senescence. It will present in depth the "white/opaque" switch of Candida albicans and the "crippled growth" degeneration of the model fungus Podospora anserina. These are two of the most thoroughly studied epigenetic phenotypic switches. I will also discuss the "sectors" presented by many filamentous ascomycetes, for which a prion-based model exists but is not demonstrated. Finally, I will also describe intriguing examples of phenotypic instability for which an explanation has yet to be provided.

Effect of Intrinsic Noise on the Phenotype of Cell Populations Featuring Solution Multiplicity: An Artificial lac Operon Network Paradigm

Heterogeneity in cell populations originates from two fundamentally different sources: the uneven distribution of intracellular content during cell division, and the stochastic fluctuations of regulatory molecules existing in small amounts. Discrete stochastic models can incorporate both sources of cell heterogeneity with sufficient accuracy in the description of an isogenic cell population; however, they lack efficiency when a systems level analysis is required, due to substantial computational requirements. In this work, we study the effect of cell heterogeneity in the behaviour of isogenic cell populations carrying the genetic network of lac operon, which exhibits solution multiplicity over a wide range of extracellular conditions. For such systems, the strategy of performing solely direct temporal solutions is a prohibitive task, since a large ensemble of initial states needs to be tested in order to drive the system--through long time simulations--to possible co-existing steady state solutions. We implement a multiscale computational framework, the so-called "equation-free" methodology, which enables the performance of numerical tasks, such as the computation of coarse steady state solutions and coarse bifurcation analysis. Dynamically stable and unstable solutions are computed and the effect of intrinsic noise on the range of bistability is efficiently investigated. The results are compared with the homogeneous model, which neglects all sources of heterogeneity, with the deterministic cell population balance model, as well as with a stochastic model neglecting the heterogeneity originating from intrinsic noise effects. We show that when the effect of intrinsic source of heterogeneity is intensified, the bistability range shifts towards higher extracellular inducer concentration values.

On the impact of discreteness and abstractions on modelling noise in gene regulatory networks

In this paper, we explore the impact of different forms of model abstraction and the role of discreteness on the dynamical behaviour of a simple model of gene regulation where a transcriptional repressor negatively regulates its own expression. We first investigate the relation between a minimal set of parameters and the system dynamics in a purely discrete stochastic framework, with the twofold purpose of providing an intuitive explanation of the different behavioural patterns exhibited and of identifying the main sources of noise. Then, we explore the effect of combining hybrid approaches and quasi-steady state approximations on model behaviour (and simulation time), to understand to what extent dynamics and quantitative features such as noise intensity can be preserved.

Determining the bistability parameter ranges of artificially induced lac operon using the root locus method

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.

What makes the lac-pathway switch: identifying the fluctuations that trigger phenotype switching in gene regulatory systems

Multistable gene regulatory systems sustain different levels of gene expression under identical external conditions. Such multistability is used to encode phenotypic states in processes including nutrient uptake and persistence in bacteria, fate selection in viral infection, cell-cycle control and development. Stochastic switching between different phenotypes can occur as the result of random fluctuations in molecular copy numbers of mRNA and proteins arising in transcription, translation, transport and binding. However, which component of a pathway triggers such a transition is generally not known. By linking single-cell experiments on the lactose-uptake pathway in E. coli to molecular simulations, we devise a general method to pinpoint the particular fluctuation driving phenotype switching and apply this method to the transition between the uninduced and induced states of the lac-genes. We find that the transition to the induced state is not caused only by the single event of lac-repressor unbinding, but depends crucially on the time period over which the repressor remains unbound from the lac-operon. We confirm this notion in strains with a high expression level of the lac-repressor (leading to shorter periods over which the lac-operon remains unbound), which show a reduced switching rate. Our techniques apply to multistable gene regulatory systems in general and allow to identify the molecular mechanisms behind stochastic transitions in gene regulatory circuits.

Coarse-grained analysis of stochastically simulated cell populations with a positive feedback genetic network architecture

Among the different computational approaches modelling the dynamics of isogenic cell populations, discrete stochastic models can describe with sufficient accuracy the evolution of small size populations. However, for a systematic and efficient study of their long-time behaviour over a wide range of parameter values, the performance of solely direct temporal simulations requires significantly high computational time. In addition, when the dynamics of the cell populations exhibit non-trivial bistable behaviour, such an analysis becomes a prohibitive task, since a large ensemble of initial states need to be tested for the quest of possibly co-existing steady state solutions. In this work, we study cell populations which carry the lac operon network exhibiting solution multiplicity over a wide range of extracellular conditions (inducer concentration). By adopting ideas from the so-called "equation-free" methodology, we perform systems-level analysis, which includes numerical tasks such as the computation of coarse steady state solutions, coarse bifurcation analysis, as well as coarse stability analysis. Dynamically stable and unstable macroscopic (population level) steady state solutions are computed by means of bifurcation analysis utilising short bursts of fine-scale simulations, and the range of bistability is determined for different sizes of cell populations. The results are compared with the deterministic cell population balance model, which is valid for large populations, and we demonstrate the increased effect of stochasticity in small size populations with asymmetric partitioning mechanisms.