Exact solution of a model DNA-inversion genetic switch with orientational control

DNA inversion is an important mechanism by which bacteria and bacteriophage switch reversibly between phenotypic states. In such switches, the orientation of a short DNA element is flipped by a site-specific recombinase enzyme. We propose a simple model for a DNA inversion switch in which recombinase production is dependent on the switch state (orientational control). Our model is inspired by the fim switch in Escherichia coli. We present an exact analytical solution of the chemical master equation for the model switch, as well as stochastic simulations. Orientational control causes the switch to deviate from Poissonian behaviour: the distribution of times in the on state shows a peak and successive flip times are correlated.

💡 Research Summary

This paper presents a rigorous analytical treatment of a DNA inversion genetic switch in which the production of the recombinase enzyme is regulated by the orientation of the switch itself—a mechanism the authors term “orientational control.” Inspired by the fim switch of Escherichia coli, the authors construct a stochastic model that captures two discrete switch states (ON and OFF) and the integer copy number of the recombinase enzyme. In the model, enzyme synthesis occurs at a state‑dependent rate (α_ON when the switch is ON, α_OFF when OFF) and degradation proceeds at a constant rate β. Crucially, the rates of switching between ON→OFF and OFF→ON are proportional to the current enzyme copy number, with proportionality constants k₁ and k₂ respectively, reflecting the biological fact that recombination requires enzyme molecules to be present.

The dynamics are described by a chemical master equation (CME) comprising coupled infinite‑dimensional linear differential equations for the joint probability P_S(n,t) of being in state S with n enzyme molecules. To solve the CME exactly, the authors introduce probability generating functions G_S(z,t)=∑_{n≥0}P_S(n,t)zⁿ for each state. By converting the CME into a pair of first‑order partial differential equations in (z,t) and applying the method of characteristics, they obtain closed‑form expressions for G_ON and G_OFF in terms of elementary and special functions (Gamma functions and confluent hypergeometric functions). Imposing normalisation and steady‑state conditions yields explicit formulas for the stationary distribution of enzyme numbers in each state and for the joint distribution of switch state and enzyme count.

From the exact solution the authors derive the dwell‑time distribution for the ON state, f_ON(τ). Unlike a simple exponential (Poisson) distribution expected for memoryless switching, f_ON(τ) exhibits an initial rise, a pronounced peak, and a long tail. This non‑Poissonian shape arises because, when the switch first turns ON, enzyme production is high, leading to a gradual accumulation of recombinase that eventually triggers a rapid OFF transition. The peak position depends on the balance between production (α_ON) and degradation (β), while the height reflects the strength of the enzyme‑dependent switching rate k₁.

The paper also computes the autocorrelation function of successive dwell times, C(Δ)=⟨τ_i τ_{i+Δ}⟩−⟨τ⟩², and finds positive correlations extending over several flips. These correlations are a direct consequence of orientational control: after a switch flips, residual enzyme molecules remain, biasing the timing of the next flip. The analytical predictions are validated by Gillespie stochastic simulations across a wide parameter space. Parameter sweeps show that increasing α_ON relative to α_OFF sharpens the peak and reduces the mean dwell time, while decreasing β (slower degradation) enhances both the peak prominence and the correlation length. Varying the ratio k₁/k₂ modulates the asymmetry between ON→OFF and OFF→ON transitions, further shaping the dwell‑time statistics.

In the discussion, the authors connect their findings to experimental observations of the fim switch, where phase variation and phenotypic heterogeneity are known to be regulated by FimE‑mediated inversion. The non‑Poissonian dwell‑time distribution and the observed temporal correlations provide a plausible mechanistic explanation for the “memory” effect seen in fim expression patterns. The model, however, abstracts away several biological complexities such as multiple recombinases, DNA‑binding accessory proteins, and external signaling cues. The authors acknowledge these limitations and propose extensions that incorporate additional regulatory layers, spatial compartmentalisation, and feedback from downstream gene products.

Overall, the study delivers an exact analytical solution to a biologically motivated stochastic switch, demonstrates how orientational control fundamentally alters switching statistics, and offers a quantitative framework that can be adapted to other phase‑variation systems in bacteria and phage. The work bridges theoretical physics, stochastic processes, and molecular microbiology, and sets the stage for future experimental tests of the predicted dwell‑time peaks and inter‑flip correlations.