Coherence of biochemical oscillations is bounded by driving force and network topology

A topological mechanism for robust and efficient global oscillations in biological networks

Long and stable timescales are often observed in complex biochemical networks, such as in emergent oscillations. How these robust dynamics persist remains unclear, given the many stochastic reactions and shorter time scales demonstrated by underlying components. We propose a topological model that produces long oscillations around the network boundary, reducing the system dynamics to a lower-dimensional current in a robust manner. Using this to model KaiC, which regulates the circadian rhythm in cyanobacteria, we compare the coherence of oscillations to that in other KaiC models. Our topological model localizes currents on the system edge, with an efficient regime of simultaneously increased precision and decreased cost. Further, we introduce a new predictor of coherence from the analysis of spectral gaps, and show that our model saturates a global thermodynamic bound. Our work presents a new mechanism and parsimonious description for robust emergent oscillations in complex biological networks.

Thermodynamic cost for precision of general counting observables

We analytically derive universal bounds that describe the tradeoff between thermodynamic cost and precision in a sequence of events related to some internal changes of an otherwise hidden physical system. The precision is quantified by the fluctuations in either the number of events counted over time or the waiting times between successive events. Our results are valid for the same broad class of nonequilibrium driven systems considered by the thermodynamic uncertainty relation, but they extend to both time-symmetric and asymmetric observables. We show how optimal precision saturating the bounds can be achieved. For waiting-time fluctuations of asymmetric observables, a phase transition in the optimal configuration arises, where higher precision can be achieved by combining several signals.

Thermodynamic bounds on the asymmetry of cross-correlations with dynamical activity and entropy production

Entropy production and dynamical activity are two complementary aspects in nonequilibrium physics. The asymmetry of cross-correlation, serving as a distinctive feature of nonequilibrium, also finds widespread utility. In this Letter, we establish two thermodynamic bounds on the normalized asymmetry of cross-correlation in terms of dynamical activity and entropy production rate. These bounds demonstrate broad applicability and offer experimental testability.

Discontinuous phase transition from ferromagnetic to oscillating states in a nonequilibrium mean-field spin model

We study a nonequilibrium ferromagnetic mean-field spin model exhibiting a phase with spontaneous temporal oscillations of the magnetization, on top of the usual paramagnetic and ferromagnetic phases. This behavior is obtained by introducing dynamic field variables coupled to the spins through nonreciprocal couplings. We determine a nonequilibrium generalization of the Landau free energy in terms of the large deviation function of the magnetization and of an appropriately defined smoothed stochastic time derivative of the magnetization. While the transition between paramagnetic and oscillating phase is continuous, the transition between ferromagnetic and oscillating phases is found to be discontinuous, with a coexistence of both phases, one being stable and the other one metastable. Depending on parameter values, the ferromagnetic points may either be inside or outside the limit cycle, leading to different transition scenarios. The stability of these steady states is determined from the large deviation function. We also show that in the coexistence region, the entropy production has a pronounced maximum as a function of system size.

Topologically constrained fluctuations and thermodynamics regulate nonequilibrium response

The limits on a system's response to external perturbations inform our understanding of how physical properties can be shaped by microscopic characteristics. Here, we derive constraints on the steady-state nonequilibrium response of physical observables in terms of the topology of the microscopic state space and the strength of thermodynamic driving. Notably, evaluation of these limits requires no kinetic information beyond the state-space structure. When applied to models of receptor binding, we find that sensitivity is bounded by the steepness of a Hill function with a Hill coefficient enhanced by the chemical driving beyond the structural equilibrium limit.

Entropy production limits all fluctuation oscillations

The oscillation of fluctuation with two state observables is investigated. Following the idea of Ohga et al. [Phys. Rev. Lett. 131, 077101 (2023)10.1103/PhysRevLett.131.077101], we find that the fluctuation oscillation relative to their autocorrelations is bounded from above by the entropy production per characteristic maximum oscillation time. Our result applies to a variety of systems including Langevin systems, chemical reaction systems, and macroscopic systems. In addition, our bound consists of experimentally tractable quantities, which enables us to examine our inequality experimentally.

