Probabilistic models of uORF-mediated ATF4 translation control.

ATF4 is a key transcription factor that activates transcription of genes needed to respond to cellular stress. Although the mRNA encoding ATF4 is present at constant levels in the cell during the initial response, translation of ATF4 increases under conditions of cellular stress while the global translation rate decreases. We study two models for the control system that regulates the translation of ATF4, both based on the Vattem-Wek hypothesis. This hypothesis is based on a race to reload, following the translation of a small upstream open reading frame (uORF), the ternary complex that brings the initiator tRNA to the ribosome as the 40S subunit scans along the mRNA, encountering first a start codon for an inhibitory uORF whose reading frame overlaps the start of the ATF4 coding sequence. We develop a pair of simple, analytic, probabilistic models, one of which assumes all nucleotide triplets have identical kinetic properties, while the other recognizes the existence of triplets at which the ternary complex loads more efficiently. We also consider two different functions representing the dependence of the rate of initiation at uORF1 on the ternary complex concentration. In keeping with the theme of this Special Issue, we studied the properties of these models in a Maple document, which can easily be modified to consider different parameters, translation rate initiation functions, and so on.

On the quasi-steady-state approximation in an open Michaelis-Menten reaction mechanism.

The conditions for the validity of the standard quasi-steady-state approximation in the Michaelis-Menten mechanism in a closed reaction vessel have been well studied, but much less so the conditions for the validity of this approximation for the system with substrate inflow. We analyze quasi-steady-state scenarios for the open system attributable to singular perturbations, as well as less restrictive conditions. For both settings we obtain distinguished invariant manifolds and time scale estimates, and we highlight the special role of singular perturbation parameters in higher order approximations of slow manifolds. We close the paper with a discussion of distinguished invariant manifolds in the global phase portrait.

Small binding-site clearance delays are not negligible in gene expression modeling.

During the templated biopolymerization processes of transcription and translation, a macromolecular machine, either an RNA polymerase or a ribosome, binds to a specific site on the template. Due to the sizes of these enzymes, there is a waiting time before one clears the binding site and another can bind. These clearance delays are relatively short, and one might think that they could be neglected. However, in the case of transcription, these clearance delays are associated with conservation laws, resulting in surprisingly large effects on the bifurcation diagrams in models of gene expression networks. We study an example of this phenomenon in a model of a gene regulated by a non-coding RNA displaying bistability. Neglecting the binding-site clearance delays in this model can only be compensated for by making ad hoc, unphysical adjustments to the model's kinetic constants.

Heineken, Tsuchiya and Aris on the mathematical status of the pseudo-steady state hypothesis: A classic from volume 1 of Mathematical Biosciences.

Volume 1, Issue 1 of Mathematical Biosciences was the venue for a now-classic paper on the application of singular perturbation theory in enzyme kinetics, "On the mathematical status of the pseudo-steady state hypothesis of biochemical kinetics" by F. G. Heineken, H. M. Tsuchiya and R. Aris. More than 50 years have passed, and yet this paper continues to be studied and mined for insights. This perspective discusses both the strengths and weaknesses of the work presented in this paper. For many, the justification of the pseudo-steady-state approximation using singular perturbation theory is the main achievement of this paper. However, there is so much more material here, which laid the foundation for a great deal of research in mathematical biochemistry in the intervening decades. The parameterization of the equations, construction of the first-order uniform singular-perturbation solution, and an attempt to apply similar principles to the pseudo-equilibrium approximation are discussed in particular detail.

Temporal metabolic partitioning of the yeast and protist cellular networks: the cell is a global scale-invariant (fractal or self-similar) multioscillator.

