Origin of cyclicity in bipolar disorders: a computational approach.

Bipolar disorders are characterized by the spontaneous, recurrent alternation of episodes of mania and depression. To investigate the type of mechanism capable of accounting for the cyclical nature of the manic depressive illness, we recently proposed a minimal model for bipolar disorders based on the assumption that the propensities to mania and depression are governed by the activities of 2 putative neural circuits that inhibit each other. When mutual inhibition is sufficiently strong, the model predicts bistability: the bipolar system is then in a stable state corresponding either to unipolar depression or mania, and can display abrupt switches between these states. To account for the cyclical nature of bipolar disorders 2 simple, additional regulations allow the model to pass from bistability to oscillations. Self-sustained oscillations provide a mechanism for the spontaneous, recurrent switching between mania and depression. The model can generate oscillations with a variety of waveforms, including periodic oscillations with comparable or unequal durations of the manic and depressive episodes, or small-amplitude oscillations around one of the 2 states preceding large-amplitude periodic changes in the propensities to mania or depression, with phases during which these propensities reach intermediate levels, a situation that could correspond to mixed bipolar states. Oscillations become irregular when fluctuations of parameter values are taken into account. The model provides a theoretical framework that covers the bipolar spectrum, i. e., cycling between the 2 poles of the disease, or evolution to a stable steady state corresponding to various degrees of unipolar depression or mania or to a "normal" state in which the -propensities to mania or depression remain low, without alternation between the 2 poles of the disease. The computational approach may help the exploration of plausible mechanisms for bipolar disorders and possible dynamic bases for clinical observations on the effect of antidepressants, which can trigger the transition to mania or increase the frequency of bipolar cycling.

Systems biology of cellular rhythms.

Rhythms abound in biological systems, particularly at the cellular level where they originate from the feedback loops present in regulatory networks. Cellular rhythms can be investigated both by experimental and modeling approaches, and thus represent a prototypic field of research for systems biology. They have also become a major topic in synthetic biology. We review advances in the study of cellular rhythms of biochemical rather than electrical origin by considering a variety of oscillatory processes such as Ca++ oscillations, circadian rhythms, the segmentation clock, oscillations in p53 and NF-κB, synthetic oscillators, and the oscillatory dynamics of cyclin-dependent kinases driving the cell cycle. Finally we discuss the coupling between cellular rhythms and their robustness with respect to molecular noise.

Modelling the dual role of Per phosphorylation and its effect on the period and phase of the mammalian circadian clock.

Circadian clocks are regulated at the post-translational level by a variety of processes among which protein phosphorylation plays a prominent, although complex, role. Thus, the phosphorylation of different sites on the clock protein PER by casein kinase I (CKI) can lead to opposite effects on the stability of the protein and on the period of circadian oscillations. Here the authors extend a computational model previously proposed for the mammalian circadian clock by incorporating two distinct phosphorylations of PER by CKI. On the basis of experimental observations the authors consider that phosphorylation at one site (denoted here PER-P1) enhances the rate of degradation of the protein and decreases the period, while phosphorylation at another site (PER-P2) stabilises the protein, enhances the transcription of the Per gene, and increases the period. The model also incorporates an additional phosphorylation of PER by the Glycogen Synthase Kinase 3 (GSK3). The authors show that the extended model incorporating the antagonistic effects of PER phosphorylations by CKI can account for observations pertaining to (i) the decrease in period in the Tau mutant, because of an increase in phosphorylation by CKI leading to PER-P1, and (ii) the familial advanced sleep phase syndrome (FASPS) in which the period is shortened and the phase of the oscillations is advanced when the rate of phosphorylation leading to PER-P2 is decreased. The model further accounts for the increase in period observed in the presence of CKI inhibitors that decrease the rate of phosphorylation leading to both PER-P1 and PER-P2. A similar increase in period results from inhibition of GSK3. [Includes supplementary material].

[Deterministic and stochastic models for circadian rhythms]. / Modèles déterministes et stochastiques pour les rythmes circadiens.

Circadian rhythms, characterized by a period of about 24h, are generated in nearly all living organisms by the negative autoregulation of clock gene expression. Deterministic models based on this genetic regulation account for circadian oscillations in constant environmental conditions (e.g., in constant darkness) and for entrainment of these rhythms by light-dark cycles. When the number of clock mRNA and protein molecules is low, it is necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. Indeed, it is possible that the autoregulatory mechanism of gene expression might not produce stable rhythms due to fluctuations if the number of molecules involved in the clock mechanism remains too low. We have compared the deterministic and stochastic approaches for a model based on the negative autoregulation of a clock gene. We show by means of stochastic simulations that robust circadian oscillations can already occur when the maximum number of mRNA and protein molecules is of the order of a few tens or hundreds, respectively. Furthermore, the results indicate that the cooperativity characterizing the repression of the transcription process strenghtens the robustness of circadian rhythms and that entrainment by light-dark cycles stabilizes the phase of the oscillations.

