Gene-gene cooperativity in small networks

We show how to construct a reduced description of interacting genes in noisy, small regulatory networks using coupled binary “spin” variables. Treating both the protein number and gene expression state variables stochastically and on equal footing we propose a mapping which connects the molecular level description of networks to the binary representation. We construct a phase diagram indicating when genes can be considered to be independent and when the coupling between them cannot be neglected leading to synchrony or correlations. We find that an appropriately mapped boolean description reproduces the probabilities of gene expression states of the full stochastic system very well and can be transfered to examples of self-regulatory systems with a larger number of gene copies.

💡 Research Summary

The paper presents a systematic method for reducing stochastic models of small gene regulatory networks to a binary “spin” representation while retaining the essential noise inherent in protein copy numbers and transcriptional states. By treating protein abundance and gene expression as coupled stochastic variables, the authors derive a mapping that translates the full master‑equation description into an effective Hamiltonian of the form

 P(s, n) ≈ Z⁻¹ exp(∑₁ⁿ h_i s_i + ∑{i<j} J{ij} s_i s_j) · ∏_i Poisson(n_i | λ_i(s)),

where s_i ∈ {0,1} denotes the on/off state of gene i, h_i is an effective field encoding basal transcription and degradation rates, and J_{ij} captures the regulatory influence of gene j on gene i. The protein numbers n_i follow Poisson distributions whose means λ_i(s) depend on the current spin configuration, thereby preserving the discrete, noisy nature of molecular counts.

The authors first examine the minimal case of two interacting genes. They show that positive J corresponds to activation, negative J to repression, and that the magnitude of J relative to the intrinsic noise (set by the average protein copy number) determines whether the genes behave independently or exhibit correlated, synchronized dynamics. By scanning protein copy number (N_p) and the ratio of transcription to degradation rates (α), they construct a phase diagram that delineates a “cooperative region” (strong coupling, low copy number) from an “independent region” (weak coupling, high copy number). In the cooperative region, the binary model reproduces the full stochastic joint distribution with errors below 5 %, confirming that the spin mapping captures the essential physics of gene‑gene interaction.

The framework is then extended to self‑regulatory motifs. Autoregulation introduces a self‑interaction term J_{ii}. Positive self‑coupling (positive feedback) yields bistability, while negative self‑coupling (negative feedback) enforces monostability. Simulations demonstrate that the spin model accurately predicts the probability of each transcriptional state and the switching rates between them, again within a few percent of the full master‑equation results.

Finally, the authors address networks with multiple identical copies of a gene (gene duplication). Because each copy shares the same h and J parameters, the system can be treated in a mean‑field approximation, dramatically reducing computational complexity. As the number of copies increases, fluctuations in the collective spin diminish, leading to more robust expression levels. The binary representation continues to match the full stochastic model’s predictions for the distribution of expression states, even when the number of copies reaches ten.

Overall, the study demonstrates that a carefully constructed spin‑based Boolean model provides a highly efficient yet quantitatively accurate description of noisy, small‑scale gene regulatory circuits. It identifies clear criteria—protein copy number, transcription/degradation ratios, and interaction strength—under which genes can be treated as independent versus when their coupling must be retained. This approach offers a valuable tool for synthetic biology design, enabling rapid exploration of circuit behavior, and for theoretical investigations of how stochasticity and cooperativity shape cellular decision‑making.