Britton Chance, electronics expert when a teenager, became an enthusiastic student of biological oscillations, passing on this enthusiasm to many students and colleagues, including one of us (DL). This historical essay traces BC's influence through the accumulated work of DL to DL's many collaborators. The overall temporal organization of mass-energy, information, and signaling networks in yeast in self-synchronized continuous cultures represents, until now, the most characterized example of in vivo elucidation of time structure. Continuous online monitoring of dissolved gases by direct measurement (membrane-inlet mass spectrometry, together with NAD(P)H and flavin fluorescence) gives strain-specific dynamic information from timescales of minutes to hours as does two-photon imaging. The predominantly oscillatory behavior of network components becomes evident, with spontaneously synchronized cellular respiration cycles between discrete periods of increased oxygen consumption (oxidative phase) and decreased oxygen consumption (reductive phase). This temperature-compensated ultradian clock provides coordination, linking temporally partitioned functions by direct feedback loops between the energetic and redox state of the cell and its growing ultrastructure. Multioscillatory outputs in dissolved gases with 13 h, 40 min, and 4 min periods gave statistical self-similarity in power spectral and relative dispersional analyses: i.e., complex nonlinear (chaotic) behavior and a functional scale-free (fractal) network operating simultaneously over several timescales.

A mathematical model of the biochemical network underlying left-right asymmetry establishment in mammals.

The expression of the TGF-ß protein Nodal on the left side of vertebrate embryos is a determining event in the development of internal-organ asymmetry. We present a mathematical model for the control of the expression of Nodal and its antagonist Lefty consisting entirely of realistic elementary reactions. We analyze the model in the absence of Lefty and find a wide range of parameters over which bistability (two stable steady states) is observed, with one stable steady state a low-Nodal state corresponding to the right-hand developmental fate, and the other a high-Nodal state corresponding to the left. We find that bistability requires a transcription factor containing two molecules of phosphorylated Smad2. A numerical survey of the full model, including Lefty, shows the effects of Lefty on the potential for bistability, and on the conditions that lead to the system reaching one or the other steady state.

Parametric Sensitivity Analysis of Oscillatory Delay Systems with an Application to Gene Regulation.

A parametric sensitivity analysis for periodic solutions of delay-differential equations is developed. Because phase shifts cause the sensitivity coefficients of a periodic orbit to diverge, we focus on sensitivities of the extrema, from which amplitude sensitivities are computed, and of the period. Delay-differential equations are often used to model gene expression networks. In these models, the parametric sensitivities of a particular genotype define the local geometry of the evolutionary landscape. Thus, sensitivities can be used to investigate directions of gradual evolutionary change. An oscillatory protein synthesis model whose properties are modulated by RNA interference is used as an example. This model consists of a set of coupled delay-differential equations involving three delays. Sensitivity analyses are carried out at several operating points. Comments on the evolutionary implications of the results are offered.

A biochemically semi-detailed model of auxin-mediated vein formation in plant leaves.

We present here a model intended to capture the biochemistry of vein formation in plant leaves. The model consists of three modules. Two of these modules, those describing auxin signaling and transport in plant cells, are biochemically detailed. We couple these modules to a simple model for PIN (auxin efflux carrier) protein localization based on an extracellular auxin sensor. We study the single-cell responses of this combined model in order to verify proper functioning of the modeled biochemical network. We then assemble a multicellular model from the single-cell building blocks. We find that the model can, under some conditions, generate files of polarized cells, but not true veins.

Turing-Hopf instability in biochemical reaction networks arising from pairs of subnetworks.

Network conditions for Turing instability in biochemical systems with two biochemical species are well known and involve autocatalysis or self-activation. On the other hand general network conditions for potential Turing instabilities in large biochemical reaction networks are not well developed. A biochemical reaction network with any number of species where only one species moves is represented by a simple digraph and is modeled by a reaction-diffusion system with non-mass action kinetics. A graph-theoretic condition for potential Turing-Hopf instability that arises when a spatially homogeneous equilibrium loses its stability via a single pair of complex eigenvalues is obtained. This novel graph-theoretic condition is closely related to the negative cycle condition for oscillations in ordinary differential equation models and its generalizations, and requires the existence of a pair of subnetworks, each containing an even number of positive cycles. The technique is illustrated with a double-cycle Goodwin type model.

