Robust model matching design methodology for a stochastic synthetic gene network

Synthetic biology has shown its potential and promising applications in the last decade. However, many synthetic gene networks cannot work properly and maintain their desired behaviors due to intrinsic parameter variations and extrinsic disturbances. In this study, the intrinsic parameter uncertainties and external disturbances are modeled in a non-linear stochastic gene network to mimic the real environment in the host cell. Then a non-linear stochastic robust matching design methodology is introduced to withstand the intrinsic parameter fluctuations and to attenuate the extrinsic disturbances in order to achieve a desired reference matching purpose. To avoid solving the Hamilton-Jacobi inequality (HJI) in the non-linear stochastic robust matching design, global linearization technique is used to simplify the design procedure by solving a set of linear matrix inequalities (LMIs). As a result, the proposed matching design methodology of the robust synthetic gene network can be efficiently designed with the help of LMI toolbox in Matlab. Finally, two in silico design examples of the robust synthetic gene network are given to illustrate the design procedure and to confirm the robust model matching performance to achieve the desired behavior in spite of stochastic parameter fluctuations and environmental disturbances in the host cell.     

Balanced input allows optimal encoding in a stochastic binary neural network model: an analytical study

Recent neurophysiological experiments have demonstrated a remarkable effect of attention on the underlying neural activity that suggests for the first time that information encoding is indeed actively influenced by attention. Single cell recordings show that attention reduces both the neural variability and correlations in the attended condition with respect to the non-attended one. This reduction of variability and redundancy enhances the information associated with the detection and further processing of the attended stimulus. Beyond the attentional paradigm, the local activity in a neural circuit can be modulated in a number of ways, leading to the general question of understanding how the activity of such circuits is sensitive to these relatively small modulations. Here, using an analytically tractable neural network model, we demonstrate how this enhancement of information emerges when excitatory and inhibitory synaptic currents are balanced. In particular, we show that the network encoding sensitivity--as measured by the Fisher information--is maximized at the exact balance. Furthermore, we find a similar result for a more realistic spiking neural network model. As the regime of balanced inputs has been experimentally observed, these results suggest that this regime is functionally important from an information encoding standpoint.     

Resonances and noise in a stochastic Hindmarsh-Rose model of thalamic neurons

Thalamic neurons exhibit subthreshold resonance when stimulated with small sine wave signals of varying frequency and stochastic resonance when noise is added to these signals. We study a stochastic Hindmarsh-Rose model using Monte-Carlo simulations to investigate how noise, in conjunction with subthreshold resonance, leads to a preferred frequency in the firing pattern. The resulting stochastic resonance (SR) exhibits a preferred firing frequency that is approximately exponential in its dependence on the noise amplitude. In similar experiments, frequency dependent SR is found in the reliability of detection of alpha-function inputs under noise, which are more realistic inputs for neurons. A mathematical analysis of the equations reveals that the frequency preference arises from the dynamics of the slow variable. Noise can then transfer the resonance over the firing threshold because of the proximity of the fast subsystem to a Hopf bifurcation point. Our results may have implications for the behavior of thalamic neurons in a network, with noise switching the membrane potential between different resonance modes.     

Floral morphogenesis: stochastic explorations of a gene network epigenetic landscape

In contrast to the classical view of development as a preprogrammed and deterministic process, recent studies have demonstrated that stochastic perturbations of highly non-linear systems may underlie the emergence and stability of biological patterns. Herein, we address the question of whether noise contributes to the generation of the stereotypical temporal pattern in gene expression during flower development. We modeled the regulatory network of organ identity genes in the Arabidopsis thaliana flower as a stochastic system. This network has previously been shown to converge to ten fixed-point attractors, each with gene expression arrays that characterize inflorescence cells and primordial cells of sepals, petals, stamens, and carpels. The network used is binary, and the logical rules that govern its dynamics are grounded in experimental evidence. We introduced different levels of uncertainty in the updating rules of the network. Interestingly, for a level of noise of around 0.5-10%, the system exhibited a sequence of transitions among attractors that mimics the sequence of gene activation configurations observed in real flowers. We also implemented the gene regulatory network as a continuous system using the Glass model of differential equations, that can be considered as a first approximation of kinetic-reaction equations, but which are not necessarily equivalent to the Boolean model. Interestingly, the Glass dynamics recover a temporal sequence of attractors, that is qualitatively similar, although not identical, to that obtained using the Boolean model. Thus, time ordering in the emergence of cell-fate patterns is not an artifact of synchronous updating in the Boolean model. Therefore, our model provides a novel explanation for the emergence and robustness of the ubiquitous temporal pattern of floral organ specification. It also constitutes a new approach to understanding morphogenesis, providing predictions on the population dynamics of cells with different genetic configurations during development.     

Improved dimensionally-reduced visual cortical network using stochastic noise modeling

In this paper, we extend our framework for constructing low-dimensional dynamical system models of large-scale neuronal networks of mammalian primary visual cortex. Our dimensional reduction procedure consists of performing a suitable linear change of variables and then systematically truncating the new set of equations. The extended framework includes modeling the effect of neglected modes as a stochastic process. By parametrizing and including stochasticity in one of two ways we show that we can improve the systems-level characterization of our dimensionally reduced neuronal network model. We examined orientation selectivity maps calculated from the firing rate distribution of large-scale simulations and stochastic dimensionally reduced models and found that by using stochastic processes to model the neglected modes, we were able to better reproduce the mean and variance of firing rates in the original large-scale simulations while still accurately predicting the orientation preference distribution.     

