Gene regulatory networks control when, where, and how much protein each gene produces, orchestrating everything from cell differentiation during development to daily circadian rhythms. Transcription factors are key proteins that bind DNA at specific promoter sequences and either activate or repress transcription of target genes. A single transcription factor can regulate hundreds of genes, and genes regulate transcription factors, creating feedback loops of remarkable complexity and precision.
The simplest regulatory motifs—network sub-graphs with defined input-output relationships—include the autoregulation loop (a transcription factor represses or activates its own gene), the toggle switch (two mutually repressing genes that produce bistability—two stable states), and the repressilator (three genes repressing each other in a ring, producing sustained oscillations analogous to a biological clock). These minimal circuits appear repeatedly across organisms from bacteria to mammals, suggesting that evolution has discovered optimal network architectures.
This simulator models the dynamics of small gene networks using ordinary differential equations for mRNA and protein concentrations, including Hill-function cooperativity terms. You can build toggle switches, repressilators, and feedforward loops, tune production/degradation rates and cooperativity coefficients, and observe how initial conditions determine which stable state is reached (bistability), or watch oscillations emerge from the repressilator when the cooperativity is high enough—demonstrating key principles of synthetic biology and developmental gene circuits.
What is a transcription factor and how does it control gene expression?
A transcription factor (TF) is a protein that binds specific DNA sequences (promoter or enhancer regions) near a target gene and influences whether RNA polymerase transcribes that gene. Activating TFs recruit polymerase and general transcription machinery, increasing mRNA production. Repressing TFs block polymerase access or recruit chromatin-compacting complexes that silence the gene. Because TFs are themselves encoded by genes, TFs can regulate each other, creating complex regulatory networks. Most development is orchestrated by combinatorial TF codes: different combinations of TFs determine cell type and gene expression programs.
What is bistability in gene networks and what causes it?
Bistability means a gene network has two stable steady states and can switch between them depending on initial conditions or transient signals. The genetic toggle switch—two mutually repressing genes A and B—is bistable: if A is high, it represses B (keeping B low), which allows A to remain high; conversely if B is high, it represses A. Both states are self-sustaining. Bistability requires positive feedback (either mutual repression or self-activation) and typically requires cooperativity (Hill coefficient n > 1) in binding. Cell fate decisions during development often exploit bistability for robust, irreversible commitment to one of two lineages.
What is the repressilator and how does it produce oscillations?
The repressilator, designed by Elowitz and Leibler (2000), is a synthetic three-gene network where gene A represses gene B, gene B represses gene C, and gene C represses gene A—a ring of repressors. When A is high, it suppresses B; low B allows C to rise; high C suppresses A; low A allows B to rise; high B suppresses C; and so on, producing sustained oscillations in protein concentrations with a period determined by production/degradation time constants. This was the first engineered synthetic gene circuit to demonstrate that simple network topology generates oscillation, and it informed our understanding of biological clocks.
The Hill function H(x) = xⁿ/(Kⁿ + xⁿ) models the sigmoidal response of gene transcription to transcription factor concentration x, where K is the concentration for half-maximal activation and n is the Hill coefficient measuring cooperativity. At n=1 (no cooperativity), the response is a hyperbola—gradual. At n=2 (e.g., TF binding as a dimer), the response is steeper. At n → ∞, the response becomes a step function (switch-like). High cooperativity (large n) is essential for bistability and oscillation in gene circuits, and is achieved biologically through TF oligomerization or multiple TF binding sites in a promoter.
Synthetic biology applies engineering principles to design new biological functions by assembling standardized genetic parts (promoters, ribosome-binding sites, coding sequences, terminators). Gene regulation models help predict whether designed circuits will function as intended, identify parameter ranges for bistability or oscillation, and optimize timing and amplitude. Synthetic circuits have been designed to detect cancer biomarkers and kill tumor cells, produce high-value biochemicals (insulin, artemisinin), create biosensors for environmental toxins, and build biological computers. Mathematical modeling of gene regulation is essential because biological parts interact in complex, nonlinear ways that intuition alone cannot navigate.
As n rises from 1 toward 6, the Hill function H(x) = xⁿ/(Kⁿ+xⁿ) becomes an increasingly sharp, switch-like response instead of a gradual curve. On the mutual-repression Toggle switch preset, a higher n makes bistability more robust: the two stable states (gene A high/B low, or B high/A low) become more sharply separated and small perturbations are less likely to flip the switch, whereas at low n the network can settle into an intermediate, less clearly-committed state.
β (1 to 40) sets the maximum rate new protein is synthesized, while γ (0.2 to 4) sets the rate it is degraded or diluted away. In the absence of repression, a gene settles to a steady-state concentration of roughly β/γ, so raising β or lowering γ both push the equilibrium level higher. The ratio β/γ relative to the Threshold K also determines whether a gene can turn fully on: if β/γ is much larger than K, the gene reaches strong activation once past its threshold.
K (0.5 to 20) is the transcription-factor concentration at which the Hill function reaches half its maximal value — it shifts the whole sigmoidal response curve left or right along the concentration axis without changing its steepness (that is controlled separately by n). A low K means the gene switches on at low TF concentrations, while a high K means much more transcription factor must accumulate before the target gene responds significantly.
No — Speed (0.5× to 8×) only scales how fast the simulation's internal clock advances relative to real time on screen, purely for viewing convenience. The actual oscillation period on the Repressilator preset is set by the production and degradation time constants (β and γ) acting on the three-gene repression ring, not by the Speed control, so changing Speed makes the same oscillation play faster or slower without altering its underlying period in simulated time units.
In a coherent feed-forward loop, an upstream regulator activates a target gene both directly and indirectly through an intermediate gene, and the target only turns on when both paths agree. A brief input pulse can activate the direct path but disappears before the slower intermediate path accumulates enough product to also cross its threshold, so the target gene filters out the transient noise; only an input that persists long enough for both paths to align switches the target gene on, making this network motif a natural noise filter.