Time delays in metabolic control systems

In this work we use mathematical models with discrete and distributed time delays to analyse the stability of metabolic pathways controlled by end product. We assume the kinetics of the intermediates of the path to be unknown, and we cover the lack of information by using a time delay. We find that above a definite substrate value, there is a critical delay Tc in which a transition from stability to instability occurs. For discrete delays, we find that even if the interaction of the end product with the first (allosteric) enzyme is not cooperative, the pathway can potentially become unstable and oscillate. We then show that the existence of cooperative inhibition extends the parametric domain of instability. The introduction of distributed delays shows, when the kernels are not monotonically decreasing, that the dispersion increases the critical delay Tc. Finally, we comment on the possibility that metabolic oscillations are physiological signals useful for triggering adaptive strategies in cell behavior.     

一类具分段常数变量的捕食-食饵系统的Neimark-Sacker分支

讨论了具时滞与分段常数变量的捕食-食饵生态模型的稳定性及Neimark-Sacker分支;通过计算得到连续模型对应的差分模型,基于特征值理论和Schur-Cohn判据得到正平衡态局部渐进稳定的充分条件;以食饵的内禀增长率为分支参数,运用分支理论和中心流形定理分析了Neimark-Sacker分支的存在性与稳定性条件;通过举例和数值模拟验证了理论的正确性。     

A juvenile-adult model with periodic vital rates

A global branch of positive cycles is shown to exist for a general discrete time, juvenile-adult model with periodically varying coefficients. The branch bifurcates from the extinction state at a critical value of the mean, inherent fertility rate. In comparison to the autonomous system with the same mean fertility rate, the critical bifurcation value can either increase or decrease with the introduction of periodicities. Thus, periodic oscillations in vital parameter can be either advantageous or deleterious. A determining factor is the phase relationship among the oscillations in the inherent fertility and survival rates.Research supported by NSF grant DMS-0414212.     

Competition in size-structured populations: mechanisms inducing cohort formation and population cycles

In this paper we investigate the consequences of size-dependent competition among the individuals of a consumer population by analyzing the dynamic properties of a physiologically structured population model. Only 2 size-classes of individuals are distinguished: juveniles and adults. Juveniles and adults both feed on one and the same resource and hence interact by means of exploitative competition. Juvenile individuals allocate all assimilated energy into development and mature on reaching a fixed developmental threshold. The combination of this fixed threshold and the resource-dependent developmental rate, implies that the juvenile delay between birth and the onset of reproduction may vary in time. Adult individuals allocate all assimilated energy to reproduction. Mortality of both juveniles and adults is assumed to be inversely proportional to the amount of energy assimilated. In this setting we study how the dynamics of the population are influenced by the relative foraging capabilities of juveniles and adults.In line with results that we previously obtained in size-structured consumer-resource models with pulsed reproduction, population cycles primarily occur when either juveniles or adults have a distinct competitive advantage. When adults have a larger per capita feeding rate and are hence competitively superior to juveniles, population oscillations occur that are primarily induced by the fact that the duration of the juvenile period changes with changing food conditions. These cycles do not occur when the juvenile delay is a fixed parameter. When juveniles are competitively superior, two different types of population fluctuations can occur: (1) rapid, low-amplitude fluctuations having a period of half the juvenile delay and (2) slow, large-amplitude fluctuations characterized by a period, which is roughly equal to the juvenile delay. The analysis of simplified versions of the structured model indicates that these two types of oscillations also occur if mortality and/or development is independent of food density, i.e. in a situation with a constant juvenile developmental delay and a constant, food-independent background mortality. Thus, the oscillations that occur when juveniles are more competitive are induced by the juvenile delay per se. When juveniles exert a larger foraging pressure on the shared resource, maturation implies an increase not only in adult density, but also in food density and consequently fecundity. Our analysis suggests that this correlation in time between adult density and fecundity is crucial for the occurrence of population cycles when juveniles are competitively superior.     

Global Hopf bifurcation analysis of an susceptible-infective-removed epidemic model incorporating media coverage with time delay

