Promoters Architecture-Based Mechanism for Noise-Induced Oscillations in a Single-Gene Circuit

It is well known that single-gene circuits with negative feedback loop can lead to oscillatory gene expression when they operate with time delay. In order to generate these oscillations many processes can contribute to properly timing such delay. Here we show that the time delay coming from the transitions between internal states of the cis-regulatory system (CRS) can drive sustained oscillations in an auto-repressive single-gene circuit operating in a small volume like a cell. We found that the cooperative binding of repressor molecules is not mandatory for a oscillatory behavior if there are enough binding sites in the CRS. These oscillations depend on an adequate balance between the CRS kinetic, and the synthesis/degradation rates of repressor molecules. This finding suggest that the multi-site CRS architecture can play a key role for oscillatory behavior of gene expression. Finally, our results can also help to synthetic biologists on the design of the promoters architecture for new genetic oscillatory circuits.     

Complex metabolic oscillations in plants forced by harmonic irradiance

Plants exposed to harmonically modulated irradiance, approximately 1 + cos(omegat), exhibit a complex periodic pattern of chlorophyll fluorescence emission that can be deconvoluted into a steady-state component, a component that is modulated with the frequency of the irradiance (omega), and into at least two upper harmonic components (2omega and 3omega). A model is proposed that accounts for the upper harmonics in fluorescence emission by nonlinear negative feedback regulation of photosynthesis. In contrast to simpler linear models, the model predicts that the steady-state fluorescence component will depend on the frequency of light modulation, and that amplitudes of all fluorescence components will exhibit resonance peak(s) when the irradiance frequency is tuned to an internal frequency of a regulatory component. The experiments confirmed that the upper harmonic components appear and exhibit distinct resonant peaks. The frequency of autonomous oscillations observed earlier upon an abrupt increase in CO(2) concentration corresponds to the sharpest of the resonant peaks of the forced oscillations. We propose that the underlying principles are general for a wide spectrum of negative-feedback regulatory mechanisms. The analysis by forced harmonic oscillations will enable us to examine internal dynamics of regulatory processes that have not been accessible to noninvasive fluorescence monitoring to date.     

Chemical oscillations arise solely from kinetic nonlinearity and hence can occur near equilibrium.

A minimal kinetic scheme for a system displaying sustained chemical oscillations is presented. The system is isothermal, and all steps in the scheme are kinetically reversible. The oscillations are analyzed and the crucial points elucidated. Both positive and negative feedback, if properly introduced, support oscillations, provided the state responsible for feedback is optimally buffered. It is shown that the requisite nonlinearity is introduced at the kinetic level because of feedback regulation and not, as is usually assumed, by large affinities that introduce nonlinearity at the thermodynamic level. Hence, sustained oscillations may occur near equilibrium.     

Internal noise-driven circadian oscillator in Drosophila

An internal noise-driven oscillator was studied in a two-variable Drosophila model, where both positive feedback and negative feedback are crucial to the circadian oscillations. It is shown that internal noise could sustain reliable oscillations for the parameter which produces a stable steady state in the deterministic system. The noise-sustained oscillations are interpreted by using phase plane analysis. The period of such oscillations fluctuates slightly around the period of deterministic oscillations and the coherence of oscillations becomes the best at an optimal internal noise intensity, indicating the occurrence of intrinsic coherence resonance. In addition, in the oscillatory region, the coherence of noisy circadian oscillations is suppressed by the internal noise, but the period is hardly affected, demonstrating the robustness of the Drosophila model for circadian rhythms to the intrinsic noise.     

Emergence of complex behaviour from simple circuit structures

