key: cord-0792641-11esz6lk
authors: Lepzelter, David; Feng, Haidong; Wang, Jin
title: Oscillation, cooperativity, and intermediates in the self-repressing gene
date: 2010-04-26
journal: Chem Phys Lett
DOI: 10.1016/j.cplett.2010.03.029
sha: 005f37ed7fb9d0a45b613c2b3bc1bdbafc4f26ed
doc_id: 792641
cord_uid: 11esz6lk

Biological oscillators are vital to living organisms, which use them as clocks for time-sensitive processes. However, much is unknown about mechanisms which can give rise to coherent oscillatory behavior, with few exceptions (e.g., explicitly delayed self-repressors and simple models of specific organisms’ circadian clocks). We present what may be the simplest possible reliable gene network oscillator, a self-repressing gene. We show that binding cooperativity, which has not been considered in detail in this context, can combine with small numbers of intermediate steps to create coherent oscillation. We also note that noise blurs the line between oscillatory and non-oscillatory behavior.

Biological oscillation and its mechanisms have recently become a subject of intense study by a number of experimental and theoretical groups. Its study is still in the early stages, and while the phenomenon itself is believed to be of great importance, knowledge of mechanisms by which it can occur is quite incomplete.

Oscillation is used in a number of biological systems, most notably in the circadian rhythms responsible for keeping organisms' biochemical processes in line with the day-night schedule of the Earth. A number of mechanisms can be used by organisms to give a reasonably coherent oscillation, including three-gene repressilator systems [1] [2] [3] , self-repressors with maturation times or other explicit, deterministic time delays (as opposed to intermediate steps which would result from a full treatment of chemical intermediates) [4] [5] [6] [7] , more complex negative feedback loops tied to population dynamics [8] combined repression and activation loops [9] , highly non-linear protein degradation [10] , very large numbers (hundreds) of intermediate steps [11] , and self-repressors whose production and gene repression involve diffusion through the nuclear membrane [12, 13] . A self-activator with one stable and one metastable state, and stochastic switching between the two, can cause what may be called incoherent oscillation [14] . Additionally, certain kinds of self-repression can cause behavior which is not coherent oscillation when one considers the deterministic average protein and mRNA concentrations, but which still appear quite oscillatory and reasonably coherent when one considers any given system's stochastic trajectory through protein and mRNA concentration space over time [15] .

This last example is of particular interest because in part of its simplicity; it is an uncomplicated system, easily modelled using straightforward Markov chain Monte Carlo methods. 1 There is, however, another reason for interest in the work. The binding mechanism suggested by the mathematics presented by McKane and Newman in Ref. [15] is in fact a combination of multiple n-mer bindings of regulatory proteins to the gene. In addition, the which n contributes most to the genetic activity here depends on the single parameter which McKane and Newman used to tweak their system and the average number of proteins present in the system. While this in itself may be non-physical, it inherently suggests that the number of bound proteins may be able to change the system from non-oscillatory to oscillatory with noise to (possibly) oscillatory even without noise. We note that the simplest possible version of this last suggestion, a deterministic model with mRNA and with n regulatory proteins binding directly and instantaneously to the gene, cannot be oscillatory. A brief proof of this is included in the section on linear stability analysis (however, the stochastic version of the system can have oscillation-like behavior, and we explore this briefly).

1 Deterministic time delays of the type t À s; the most common method used to get a simple self-repressor to oscillate, can reproduce most results from other methods. However, for stochastics they require non-Markov simulation. Furthermore, while the equations are straightforward, the authors know of no way to determine a priori the numerical values of time delays that substitute for the delaying steps. The methods outlined in this Letter, in contrast, are purely Markov and involve only the actual chemical rates of the represented processes.

Chemical Physics Letters 490 (2010) 216-220

Chemical Physics Letters j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / c p l e t t mRNA lifetimes (also monomer protein lifetimes). The following were constant for all calculations: the degradation rate of mRNA k m ¼ 1 (hence, s m ¼ 1); the degradation rate of monomer protein k p ¼ 1 (hence, s p P 1); the rates of gene binding and unbinding x ) 1 (and so the system is very adiabatic); the ratio of the nmer dissolution constant f to the formation constant h is