An susceptible-infective-removed epidemic model incorporating media coverage with time delay is proposed. The stability of the disease-free equilibrium and endemic equilibrium is studied. And then, the conditions which guarantee the existence of local Hopf bifurcation are given. Furthermore, we show that the local Hopf bifurcation implies the global Hopf bifurcation after the second critical value of delay. The obtained results show that the time delay in media coverage can not affect the stability of the disease-free equilibrium when the basic reproduction number is less than unity. However, the time delay affects the stability of the endemic equilibrium and produces limit cycle oscillations while the basic reproduction number is greater than unity. Finally, some examples for numerical simulations are included to support the theoretical prediction.     

Distinguishing intrinsic limit cycles from forced oscillations in ecological time series

Ecological cycles are ubiquitous in nature and have triggered ecologists’ interests for decades. Deciding whether a cyclic ecological variable, such as population density, is part of an intrinsically emerging limit cycle or simply driven by a varying environment is still an unresolved issue, particularly when the only available information is in the form of a recorded time series. We investigate the possibility of discerning intrinsic limit cycles from oscillations forced by a cyclic environment based on a single time series. We argue that such a distinction is possible because of the fundamentally different effects that perturbations have on the focal system in these two cases. Using a set of generic mathematical models, we show that random perturbations leave characteristic signatures on the power spectrum and autocovariance that differ between limit cycles and forced oscillations. We quantify these differences through two summary variables and demonstrate their predictive power using numerical simulations. Our work demonstrates that random perturbations of ecological cycles can give valuable insight into the underlying deterministic dynamics.     

A two variable delay model for the circadian rhythm of Neurospora crassa

A two variable model with delay in both the variables, is proposed for the circadian oscillations of protein concentrations in the fungal species Neurospora crassa. The dynamical variables chosen are the concentrations of FRQ and WC-1 proteins. Our model is a two variable simplification of the detailed model of Smolen et al. (J. Neurosci. 21 (2001) 6644) modeling circadian oscillations with interlocking positive and negative feedback loops, containing 23 variables. In our model, as in the case of Smolen's model, a sustained limit cycle oscillation takes place in both FRQ and WC-1 protein in continuous darkness, and WC-1 is anti-phase to FRQ protein, as observed in experiments. The model accounts for various characteristic features of circadian rhythms such as entrainment to light dark cycles, phase response curves and robustness to parameter variation and molecular fluctuations. Simulations are carried out to study the effect of periodic forcing of circadian oscillations by light-dark cycles. The periodic forcing resulted in a rich bifurcation diagram that includes quasiperiodicity and chaotic oscillations, depending on the magnitude of the periodic changes in the light controlled parameter. When positive feedback is eliminated, our model reduces to the generic one dimensional delay model of Lema et al. (J. Theor. Biol. 204 (2000) 565), delay model of the circadian pace maker with FRQ protein as the dynamical variable which represses its own production. This one-dimensional model also exhibits all characteristic features of circadian oscillations and gives rise to circadian oscillations which are reasonably robust to parameter variations and molecular noise.     

Oscillatory dynamics in rock-paper-scissors games with mutations

We study the oscillatory dynamics in the generic three-species rock-paper-scissors games with mutations. In the mean-field limit, different behaviors are found: (a) for high mutation rate, there is a stable interior fixed point with coexistence of all species; (b) for low mutation rates, there is a region of the parameter space characterized by a limit cycle resulting from a Hopf bifurcation; (c) in the absence of mutations, there is a region where heteroclinic cycles yield oscillations of large amplitude (not robust against noise). After a discussion on the main properties of the mean-field dynamics, we investigate the stochastic version of the model within an individual-based formulation. Demographic fluctuations are therefore naturally accounted and their effects are studied using a diffusion theory complemented by numerical simulations. It is thus shown that persistent erratic oscillations (quasi-cycles) of large amplitude emerge from a noise-induced resonance phenomenon. We also analytically and numerically compute the average escape time necessary to reach a (quasi-)cycle on which the system oscillates at a given amplitude.     

Control Analysis for Autonomously Oscillating Biochemical Networks

