Dissipation in noisy chemical networks: The role of deficiency

Milestoning estimators of dissipation in systems observed at a coarse resolution

Many nonequilibrium, active processes are observed at a coarse-grained level, where different microscopic configurations are projected onto the same observable state. Such "lumped" observables display memory, and in many cases, the irreversible character of the underlying microscopic dynamics becomes blurred, e.g., when the projection hides dissipative cycles. As a result, the observations appear less irreversible, and it is very challenging to infer the degree of broken time-reversal symmetry. Here we show, contrary to intuition, that by ignoring parts of the already coarse-grained state space we may-via a process called milestoning-improve entropy-production estimates. We present diverse examples where milestoning systematically renders observations "closer to underlying microscopic dynamics" and thereby improves thermodynamic inference from lumped data assuming a given range of memory, and we hypothesize that this effect is quite general. Moreover, whereas the correct general physical definition of time reversal in the presence of memory remains unknown, we here show by means of physically relevant examples that at least for semi-Markov processes of first and second order, waiting-time contributions arising from adopting a naive Markovian definition of time reversal generally must be discarded.

Trade-offs between number fluctuations and response in nonequilibrium chemical reaction networks

We study the response of chemical reaction networks driven far from equilibrium to logarithmic perturbations of reaction rates. The response of the mean number of a chemical species is observed to be quantitively limited by number fluctuations and the maximum thermodynamic driving force. We prove these trade-offs for linear chemical reaction networks and a class of nonlinear chemical reaction networks with a single chemical species. Numerical results for several model systems support the conclusion that these trade-offs continue to hold for a broad class of chemical reaction networks, though their precise form appears to sensitively depend on the deficiency of the network.

Long-term behaviours of Autocatalytic Sets

Autocatalytic Sets are reaction networks theorised as networks at the basis of life. Their main feature is the ability of spontaneously emerging and self-reproducing. The Reflexively and Food-generated theory provides a formal definition of Autocatalytic Sets in terms of graphs with peculiar topological properties. This formalisation has been proved to be a powerful tool for the study of the chemical networks underlying life, and it was able to identify autocatalytic structures in real metabolic networks. However, the dynamical behaviour of such networks has not been yet complitely clarified. In this work, I present a first attempt to connect the topology of an Autocatalytic Set with its dynamics. For this purpose, I represent Autocatalytic Sets in terms of Chemical Reaction Networks, and I use the Chemical Reaction Network theory to detect motifs in the networks'structure, that allow to determine the long-term behaviour of the system.

Stochastic approach to entropy production in chemical chaos

Methods are presented to evaluate the entropy production rate in stochastic reactive systems. These methods are shown to be consistent with known results from nonequilibrium chemical thermodynamics. Moreover, it is proved that the time average of the entropy production rate can be decomposed into the contributions of the cycles obtained from the stoichiometric matrix in both stochastic processes and deterministic systems. These methods are applied to a complex reaction network constructed on the basis of Rössler's reinjection principle and featuring chemical chaos.

Deciphering noise amplification and reduction in open chemical reaction networks

The impact of fluctuations on the dynamical behaviour of complex biological systems is a longstanding issue, whose understanding would elucidate how evolutionary pressure tends to modulate intrinsic noise. Using the Itō stochastic differential equation formalism, we performed analytic and numerical analyses of model systems containing different molecular species in contact with the environment and interacting with each other through mass-action kinetics. For networks of zero deficiency, which admit a detailed- or complex-balanced steady state, all molecular species are uncorrelated and their Fano factors are Poissonian. Systems of higher deficiency have non-equilibrium steady states and non-zero reaction fluxes flowing between the complexes. When they model homo-oligomerization, the noise on each species is reduced when the flux flows from the oligomers of lowest to highest degree, and amplified otherwise. In the case of hetero-oligomerization systems, only the noise on the highest-degree species shows this behaviour.

Insights into noise modulation in oligomerization systems of increasing complexity

Understanding under which conditions the increase of systems complexity is evolutionarily advantageous, and how this trend is related to the modulation of the intrinsic noise, are fascinating issues of utmost importance for synthetic and systems biology. To get insights into these matters, we analyzed a series of chemical reaction networks with different topologies and complexity, described by mass-action kinetics. We showed, analytically and numerically, that the global level of fluctuations at the steady state, measured by the sum over all species of the Fano factors of the number of molecules, is directly related to the network's deficiency. For zero-deficiency systems, this sum is constant and equal to the rank of the network. For higher deficiencies, additional terms appear in the Fano factor sum, which are proportional to the net reaction fluxes between the molecular complexes. We showed that the system's global intrinsic noise is reduced when all fluxes flow from lower to higher degree oligomers, or equivalently, towards the species of higher complexity, whereas it is amplified when the fluxes are directed towards lower complexity species.

Effective Thermodynamics for a Marginal Observer