The follicular automaton model: effect of stochasticity and of synchronization of hair cycles.

Human scalp hair consists of a set of about 10(5)follicles which progress independently through developmental cycles. Each hair follicle successively goes through the anagen (A), catagen (C), telogen (T) and latency (L) phases that correspond, respectively, to growth, arrest and hair shedding before a new anagen phase is initiated. Long-term experimental observations in a group of ten male, alopecic and non-alopecic volunteers allowed determination of the characteristics of hair follicle cycles. On the basis of these observations, we previously proposed a follicular automaton model to simulate the dynamics of human hair cycles and the development of different patterns of alopecia [Halloy et al. (2000) Proc. Natl Acad. Sci. U.S.A.97, 8328-8333]. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states A, T, L and the subsequent return to A. These transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a standard deviation. The follicular automaton model was shown to account both for the dynamical transitions observed in a single follicle, and for the behaviour of an ensemble of independently cycling follicles. Here, we extend these results and investigate additional properties of the model. We present a deterministic version of the follicular automaton. We show that numerical simulations of the stochastic version of the automaton yield steady-state level of follicles in the different phases which approach the levels predicted by the deterministic equations as the number of follicles progressively increases. Only the stochastic version can successfully reproduce the fluctuations of the fractions of follicles in each of the three phases, observed in small follicle populations. When the standard deviation is reduced or when the follicles become otherwise synchronized, e.g. by a periodic external signal inducing the transition of anagen follicles into telogen phase, large-amplitude oscillations occur in the fractions of follicles in the three phases. These oscillations are not observed in humans but are reminiscent of the phenomenon of moulting observed in a number of mammalian species.

Deterministic versus stochastic models for circadian rhythms.

Circadian rhythms which occur with a period close to 24 h in nearly all living organisms originate from the negative autoregulation of gene expression.Deterministic models based on genetic regulatory processes account for theoccurrence 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. At low numbers of protein and mRNA molecules, 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 two stochastic versions of a core model for circadian rhythms. The deterministic version of this core modelwas previously proposed for circadian oscillations of the PER protein in Drosophila and of the FRQ protein in Neurospora. In the first, non-developed version of the stochastic model, we introduce molecular noise without decomposing the deterministic mechanism into detailed reaction steps while in the second, developed version we carry out such a detailed decomposition. Numerical simulations of the two stochastic versions of the model are performed by means of the Gillespie method. We compare the predictions of the deterministic approach with those of the two stochastic models, with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity of a bifurcation point beyond which the system evolves to a stable steady state. The results indicate that robust circadian oscillations can occur even when the numbers of mRNA and nuclear protein involved in the oscillatory mechanism are reduced to a few tens orhundreds, respectively. The non-developed and developed versions of the stochastic model yield largely similar results and provide good agreement with the predictions of the deterministic model for circadian rhythms.

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.

A molecular explanation for the long-term suppression of circadian rhythms by a single light pulse.

With the use of a molecular model for circadian rhythms in Drosophila based on transcriptional regulation, we show how a single, critical pulse of light can permanently suppress circadian rhythmicity, whereas a second light pulse can restore the abolished rhythm. The phenomena occur via the pulsatile induction of either protein degradation or gene expression in conditions in which a stable steady state coexists with stable circadian oscillations of the limit cycle type. The model indicates that suppression by a light pulse can only be accounted for by assuming that the biochemical effects of such a pulse much outlast its actual duration. We determine the characteristics of critical pulses suppressing the oscillations as a function of the phase at which the rhythm is perturbed. The model predicts how the amplitude and duration of the biochemical changes induced by critical pulses vary with this phase. The results provide a molecular, dynamic explanation for the long-term suppression of circadian rhythms observed in a variety of organisms in response to a single light pulse and for the subsequent restoration of the rhythms by a second light pulse.

Modeling the dynamics of human hair cycles by a follicular automaton.