Feedforward non-Michaelis-Menten mechanism for CO(2) uptake by Rubisco: contribution of carbonic anhydrases and photorespiration to optimization of photosynthetic carbon assimilation.

Rubisco, the most abundant protein serving as the primary engine generating organic biomass on Earth, is characterized by a low catalytic constant (in higher plants approx. 3s(-1)) and low specificity for CO(2) leading to photorespiration. We analyze here why this enzyme evolved as the main carbon fixation engine. The high concentration of Rubisco exceeding the concentration of its substrate CO(2) by 2-3 orders of magnitude makes application of Michaelis-Menten kinetics invalid and requires alternative kinetic approaches to describe photosynthetic CO(2) assimilation. Efficient operation of Rubisco is supported by a strong flux of CO(2) to the chloroplast stroma provided by fast equilibration of bicarbonate and CO(2) and forwarding the latter to Rubisco reaction centers. The main part of this feedforward mechanism is a thylakoidal carbonic anhydrase associated with photosystem II and pumping CO(2) from the thylakoid lumen in coordination with the rate of electron transport, water splitting and proton gradient across the thylakoid membrane. This steady flux of CO(2) limits photosynthesis at saturating CO(2) concentrations. At low ambient CO(2) and correspondingly limited capacity of the bicarbonate pool in the stroma, its depletion at the sites of Rubisco is relieved by utilizing O(2) instead of CO(2), i.e. by photorespiration, a process which supplies CO(2) back to Rubisco and buffers the redox state and energy level in the chloroplast. Thus, the regulation of Rubisco function aims to keep steady non-equilibrium levels of CO(2), NADPH/NADP and ATP/ADP in the chloroplast stroma and to optimize the condition of homeostatic photosynthetic flux of matter and energy.

Dynamics and mechanisms of oscillatory photosynthesis.

We classify mathematical models that can be used to describe photosynthetic oscillations using ideas from nonlinear dynamics, and discuss potential mechanisms for photosynthetic oscillations in the context of this classification. We then turn our attention to recent experiments with leaves transferred to a low CO2 atmosphere which revealed stochastic oscillations with a period of a few seconds. Rubisco is the enzyme that takes both CO2 and O2 as substrates correspondingly for photosynthetic assimilation and for photorespiration. Photosynthesis depletes CO2 and produces O2 while respiration and photorespiration work in the opposite direction, so the product of one process becomes the reactant of the other coupled process. We examine the possibility of oscillations of CO2 and O2 in the leaf in relation to photorespiration. We suggest that in the cell, oscillations with a period of a few seconds, corresponding to the time between photosynthetic CO2 fixation and photorespiratory CO2 release, underlie the dynamics of metabolism in C3 plants.

The scale-free dynamics of eukaryotic cells.

Temporal organization of biological processes requires massively parallel processing on a synchronized time-base. We analyzed time-series data obtained from the bioenergetic oscillatory outputs of Saccharomyces cerevisiae and isolated cardiomyocytes utilizing Relative Dispersional (RDA) and Power Spectral (PSA) analyses. These analyses revealed broad frequency distributions and evidence for long-term memory in the observed dynamics. Moreover RDA and PSA showed that the bioenergetic dynamics in both systems show fractal scaling over at least 3 orders of magnitude, and that this scaling obeys an inverse power law. Therefore we conclude that in S. cerevisiae and cardiomyocytes the dynamics are scale-free in vivo. Applying RDA and PSA to data generated from an in silico model of mitochondrial function indicated that in yeast and cardiomyocytes the underlying mechanisms regulating the scale-free behavior are similar. We validated this finding in vivo using single cells, and attenuating the activity of the mitochondrial inner membrane anion channel with 4-chlorodiazepam to show that the oscillation of NAD(P)H and reactive oxygen species (ROS) can be abated in these two evolutionarily distant species. Taken together these data strongly support our hypothesis that the generation of ROS, coupled to redox cycling, driven by cytoplasmic and mitochondrial processes, are at the core of the observed rhythmicity and scale-free dynamics. We argue that the operation of scale-free bioenergetic dynamics plays a fundamental role to integrate cellular function, while providing a framework for robust, yet flexible, responses to the environment.