It has hitherto not been possible to analyze the control of oscillatory dynamic cellular processes in other than qualitative ways. The control coefficients, used in metabolic control analyses of steady states, cannot be applied directly to dynamic systems. We here illustrate a way out of this limitation that uses Fourier transforms to convert the time domain into the stationary frequency domain, and then analyses the control of limit cycle oscillations. In addition to the already known summation theorems for frequency and amplitude, we reveal summation theorems that apply to the control of average value, waveform, and phase differences of the oscillations. The approach is made fully operational in an analysis of yeast glycolytic oscillations. It follows an experimental approach, sampling from the model output and using discrete Fourier transforms of this data set. It quantifies the control of various aspects of the oscillations by the external glucose concentration and by various internal molecular processes. We show that the control of various oscillatory properties is distributed over the system enzymes in ways that differ among those properties. The models that are described in this paper can be accessed on http://jjj.biochem.sun.ac.za.     

A size dependent predator-prey interaction: who pursues whom?

We investigate the properties of an (age, size) -structured model for a population of Daphnia that feeds on a dynamical algal food source. The stability of the internal equilibrium is studied in detail and combined with numerical studies on the dynamics of the model to obtain insight in the relation between individual behaviour and population dynamical phenomena. Particularly the change in the (age, size)-relation with a change in the food availability seems to be an important behavioural mechanism that strongly influences the dynamics. This influence is partly stabilizing and partly destabilizing and leads to the coexistence of a stable equilibrium and a stable limit cycle or even coexistence of two stable limit cycles for the same parameter values. The oscillations in this case are characterized by drastic changes in the size-structure of the population during a cycle. In addition the model exhibits the usual predator-prey oscillations that characterize Lotka-Volterra models.     

Modeling the glucose-insulin regulatory system and ultradian insulin secretory oscillations with two explicit time delays

In the glucose-insulin regulatory system, ultradian insulin secretory oscillations are observed to have a period of 50-150 min. After pioneering work traced back to the 1960s, several mathematical models have been proposed during the last decade to model these ultradian oscillations as well as the metabolic system producing them. These currently existing models still lack some of the key physiological aspects of the glucose-insulin system. Applying the mass conservation law, we introduce two explicit time delays and propose a more robust alternative model for better understanding the glucose-insulin endocrine metabolic regulatory system and the ultradian insulin secretory oscillations for the cases of continuous enteral nutrition and constant glucose infusion. We compare the simulation profiles obtained from this two time delay model with those from the other existing models. As a result, we notice many unique features of this two delay model. Based on our intensive simulations, we suspect that one of the possibly many causes of ultradian insulin secretion oscillations is the time delay of the insulin secretion stimulated by the elevated glucose concentration.     

A model for a network of phosphorylation-dephosphorylation cycles displaying the dynamics of dominoes and clocks

We consider a model for a network of phosphorylation-dephosphorylation cycles coupled through forward and backward regulatory interactions, such that a protein phosphorylated in a given cycle activates the phosphorylation of a protein by a kinase in the next cycle as well as the dephosphorylation of a protein by a phosphatase in a preceding cycle. The network is cyclically organized in such a way that the protein phosphorylated in the last cycle activates the kinase in the first cycle. We study the dynamics of the network in the presence of both forward and backward coupling, in conditions where a threshold exists in each cycle in the amount of protein phosphorylated as a function of the ratio of kinase to phosphatase maximum rates. We show that this system can display sustained (limit-cycle) oscillations in which each cycle in the pathway is successively turned on and off, in a sequence resembling the fall of a series of dominoes. The model thus provides an example of a biochemical system displaying the dynamics of dominoes and clocks (Murray & Kirschner, 1989). It also shows that a continuum of clock waveforms exists of which the fall of dominoes represents a limit. When the cycles in the network are linked through only forward (positive) coupling, bistability is observed, while in the presence of only backward (negative) coupling, the system can display multistability or oscillations, depending on the number of cycles in the network. Inhibition or activation of any kinase or phosphatase in the network immediately stops the oscillations by bringing the system into a stable steady state; oscillations resume when the initial value of the kinase or phosphatase rate is restored. The progression of the system on the limit cycle can thus be temporarily halted as long as an inhibitor is present, much as when a domino is held in place. These results suggest that the eukaryotic cell cycle, governed by a network of phosphorylation-dephosphorylation reactions in which the negative control of cyclin-dependent kinases plays a prominent role, behaves as a limit-cycle oscillator impeded in the presence of inhibitors. We contrast the case where the sequence of domino-like transitions constitutes the clock with the case where the sequence of transitions is passively coupled to a biochemical oscillator operating as an independent clock.     