The set of (feedback) circuits of a complex system is the machinery that allows the system to be aware of the levels of its crucial constituents. Circuits can be identified without ambiguity from the elements of the Jacobian matrix of the system. There are two types of circuits: positive if they comprise an even number of negative interactions, negative if this number is odd. The two types of circuits play deeply different roles: negative circuits are required for homeostasis, with or without oscillations, positive circuits are required for multistationarity, and hence, in biology, for differentiation and memory. In non-linear systems, a circuit can positive or negative (an 'ambiguous circuit', depending on the location in phase space. Full circuits are those circuits (or unions of disjoint circuits) that imply all the variables of the system. There is a tight relation between circuits and steady states. Each full circuit, if isolated, generates steady state(s) whose nature (eigenvalues) is determined by the structure of the circuit. Multistationarity requires the presence of at least two full circuits of opposite Eisenfeld signs, or else, an ambiguous circuit. We show how a significant part of the dynamical behaviour of a system can be predicted by a mere examination of its Jacobian matrix. We also show how extremely complex dynamics can be generated by such simple logical structures as a single (full and ambiguous) circuit.     

Stable linear reaction oscillator

The linear chemical reaction system, described in a preceding paper, which showed undamped oscillations under certain conditions is reviewed in the light of the complete set of eigenvalues. Since one eigenvalue is found to be positive real, the system is extremely sensitive to the correct choice of initial conditions. A minute change in the “circuitry” of the system has the effect that all eigenvalues can be moved into the left complex half-plane (real part negative) including the imaginary axis. Thus, within a wide range, an arbitrarily chosen set of initial conditions will yield undamped sinusoidal oscillations after a short period of transient behaviour. Besides, the minimum number of steps in the feedback path is now three instead of five, making the original analogy to the electronic RC phase-shift oscillator even closer.     

Positive feedback circuits and memory

The concept of regulatory feedback circuit refers to oriented cyclic interactions between elements of a system. There are two classes of circuits, positive and negative, whose properties are in striking contrast. Positive circuits are a prerequisite for the occurrence of multiple steady states (multistationarity), and hence, they are involved in all processes showing hysteresis or memory. Endogenous or exogenous perturbations can lead the system to exhibit or to evoke one particular stable regime. The role of positive circuits in cell differentiation and in immunology is well documented. Negative circuits are involved in homeostatic regulation, with or without oscillations. The aim of this paper is to show: a) that positive circuits account for many features of memory stricto sensu (i.e., neural memory and mnesic evocation) as well as largo sensu (e.g. differentiation or immunological memory); and b) that simple combinations of positive and negative circuits provide powerful regulatory modules, which can also be associated in batteries. These entities have vast dynamical possibilities in the field of neurobiology, as well as in the fields of differentiation and immunology. Here we consider a universal minimal regulatory module, for which we suggest to adopt the term 'logical regulon', which can be considered as an atom of Jacob's integron. It comprises a positive and a negative circuit in its interaction matrix, and we recall the main results related to the simultaneous presence of these circuits. Finally, we give three applications of this type of interaction matrix. The first two deal with the coexistence of multiple stable steady states and periodicity in differentiation and in an immunological system showing hysteretic properties. The third deals with the dual problems of synchronization and desynchronization of a neural model for hippocampus memory evocation processes.     

The modularity of biological regulatory networks

A useful approach to complex regulatory networks consists of modeling their elements and interactions by Boolean equations. In this context, feedback circuits (i.e. circular sequences of interactions) have been shown to play key dynamical roles: whereas positive circuits are able to generate multistationarity, negative circuits may generate oscillatory behavior. In this paper, we principally focus on the case of gene networks. These are represented by fully connected Boolean networks where each element interacts with all elements including itself. Flexibility in network design is introduced by the use of Boolean parameters, one associated with each interaction or group of interactions affecting a given element. Within this formalism, a feedback circuit will generate its typical dynamical behavior (i.e. multistationarity or oscillations) only for appropriate values of some of the logical parameters. Whenever it does, we say that the circuit is 'functional'. More interestingly, this formalism allows the computation of the constraints on the logical parameters to have any feedback circuit functional in a network. Using this methodology, we found that the fraction of the total number of consistent combinations of parameter values that make a circuit functional decreases geometrically with the circuit length. From a biological point of view, this suggests that regulatory networks could be decomposed into small and relatively independent feedback circuits or 'regulatory modules'.     

Modelling receptor-controlled intracellular calcium oscillators.

This paper presents mathematical models for the hepatocyte calcium oscillator which follow the concepts in a class of informal models developed to account for the striking dependence on the receptor type of several features of the calcium oscillations, in particular the shape and duration of the free calcium transients. The essence of these models is that the transients should be timed by a build-up of activated GTP-binding proteins, which, combined with positive feedback processes and perhaps with cooperative effects, leads to a sudden activation of phospholipase C (PLC), followed by negative feedback processes which switch off the calcium rise and lead to a fall in free calcium back to resting levels. These models predict pulsatile oscillations in inositol (1,4,5)P3 as well as in free calcium. We show that receptor-controlled intracellular calcium oscillators involving an unknown positive feedback pathway onto PLC and negative feedback from protein kinase C (PKC) onto G-proteins and receptors, or negative feedback by stimulation of GTPase activity can simulate many of the features of observed intracellular calcium oscillations. These oscillators exhibit a dependence of frequency on agonist concentration and a dependence of transient duration on receptor and G-protein type. We also show that a PLC-dependent GTPase activating factor (GAF) could provide explanations for some otherwise puzzling features of intracellular calcium oscillations.     