; the rate of n-mer dissolution f ) 1 in the simplest case, Fig. 1 , and f ¼ 1 in all other calculations; the number of complete n-mers which will cause the gene to be repressed by a factor of 1=e, X eq ¼ 10; the protein synthesis rate from a single strand of mRNA g p ¼ 3; and the ratio of mRNA synthesis in the repressed versus the unrepressed gene state g 0 ¼ g 1 =100. It should be noted that g p and X eq can be rescaled in systems without internal noise; only when the actual number becomes important, as opposed to relative concentrations, does do these quantities have any significant effect other than a rescaling. It should also be noted that the ratio s p =s m is in rough agreement with the average for this sort of gene in yeast, given as approximately 2 in Ref. [16] . Deterministic calculations (involving the average numbers of mRNA, proteins, and protein aggregates) used the simple set of equations

where m is the number of mRNA molecules, P off ðP on Þ is the probability or fraction of the time that the gene is repressed (unrepressed), p is the number of monomer proteins, and c is the number of n-mer proteins in the system at time t. Simple modifications made additional intermediate steps possible (c becoming c 1 , n becoming n 1 , additional terms in the last equation for the n 2 c 1 ! c 2 step, the extra equations for c 2 in terms of c 1 and c 3 , etc.). We should also note that we assume all binding steps involved are fast compared to the other processes involved in the system, though binding events are relatively rare, and binding is highly cooperative. Slower binding, or binding that is less cooperative, can introduce very different elements to a system, as shown in Ref. [17] .

The authors used time-series data to determine the oscillatory nature of the systems. A general solution to the linear stability analysis of the system would involve analytically solving an n þ 1-degree polynomial with non-integer parameter-dependent coefficients, which is infeasible [18] . Even numerical solutions to the linear stability analysis problem alone could have glossed over complex behavior less useful to an organism than regular oscillation (i.e., chaos). However, individual points in different regions identified by the time-series data have been checked numerically using linear stability analysis. In the regime noted as 'oscillatory' via the time-series data, the Jacobian at the equilibrium point has one real negative eigenvalue and two complex eigenvalues with positive real parts, confirming that the point can be oscillatory. Outside that region, all three eigenvalues' real parts are negative, confirming the non-oscillatory classification. We note that the line between the zones labeled 'decaying oscillation' and 'no oscillation' is much less well-defined than the line between those labeled 'oscillation' and 'decaying oscillation', as a system which displays decaying oscillation given one set of initial conditions can have no oscillations with another set of initial conditions. However, we believe the general trend is correct and useful, and so we indicate the distinction in the figure.

Stochastic calculations were straightforward, using the following equations:

R monomer protein degradation ¼ k p p:

In the calculations for Fig. 1 (in which binding to the gene is coupled with n-mer formation), we used

derived from the standard master equation rates of binding hpðp À 1Þ Á Á Á ðp À n þ 1Þ and unbinding f. For other stochastic calculations, in which the proteins bind to each other before binding to the gene, we used instead R nÀmer formation ¼ hpðp À 1Þ Á Á Á ðp À n þ 1Þ;

R nÀmer breaking ¼ fc:

A time-step dt was then calculated using these rates, ensuring that events generally occurred one at a time by using

(this is subtly different from the traditional Gillespie simulation in that multiple events were in theory allowed, if quite unlikely; however, the difference is very small and may reasonably be considered to be an advantage of our algorithm.) For practicality's sake, at very high levels of synthesis, multiple synthesis and decay events were allowed and occurred in the background of other events, but were kept to at most a 0.01% mean change in the number of molecules.

Each process was then treated as a Poisson process using the same dt, with the probability of a events of type b occurring being e ÀR b dt ðR b dtÞ a a! . Fig. 1. (a) Stochastic calculation of coherence in a system which has mRNA and proteins cooperatively binding to the gene, with coherence given by 2 P Hðd/Þ P jd/j À 1 where d/ is the difference in angle in mRNA-protein space. (b) Stochastic calculation of the standard deviation of the period distribution divided by its mean. Both colormaps are on the same scale as Fig. 3 .

