Qualitative simulation of the carbon starvation response in Escherichia coli

In case of nutritional stress, like carbon starvation, Escherichia coli cells abandon their exponential-growth state to enter a more resistant, non-growth state called stationary phase. This growth-phase transition is controlled by a genetic regulatory network integrating various environmental signals. Although E. coli is a paradigm of the bacterial world, it is little understood how its response to carbon starvation conditions emerges from the interactions between the different components of the regulatory network. Using a qualitative method that is able to overcome the current lack of quantitative data on kinetic parameters and molecular concentrations, we model the carbon starvation response network and simulate the response of E. coli cells to carbon deprivation. This allows us to identify essential features of the transition between exponential and stationary phase and to make new predictions on the qualitative system behavior following a carbon upshift.     

Starvation proteins in Escherichia coli: kinetics of synthesis and role in starvation survival

Starvation proteins synthesized by Escherichia coli at the onset of carbon starvation (R. G. Groat and A. Matin, J. Indust. Microbiol. 1:69-73, 1986) exhibited four temporal classes of synthesis in response to glucose or succinate starvation, indicating sequential expression of carbon starvation response (cst) genes. A cst mutant of E. coli showed greatly impaired carbon starvation survival. Thus, it appears that E. coli undergoes a significant molecular realignment in response to starvation, which increases its resistance to this stress. New polypeptides were also synthesized by E. coli in response to phosphate or nitrogen starvation. Some of these polypeptides were unique to a given starvation regimen, but at least 13 appeared to be synthesized regardless of the nutrient deprivation causing the starvation.     

A versatile ODE approximation to a network model for the spread of sexually transmitted diseases

 We develop a moment closure approximation (MCA) to a network model of sexually transmitted disease (STD) spread through a steady/casual partnership network. MCA has been used previously to approximate static, regular lattices, whereas application to dynamic, irregular networks is a new endeavour, and application to sociologically-motivated network models has not been attempted. Our goals are 1) to investigate issues relating to the application of moment closure approximations to dynamic and irregular networks, and 2) to understand the impact of concurrent casual partnerships on STD transmission through a population of predominantly steady monogamous partnerships. We are able to derive a moment closure approximation for a dynamic irregular network representing sexual partnership dynamics, however, we are forced to use a triple approximation due to the large error of the standard pair approximation. This example underscores the importance of doing error analysis for moment closure approximations. We also find that a small number of casual partnerships drastically increases the prevalence and rate of spread of the epidemic. Finally, although the approximation is derived for a specific network model, we can recover approximations to a broad range of network models simply by varying model parameters which control the structure of the dynamic network. Thus our moment closure approximation is very flexible in the kinds of network models it can approximate. Received: 26 August 2001 / Revised version: 15 March 2002 / Published online: 23 August 2002 C.T.B. was supported by the NSF. Key words or phrases: Moment closure approximation – Network model – Pair approximation – Sexually transmitted diseases – Steady/casual partnership network     

Gene expression analysis on biochemical networks using the Potts spin model

MOTIVATION: Microarray technology allows us to profile the expression of a large subset or all genes of a cell. Biochemical research over the last three decades has elucidated an increasingly complete image of the metabolic architecture. For less complex organisms, such as Escherichia coli, the biochemical network has been described in much detail. Here, we investigate the clustering of such networks by applying gene expression data that define edge lengths in the network. RESULTS: The Potts spin model is used as a nearest neighbour based clustering algorithm to discover fragmentation of the network in mutants or in biological samples when treated with drugs. As an example, we tested our method with gene expression data from E.coli treated with tryptophan excess, starvation and trpyptophan repressor mutants. We observed fragmentation of the tryptophan biosynthesis pathway, which corresponds well to the commonly known regulatory response of the cells.     

Comparing Boolean and Piecewise Affine Differential Models for Genetic Networks