Positive Feedback Promotes Oscillations in Negative Feedback Loops

A simple three-component negative feedback loop is a recurring motif in biochemical oscillators. This motif oscillates as it has the three necessary ingredients for oscillations: a three-step delay, negative feedback, and nonlinearity in the loop. However, to oscillate, this motif under the common Goodwin formulation requires a high degree of cooperativity (a measure of nonlinearity) in the feedback that is biologically “unlikely.” Moreover, this recurring negative feedback motif is commonly observed augmented by positive feedback interactions. Here we show that these positive feedback interactions promote oscillation at lower degrees of cooperativity, and we can thus unify several common kinetic mechanisms that facilitate oscillations, such as self-activation and Michaelis-Menten degradation. The positive feedback loops are most beneficial when acting on the shortest lived component, where they function by balancing the lifetimes of the different components. The benefits of multiple positive feedback interactions are cumulative for a majority of situations considered, when benefits are measured by the reduction in the cooperativity required to oscillate. These positive feedback motifs also allow oscillations with longer periods than that determined by the lifetimes of the components alone. We can therefore conjecture that these positive feedback loops have evolved to facilitate oscillations at lower, kinetically achievable, degrees of cooperativity. Finally, we discuss the implications of our conclusions on the mammalian molecular clock, a system modeled extensively based on the three-component negative feedback loop.     

Slow oscillations in blood pressure via a nonlinear feedback model

Blood pressure is well established to contain a potential oscillation between 0.1 and 0.4 Hz, which is proposed to reflect resonant feedback in the baroreflex loop. A linear feedback model, comprising delay and lag terms for the vasculature, and a linear proportional derivative controller have been proposed to account for the 0.4-Hz oscillation in blood pressure in rats. However, although this model can produce oscillations at the required frequency, some strict relationships between the controller and vasculature parameters must be true for the oscillations to be stable. We developed a nonlinear model, containing an amplitude-limiting nonlinearity that allows for similar oscillations under a very mild set of assumptions. Models constructed from arterial pressure and sympathetic nerve activity recordings obtained from conscious rabbits under resting conditions suggest that the nonlinearity in the feedback loop is not contained within the vasculature, but rather is confined to the central nervous system. The advantage of the model is that it provides for sustained stable oscillations under a wide variety of situations even where gain at various points along the feedback loop may be altered, a situation that is not possible with a linear feedback model. Our model shows how variations in some of the nonlinearity characteristics can account for growth or decay in the oscillations and situations where the oscillations can disappear altogether. Such variations are shown to accord well with observed experimental data. Additionally, using a nonlinear feedback model, it is straightforward to show that the variation in frequency of the oscillations in blood pressure in rats (0.4 Hz), rabbits (0.3 Hz), and humans (0.1 Hz) is primarily due to scaling effects of conduction times between species.     

Nonrandom variability in respiratory cycle parameters of humans during stage 2 sleep

We analyzed breath-to-breath inspiratory time (TI), expiratory time (TE), inspiratory volume (VI), and minute ventilation (Vm) from 11 normal subjects during stage 2 sleep. The analysis consisted of 1) fitting first- and second-order autoregressive models (AR1 and AR2) and 2) obtaining the power spectra of the data by fast-Fourier transform. For the AR2 model, the only coefficients that were statistically different from zero were the average alpha 1 (a1) for TI, VI, and Vm (a1 = 0.19, 0.29, and 0.15, respectively). However, the power spectra of all parameters often exhibited peaks at low frequency (less than 0.2 cycles/breath) and/or at high frequency (greater than 0.2 cycles/breath), indicative of periodic oscillations. After accounting for the corrupting effects of added oscillations on the a1 estimates, we conclude that 1) breath-to-breath fluctuations of VI, and to a lesser extent TI and Vm, exhibit a first-order autoregressive structure such that fluctuations of each breath are positively correlated with those of immediately preceding breaths and 2) the correlated components of variability in TE are mostly due to discrete high- and/or low-frequency oscillations with no underlying autoregressive structure. We propose that the autoregressive structure of VI, TI, and Vm during spontaneous breathing in stage 2 sleep may reflect either a central neural mechanism or the effects of noise in respiratory chemical feedback loops; the presence of low-frequency oscillations, seen more often in Vm, suggests possible instability in the chemical feedback loops. Mechanisms of high-frequency periodicities, seen more often in TE, are unknown.     

Oscillations in MAPK cascade triggered by two distinct designs of coupled positive and negative feedback loops