Thermodynamic Bound on the Asymmetry of Cross-Correlations

The principle of microscopic reversibility says that, in equilibrium, two-time cross-correlations are symmetric under the exchange of observables. Thus, the asymmetry of cross-correlations is a fundamental, measurable, and often-used statistical signature of deviation from equilibrium. Here we find a simple and universal inequality that bounds the magnitude of asymmetry by the cycle affinity, i.e., the strength of thermodynamic driving. Our result applies to a large class of systems and all state observables, and it suggests a fundamental thermodynamic cost for various nonequilibrium functions quantified by the asymmetry. It also provides a powerful tool to infer affinity from measured cross-correlations, in a different and complementary way to the thermodynamic uncertainty relations. As an application, we prove a thermodynamic bound on the coherence of noisy oscillations, which was previously conjectured by Barato and Seifert [Phys. Rev. E 95, 062409 (2017)PRESCM2470-004510.1103/PhysRevE.95.062409]. We also derive a thermodynamic bound on directed information flow in a biochemical signal transduction model.

Nonequilibrium Phase Transition to Temporal Oscillations in Mean-Field Spin Models

We propose a mean-field theory to describe the nonequilibrium phase transition to a spontaneously oscillating state in spin models. A nonequilibrium generalization of the Landau free energy is obtained from the joint distribution of the magnetization and its smoothed stochastic time derivative. The order parameter of the transition is a Hamiltonian, whose nonzero value signals the onset of oscillations. The Hamiltonian and the nonequilibrium Landau free energy are determined explicitly from the stochastic spin dynamics. The oscillating phase is also characterized by a nontrivial overlap distribution reminiscent of a continuous replica symmetry breaking, in spite of the absence of disorder. An illustration is given on an explicit kinetic mean-field spin model.

Thermodynamic bounds on ultrasensitivity in covalent switching

Switch-like motifs are among the basic building blocks of biochemical networks. A common motif that can serve as an ultrasensitive switch consists of two enzymes acting antagonistically on a substrate, one making and the other removing a covalent modification. To work as a switch, such covalent modification cycles must be held out of thermodynamic equilibrium by continuous expenditure of energy. Here, we exploit the linear framework for timescale separation to establish tight bounds on the performance of any covalent-modification switch in terms of the chemical potential difference driving the cycle. The bounds apply to arbitrary enzyme mechanisms, not just Michaelis-Menten, with arbitrary rate constants and thereby reflect fundamental physical constraints on covalent switching.

Universal minimal cost of coherent biochemical oscillations

Biochemical clocks are essential for virtually all living systems. A biochemical clock that is isolated from an external periodic signal and subjected to fluctuations can oscillate coherently only for a finite number of oscillations. Furthermore, such an autonomous clock can oscillate only if it consumes free energy. What is the minimum amount of free-energy consumption required for a certain number of coherent oscillations? We conjecture a universal bound that answers this question. A system that oscillates coherently for N oscillations has a minimal free-energy cost per oscillation of 4π^{2}Nk_{B}T. Our bound is valid for general finite Markov processes, is conjectured based on extensive numerical evidence, is illustrated with numerical simulations of a known model for a biochemical oscillator, and applies to existing experimental data.

Coherence of oscillations in the weak-noise limit

In a noisy environment, oscillations lose their coherence, which can be characterized by a quality factor. We determine this quality factor for oscillations arising from a driven Fokker-Planck dynamics along a periodic one-dimensional potential analytically in the weak-noise limit. With this expression, we can prove for this continuum model the analog of an upper bound that has been conjectured for the coherence of oscillations in discrete Markov network models. We show that our approach can also be adapted to motion along a noisy two-dimensional limit cycle. Specifically, we apply our scheme to the noisy Stuart-Landau oscillator and the thermodynamically consistent Brusselator as a simple model for a chemical clock. Our approach thus complements the fairly sophisticated extant general framework based on techniques from Hamilton-Jacobi theory with which we compare our results numerically.