Modeling a hox gene network in silico using a stochastic simulation algorithm

The amount of molecular information that has been gathered about Hox cis-regulatory mechanisms allows us to take the next important step: integrating the results and constructing a higher-level model for the interaction and regulation of the Hox genes. Here, we present the results of our investigation into a cis-regulatory network for the early Hox genes. Instead of using conventional differential equation approaches for analyzing the system, we have adopted the use of a stochastic simulation algorithm (SSA) to model the network. The model allows us to track in detail the behavior of each component of a biochemical pathway and to produce computerized movies of the time evolution of the system that is a result of the dynamic interplay of these various components. The simulation is able to reproduce key features of the wild-type pattern of gene expression, and in silico experiments yield results similar to their corresponding in vivo experiments. This analysis shows the utility of using stochastic methods to model biochemical networks. In addition, the model has suggested several intriguing new results that are currently being investigated in vivo.     

Intrinsic noise, gene regulation and steady-state statistics in a two-gene network

The intrinsic noise in a two-gene network model is analysed. The technique of the Fokker-Planck approximation is used to investigate the statistics of noise when the system state is near a stable equilibrium. This is called also the steady-state statistics. The relative size of noise is measured by the Fano factor that is defined as the ratio of the variance to the mean. Our main result shows that in general, the noise control in a two-gene network might be a very complicated process, but for the repressor-repressor system that is a very important case in investigating the genetic switch, the relative size of noise, i.e. the Fano factor, must be bigger than one for both the repressor proteins.     

A simple stochastic gene substitution model

If the fitnesses of n haploid alleles in a finite population are assigned at random and if the alleles can mutate to one another, and if the population is initially fixed for the kth most fit allele, then the mean number of substitutions that will occur before the most fit allele is fixed is shown to be

$\text{1}\text{2}\text{+}\text{1}\text{k}\text{+}\text{∑}\text{i=2}\text{k?1}\text{(i+3)}\text{(2i(i+1))}$

when selection is strong and mutation is weak. This result is independent of the parameters that went into the model. The result is used to provide a partial explanation for the large variance observed in the rates of molecular evolution.     

Metabolic robustness and network modularity: a model study

Background

Several studies have mentioned network modularity—that a network can easily be decomposed into subgraphs that are densely connected within and weakly connected between each other—as a factor affecting metabolic robustness. In this paper we measure the relation between network modularity and several aspects of robustness directly in a model system of metabolism.

Methodology/Principal Findings

By using a model for generating chemical reaction systems where one can tune the network modularity, we find that robustness increases with modularity for changes in the concentrations of metabolites, whereas it decreases with changes in the expression of enzymes. The same modularity scaling is true for the speed of relaxation after the perturbations.

Conclusions/Significance

Modularity is not a general principle for making metabolism either more or less robust; this question needs to be addressed specifically for different types of perturbations of the system.     

The IUD follow-up study: A stochastic model

Recurrence equations are derived for the change in genotype frequencies for general fecundity functions. Statistical models that will result in fecundity functions factorable into male and female contributions are developed. Some attention is given to the genetic structure and equilibrium behavior of single-locus genetic systems with additive fecundities.     

And the noise played on: stochastic gene expression and HIV-1 infection

Stochastic gene expression has been implicated in a variety of cellular processes, including cell differentiation and disease. In this issue of Cell, take an integrated computational-experimental approach to study the Tat transactivation feedback loop of HIV-1. They show that fluctuations in a key regulator, Tat, in an isogenic population of infected cells result in two distinct expression states corresponding to latent and productive HIV-1 infection. These findings demonstrate the importance of stochastic gene expression in molecular "decision-making."     

Listening to the noise: random fluctuations reveal gene network parameters

The cellular environment is abuzz with noise originating from the inherent random motion of reacting molecules in the living cell. In this noisy environment, clonal cell populations show cell‐to‐cell variability that can manifest significant phenotypic differences. Noise‐induced stochastic fluctuations in cellular constituents can be measured and their statistics quantified. We show that these random fluctuations carry within them valuable information about the underlying genetic network. Far from being a nuisance, the ever‐present cellular noise acts as a rich source of excitation that, when processed through a gene network, carries its distinctive fingerprint that encodes a wealth of information about that network. We show that in some cases the analysis of these random fluctuations enables the full identification of network parameters, including those that may otherwise be difficult to measure. This establishes a potentially powerful approach for the identification of gene networks and offers a new window into the workings of these networks.     

Mixture distributions in a stochastic gene expression model with delayed feedback: a WKB approximation approach

Journal of Mathematical Biology - Noise in gene expression can be substantively affected by the presence of production delay. Here we consider a mathematical model with bursty production of...     

Parameter identification for a stochastic SEIRS epidemic model: case study influenza

Journal of Mathematical Biology - A recent parameter identification technique, the local lagged adapted generalized method of moments, is used to identify the time-dependent disease transmission...     

Varieties of stochastic model: a comparative study of the Gompertz effect

A comparative study is made of various models for the Gompertz phenomenon, which is a form of growth rate limitation in population dynamics. Deterministic, Markov birth-death, diffusion and stochastic differential equation models are studied, with a view to assessing their advantages and limitations.