Search

Maintaining stability in an uncertain environment is essential for proper functioning of living systems. Robust perfect adaptation (RPA) is a property of a system that generates an output at a fixed level even after fluctuations in input stimulus without fine-tuning parameters, and it is important to understand how this feature is implemented through biochemical networks. The existing literature has mainly focused on RPA of the concentration of a chosen chemical species, and no generic analysis has been made on RPA of reaction fluxes, that play an equally important role. Here, we identify structural conditions on reaction networks under which all the reaction fluxes exhibit RPA against the perturbation of the parameters inside a subnetwork. Based on this understanding, we give a recipe for obtaining a simpler reaction network, from which we can fully recover the steady-state reaction fluxes of the original system. This helps us identify key parameters that determine the fluxes and study the properties of complex reaction networks using a smaller one without losing any information about steady-state reaction fluxes., Comment: 5 pages, 2 figures (Supplemental Material: 14 pages, 1 figure)

Long-term behaviors of biochemical systems are described by steady states in deterministic models and stationary distributions in stochastic models. Obtaining their analytic solutions can be done for limited cases, such as linear or finite-state systems, as it generally requires solving many coupled equations. Interestingly, analytic solutions can be easily obtained when underlying networks have special topologies, called weak reversibility (WR) and zero deficiency (ZD), and the kinetic law follows a generalized form of mass-action kinetics. However, such desired topological conditions do not hold for the majority of cases. Thus, translating networks to have WR and ZD while preserving the original dynamics was proposed. Yet, this approach is limited because manually obtaining the desired network translation among the large number of candidates is challenging. Here, we prove necessary conditions for having WR and ZD after translation, and based on these conditions, we develop a user-friendly computational package, TOWARDZ, that automatically and efficiently identifies translated networks with WR and ZD. This allows us to quantitatively examine how likely it is to obtain WR and ZD after translation depending on the number of species and reactions. Importantly, we also describe how our package can be used to analytically derive steady states of deterministic models and stationary distributions of stochastic models. TOWARDZ provides an effective tool to analyze biochemical systems., Comment: 24 pages, 7 figures

Circadian clocks play a pivotal role in orchestrating numerous physiological and developmental events. Waveform shapes of the oscillations of protein abundances can be informative about the underlying biochemical processes of circadian clocks. We derive a mathematical framework where waveforms do reveal hidden biochemical mechanisms of circadian timekeeping. We find that the cost of synthesizing proteins with particular waveforms can be substantially reduced by rhythmic protein half-lives over time, as supported by previous plant and mammalian data, as well as our own seedling experiment. We also find that previously-enigmatic, cyclic expression of positive arm components within the mammalian and insect clocks allows both a broad range of peak time differences between protein waveforms and the symmetries of the waveforms about the peak times. Such various peak-time differences may facilitate tissue-specific or developmental stage-specific multicellular processes. Our waveform-guided approach can be extended to various biological oscillators, including cell-cycle and synthetic genetic oscillators., Comment: Supplementary material is available at the journal website

Biochemical reaction networks frequently consist of species evolving on multiple timescales. Stochastic simulations of such networks are often computationally challenging and therefore various methods have been developed to obtain sensible stochastic approximations on the timescale of interest. One of the rigorous and popular approaches is the multiscale approximation method for continuous time Markov processes. In this approach, by scaling species abundances and reaction rates, a family of processes parameterized by a scaling parameter is defined. The limiting process of this family is then used to approximate the original process. However, we find that such approximations become inaccurate when combinations of species with disparate abundances either constitute conservation laws or form virtual slow auxiliary species. To obtain more accurate approximation in such cases, we propose here an appropriate modification of the original method., Comment: 27 pages, 5 figures, This pre-print has been accepted for publication in SIAM Multiscale Modeling & Simulation. The final copyedited version of this paper will be available at https://www.siam.org/journals/mms.php

Circadian (~24hr) clocks are self-sustained endogenous oscillators with which organisms keep track of daily and seasonal time. Circadian clocks frequently rely on interlocked transcriptional- translational feedback loops to generate rhythms that are robust against intrinsic and extrinsic perturbations. To investigate the dynamics and mechanisms of the intracellular feedback loops in circadian clocks, a number of mathematical models have been developed. The majority of the models use Hill functions to describe transcriptional repression in a way that is similar to the Goodwin model. Recently, a new class of models with protein sequestration-based repression has been introduced. Here, we discuss how this new class of models differs dramatically from those based on Hill-type repression in several fundamental aspects: conditions for rhythm generation, robust network designs and the periods of coupled oscillators. Consistently, these fundamental properties of circadian clocks also differ among Neurospora, Drosophila, and mammals depending on their key transcriptional repression mechanisms (Hill-type repression or protein sequestration). Based on both theoretical and experimental studies, this review highlights the importance of careful modeling of transcriptional repression mechanisms in molecular circadian clocks., Comment: 31 pages, 5 figures