The hair follicle cycle successively goes through the anagen, catagen, telogen, and latency phases, which correspond, respectively, to hair growth, arrest, shedding, and absence before a new anagen phase is initiated. Experimental observations collected over a period of 14 years in a group of 10 male volunteers, alopecic and nonalopecic, allowed us to determine the characteristics of scalp hair follicle cycles. On the basis of these observations, we propose a follicular automaton model to simulate the dynamics of human hair cycles. The automaton model is defined by a set of rules that govern the stochastic transitions of each follicle between the successive states anagen, telogen, and latency, and the subsequent return to anagen. The transitions occur independently for each follicle, after time intervals given stochastically by a distribution characterized by a mean and a variance. The follicular automaton model accounts both for the dynamical transitions observed in a single follicle and for the behavior of an ensemble of independently cycling follicles. Thus, the model successfully reproduces the evolution of the fractions of follicle populations in each of the three phases, which fluctuate around steady-state or slowly drifting values. We apply the follicular automaton model to the study of spatial patterns of follicular growth that result from a spatially heterogeneous distribution of parameters such as the mean duration of anagen phase. When considering that follicles die or miniaturize after going through a critical number of successive cycles, the model can reproduce the evolution to hair patterns similar to well known types of diffuse or androgenetic alopecia.

Modeling the differential fitness of cyanobacterial strains whose circadian oscillators have different free-running periods: comparing the mutual inhibition and substrate depletion hypotheses.

In a recent experimental study, Ouyang et al. (1998, Proc. Natl. Acad. Sci. U.S.A.95, 8660-8664) have shown that, in direct competition, cyanobacterial strains whose circadian clocks have free-running periods (FRPs) which match the period of an imposed light/dark (LD) cycle exclude strains whose FRPs are out of resonance with the LD cycle. These differences in competitive fitness are observed despite the lack of measurable differences in monoculture growth rates between the strains. Here we show that the experimental results are consistent with a mathematical model in which cells rhythmically produce a metabolic inhibitor to which they display a sensitivity modulated by their circadian rhythm. We argue that models in which there is a circadian modulation of nutrient uptake kinetics cannot account for the results of these experiments. We discuss possible experiments to further characterize this phenomenon. The experimental protocol we propose can be used to distinguish between mutual inhibition and substrate depletion as underlying causes of the competitive advantage of circadian resonance.

Theoretical models for circadian rhythms in Neurospora and Drosophila.

We examine theoretical models proposed for the molecular mechanism of circadian rhythms in Drosophila. The models are based on the negative feedback exerted by a complex between the PER and TIM proteins on the expression of the per and tim genes. We show that a similar model can account for circadian oscillations in Neurospora, where the protein FRQ negatively regulates the expression of the frq gene. The effect of light on the circadian rhythms is included by considering that it elicits a rise in the rate of TIM degradation in Drosophila, whereas in Neurospora it enhances the rate of frq transcription. The models account for the occurrence of sustained circadian oscillations in continuous darkness in Drosophila and Neurospora. Numerical simulations further indicate that the periodic forcing of circadian oscillations by light-dark cycles can result either in the entrainment to the external periodicity or in aperiodic oscillations (i.e. chaos), depending on the magnitude of the periodic changes in the light-controlled parameter.

The frequency encoding of pulsatility.

Examples of pulsatile signalling abound in intercellular communication, suggesting that this phenomenon represents a major function of biological rhythms. Pulsatile signals can be encoded in terms of their frequency and prove more efficient than monotonous ones whenever constant stimulation induces desensitization of target cells. We address the main examples of frequency encoding of pulsatility, besides those of neuronal nature. Considered in turn are cAMP oscillations in the slime mould Dictyostelium discoideum, the pulsatile secretion of hormones such as gonadotropin-releasing hormone or growth hormone, intracellular Ca2+ oscillations, and circadian rhythms. Models based on receptor desensitization show the possibility of optimizing cellular responses to cAMP signals in Dictyostelium or to pulsatile hormonal stimulation. The models indicate how the optimal duration of the pulsatile signal and the optimal interval between successive pulses vary as a function of the rates or receptor desensitization and resensitization and of the maximum ligand level during stimulation. The frequency encoding of intracellular Ca2+ oscillations appears to rely on another molecular mechanism. Models based on protein phosphorylation by a Ca(2+)-calmodulin activated kinase show that the mean level of phosphorylated protein increases with the frequency of calcium spikes--which itself rises with the degree of stimulation--and that distinct levels of different phosphorylated proteins can be reached for a Ca2+ signal of given frequency.

Oscillations and bistability predicted by a model for a cyclical bienzymatic system involving the regulated isocitrate dehydrogenase reaction.