Kinetic uncertainty relation on first-passage time for accumulated current

The kinetic uncertainty relation (KUR) is a trade-off relation between the precision of an observable and the mean dynamical activity in a fixed time interval for a time-homogeneous and continuous-time Markov chain. In this Letter, we derive the KUR on the first passage time for the time-integrated current from the information inequality at stopping times. The relation shows that the precision of the first passage time is bounded from above by the mean number of jumps up to that time. We apply our result to simple systems and demonstrate that the activity constraint gives a tighter bound than the thermodynamic uncertainty relation in the regime far from equilibrium.

Precision and dissipation of a stochastic Turing pattern

Spontaneous pattern formation is a fundamental scientific problem that has received much attention since the seminal theoretical work of Turing on reaction-diffusion systems. In molecular biophysics, this phenomenon often takes place under the influence of large fluctuations. It is then natural to inquire about the precision of such pattern. In particular, spontaneous pattern formation is a nonequilibrium phenomenon, and the relation between the precision of a pattern and the thermodynamic cost associated with it remains largely unexplored. Here, we analyze this relation with a paradigmatic stochastic reaction-diffusion model, i.e., the Brusselator in one spatial dimension. We find that the precision of the pattern is maximized for an intermediate thermodynamic cost, i.e., increasing the thermodynamic cost beyond this value makes the pattern less precise. Even though fluctuations get less pronounced with an increase in thermodynamic cost, we argue that larger fluctuations can also have a positive effect on the precision of the pattern.

Robust oscillations in multi-cyclic Markov state models of biochemical clocks

Organisms often use cyclic changes in the concentrations of chemical species to precisely time biological functions. Underlying these biochemical clocks are chemical reactions and transport processes, which are inherently stochastic. Understanding the physical basis for robust biochemical oscillations in the presence of fluctuations has thus emerged as an important problem. In a previous paper [C. del Junco and S. Vaikuntanathan, Phys. Rev. E 101, 012410 (2020)], we explored this question using the non-equilibrium statistical mechanics of single-ring Markov state models of biochemical networks that support oscillations. Our finding was that they can exploit non-equilibrium driving to robustly maintain the period and coherence of oscillations in the presence of randomness in the rates. Here, we extend our work to Markov state models consisting of a large cycle decorated with multiple small cycles. These additional cycles are intended to represent alternate pathways that the oscillator may take as it fluctuates about its average path. Combining a mapping to single-cycle networks based on first passage time distributions with our previously developed theory, we are able to make analytical predictions for the period and coherence of oscillations in these networks. One implication of our predictions is that a high energy budget can make different network topologies and arrangements of rates degenerate as far as the period and coherence of oscillations are concerned. Excellent agreement between analytical and numerical results confirms that this is the case. Our results suggest that biochemical oscillators can be more robust to fluctuations in the path of the oscillator when they have a high energy budget.

Reaching and violating thermodynamic uncertainty bounds in information engines

Thermodynamic uncertainty relations (TURs) set fundamental bounds on the fluctuation and dissipation of stochastic systems. Here, we examine these bounds, in experiment and theory, by exploring the entire phase space of a cyclic information engine operating in a nonequilibrium steady state. Close to its maximal efficiency, we find that the engine violates the original TUR. This experimental demonstration of TUR violation agrees with recently proposed softer bounds: The engine satisfies two generalized TUR bounds derived from the detailed fluctuation theorem with feedback control and another bound linking fluctuation and dissipation to mutual information and Renyi divergence. We examine how the interplay of work fluctuation and dissipation shapes the information conversion efficiency of the engine, and find that dissipation is minimal at a finite noise level, where the original TUR is violated.