Sustained oscillations and threshold phenomena in an operon control circuit

A genetic regulatory model involving a positive feedback (via induction) and a negative feedback (via catabolite repression) is analyzed and applied to the problem of the lac operon regulation in E. coli. Damped and sustained oscillations of the limit cycle type are found along with threshold phenomena corresponding to multiple limit cycles or to multiple steady states, for values of the parameters compatible with experimental data. A comparison With the observations of Knorre and Goodwin is outlined.     

Oscillations in a model of repression with external control

A mathematical model for control by repression by an extracellular substance is developed, including diffusion and time delays. The model examines how active transport of a nutrient can produce either oscillatory or stable responses depending on a variety of parameters, such as diffusivity, cell size, or nutrient concentration. The system of equations for the mathematical model is reduced to a system of delay differential equations and linear Volterra equations. After linearizing these equations and forming the limiting Volterra equations, the resulting linear system no longer has any spatial dependence. Local stability analysis of the radially symmetric model shows that the system of equations can undergo Hopf bifurcations for certain parameter values, while other ranges of the parameters guarantee asymptotic stability. One numerical study shows that the model can exhibit intracellular biochemical oscillations with increasing extracellular concentrations of the nutrient, which suggests a possible trigger mechanism for morphogenesis.The work of this author was supported in part by NSF grants DMS-8603787 and DMS-8807360The work of this author was supported under the REU program of NSF by grant DMS-8807360The work of this author was supported under the REU program of NSF by grant DMS-8807360     

Stochastic models for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour

Circadian rhythms are endogenous oscillations that occur with a period close to 24 h in nearly all living organisms. These rhythms originate from the negative autoregulation of gene expression. Deterministic models based on such genetic regulatory processes account for the occurrence of circadian rhythms in constant environmental conditions (e.g., constant darkness), for entrainment of these rhythms by light-dark cycles, and for their phase-shifting by light pulses. When the numbers of protein and mRNA molecules involved in the oscillations are small, as may occur in cellular conditions, it becomes necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. We address the effect of molecular noise by considering the stochastic version of a deterministic model previously proposed for circadian oscillations of the PER and TIM proteins and their mRNAs in Drosophila. The model is based on repression of the per and tim genes by a complex between the PER and TIM proteins. Numerical simulations of the stochastic version of the model are performed by means of the Gillespie method. The predictions of the stochastic approach compare well with those of the deterministic model with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity from a bifurcation point beyond which the system evolves to stable steady state. Stochastic simulations indicate that robust circadian oscillations can emerge at the cellular level even when the maximum numbers of mRNA and protein molecules involved in the oscillations are of the order of only a few tens or hundreds. The stochastic model also reproduces the evolution to a strange attractor in conditions where the deterministic PER-TIM model admits chaotic behaviour. The difference between periodic oscillations of the limit cycle type and aperiodic oscillations (i.e. chaos) persists in the presence of molecular noise, as shown by means of Poincaré sections. The progressive obliteration of periodicity observed as the number of molecules decreases can thus be distinguished from the aperiodicity originating from chaotic dynamics. As long as the numbers of molecules involved in the oscillations remain sufficiently large (of the order of a few tens or hundreds, or more), stochastic models therefore provide good agreement with the predictions of the deterministic model for circadian rhythms.     

Dissipative Structures for an Allosteric Model: Application to Glycolytic Oscillations

An allosteric model of an open monosubstrate enzyme reaction is analyzed for the case where the enzyme, containing two protomers, is activated by the product. It is shown that this system can lead to instabilities beyond which a new state organized in time or in space (dissipative structure) can be reached. The conditions for both types of instabilities are presented and the occurrence of a temporal structure, consisting of a limit cycle behavior, is determined numerically as a function of the important parameters involved in the system. Sustained oscillations in the product and substrate concentrations are shown to occur for acceptable values of the allosteric and kinetic constants; moreover, they seem to be favored by substrate activation. The model is applied to phosphofructokinase, which is the enzyme chiefly responsible for glycolytic oscillations and which presents the same pattern of regulation as the allosteric enzyme appearing in the model. A qualitative and quantitative agreement is obtained with the experimental observations concerning glycolytic self-oscillations.     