Thermodynamics is usually formulated on the presumption that the observer has complete information about the system he or she deals with: no parasitic current, exact evaluation of the forces that drive the system. For example, the acclaimed fluctuation relation (FR), relating the probability of time-forward and time-reversed trajectories, assumes that the measurable transitions suffice to characterize the process as Markovian (in our case, a continuous-time jump process). However, most often the observer only measures a marginal current. We show that he or she will nonetheless produce an effective description that does not dispense with the fundamentals of thermodynamics, including the FR and the 2nd law. Our results stand on the mathematical construction of a hidden time reversal of the dynamics, and on the physical requirement that the observed current only accounts for a single transition in the configuration space of the system. We employ a simple abstract example to illustrate our results and to discuss the feasibility of generalizations.

Maximum Entropy Prediction of Non-Equilibrium Stationary Distributions for Stochastic Reaction Networks with Oscillatory Dynamics

Many chemical reaction networks in biological systems present complex oscillatory dynamics. In systems such as regulatory gene networks, cell cycle, and enzymatic processes, the number of molecules involved is often far from the thermodynamic limit. Although stochastic models based on the probabilistic approach of the Chemical Master Equation (CME) have been proposed, studies in the literature have been limited by the challenges of solving the CME and the lack of computational power to perform large-scale stochastic simulations. In this paper, we show that the infinite set of stationary moment equations describing the stochastic Brusselator and Schnakenberg oscillatory reactions networks can be truncated and solved using maximization of the entropy of the distributions. The results from our numerical experiments compare with the distributions obtained from well-established kinetic Monte Carlo methods and suggest that the accuracy of the prediction increases exponentially with the closure order chosen for the system. We conclude that maximum entropy models can be used as an efficient closure scheme alternative for moment equations to predict the non-equilibrium stationary distributions of stochastic chemical reactions with oscillatory dynamics. This prediction is accomplished without any prior knowledge of the system dynamics and without imposing any biased assumptions on the mathematical relations among species involved.

Solving moment hierarchies for chemical reaction networks

The study of chemical reaction networks (CRN's) is a very active field. Earlier well-known results (Feinberg 1987 Chem. Enc. Sci. 42 2229, Anderson et al 2010 Bull. Math. Biol. 72 1947) identify a topological quantity called deficiency, for any CRN, which, when exactly equal to zero, leads to a unique factorized steady-state for these networks. No results exist however for the steady states of non-zero-deficiency networks. In this paper, we show how to write the full moment-hierarchy for any non-zero-deficiency CRN obeying mass-action kinetics, in terms of equations for the factorial moments. Using these, we can recursively predict values for lower moments from higher moments, reversing the procedure usually used to solve moment hierarchies. We show, for nontrivial examples, that in this manner we can predict any moment of interest, for CRN's with non-zero deficiency and non-factorizable steady states.

Optimal inference strategies and their implications for the linear noise approximation

We study the information loss of a class of inference strategies that is solely based on time averaging. For an array of independent binary sensors (e.g., receptors, single electron transistors) measuring a weak random signal (e.g., ligand concentration, gate voltage) this information loss is up to 0.5  bit per measurement irrespective of the number of sensors. We derive a condition related to the local detailed balance relation that determines whether or not such a loss of information occurs. Specifically, if the free-energy difference arising from the signal is symmetrically distributed among the forward and backward rates, time integration mechanisms will capture the full information about the signal. As an implication, for the linear noise approximation, we can identify the same loss of information, arising from its inherent simplification of the dynamics.

Conservation laws and symmetries in stochastic thermodynamics

Phenomenological nonequilibrium thermodynamics describes how fluxes of conserved quantities, such as matter, energy, and charge, flow from outer reservoirs across a system and how they irreversibly degrade from one form to another. Stochastic thermodynamics is formulated in terms of probability fluxes circulating in the system's configuration space. The consistency of the two frameworks is granted by the condition of local detailed balance, which specifies the amount of physical quantities exchanged with the reservoirs during single transitions between configurations. We demonstrate that the topology of the configuration space crucially determines the number of independent thermodynamic affinities (forces) that the reservoirs generate across the system and provides a general algorithm that produces the fundamental affinities and their conjugate currents contributing to the total dissipation, based on the interplay between macroscopic conservations laws for the currents and microscopic symmetries of the affinities.

Fluctuation-Dissipation Relations Far from Equilibrium

Near equilibrium, where all currents of a system vanish on average, the fluctuation-dissipation relation (FDR) connects a current's spontaneous fluctuations with its response to perturbations of the conjugate thermodynamic force. Out of equilibrium, fluctuation-response relations generally involve additional nondissipative contributions. Here, in the framework of stochastic thermodynamics, we show that an equilibriumlike FDR holds for internally equilibrated currents, if the perturbing conjugate force only affects the microscopic transitions that contribute to the current. We discuss the physical requirements for the validity of our result and apply it to nanosized electronic devices.