Stochastic thermodynamics of chemical reactions coupled to finite reservoirs: A case study for the Brusselator

Biomolecular processes are typically modeled using chemical reaction networks coupled to infinitely large chemical reservoirs. A difference in chemical potential between these reservoirs can drive the system into a non-equilibrium steady-state (NESS). In reality, these processes take place in finite systems containing a finite number of molecules. In such systems, a NESS can be reached with the help of an externally driven pump for which we introduce a simple model. The crucial parameters are the pumping rate and the finite size of the chemical reservoir. We apply this model to a simple biochemical oscillator, the Brusselator, and quantify the performance using the number of coherent oscillations. As a surprising result, we find that higher precision can be achieved with finite-size reservoirs even though the corresponding current fluctuations are larger than in the ideal infinite case.

The thermodynamic uncertainty relation in biochemical oscillations

Living systems regulate many aspects of their behaviour through periodic oscillations of molecular concentrations, which function as 'biochemical clocks.' The chemical reactions that drive these clocks are intrinsically stochastic at the molecular level, so that the duration of a full oscillation cycle is subject to random fluctuations. Their success in carrying out their biological function is thought to depend on the degree to which these fluctuations in the cycle period can be suppressed. Biochemical oscillators also require a constant supply of free energy in order to break detailed balance and maintain their cyclic dynamics. For a given free energy budget, the recently discovered 'thermodynamic uncertainty relation' yields the magnitude of period fluctuations in the most precise conceivable free-running clock. In this paper, we show that computational models of real biochemical clocks severely underperform this optimum, with fluctuations several orders of magnitude larger than the theoretical minimum. We argue that this suboptimal performance is due to the small number of internal states per molecule in these models, combined with the high level of thermodynamic force required to maintain the system in the oscillatory phase. We introduce a new model with a tunable number of internal states per molecule and confirm that it approaches the optimal precision as this number increases.

Efficiency fluctuations and noise induced refrigerator-to-heater transition in information engines

Understanding noisy information engines is a fundamental problem of non-equilibrium physics, particularly in biomolecular systems agitated by thermal and active fluctuations in the cell. By the generalized second law of thermodynamics, the efficiency of these engines is bounded by the mutual information passing through their noisy feedback loop. Yet, direct measurement of the interplay between mutual information and energy has so far been elusive. To allow such examination, we explore here the entire phase-space of a noisy colloidal information engine, and study efficiency fluctuations due to the stochasticity of the mutual information and extracted work. We find that the average efficiency is maximal for non-zero noise level, at which the distribution of efficiency switches from bimodal to unimodal, and the stochastic efficiency often exceeds unity. We identify a line of anomalous, noise-driven equilibrium states that defines a refrigerator-to-heater transition, and test the generalized integral fluctuation theorem for continuous engines.

High chemical affinity increases the robustness of biochemical oscillations

Biochemical oscillations are ubiquitous in nature and allow organisms to properly time their biological functions. In this paper, we consider minimal Markov state models of nonequilibrium biochemical networks that support oscillations. We obtain analytical expressions for the coherence and period of oscillations in these networks. These quantities are expected to depend on all details of the transition rates in the Markov state model. However, our analytical calculations reveal that driving the system out of equilibrium makes many of these details-specifically, the location and arrangement of the transition rates-irrelevant to the coherence and period of oscillations. This theoretical prediction is confirmed by excellent agreement with numerical results. As a consequence, the coherence and period of oscillations can be robustly maintained in the presence of fluctuations in the irrelevant variables. While recent work has established that increasing energy consumption improves the coherence of oscillations, our findings suggest that it plays the additional role of making the coherence and the average period of oscillations robust to fluctuations in rates that can result from the noisy environment of the cell.

Energy dissipation bounds for autonomous thermodynamic cycles

