Testing the transition state theory in stochastic dynamics of a genetic switch

We present computer-assisted methods for analyzing stochastic models of gene regulatory networks. The main idea that underlies this equation-free analysis is the design and execution of appropriately-initialized short bursts of stochastic simulations; the results of these are processed to estimate coarse-grained quantities of interest, such as mesoscopic transport coefficients. In particular, using a simple model of a genetic toggle switch, we illustrate the computation of an...

Many processes in biology, from the regulation of gene expression in bacteria to memory in the brain, involve switches constructed from networks of biochemical reactions. Crucial molecules are present in small numbers, raising questions about noise and stability. Analysis of noise in simple reaction schemes indicates that switches stable for years and switchable in milliseconds can be built from fewer than one hundred molecules. Prospects for direct tests of this prediction, ...

From a coarse-grained perspective the motif of a self-activating species, activating a second species which acts as its own repressor, is widely found in biological systems, in particular in genetic systems with inherent oscillatory behavior. Here we consider a specific realization of this motif as a genetic circuit, in which genes are described as directly producing proteins, leaving out the intermediate step of mRNA production. We focus on the effect that inherent time scal...

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 transcriptio...

We present a detailed analysis, based on the Forward Flux Sampling (FFS) simulation method, of the switching dynamics and stability of two models of genetic toggle switches, consisting of two mutually-repressing genes encoding transcription factors (TFs); in one model (the exclusive switch), they mutually exclude each other's binding, while in the other model (general switch) the two transcription factors can bind simultaneously to the shared operator region. We assess the ro...

The nonequilibrium stationary state of an exclusive genetic switch is considered. The model comprises two competing species and a single binding site which, when bound to by a protein of one species, causes the other species to be repressed. The model may be thought of as a minimal model of the power struggle between two competing parties. Exact solutions are given for the limits of vanishing binding/unbinding rates and infinite binding/unbinding rates. A mean field theory is...

Two synthetic gene circuits -- the genetic toggle switch and the repressilator -- are analyzed quantitatively and discussed in the context of an educational module on gene circuits and feedback that constitutes the final topic of a year-long introductory physics sequence, aimed at biology and premedical undergraduate students. The genetic toggle switch consists of two genes, each of whose protein product represses the other's expression, while the repressilator consists of th...

The lactose uptake-pathway of E. coli is a paradigmatic example of multistability in gene-regulatory circuits. In the induced state of the lac-pathway, the genes comprising the lac-operon are transcribed, leading to the production of proteins which import and metabolize lactose. In the uninduced state, a stable repressor-DNA loop frequently blocks the transcription of the lac-genes. Transitions from one phenotypic state to the other are driven by fluctuations, which arise fro...

We present an analytical treatment of a genetic switch model consisting of two mutually inhibiting genes operating without cooperative binding of the corresponding transcription factors. Previous studies have numerically shown that these systems can exhibit bimodal dynamics without possessing two stable fixed points at the deterministic level. We analytically show that bimodality is induced by the noise and find the critical repression strength that controls a transition betw...

Molecular noise in gene regulatory networks has two intrinsic components, one part being due to fluctuations caused by the birth and death of protein or mRNA molecules which are often present in small numbers and the other part arising from gene state switching, a single molecule event. Stochastic dynamics of gene regulatory circuits appears to be largely responsible for bifurcations into a set of multi-attractor states that encode different cell phenotypes. The interplay of ...