Global entrainment of transcriptional systems to periodic inputs

This paper addresses the problem of giving conditions for transcriptional systems to be globally entrained to external periodic inputs. By using contraction theory, a powerful tool from dynamical systems theory, it is shown that certain systems driven by external periodic signals have the property that all solutions converge to a fixed limit cycle. General results are proved, and the properties are verified in the specific case of some models of transcriptional systems. The basic mathematical results needed from contraction theory are proved in the paper, making it self-contained.

💡 Research Summary

The paper tackles the fundamental question of when a transcriptional regulatory network will globally entrain to an external periodic signal. Using contraction theory—a framework that guarantees exponential convergence of any two trajectories in a suitably defined metric—the authors derive sufficient conditions under which all solutions of a nonlinear dynamical system converge to a unique T‑periodic orbit, regardless of initial conditions.

The theoretical development begins with a general time‑varying system (\dot{x}=f(x,t)). By introducing a uniformly positive definite Riemannian metric (M(x,t)), the authors show that if the symmetric part of the generalized Jacobian, (\frac{1}{2}(MJ+J^{!T}M)), is uniformly negative definite for all states and times, the system is contracting. They provide a self‑contained proof of this result, making the paper accessible to readers without prior expertise in contraction analysis.

Next, the authors consider the case where the system is driven by a periodic input (u(t)=u(t+T)). Because contraction is input‑independent, a contracting system necessarily possesses a unique T‑periodic solution, and every trajectory converges exponentially to this limit cycle. This property is termed global entrainment. The analysis accommodates nonlinear input functions, as the input appears explicitly in the Jacobian and therefore does not break the contraction condition.

To demonstrate the practical relevance of the theory, three canonical transcriptional motifs are examined:

1. Single‑gene negative feedback – a one‑dimensional system where a gene represses its own expression. The Jacobian reduces to a scalar; the contraction condition translates into a simple inequality relating the maximal transcription rate, degradation rate, and Hill coefficient.

2. Toggle switch – a two‑gene mutually repressive circuit. By constructing a diagonal metric, the authors derive explicit bounds on the repression strengths and degradation rates that guarantee both eigenvalues of the symmetric Jacobian are negative, ensuring contraction.

3. Repressilator – a three‑gene cyclic repression network. Although the 3‑dimensional Jacobian is more complex, the authors employ a block‑diagonal metric and a Lyapunov‑like auxiliary function to obtain tractable sufficient conditions. In particular, sufficiently steep Hill functions (Hill exponent ≥ 2) and high protein degradation rates drive the system into a contracting regime.

For each motif, the admissible parameter region is identified analytically and visualized numerically. The authors then subject the models to two types of periodic forcing: a sinusoidal input mimicking light‑induced promoter activation, and a rectangular pulse representing timed drug delivery. Simulations across a wide range of initial conditions consistently converge to the same T‑periodic orbit, confirming the theoretical prediction of global entrainment. Notably, even when the forcing period differs markedly from the intrinsic oscillation period of the unforced system, the trajectories still lock onto the external rhythm provided the contraction conditions hold.

The discussion emphasizes the engineering implications. In synthetic biology, designing a gene circuit that is guaranteed to entrain to an external clock (e.g., a circadian light‑dark cycle) is highly desirable. The contraction‑based criteria give designers a clear, quantitative recipe: select promoter strengths, degradation tags, and Hill coefficients such that the derived matrix inequality is satisfied. This ensures that the circuit will be robust to variations in initial protein concentrations, stochastic fluctuations, and moderate parameter drift. Compared with traditional linearization‑based local stability analysis, contraction offers a truly global guarantee, making it especially valuable for highly nonlinear feedback architectures.

Finally, the paper outlines future directions, including extensions to networks with multiple periodic inputs, explicit time delays, and stochastic noise, as well as experimental validation in living cells. By bridging rigorous dynamical‑systems theory with concrete models of transcriptional regulation, the work provides a solid foundation for the systematic design of globally entrainable synthetic gene oscillators.