The programs used were written in C, and run using FEDORA 10 LINUX on a Dell desktop computer.

We present, as an aside, a brief proof of the non-oscillation of the system described without delays. This system is described by the equations dm dt ¼ Àk m m þ g 0 þ g 1 À g 0 1 þ p n =X eq ; dp dt ¼ Àk p p þ g p m:

The Jacobian of the system is given by

with eigenvalues that add to a real negative number and multiply to a real positive number. If the eigenvalues are complex, due to the result of multiplication they have the same real part. Due to this and the result of the addition, both eigenvalues must have negative real parts, meaning oscillation is impossible. If the eigenvalues are real, both must be negative, again making oscillation impossible.

It should be noted that we made these conclusions based on the structure of the Jacobian without finding the fixed point, based on the fact that there are only two eigenvalues. The more complex cases described in this work would have three or more eigenvalues, which would make finding the fixed point necessary. They also add to the complexity of the equations for finding the fixed point, which are already of degree n þ 1 with multiple terms and noninteger coefficients.

We begin with the stochastic system described first, in which binding to the gene is coupled with n-mer formation (this could be either because the binding to the gene is itself cooperative or because it is very fast once the n-mer is formed). In the described regime, whose deterministic solutions yield at best decaying oscillation, we now explore the possibility of noise-induced oscillation. We note that the system in this case is between truly oscillatory and simply two-state with reasonably frequent switching. Simple two-state switching would lead to a 'period distribution' (time between maxima) with a normalized standard deviation of that period rt t of 1 ffiffi 2 p . Fig. 1b shows rt t , and Fig. 1a shows the coherence (generalized from the definition used in Ref. [2] ). Coherence is defined as 2 P t Hðd/ðtÞÞ P t jd/ðtÞj À 1 where / is the angle in mRNA-protein space from the average concentration, and d/ðtÞ ¼ /ðtÞ À /ðt À 0:1Þ. It is, in essence, a measure of how often a system goes forward (clockwise) instead of backward (counterclockwise) in its oscillation. Both the coherence and rt t strongly imply that increased cooperativity yields a steadier oscillatory behavior, but that in the regions explored there is no coherent oscillation (which would require a coherence close to 1 and a value of rt t significantly less than 1). So far, the only pieces of the puzzle considered have been multimer binding of proteins to genes, mRNA, and noise. We now add the final piece, the possibility of proteins binding together before they bind to the gene. The number of steps involved here keeps oscillation from occurring at a cooperativity of 1, but increased n causes oscillatory behavior when combined with the rest of the system.

Specifically, Fig. 2 shows a stochastic system with n ¼ 16 and the distribution of periods, defined here as time between local maxima in protein number which are at least 0:4 Á s mRNA (in order to remove less significant fluctuations from consideration). This system is clearly oscillatory, with a reasonably sharp period distribution.

It should be noted that 16 is at best a marginally reasonable value for cooperativity in a simple genetic system. However, it is clear now that simple intermediate steps can interact with cooperativity, and between them can produce oscillations in which neither is the primary factor in the behavior. Furthermore, examining the parameter space in more detail, we find telling behavior in both deterministic and stochastic cases. Fig. 3a shows the deterministic phase diagram of the system, made using timeseries calculations. Region I is oscillatory, region II displays decaying oscillation, and region III is non-oscillatory. High cooperativity and synthesis rate are both clearly necessary for this; also, it should again be noted that the additional intermediate step as compared with the system from Fig. 1 is necessary for deterministic oscillation.

In Fig. 3b , we see, in the same region, a graph of the coherence. It corresponds well with the deterministic calculations: the system reliably goes in the correct direction (associated with time going forward) in the stochastic case when its parameters would predict that it is deterministically oscillatory. Fig. 3c shows rt t , which again corresponds well, although not perfectly; some loss of specificity in period at very high cooperativity and synthesis implies that the oscillation may be imperfect or require that other parameters be very specific at some level. However, the increase in rt t at these points is slight.