Multi-level discrete models of genetic networks, or the more general piecewise affine differential models, provide qualitative information on the dynamics of the system, based on a small number of parameters (such as synthesis and degradation rates). Boolean models also provide qualitative information, but are based simply on the structure of interconnections. To explore the relationship between the two formalisms, a piecewise affine differential model and a Boolean model are compared, for the carbon starvation response network in E. coli. The asymptotic dynamics of both models are shown to be quite similar. This study suggests new tools for analysis and reduction of biological networks.     

Comparing different ODE modelling approaches for gene regulatory networks

A fundamental step in synthetic biology and systems biology is to derive appropriate mathematical models for the purposes of analysis and design. For example, to synthesize a gene regulatory network, the derivation of a mathematical model is important in order to carry out in silico investigations of the network dynamics and to investigate parameter variations and robustness issues. Different mathematical frameworks have been proposed to derive such models. In particular, the use of sets of nonlinear ordinary differential equations (ODEs) has been proposed to model the dynamics of the concentrations of mRNAs and proteins. These models are usually characterized by the presence of highly nonlinear Hill function terms. A typical simplification is to reduce the number of equations by means of a quasi-steady-state assumption on the mRNA concentrations. This yields a class of simplified ODE models. A radically different approach is to replace the Hill functions by piecewise-linear approximations [Casey, R., de Jong, H., Gouz, J.-L., 2006. Piecewise-linear models of genetic regulatory networks: equilibria and their stability. J. Math. Biol. 52 (1), 27-56]. A further modelling approach is the use of discrete-time maps [Coutinho, R., Fernandez, B., Lima, R., Meyroneinc, A., 2006. Discrete time piecewise affine models of genetic regulatory networks. J. Math. Biol. 52, 524-570] where the evolution of the system is modelled in discrete, rather than continuous, time. The aim of this paper is to discuss and compare these different modelling approaches, using a representative gene regulatory network. We will show that different models often lead to conflicting conclusions concerning the existence and stability of equilibria and stable oscillatory behaviours. Moreover, we shall discuss, where possible, the viability of making certain modelling approximations (e.g. quasi-steady-state mRNA dynamics or piecewise-linear approximations of Hill functions) and their effects on the overall system dynamics.     

Firing-rate models for neurons with a broad repertoire of spiking behaviors

Capturing the response behavior of spiking neuron models with rate-based models facilitates the investigation of neuronal networks using powerful methods for rate-based network dynamics. To this end, we investigate the responses of two widely used neuron model types, the Izhikevich and augmented multi-adapative threshold (AMAT) models, to a range of spiking inputs ranging from step responses to natural spike data. We find (i) that linear-nonlinear firing rate models fitted to test data can be used to describe the firing-rate responses of AMAT and Izhikevich spiking neuron models in many cases; (ii) that firing-rate responses are generally too complex to be captured by first-order low-pass filters but require bandpass filters instead; (iii) that linear-nonlinear models capture the response of AMAT models better than of Izhikevich models; (iv) that the wide range of response types evoked by current-injection experiments collapses to few response types when neurons are driven by stationary or sinusoidally modulated Poisson input; and (v) that AMAT and Izhikevich models show different responses to spike input despite identical responses to current injections. Together, these findings suggest that rate-based models of network dynamics may capture a wider range of neuronal response properties by incorporating second-order bandpass filters fitted to responses of spiking model neurons. These models may contribute to bringing rate-based network modeling closer to the reality of biological neuronal networks.     

Deterministic epidemiological models at the individual level

In many fields of science including population dynamics, the vast state spaces inhabited by all but the very simplest of systems can preclude a deterministic analysis. Here, a class of approximate deterministic models is introduced into the field of epidemiology that reduces this state space to one that is numerically feasible. However, these reduced state space master equations do not in general form a closed set. To resolve this, the equations are approximated using closure approximations. This process results in a method for constructing deterministic differential equation models with a potentially large scope of application including dynamic directed contact networks and heterogeneous systems using time dependent parameters. The method is exemplified in the case of an SIR (susceptible-infectious-removed) epidemiological model and is numerically evaluated on a range of networks from spatially local to random. In the context of epidemics propagated on contact networks, this work assists in clarifying the link between stochastic simulation and traditional population level deterministic models.     