Graph-theoretic methods for the analysis of chemical and biochemical networks. II. Oscillations in networks with delays.

Delay-differential equations are commonly used to model genetic regulatory systems with the delays representing transcription and translation times. Equations with delayed terms can also be used to represent other types of chemical processes. Here we analyze delayed mass-action systems, i.e. systems in which the rates of reaction are given by mass-action kinetics, but where the appearance of products may be delayed. Necessary conditions for delay-induced instability are presented in terms both of a directed graph (digraph) constructed from the Jacobian matrix of the corresponding ODE model and of a species-reaction bipartite graph which directly represents a chemical mechanism. Methods based on the bipartite graph are particularly convenient and powerful. The condition for a delay-induced instability in this case is the existence of a subgraph of the bipartite graph containing an odd number of cycles of which an odd number are negative.

Graph-theoretic methods for the analysis of chemical and biochemical networks. I. Multistability and oscillations in ordinary differential equation models.

A chemical mechanism is a model of a chemical reaction network consisting of a set of elementary reactions that express how molecules react with each other. In classical mass-action kinetics, a mechanism implies a set of ordinary differential equations (ODEs) which govern the time evolution of the concentrations. In this article, ODE models of chemical kinetics that have the potential for multiple positive equilibria or oscillations are studied. We begin by considering some methods of stability analysis based on the digraph of the Jacobian matrix. We then prove two theorems originally given by A. N. Ivanova which correlate the bifurcation structure of a mass-action model to the properties of a bipartite graph with nodes representing chemical species and reactions. We provide several examples of the application of these theorems.

Observation of a chaotic multioscillatory metabolic attractor by real-time monitoring of a yeast continuous culture.

We monitored a continuous culture of the yeast Saccharomyces cerevisiae by membrane-inlet mass spectrometry. This technique allows very rapid simultaneous measurements (one point every 12 s) of several dissolved gases. During our experiment, the culture exhibited a multioscillatory mode in which the dissolved oxygen and carbon dioxide records displayed periodicities of 13 h, 36 min and 4 min. The 36- and 4-min modes were not visible at all times, but returned at regular intervals during the 13-h cycle. The 4-min mode, which has not previously been described in continuous culture, can also be seen when the culture displays simpler oscillatory behavior. The data can be used to visualize a metabolic attractor of this system, i.e. the set of dissolved gas concentrations which are consistent with the multioscillatory state. Computation of the leading Lyapunov exponent reveals the dynamics on this attractor to be chaotic.

Oscillations of the internal CO(2) concentration in tobacco leaves transferred to low CO(2).

Measurement of the internal CO(2) concentration (Ci) in tobacco leaves using a fast-response CO(2) exchange system showed that in the light, switching from 350 microLL(-1) to a low CO(2) concentration of 36.5 microLL(-1) (promoting high photorespiration) resulted in the Ci oscillating near the value of CO(2) compensation point (Gamma*). The oscillations are highly irregular, the range of Ci varying by 2-4 microLL(-1) in substomatal cavities with a period of a few seconds. The statistical properties of the time series became stationary after a transient of approximately 100s following transfer to low CO(2). Attractor reconstruction shows that the observed oscillations are not chaotic but exhibit stochastic behavior. The period of oscillations is consistent with the duration of photorespiratory post-illumination burst (PIB). We suggest that the observed oscillations may be due to a similar mechanism to that which leads to PIB, and may play a role in switching mitochondrial operation between oxidation of the photorespiratory glycine and of the tricarboxylic acid cycle substrates.

The functional equation truncation method for approximating slow invariant manifolds: a rapid method for computing intrinsic low-dimensional manifolds.