From these figures, it is apparent that coherence is possible with higher rates of synthesis, and that there is a tendency towards higher coherence with higher cooperativity. Additionally, we note that in the stochastic case there is no jump in coherence or rt t as we would find if the system exhibited the phase-transition behavior from the deterministic system, Fig. 3a ; in the stochastic case ( Fig. 3b and c) , the line between oscillatory and non-oscillatory is blurred. It is important to note that noise alone is not sufficient for reliable or even semi-reliable oscillation. The deterministically oscillating region clearly has more stable oscillation in the stochastic regime, and systems close to it are more capable of reliable behavior than those far away from it.

For comparison, we now choose a point in the region of deterministic decaying oscillations, and make plots comparable to Fig. 2. In Fig. 4 we see a system whose behavior seems similar in many ways; the main difference in this case is a slightly widened period distribution. While this system is definitely in region II of Fig. 3a (with a cooperativity of 8) , and the period distribution may be too wide for a truly accurate clock, it is obvious that the line between coherent oscillation and incoherent oscillation-like behavior is blurred in this case by noise.

We see in these figures that the optimal cooperativity for this number of intermediate steps in binding is still, at best, at the very high end of cooperativity in reasonable biological systems. Therefore, we have attempted a less complete search of parameters considering additional steps (instead of monomer to n-mer, monomer to dimer to tetramer, etc.). While our searches were not exhaustive, we found at least one region in which such a system can oscillate at quite reasonable values for n, as low as 8 (octamer). For reference, this is not an uncommon amount of cooperativity for these sorts of proteins. SARS-CoV nsp10, for instance, forms a dodecameric structure (n ¼ 12) [19] ; Rhizobium leguminosarum NodD has been shown to bind to DNA preferentially in octamer form [20] ; and k phage CI, while it has a small autoregulatory effect with tetramer and even dimer forms, only fully represses itself when an octamer form of the protein links distant sites on DNA together [21] .

Our data show that oscillation occurs when cooperativity and intermediate steps are both present, when neither cooperativity nor intermediates would otherwise be sufficient. More, it implies that, with many intermediate steps, lower cooperativity can yield oscillation, while with higher cooperativity fewer intermediate steps are necessary. Oscillation behavior in deterministic calculations gives rise to more stable oscillation in the stochastic case.

Additionally, it should be noted that, even though we found no coherent oscillation for lower values of cooperativity, moderate coherence in an oscillation-like behavior may not necessarily be detrimental to a biological system which relies on it as a clock for functions which are not truly vital. Many biological clocks receive input of some kind from outside sources, whether light and food for 24-h clocks or some other form [22, 23] . These inputs can reset a biological clock, and one may argue that they might be able to do so more easily in a fundamentally inexact clock than a single-mode oscillator.

In summary, we have discovered a relationship between generally realistic intermediate steps (as opposed to set-time time delays of the form t À s, which are realistic for a limited set of systems), cooperativity, and oscillation. Intermediate steps result- region I is oscillatory, region II has decaying oscillations, and region III is non-oscillatory. Due to the fact that the distinction betweeen II and III is somewhat dependent on initial conditions, there may be some overlap between those two regions as discussed in the Methods section. (b) Stochastic calculation of coherence, given by 2 P Hðd/Þ P jd/j À 1 where d/ is the difference in angle in mRNA-protein space. (c) Stochastic calculation of the standard deviation of the period distribution divided by its mean. ing from slow cooperative binding can cause oscillation at biologically relevant cooperativity. We have also explored the noise effect, in which the line between oscillatory and non-oscillatory (relatively clear in the deterministic case) is blurred in the stochastic case.

A First Course in Abstract Algebra

The Entrainment of Circadian Systems

Synchronization of the molecular clockwork by light-and food-related cues in mammals

We performed both deterministic and stochastic calculations.