ABSTRACT: BACKGROUND: Feedback loops, both positive and negative are embedded in the Mitogen Activated Protein Kinase (MAPK) cascade. In the three layer MAPK cascade, both feedback loops originate from the terminal layer and their sites of action are either of the two upstream layers. Recent studies have shown that the cascade uses coupled positive and negative feedback loops in generating oscillations. Two plausible designs of coupled positive and negative feedback loops can be elucidated from the literature; in one design the positive feedback precedes the negative feedback in the direction of signal flow and vice-versa in another. But it remains unexplored how the two designs contribute towards triggering oscillations in MAPK cascade. Thus it is also not known how amplitude, frequency, robustness or nature (analogous/digital) of the oscillations would be shaped by these two designs. RESULTS: We built two models of MAPK cascade that exhibited oscillations as function of two underlying designs of coupled positive and negative feedback loops. Frequency, amplitude and nature (digital/analogous) of oscillations were found to be differentially determined by each design. It was observed that the positive feedback emerging from an oscillating MAPK cascade and functional in an external signal processing module can trigger oscillations in the target module, provided that the target module satisfy certain parametric requirements. The augmentation of the two models was done to incorporate the nuclear-cytoplasmic shuttling of cascade components followed by induction of a nuclear phosphatase. It revealed that the fate of oscillations in the MAPK cascade is governed by the feedback designs. Oscillations were unaffected due to nuclear compartmentalization owing to one design but were completely abolished in the other case. CONCLUSION: The MAPK cascade can utilize two distinct designs of coupled positive and negative feedback loops to trigger oscillations. The amplitude, frequency and robustness of the oscillations in presence or absence of nuclear compartmentalization were differentially determined by two designs of coupled positive and negative feedback loops. A positive feedback from an oscillating MAPK cascade was shown to induce oscillations in an external signal processing module, uncovering a novel regulatory aspect of MAPK signal processing.     

Modelling the effects of a negative feedback circuit from horizontal cells to cones on the impulse responses of cones and horizontal cells in the catfish retina

A model of the cone-L-HC circuit for the catfish retina is presented with the following features: the outer segment consists of a compression factor and 7 low-pass filters in tandem; the cone pedicle consists of an internal negative feedback circuit in series with a low-pass filter; and the L-HC consists of a low-pass filter and forms a negative feedback circuit with the cone pedicle. By proper adjustment of the various time constants of the low-pass filters and the gain factors, the impulse responses for cones and L-HCs of the catfish retina (and turtle) can be duplicated. The negative feedback gain increases with increasing levels of mean illuminance which causes the monophasic impulse responses to become faster, biphasic and decrease in amplitude, i.e. in gain. This is an expression of the Weber-Fechner law.     

Winner-take-all selection in a neural system with delayed feedback

We consider the effects of temporal delay in a neural feedback system with excitation and inhibition. The topology of our model system reflects the anatomy of the avian isthmic circuitry, a feedback structure found in all classes of vertebrates. We show that the system is capable of performing a 'winner-take-all' selection rule for certain combinations of excitatory and inhibitory feedback. In particular, we show that when the time delays are sufficiently large a system with local inhibition and global excitation can function as a 'winner-take-all' network and exhibit oscillatory dynamics. We demonstrate how the origin of the oscillations can be attributed to the finite delays through a linear stability analysis.     

Roles of protein ubiquitination and degradation kinetics in biological oscillations

Protein ubiquitination and degradation play important roles in many biological functions and are associated with many human diseases. It is well known that for biochemical oscillations to occur, proper degradation rates of the participating proteins are needed. In most mathematical models of biochemical reactions, linear degradation kinetics has been used. However, the degradation kinetics in real systems may be nonlinear, and how nonlinear degradation kinetics affects biological oscillations are not well understood. In this study, we first develop a biochemical reaction model of protein ubiquitination and degradation and calculate the degradation rate against the concentration of the free substrate. We show that the protein degradation kinetics mainly follows the Michaelis-Menten formulation with a time delay caused by ubiquitination and deubiquitination. We then study analytically how the Michaelis-Menten degradation kinetics affects the instabilities that lead to oscillations using three generic oscillation models: 1) a positive feedback mediated oscillator; 2) a positive-plus-negative feedback mediated oscillator; and 3) a negative feedback mediated oscillator. In all three cases, nonlinear degradation kinetics promotes oscillations, especially for the negative feedback mediated oscillator, resulting in much larger oscillation amplitudes and slower frequencies than those observed with linear kinetics. However, the time delay due to protein ubiquitination and deubiquitination generally suppresses oscillations, reducing the amplitude and increasing the frequency of the oscillations. These theoretical analyses provide mechanistic insights into the effects of specific proteins in the ubiquitination-proteasome system on biological oscillations.