How much free energy is irreversibly lost during a thermodynamic process? For deterministic protocols, lower bounds on energy dissipation arise from the thermodynamic friction associated with pushing a system out of equilibrium in finite time. Recent work has also bounded the cost of precisely moving a single degree of freedom. Using stochastic thermodynamics, we compute the total energy cost of an autonomously controlled system by considering both thermodynamic friction and the entropic cost of precisely directing a single control parameter. Our result suggests a challenge to the usual understanding of the adiabatic limit: Here, even infinitely slow protocols are energetically irreversible.

Oscillations in feedback-driven systems: Thermodynamics and noise

Oscillations in nonequilibrium noisy systems are important physical phenomena. These oscillations can happen in autonomous biochemical oscillators such as circadian clocks. They can also manifest as subharmonic oscillations in periodically driven systems such as time crystals. Oscillations in autonomous systems and, to a lesser degree, subharmonic oscillations in periodically driven systems have been both thoroughly investigated, including their relation with thermodynamic cost and noise. We perform a systematic study of oscillations in a third class of nonequilibrium systems: feedback-driven systems. In particular, we use the apparatus of stochastic thermodynamics to investigate the role of noise and thermodynamic cost in feedback-driven oscillations. For a simple two-state model that displays oscillations, we analyze the relation between precision and dissipation, revealing that oscillations can remain coherent for an indefinite time in a finite system with thermal fluctuations in a limit of diverging thermodynamic cost. We consider oscillations in a more complex system with several degrees of freedom, an Ising model driven by feedback between the magnetization and the external field. This feedback-driven system can display subharmonic oscillations similar to the ones observed in time crystals. We illustrate the second law for feedback-driven systems that display oscillations. For the Ising model, the oscillating dissipated heat can be negative. However, when we consider the total entropy that also includes an informational term related to measurements, the oscillating total entropy change is always positive. We also study the finite-size scaling of the dissipated heat, providing evidence for the existence of a first-order phase transition for certain parameter regimes.

The energy cost and optimal design for synchronization of coupled molecular oscillators

A model of coupled molecular biochemical oscillators is proposed to study nonequilibrium thermodynamics of synchronization. We find that synchronization of nonequilibrium oscillators costs addition energy to drive the exchange reaction (chemical interaction) between individual oscillators. By solving the steady state of the many-body system analytically, we show that the system goes through a nonequilibrium phase transition driven by energy dissipation, and the critical energy dissipation depends on both the frequency and strength of the exchange reaction. Moreover, our study reveals the optimal design for achieving maximum synchronization with a fixed energy budget. We apply our general theory to the Kai system in Cyanobacteria circadian clock and predict a relationship between the KaiC ATPase activity and synchronization of the KaiC hexamers. The theoretical framework can be extended to study thermodynamics of collective behaviors in other extended nonequilibrium active systems.

Interlinked GTPase cascades provide a motif for both robust switches and oscillators

GTPases regulate a wide range of cellular processes, such as intracellular vesicular transport, signal transduction and protein translation. These hydrolase enzymes operate as biochemical switches by toggling between an active guanosine triphosphate (GTP)-bound state and an inactive guanosine diphosphate (GDP)-bound state. We compare two network motifs, a single-species switch and an interlinked cascade that consists of two species coupled through positive and negative feedback loops. We find that interlinked cascades are closer to the ideal all-or-none switch and are more robust against fluctuating signals. While the single-species switch can only achieve bistability, interlinked cascades can be converted into oscillators by tuning the cofactor concentrations, which catalyse the activity of the cascade. These regimes can only be achieved with sufficient chemical driving provided by GTP hydrolysis. In this study, we present a thermodynamically consistent model that can achieve bistability and oscillations with the same feedback motif.

Subharmonic oscillations in stochastic systems under periodic driving