A slow manifold is a low-dimensional invariant manifold to which trajectories nearby are rapidly attracted on the way to the equilibrium point. The exact computation of the slow manifold simplifies the model without sacrificing accuracy on the slow time scales of the system. The Maas-Pope intrinsic low-dimensional manifold (ILDM) [Combust. Flame 88, 239 (1992)] is frequently used as an approximation to the slow manifold. This approximation is based on a linearized analysis of the differential equations and thus neglects curvature. We present here an efficient way to calculate an approximation equivalent to the ILDM. Our method, called functional equation truncation (FET), first develops a hierarchy of functional equations involving higher derivatives which can then be truncated at second-derivative terms to explicitly neglect the curvature. We prove that the ILDM and FET-approximated (FETA) manifolds are identical for the one-dimensional slow manifold of any planar system. In higher-dimensional spaces, the ILDM and FETA manifolds agree to numerical accuracy almost everywhere. Solution of the FET equations is, however, expected to generally be faster than the ILDM method.

A graph-theoretic method for detecting potential Turing bifurcations.

The conditions for diffusion-driven (Turing) instabilities in systems with two reactive species are well known. General methods for detecting potential Turing bifurcations in larger reaction schemes are, on the other hand, not well developed. We prove a theorem for a graph-theoretic condition originally given by Volpert and Ivanova [Mathematical Modeling (Nauka, Moscow, 1987) (in Russian), p. 57] for Turing instabilities in a mass-action reaction-diffusion system involving n substances. The method is based on the representation of a reaction mechanism as a bipartite graph with two types of nodes representing chemical species and reactions, respectively. The condition for diffusion-driven instability is related to the existence of a structure in the graph known as a critical fragment. The technique is illustrated using a substrate-inhibited bifunctional enzyme mechanism which involves seven chemical species.

Stochastic kinetics description of a simple transcription model.

We study a stochastic model of transcription kinetics in order to characterize the distributions of transcriptional delay and of elongation rates. Transcriptional delay is the time which elapses between the binding of RNA polymerase to a promoter sequence and its dissociation from the DNA template strand with consequent release of the transcript. Transcription elongation is the process by which the RNA polymerase slides along the template strand. The model considers a DNA template strand with one promoter site and n nucleotide sites, and five types of reaction processes, which we think are key ones in transcription. The chemical master equation is a set of ordinary differential equations in 3(n) variables, where n is the number of bases in the template. This model is too huge to be handled if n is large. We manage to get a reduced Markov model which has only 2n independent variables and can well approximate the original dynamics. We obtain a number of analytical and numerical results for this model, including delay and transcript elongation rate distributions. Recent studies of single-RNA polymerase transcription by using optical-trapping techniques raise an issue of whether the elongation rates measured in a population are heterogeneous or not. Our model implies that in the cases studied, different RNA polymerase molecules move at different characteristic rates along the template strand. We also discuss the implications of this work for the mathematical modeling of genetic regulatory circuits.

The mechanism distinguishability problem in biochemical kinetics: the single-enzyme, single-substrate reaction as a case study.

A theoretical analysis of the distinguishability problem of two rival models of the single enzyme-single substrate reaction, the Michaelis-Menten and Henri mechanisms, is presented. We also outline a general approach for analysing the structural indistinguishability between two mechanisms. The approach involves constructing, if possible, a smooth mapping between the two candidate models. Evans et al. [N.D. Evans, M.J. Chappell, M.J. Chapman, K.R. Godfrey, Structural indistinguishability between uncontrolled (autonomous) nonlinear analytic systems, Automatica 40 (2004) 1947-1953] have shown that if, in addition, either of the mechanisms satisfies a particular criterion then such a transformation always exists when the models are indistinguishable from their experimentally observable outputs. The approach is applied to the single enzyme-single substrate reaction mechanism. In principle, mechanisms can be distinguished using this analysis, but we show that our ability to distinguish mechanistic models depends both on the precise measurements made, and on our knowledge of the system prior to performing the kinetics experiments.