Cellular Adaptation as Learning: Gene Self-Regulation During Signal Fluctuations ✛ ✦

Original: "Cellular Adaptation to Signal Fluctuations as Learning"
· Tuan Minh Pham, David Saad
arXiv:2607.03545 · 2026-07-03 · CC BY 4.0 · 4 min
A mathematical model shows that cellular self-regulation under fluctuating external signals is equivalent to feedback learning, leading to stable collective correlations.

Context

In development, cells must make fateful decisions based on external chemical signals—morphogens. However, morphogen concentration in tissue is subject to temporal fluctuations, threatening the stability of cellular identity. Until now, it has remained an open question how a cell can stably interpret a noisy signal without resorting to rigid noise suppression. New research tackles this problem through the lens of machine learning: an intracellular signaling cascade with feedback acts as an adaptive controller, tuning sensitivity to the signal. This allows an analogy with neural network training and the application of statistical physics methods originally introduced by Ludwig Boltzmann.

Methods

Using the Shh signaling pathway in the vertebrate neural tube as an example, the authors built a two-level model. Fast variables—concentrations of target genes—evolve according to stochastic equations with sigmoidal activation, forming a multistable gene regulatory network. Slow variables—the effective sensitivity of effectors—update via gradient descent, minimizing the squared mismatch between total gene activity and effectors. In the limit of a large number of signaling molecules, a closed system of equations for macroscopic correlations (covariances) between pathway components is derived. This numerical simulation allowed the exploration of the system's phase behavior depending on the adaptation rate η and feedback strength γ. The approach relies on statistical physics methods of learning and the concept of entropy as a measure of mismatch.

Results

Analysis revealed three regimes. With weak adaptation (η < 0.9), all correlations decay, and the cell remains unstable to fluctuations. As η increases, a symmetry-broken phase appears: each effector molecule becomes strongly correlated with one specific target gene and anticorrelated with others. This is reminiscent of spontaneous ordering in a Bose-Einstein condensate, where a macroscopic number of particles occupy a single quantum state. Diagonal elements of covariance matrices T and Q significantly exceed off-diagonal ones, and the cross-covariance matrix R indicates selective affinity. At even higher η and γ, undamped oscillations arise—an analog of plasma oscillations due to excessive feedback. The critical threshold also depends on network size: with increasing K, the transition to the structured regime occurs at lower η. Importantly, the stability phase corresponds not to reduced noise variance but to the formation of structured correlations—similar to how in quantum information entanglement generates order. Signal noise was modeled as an uncorrelated stochastic process, reminiscent of the cosmic microwave background—a random backdrop against which meaningful structure stands out.

Implications

The work offers a new perspective on cellular stability: instead of passive noise filtering, the cell actively shapes internal correlations, aligning the states of effectors and target genes. This self-organization mechanism, interpreted as learning, explains how a cell can maintain identity under continuously changing external conditions. Conceptually, the results echo Ludwig Boltzmann's idea of a system's tendency towards the most probable state, as well as Waddington's epigenetic landscape where cell types correspond to attractors. Practically, the predicted collective modes can be tested in single-cell sequencing experiments, opening the way to quantitative verification of the theory.

Future development

A natural extension will include more realistic details: time-dependent adaptation rate, epigenetic modifications, diverse network motifs, and consideration of evolving gene networks. It would be interesting to investigate an optimal schedule of adaptation η(τ) for the fastest mismatch suppression—an optimal control problem akin to training neural networks with variable step size. Furthermore, scaling the theory to hundreds of genes and calibrating with experimental data will allow predicting cell-type distributions in developing tissues.

Impact

The results will impact developmental biology (understanding tissue formation), synthetic biology (designing robust gene circuits), and bioinformatics (interpretation of sequencing data).

Next steps

Experimental verification of predictions is needed through analysis of transcription factor covariances in cell populations, as well as calibration of parameters η and γ using single-cell RNA-seq data.

Key open problems

The work touches on the fundamental problem of linking microscopic interactions to macroscopic order in living systems, as well as the question of how a cell maintains identity under constantly fluctuating signals—analogous to the problem of attractor stability in nonlinear dynamical systems.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors built a two-level model of the Shh signaling pathway in the neural tube, where fast variables (gene concentrations) evolve stochastically, while slow ones (effector sensitivity) are updated via gradient descent, minimizing mismatch. In the limit of a large number of molecules, equations for macroscopic correlations are derived, and numerically three regimes are found: decay of correlations under weak adaptation, a symmetry-broken phase with strong selective correlation at intermediate adaptation speed, and oscillations under excessive feedback.

Where its strength is

This work transplants the concept of learning from machine learning into intracellular regulation, showing that robustness to noise can be achieved not by suppressing noise but by shaping structured correlations through sensitivity adaptation. It provides a concrete mechanism and a phase diagram that can be tested and generalized.

Where it could go next

Worth trying

  • Test whether the symmetry-broken phase arises with other activation functions, e.g., a steeper sigmoid.
  • Apply the covariance matrix formalism to experimental data on gene expression variability in the neural tube.
  • Investigate the effect of heterogeneous adaptation rates of different effectors on the phase diagram.

Why this may matter

If the proposed adaptive self-regulation mechanism is confirmed in living systems, it will change our understanding of how cells maintain identity under noise and provide principles for building robust synthetic gene circuits.

Original abstract

as written by the authors arXiv:2607.03545

Cells represent one of the most fundamental units of life. Underlying their robust performance against environmental variability, such as temporal fluctuations of chemical signals, between different cell types, is a dynamical interrelation between the two components of an intracellular pathway: a gene-regulatory network and its upstream signal transducers. To understand how a single cell utilizes this feedback to self-regulate its gene-expressions, we develop a multiscale model of the pathway's components, in which the adaptive variables responsible for signal interpretation follow a feedback-induced learning process. We then derive a macroscopic theory capturing the covariations between these components - so-called collective modes. Our theory shows how cells can achieve robust output against signal fluctuations via self-regulation rather than simple noise suppression. Such robustness corresponds to a transition from random- to structured collective modes beyond a critical adaptation rate.

Vertebrate neural tube cells interpret the same Shh morphogen gradient differently depending on context: the same signal can activate different genetic programs, much like how a single note sounds different in different musical keys.

J_{k\alpha}(\tau_{\mu+1}) = J_{k\alpha}(\tau_\mu) - \frac{\eta}{N} \frac{\partial \epsilon(\tau_\mu)}{\partial J_{k\alpha}(\tau_\mu)}
Updating the sensitivity of effector k to signal component α is proportional to the learning rate η and the gradient of the mismatch ε.

Key numbers

  • Number of target genes in the model: 4
  • Critical adaptation rate η: ≈ 0.9 at γ=0.5
  • Number of averaged random matrices for phase diagram: 100
  • Noise level in fast dynamics (σ): 0.01
  • Network size (number of genes and effectors): K=4, tested up to K=16
Scientists
Andrey Markov
Concepts
cell neuron living cells collective cell migration cellular neural interface control theory neural language acquisition biological stochastic dynamics stimulus hormone feedback loop biological rhythm modeling real time control learning gene expression peptide gene species Reinforcement Learning modulation