RpoS proteolysis is regulated by a mechanism that does not require the SprE (RssB) response regulator phosphorylation site

In Escherichia coli the response regulator SprE (RssB) facilitates degradation of the sigma factor RpoS by delivering it to the ClpXP protease. This process is regulated: RpoS is degraded in logarithmic phase but becomes stable upon carbon starvation, resulting in its accumulation. Because SprE contains a CheY domain with a conserved phosphorylation site (D58), the prevailing model posits that this control is mediated by phosphorylation. To test this model, we mutated the conserved response regulator phosphorylation site (D58A) of the chromosomal allele of sprE and monitored RpoS levels in response to carbon starvation. Though phosphorylation contributed to the SprE basal activity, we found that RpoS proteolysis was still regulated upon carbon starvation. Furthermore, our results indicate that phosphorylation of wild-type SprE occurs by a mechanism that is independent of acetyl phosphate.     

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems-oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta-ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations.     

Quantifying the dynamics of coupled networks of switches and oscillators

Complex network dynamics have been analyzed with models of systems of coupled switches or systems of coupled oscillators. However, many complex systems are composed of components with diverse dynamics whose interactions drive the system's evolution. We, therefore, introduce a new modeling framework that describes the dynamics of networks composed of both oscillators and switches. Both oscillator synchronization and switch stability are preserved in these heterogeneous, coupled networks. Furthermore, this model recapitulates the qualitative dynamics for the yeast cell cycle consistent with the hypothesized dynamics resulting from decomposition of the regulatory network into dynamic motifs. Introducing feedback into the cell-cycle network induces qualitative dynamics analogous to limitless replicative potential that is a hallmark of cancer. As a result, the proposed model of switch and oscillator coupling provides the ability to incorporate mechanisms that underlie the synchronized stimulus response ubiquitous in biochemical systems.     

A Semiquantitative Framework for Gene Regulatory Networks: Increasing the Time and Quantitative Resolution of Boolean Networks

Boolean models have been instrumental in predicting general features of gene networks and more recently also as explorative tools in specific biological applications. In this study we introduce a basic quantitative and a limited time resolution to a discrete (Boolean) framework. Quantitative resolution is improved through the employ of normalized variables in unison with an additive approach. Increased time resolution stems from the introduction of two distinct priority classes. Through the implementation of a previously published chondrocyte network and T helper cell network, we show that this addition of quantitative and time resolution broadens the scope of biological behaviour that can be captured by the models. Specifically, the quantitative resolution readily allows models to discern qualitative differences in dosage response to growth factors. The limited time resolution, in turn, can influence the reachability of attractors, delineating the likely long term system behaviour. Importantly, the information required for implementation of these features, such as the nature of an interaction, is typically obtainable from the literature. Nonetheless, a trade-off is always present between additional computational cost of this approach and the likelihood of extending the model’s scope. Indeed, in some cases the inclusion of these features does not yield additional insight. This framework, incorporating increased and readily available time and semi-quantitative resolution, can help in substantiating the litmus test of dynamics for gene networks, firstly by excluding unlikely dynamics and secondly by refining falsifiable predictions on qualitative behaviour.     

A discrete mathematical model applied to genetic regulation and metabolic networks