Complex dynamics in biological systems arising from multiple limit cycle bifurcation

In this paper, we study complex dynamical behaviour in biological systems due to multiple limit cycles bifurcation. We use simple epidemic and predator–prey models to show exact routes to new types of bistability, that is, bistability between equilibrium and periodic oscillation, and bistability between two oscillations, which may more realistically describe the real situations. Bifurcation theory and normal form theory are applied to investigate the multiple limit cycles bifurcating from Hopf critical point.     

Coupling oscillations and switches in genetic networks

Switches (bistability) and oscillations (limit cycle) are omnipresent in biological networks. Synthetic genetic networks producing bistability and oscillations have been designed and constructed experimentally. However, in real biological systems, regulatory circuits are usually interconnected and the dynamics of those complex networks is often richer than the dynamics of simple modules. Here we couple the genetic Toggle switch and the Repressilator, two prototypic systems exhibiting bistability and oscillations, respectively. We study two types of coupling. In the first type, the bistable switch is under the control of the oscillator. Numerical simulation of this system allows us to determine the conditions under which a periodic switch between the two stable steady states of the Toggle switch occurs. In addition we show how birhythmicity characterized by the coexistence of two stable small-amplitude limit cycles, can easily be obtained in the system. In the second type of coupling, the oscillator is placed under the control of the Toggleswitch. Numerical simulation of this system shows that this construction could for example be exploited to generate a permanent transition from a stable steady state to self-sustained oscillations (and vice versa) after a transient external perturbation. Those results thus describe qualitative dynamical behaviors that can be generated through the coupling of two simple network modules. These results differ from the dynamical properties resulting from interlocked feedback loops systems in which a given variable is involved at the same time in both positive and negative feedbacks. Finally the models described here may be of interest in synthetic biology, as they give hints on how the coupling should be designed to get the required properties.     

A model for the dynamics of human weight cycling

The resolution to lose weight by cognitive restraint of nutritional intake often leads to repeated bouts of weight loss and regain, a phenomenon known as weight cycling or “yo-yo dieting”. A simple mathematical model for weight cycling is presented. The model is based on a feedback of psychological nature by which a subject decides to reduce dietary intake once a threshold weight is exceeded. The analysis of the model indicates that sustained oscillations in body weight occur in a parameter range bounded by critical values. Only outside this range can body weight reach a stable steady state. The model provides a theoretical framework that captures key facets of weight cycling and suggests ways to control the phenomenon. The view that weight cycling represents self-sustained oscillations has indeed specific implications. In dynamical terms, to bring weight cycling to an end, parameter values should change in such a way as to induce the transition of body weight from sustained oscillations around an unstable steady state to a stable steady state. Maintaining weight under a critical value should prevent weight cycling and allow body weight to stabilize below the oscillatory range.     

Effects of conduction delays on the existence and stability of one to one phase locking between two pulse-coupled oscillators

Gamma oscillations can synchronize with near zero phase lag over multiple cortical regions and between hemispheres, and between two distal sites in hippocampal slices. How synchronization can take place over long distances in a stable manner is considered an open question. The phase resetting curve (PRC) keeps track of how much an input advances or delays the next spike, depending upon where in the cycle it is received. We use PRCs under the assumption of pulsatile coupling to derive existence and stability criteria for 1:1 phase-locking that arises via bidirectional pulse coupling of two limit cycle oscillators with a conduction delay of any duration for any 1:1 firing pattern. The coupling can be strong as long as the effect of one input dissipates before the next input is received. We show the form that the generic synchronous and anti-phase solutions take in a system of two identical, identically pulse-coupled oscillators with identical delays. The stability criterion has a simple form that depends only on the slopes of the PRCs at the phases at which inputs are received and on the number of cycles required to complete the delayed feedback loop. The number of cycles required to complete the delayed feedback loop depends upon both the value of the delay and the firing pattern. We successfully tested the predictions of our methods on networks of model neurons. The criteria can easily be extended to include the effect of an input on the cycle after the one in which it is received.