Abstract:Background: Bacteria have developed a repertoire of signalling mechanisms that enable adaptive responses to fluctuating environmental conditions. The formation of biofilm, for example, allows persisting in times of external stresses, e.g. induced by antibiotics or a lack of nutrients. Adhesive curli fibers, the major extracellular matrix components in Escherichia coli biofilms, exhibit heterogeneous expression in isogenic cells exposed to identical external conditions. The dynamical mechanisms underlying this heterogeneity remain poorly understood. In this work, we elucidate the potential role of post-translational bistability as a source for this heterogeneity. Results: We introduce a structured modelling workflow combining logical network topology analysis with time-continuous deterministic and stochastic modelling. The aim is to evaluate the topological structure of the underlying signalling network and to identify and analyse model parameterisations that satisfy observations from a set of genetic knockout experiments. Our work supports the hypothesis that the phenotypic heterogeneity of curli expression in biofilm cells is induced by bistable regulation at the post-translational level. Stochastic modelling suggests diverse noise-induced switching behaviours between the stable states, depending on the expression levels of the c-di-GMP-producing (diguanylate cyclases, DGCs) and -degrading (phosphodiesterases, PDEs) enzymes and reveals the quantitative difference in stable c-di-GMP levels between distinct phenotypes. The most dominant type of behaviour is characterised by a fast switching from curli-off to curli-on with a slow switching in the reverse direction and the second most dominant type is a long-term differentiation into curli-on or curli-off cells. This behaviour may implicate an intrinsic feature of the system allowing for a fast adaptive response (curli-on) versus a slow transition to the curli-off state, in line with experimental observations. Conclusion: The combination of logical and continuous modelling enables a thorough analysis of different determinants of bistable regulation, i.e. network topology and biochemical kinetics, and allows for an incorporation of experimental data from heterogeneous sources. Our approach yields a mechanistic explanation for the phenotypic heterogeneity of curli fiber expression. Furthermore, the presented work provides a detailed insight into the interactions between the multiple DGC- and PDE-type enzymes and the role of c-di-GMP in dynamical regulation of cellular decisions.

The regulatory network of E. coli metabolism as a Boolean dynamical system exhibits both homeostasis and flexibility of response

The interdependent network of gene regulation and metabolism is robust where it needs to be

Metabolic regulation is sufficient for global and robust coordination of glucose uptake, catabolism, energy production and growth in Escherichia coli

Dynamics and Design Principles of a Basic Regulatory Architecture Controlling Metabolic Pathways

Feedbacks from the metabolic network to the genetic network reveal regulatory modules in E. coli and B. subtilis

A system-wide network reconstruction of gene regulation and metabolism in Escherichia coli

A simple regulatory network coordinates a bacterial stress response in space and time

A simple regulatory architecture allows learning the statistical structure of a changing environment

Stability of gene regulatory networks

Metabolism Regulated by Network Lyapunov Function

Engineering regulatory networks for complex phenotypes in E. coli

Perturbation-response analysis of in silico metabolic dynamics in nonlinear regime: Hard-coded responsiveness in the cofactors and network sparsity

On the basic computational structure of gene regulatory networks

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism

Synchrony in a Boolean network model of the L-arabinose operon in Escherichia coli

Biological Design Principles of Complex Feedback Modules in the E. coli Ammonia Assimilation System

Bistability and Asynchrony in a Boolean Model of the l-arabinose Operon in Escherichia coli

Two distinct logical types of network control in gene expression profiles

Regulatory Architecture Determines Optimal Regulation of Gene Expression in Metabolic Pathways.

Logical-continuous modelling of post-translationally regulated bistability of curli fiber expression in Escherichia coli