We analyze the dynamics of a bienzymatic system consisting of isocitrate dehydrogenase (IDH, EC. 1.1.1.42), which transforms NADP+ into NADPH, and of diaphorase (DIA, EC 1.8.1.4), which catalyzes the reverse reaction. Experimental evidence as well as a theoretical model showed the possibility of a coexistence between two stable steady states in this reaction system G.M. Guidi et al. Biophys. J. 74 (1998) 1229-1240[, owing to the regulatory properties of IDH. Here we extend this analysis by considering the behavior of the model proposed for the IDH-DIA bienzymatic system in conditions where the system is open to an influx of its substrates isocitrate and NADP+ and to an efflux of all metabolic species. The analysis indicates that in addition to different modes of bistability (including mushrooms and isolas), sustained oscillations can be observed in such conditions. These results point to the isocitrate dehydrogenase reaction coupled to diaphorase as a suitable candidate for further experimental and theoretical studies of bistability and oscillations in biochemical systems. The results obtained in this particular bienzymatic system bear on other enzymatic systems possessing a cyclical nature, which are known to play significant roles in a variety of metabolic and cellular regulatory processes.

Modeling the molecular regulatory mechanism of circadian rhythms in Drosophila.

Thanks to genetic and biochemical advances on the molecular mechanism of circadian rhythms in Drosophila, theoretical models closely related to experimental observations can be considered for the regulatory mechanism of the circadian clock in this organism. Modeling is based on the autoregulatory negative feedback exerted by a complex between PER and TIM proteins on the expression of per and tim genes. The model predicts the occurrence of sustained circadian oscillations in continuous darkness. When incorporating light-induced TIM degradation, the model accounts for damping of oscillations in constant light, entrainment of the rhythm by light-dark cycles of varying period or photoperiod, and phase shifting by light pulses. The model further provides a molecular dynamical explanation for the permanent or transient suppression of circadian rhythmicity triggered in a variety of organisms by a critical pulse of light. Finally, the model shows that to produce a robust rhythm the various clock genes must be expressed at the appropriate levels since sustained oscillations only occur in a precise range of parameter values. BioEssays 22:84-93, 2000.

Alternating oscillations and chaos in a model of two coupled biochemical oscillators driving successive phases of the cell cycle.

The animal cell cycle is controlled by the periodic variation of two cyclin-dependent protein kinases, cdk1 and cdk2, which govern the entry into the M (mitosis) and S (DNA replication) phases, respectively. The ordered progression between these phases is achieved thanks to the existence of checkpoint mechanisms based on mutual inhibition of these processes. Here we study a simple theoretical model for oscillations in cdk1 and cdk2 activity, involving mutual inhibition of the two oscillators. Each minimal oscillator is described by a three-variable cascade involving a cdk, together with the associated cyclin and cyclin-degrading enzyme. The dynamics of this skeleton model of coupled oscillators is determined as a function of the strength of their mutual inhibition. The most common mode of dynamic behavior, obtained under conditions of strong mutual inhibition, is that of alternating oscillations in cdk1 and cdk2, which correspond to the physiological situation of the ordered recurrence of the M and S phases. In addition, for weaker inhibition we obtain evidence for a variety of dynamic phenomena such as complex periodic oscillations, chaos, and the coexistence between multiple periodic or chaotic attractors. We discuss the conditions of occurrence of these various modes of oscillatory behavior, as well as their possible physiological significance.

Chaos and birhythmicity in a model for circadian oscillations of the PER and TIM proteins in drosophila

In Drosophila, circadian oscillations in the levels of two proteins, PER and TIM, result from the negative feedback exerted by a PER-TIM complex on the expression of the per and tim genes which code for these two proteins. On the basis of these experimental observations, we have recently proposed a theoretical model for circadian oscillations of the PER and TIM proteins in Drosophila. Here we show that for constant environmental conditions this model is capable of generating autonomous chaotic oscillations. For other parameter values, the model can also display birhythmicity, i.e. the coexistence between two stable regimes of limit cycle oscillations. We analyse the occurrence of chaos and birhythmicity by means of bifurcation diagrams and locate the different domains of complex oscillatory behavior in parameter space. The relative smallness of these domains raises doubts as to the possible physiological significance of chaos and birhythmicity in regard to circadian rhythm generation. Beyond the particular context of circadian rhythms we discuss the results in the light of other mechanisms underlying chaos and birhythmicity in regulated biological systems. Copyright 1999 Academic Press.

Limit cycle models for circadian rhythms based on transcriptional regulation in Drosophila and Neurospora.

We examine theoretical models for circadian oscillations based on transcriptional regulation in Drosophila and Neurospora. For Drosophila, the molecular model is based on the negative feedback exerted on the expression of the per and tim genes by the complex formed between the PER and TIM proteins. For Neurospora, similarly, the model relies on the feedback exerted on the expression of the frq gene by its protein product FRQ. In both models, sustained rhythmic variations in protein and mRNA levels occur in continuous darkness, in the form of limit cycle oscillations. The effect of light on circadian rhythms is taken into account in the models by considering that it triggers degradation of the TIM protein in Drosophila, and frq transcription in Neurospora. When incorporating the control exerted by light at the molecular level, we show that the models can account for the entrainment of circadian rhythms by light-dark cycles and for the damping of the oscillations in constant light, though such damping occurs more readily in the Drosophila model. The models account for the phase shifts induced by light pulses and allow the construction of phase response curves. These compare well with experimental results obtained in Drosophila. The model for Drosophila shows that when applied at the appropriate phase, light pulses of appropriate duration and magnitude can permanently or transiently suppress circadian rhythmicity. We investigate the effects of the magnitude of light-induced changes on oscillatory behavior. Finally, we discuss the common and distinctive features of circadian oscillations in the two organisms.

Bursting, chaos and birhythmicity originating from self-modulation of the inositol 1,4,5-trisphosphate signal in a model for intracellular Ca2+ oscillations.

We investigate the various types of complex Ca2+ oscillations which can arise in a model based on the mechanism of Ca2+-induced Ca2+ release (CICR), that takes into account the Ca2+-stimulated degradation of inositol 1,4,5-trisphosphate (InsP3) by a 3-kinase. This model was previously proposed in the course of an investigation of plausible mechanisms capable of generating complex Ca2+ oscillations. Besides simple periodic behavior, this model for cytosolic Ca2+ oscillations in nonexcitable cells shows complex oscillatory phenomena like bursting or chaos. We show that the model also admits a coexistence between two stable regimes of sustained oscillations (birhythmicity). The occurrence of these various modes of oscillatory behavior is analysed by means of bifurcation diagrams. Complex oscillations are characterized by means of Poincaré sections, power spectra and Lyapounov exponents. The results point to the role of self-modulation of the InsP3 signal by 3-kinase as a possible source for complex temporal patterns in Ca2+ signaling.

CaM kinase II as frequency decoder of Ca2+ oscillations.

In many cell types, Ca2+ signals are organized in the form of repetitive spikes. The frequency of these intracellular Ca2+ oscillations increases with the level of stimulation, suggesting the existence of a frequency encoding phenomenon. The question arises as to how the frequency of Ca2+ oscillations can be decoded inside the cell. Ca2+/calmodulin kinase II has long been proposed as an attractive candidate, as it is a key target of Ca2+ signals. By immobilizing the Ca2+/calmodulin kinase II and subjecting it to pulses of Ca2+ of variable amplitude, duration, and frequency, De Koninck and Schulman have shown for the first time that the autonomous activity of Ca2+/calmodulin kinase II is highly sensitive to the temporal pattern of Ca2+ oscillations.

Modeling oscillations and waves of cAMP in Dictyostelium discoideum cells.

We examine the theoretical aspects of temporal and spatiotemporal organization in the cAMP signaling system of Dictyostelium discoideum amoebae which aggregate in a wavelike manner after starvation, in response to pulses of cAMP emitted with a periodicity of several minutes by cells behaving as aggregation centers. We first extend the model based on receptor desensitization, previously proposed by Martiel and Goldbeter, by incorporating the role of G proteins in signal transduction. The extended model accounts for observations on the response of the signaling system to successive step increases in extracellular cAMP. In the presence of the positive feedback loop in cAMP synthesis, this model generates sustained oscillations in cAMP and in the fraction of active cAMP receptor, similar to those obtained in the simpler model where the role of the G proteins is not taken into account explicitly. We use the latter model to address the formation of concentric and spiral waves of cAMP in the course of D. discoideum aggregation. Previous analyses of the model showed that a progressive increase in the activity of adenylate cyclase and phosphodiesterase can account for the transitions no relay-relay-oscillations-relay observed in the experiments. We show that the degree of cellular synchronization on such a developmental path in parameter space markedly affects the nature of the spatial patterns generated by the model. These patterns range from concentric waves to a small number of large spirals, and finally to a large number of smaller spirals, as the degree of developmental desynchronization between cells increases.