Subharmonic response is a well-known phenomenon in, e.g., deterministic nonlinear dynamical systems. We investigate the conditions under which such subharmonic oscillations can persist for a long time in open systems with stochastic dynamics due to thermal fluctuations. In contrast to stochastic autonomous systems in a stationary state, for which the number of coherent oscillations is fundamentally bounded by the number of states in the underlying network, we demonstrate that in periodically driven systems, subharmonic oscillations can in principle remain coherent forever, even in networks with a small number of states. We also show that, inter alia, the thermodynamic cost rises only logarithmically with the number of coherent oscillations in a model calculation and that the possible periods of the persistent subharmonic response grow linearly with the number of states. We argue that our results can be relevant for biochemical oscillations and for stochastic models of time crystals.

Phase transition in thermodynamically consistent biochemical oscillators

Biochemical oscillations are ubiquitous in living organisms. In an autonomous system, not influenced by an external signal, they can only occur out of equilibrium. We show that they emerge through a generic nonequilibrium phase transition, with a characteristic qualitative behavior at criticality. The control parameter is the thermodynamic force which must be above a certain threshold for the onset of biochemical oscillations. This critical behavior is characterized by the thermodynamic flux associated with the thermodynamic force, its diffusion coefficient, and the stationary distribution of the oscillating chemical species. We discuss metrics for the precision of biochemical oscillations by comparing two observables, the Fano factor associated with the thermodynamic flux and the number of coherent oscillations. Since the Fano factor can be small even when there are no biochemical oscillations, we argue that the number of coherent oscillations is more appropriate to quantify the precision of biochemical oscillations. Our results are obtained with three thermodynamically consistent versions of known models: the Brusselator, the activator-inhibitor model, and a model for KaiC oscillations.

Quantifying fluctuations in reversible enzymatic cycles and clocks

Biochemical reactions are fundamentally noisy at a molecular scale. This limits the precision of reaction networks, but it also allows fluctuation measurements that may reveal the structure and dynamics of the underlying biochemical network. Here, we study nonequilibrium reaction cycles, such as the mechanochemical cycle of molecular motors, the phosphorylation cycle of circadian clock proteins, or the transition state cycle of enzymes. Fluctuations in such cycles may be measured using either of two classical definitions of the randomness parameter, which we show to be equivalent in general microscopically reversible cycles. We define a stochastic period for reversible cycles and present analytical solutions for its moments. Furthermore, we associate the two forms of the randomness parameter with the thermodynamic uncertainty relation, which sets limits on the timing precision of the cycle in terms of thermodynamic quantities. Our results should prove useful also for the study of temporal fluctuations in more general networks.

Design principles for enhancing phase sensitivity and suppressing phase fluctuations simultaneously in biochemical oscillatory systems

Biological systems need to function accurately in the presence of strong noise and at the same time respond sensitively to subtle external cues. Here we study design principles in biochemical oscillatory circuits to achieve these two seemingly incompatible goals. We show that energy dissipation can enhance phase sensitivity linearly by driving the phase-amplitude coupling and increase timing accuracy by suppressing phase diffusion. Two general design principles in the key underlying reaction loop formed by two antiparallel pathways are found to optimize oscillation performance with a given energy budget: balancing the forward-to-backward flux ratio between the two pathways to reduce phase diffusion and maximizing the net flux of the phase-advancing pathway relative to that of the phase-retreating pathway to enhance phase sensitivity. Experimental evidences consistent with these design principles are found in the circadian clock of cyanobacteria. Future experiments to test the predicted dependence of phase sensitivity on energy dissipation are proposed.

Fundamental Bounds on First Passage Time Fluctuations for Currents

Current is a characteristic feature of nonequilibrium systems. In stochastic systems, these currents exhibit fluctuations constrained by the rate of dissipation in accordance with the recently discovered thermodynamic uncertainty relation. Here, we derive a conjugate uncertainty relationship for the first passage time to accumulate a fixed net current. More generally, we use the tools of large-deviation theory to simply connect current fluctuations and first passage time fluctuations in the limit of long times and large currents. With this connection, previously discovered symmetries and bounds on the large-deviation function for currents are readily transferred to first passage times.