This paper describes the use of a discrete mathematical model to represent the basic mechanisms of regulation of the bacteria E. coli in batch fermentation. The specific phenomena studied were the changes in metabolism and genetic regulation when the bacteria use three different carbon substrates (glucose, glycerol, and acetate). The model correctly predicts the behavior of E. coli vis-à-vis substrate mixtures. In a mixture of glucose, glycerol, and acetate, it prefers glucose, then glycerol, and finally acetate. The model included 67 nodes; 28 were genes, 20 enzymes, and 19 regulators/biochemical compounds. The model represents both the genetic regulation and metabolic networks in an inrtegrated form, which is how they function biologically. This is one of the first attempts to include both of these networks in one model. Previously, discrete mathematical models were used only to describe genetic regulation networks. The study of the network dynamics generated 8 (2(3)) fixed points, one for each nutrient configuration (substrate mixture) in the medium. The fixed points of the discrete model reflect the phenotypes described. Gene expression and the patterns of the metabolic fluxes generated are described accurately. The activation of the gene regulation network depends basically on the presence of glucose and glycerol. The model predicts the behavior when mixed carbon sources are utilized as well as when there is no carbon source present. Fictitious jokers (Joker1, Joker2, and Repressor SdhC) had to be created to control 12 genes whose regulation mechanism is unknown, since glycerol and glucose do not act directly on the genes. The approach presented in this paper is particularly useful to investigate potential unknown gene regulation mechanisms; such a novel approach can also be used to describe other gene regulation situations such as the comparison between non-recombinant and recombinant yeast strain, producing recombinant proteins, presently under investigation in our group.     

WebCell: a web-based environment for kinetic modeling and dynamic simulation of cellular networks

SUMMARY: WebCell is a web-based environment for managing quantitative and qualitative information on cellular networks and for interactively exploring their steady-state and dynamic behaviors in response to systemic perturbations. It is designed as a user-friendly web interface, allowing users to efficiently construct, visualize, analyze and store reaction network models, thereby facilitating kinetic modeling and in silico simulation of biological systems of interest. A collected model library is also available to provide comprehensive implications for cellular dynamics of the published models.     

Oscillations and Synchrony in Large-scale Cortical Network Models

Intrinsic neuronal and circuit properties control the responses of large ensembles of neurons by creating spatiotemporal patterns of activity that are used for sensory processing, memory formation, and other cognitive tasks. The modeling of such systems requires computationally efficient single-neuron models capable of displaying realistic response properties. We developed a set of reduced models based on difference equations (map-based models) to simulate the intrinsic dynamics of biological neurons. These phenomenological models were designed to capture the main response properties of specific types of neurons while ensuring realistic model behavior across a sufficient dynamic range of inputs. This approach allows for fast simulations and efficient parameter space analysis of networks containing hundreds of thousands of neurons of different types using a conventional workstation. Drawing on results obtained using large-scale networks of map-based neurons, we discuss spatiotemporal cortical network dynamics as a function of parameters that affect synaptic interactions and intrinsic states of the neurons.     

Maximum Entropy Reconstructions of Dynamic Signaling Networks from Quantitative Proteomics Data

Advances in mass spectrometry among other technologies have allowed for quantitative, reproducible, proteome-wide measurements of levels of phosphorylation as signals propagate through complex networks in response to external stimuli under different conditions. However, computational approaches to infer elements of the signaling network strictly from the quantitative aspects of proteomics data are not well established. We considered a method using the principle of maximum entropy to infer a network of interacting phosphotyrosine sites from pairwise correlations in a mass spectrometry data set and derive a phosphorylation-dependent interaction network solely from quantitative proteomics data. We first investigated the applicability of this approach by using a simulation of a model biochemical signaling network whose dynamics are governed by a large set of coupled differential equations. We found that in a simulated signaling system, the method detects interactions with significant accuracy. We then analyzed a growth factor mediated signaling network in a human mammary epithelial cell line that we inferred from mass spectrometry data and observe a biologically interpretable, small-world structure of signaling nodes, as well as a catalog of predictions regarding the interactions among previously uncharacterized phosphotyrosine sites. For example, the calculation places a recently identified tumor suppressor pathway through ARHGEF7 and Scribble, in the context of growth factor signaling. Our findings suggest that maximum entropy derived network models are an important tool for interpreting quantitative proteomics data.