Modeling the circadian clock: From molecular mechanism to physiological disorders

Emergent robust oscillatory dynamics in the interlocked feedback-feedforward loops

One of the challenges that beset modelling complex biological networks is to relate networks to function to dynamics. A further challenge is deciphering the cellular function and dynamics that can change drastically when the network edge is tinkered with by adding or removing it. To illustrate this, the authors took a well-studied three-variable Goodwin oscillatory motif with only a negative feedback loop. To this motif, an edge was added that results in an emergent structure consisting of new feedforward and feedback loops while retaining Goodwin's original negative feedback loop. To relate emergent structure to oscillatory dynamics, the authors took all the combinations of edge signs in the interlocked motif. Bifurcation analysis reveals that all the structural combinations can be grouped into two categories based on their unique dynamics. These two groups also exhibit unique amplitude-frequency (amp-freq) plots. These two categories are attributed to the emergence of interlocked motifs with specific edge signs. To support the ideas, a well-studied plant circadian model of Arabidopsis thaliana was taken to illustrate the importance of interlocked motifs in fine-tuning amplitude and frequency in circadian oscillators. The authors briefly discuss its implications for central oscillators' adaptation to different environmental cues.

Limits of entrainment of circadian neuronal networks

Circadian rhythmicity lies at the center of various important physiological and behavioral processes in mammals, such as sleep, metabolism, homeostasis, mood changes, and more. Misalignment of intrinsic neuronal oscillations with the external day-night cycle can disrupt such processes and lead to numerous disorders. In this work, we computationally determine the limits of circadian synchronization to external light signals of different frequency, duty cycle, and simulated amplitude. Instead of modeling circadian dynamics with generic oscillator models (e.g., Kuramoto-type), we use a detailed computational neuroscience model, which integrates biomolecular dynamics, neuronal electrophysiology, and network effects. This allows us to investigate the effect of small drug molecules, such as Longdaysin, and connect our results with experimental findings. To combat the high dimensionality of such a detailed model, we employ a matrix-free approach, while our entire algorithmic pipeline enables numerical continuation and construction of bifurcation diagrams using only direct simulation. We, thus, computationally explore the effect of heterogeneity in the circadian neuronal network, as well as the effect of the corrective therapeutic intervention of Longdaysin. Last, we employ unsupervised learning to construct a data-driven embedding space for representing neuronal heterogeneity.

Cycles in epilepsy

Epilepsy is among the most dynamic disorders in neurology. A canonical view holds that seizures, the characteristic sign of epilepsy, occur at random, but, for centuries, humans have looked for patterns of temporal organization in seizure occurrence. Observations that seizures are cyclical date back to antiquity, but recent technological advances have, for the first time, enabled cycles of seizure occurrence to be quantitatively characterized with direct brain recordings. Chronic recordings of brain activity in humans and in animals have yielded converging evidence for the existence of cycles of epileptic brain activity that operate over diverse timescales: daily (circadian), multi-day (multidien) and yearly (circannual). Here, we review this evidence, synthesizing data from historical observational studies, modern implanted devices, electronic seizure diaries and laboratory-based animal neurophysiology. We discuss advances in our understanding of the mechanistic underpinnings of these cycles and highlight the knowledge gaps that remain. The potential clinical applications of a knowledge of cycles in epilepsy, including seizure forecasting and chronotherapy, are discussed in the context of the emerging concept of seizure risk. In essence, this Review addresses the broad question of why seizures occur when they occur.

Brain aging and garbage cleaning : Modelling the role of sleep, glymphatic system, and microglia senescence in the propagation of inflammaging

Brain aging is a complex process involving many functions of our body and described by the interplay of a sleep pattern and changes in the metabolic waste concentration regulated by the microglial function and the glymphatic system. We review the existing modelling approaches to this topic and derive a novel mathematical model to describe the crosstalk between these components within the conceptual framework of inflammaging. Analysis of the model gives insight into the dynamics of garbage concentration and linked microglial senescence process resulting from a normal or disrupted sleep pattern, hence, explaining an underlying mechanism behind healthy or unhealthy brain aging. The model incorporates accumulation and elimination of garbage, induction of glial activation by garbage, and glial senescence by over-activation, as well as the production of pro-inflammatory molecules by their senescence-associated secretory phenotype (SASP). Assuming that insufficient sleep leads to the increase of garbage concentration and promotes senescence, the model predicts that if the accumulation of senescent glia overcomes an inflammaging threshold, further progression of senescence becomes unstoppable even if a normal sleep pattern is restored. Inverting this process by "rejuvenating the brain" is only possible via a reset of concentration of senescent glia below this threshold. Our model approach enables analysis of space-time dynamics of senescence, and in this way, we show that heterogeneous patterns of inflammation will accelerate the propagation of senescence profile through a network, confirming a negative effect of heterogeneity.

Astrocytic Modulation of Neuronal Activity in the Suprachiasmatic Nucleus: Insights from Mathematical Modeling

The suprachiasmatic nucleus (SCN) of the hypothalamus consists of a highly heterogeneous neuronal population networked together to allow precise and robust circadian timekeeping in mammals. While the critical importance of SCN neurons in regulating circadian rhythms has been extensively studied, the roles of SCN astrocytes in circadian system function are not well understood. Recent experiments have demonstrated that SCN astrocytes are circadian oscillators with the same functional clock genes as SCN neurons. Astrocytes generate rhythmic outputs that are thought to modulate neuronal activity through pre- and postsynaptic interactions. In this study, we developed an in silico multicellular model of the SCN clock to investigate the impact of astrocytes in modulating neuronal activity and affecting key clock properties such as circadian rhythmicity, period, and synchronization. The model predicted that astrocytes could alter the rhythmic activity of neurons via bidirectional interactions at tripartite synapses. Specifically, astrocyte-regulated extracellular glutamate was predicted to increase neuropeptide signaling from neurons. Consistent with experimental results, we found that astrocytes could increase the circadian period and enhance neural synchronization according to their endogenous circadian period. The impact of astrocytic modulation of circadian rhythm amplitude, period, and synchronization was predicted to be strongest when astrocytes had periods between 0 and 2 h longer than neurons. Increasing the number of neurons coupled to the astrocyte also increased its impact on period modulation and synchrony. These computational results suggest that signals that modulate astrocytic rhythms or signaling (e.g., as a function of season, age, or treatment) could cause disruptions in circadian rhythm or serve as putative therapeutic targets.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality, i.e., whether these networks can perform their function or not, depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease.

Dissipative structures and biological rhythms

Sustained oscillations abound in biological systems. They occur at all levels of biological organization over a wide range of periods, from a fraction of a second to years, and with a variety of underlying mechanisms. They control major physiological functions, and their dysfunction is associated with a variety of physiological disorders. The goal of this review is (i) to give an overview of the main rhythms observed at the cellular and supracellular levels, (ii) to briefly describe how the study of biological rhythms unfolded in the course of time, in parallel with studies on chemical oscillations, (iii) to present the major roles of biological rhythms in the control of physiological functions, and (iv) the pathologies associated with the alteration, disappearance, or spurious occurrence of biological rhythms. Two tables present the main examples of cellular and supracellular rhythms ordered according to their period, and their role in physiology and pathophysiology. Among the rhythms discussed are neural and cardiac rhythms, metabolic oscillations such as those occurring in glycolysis in yeast, intracellular Ca⁺⁺ oscillations, cyclic AMP oscillations in Dictyostelium amoebae, the segmentation clock that controls somitogenesis, pulsatile hormone secretion, circadian rhythms which occur in all eukaryotes and some bacteria with a period close to 24 h, the oscillatory dynamics of the enzymatic network driving the cell cycle, and oscillations in transcription factors such as NF-ΚB and tumor suppressors such as p53. Ilya Prigogine's concept of dissipative structures applies to temporal oscillations and allows us to unify within a common framework the various rhythms observed at different levels of biological organization, regardless of their period and underlying mechanism.

Computational modeling of the cell-autonomous mammalian circadian oscillator

This review summarizes various mathematical models of cell-autonomous mammalian circadian clock. We present the basics necessary for understanding of the cell-autonomous mammalian circadian oscillator, modern experimental data essential for its reconstruction and some special problems related to the validation of mathematical circadian oscillator models. This work compares existing mathematical models of circadian oscillator and the results of the computational studies of the oscillating systems. Finally, we discuss applications of the mathematical models of mammalian circadian oscillator for solving specific problems in circadian rhythm biology.

Mechanists Must be Holists Too! Perspectives from Circadian Biology

The pursuit of mechanistic explanations in biology has produced a great deal of knowledge about the parts, operations, and organization of mechanisms taken to be responsible for biological phenomena. Holist critics have often raised important criticisms of proposed mechanistic explanations, but until recently holists have not had alternative research strategies through which to advance explanations. This paper argues both that the results of mechanistic strategies has forced mechanists to confront ways in which whole systems affect their components and that new representational and modeling strategies are providing tools for understanding these effects of whole systems upon components. Drawing from research on the mechanism responsible for circadian rhythms in mammals, I develop two examples in which mechanistic analysis is being integrated into a more holist perspective: research revealing intercellular integration of circadian mechanisms with those involved in cell metabolism and research revealing that stable␣rhythms are dependent on how individual cells in the suprachiasmatic nucleus synchronize with each other to generate regular rhythms. Tools such as network diagramming and computational modeling are providing means to integrate mechanistic models into accounts of whole systems.

Modulation of E-cadherin expression promotes migration ability of esophageal cancer cells

Losing the E-cadherin plays an important role in the metastasis of cancer. The regulation of the expression of E-cadherin is unclear. Circadian rhythm alteration is associated with the pathogenesis of a number of cancers. This study aims to investigate the role of one of the circadian proteins, period-2 (Per2) in repressing the expression of E-cadherin in esophageal cancer (esophageal cancer). We observed that the levels of circadian protein Per2 were significantly increased and E-cadherin was significantly decreased in the tissue of human esophageal cancer with metastasis as compared with non-metastatic esophageal cancer. Overexpression of Per2 in the esophageal cancer cells markedly repressed the expression of E-cadherin. The pHDAC1 was detected in human esophageal cancer with metastasis, which was much less in the esophageal cancer tissue without metastasis. Overexpression of Per2 increased the levels of pHDAC1 as well as the E-cadherin repressors at the E-cadherin promoter locus. Overexpression of Per2 markedly increased the migratory capacity of esophageal cancer cells, which was abolished by the inhibition of HDAC1. We conclude that Per-2 plays an important role in the esophageal cancer cell metastasis, which may be a novel therapeutic target for the treatment of esophageal cancer.

Analytical study of robustness of a negative feedback oscillator by multiparameter sensitivity

Background

One of the distinctive features of biological oscillators such as circadian clocks and cell cycles is robustness which is the ability to resume reliable operation in the face of different types of perturbations. In the previous study, we proposed multiparameter sensitivity (MPS) as an intelligible measure for robustness to fluctuations in kinetic parameters. Analytical solutions directly connect the mechanisms and kinetic parameters to dynamic properties such as period, amplitude and their associated MPSs. Although negative feedback loops are known as common structures to biological oscillators, the analytical solutions have not been presented for a general model of negative feedback oscillators.

Results

We present the analytical expressions for the period, amplitude and their associated MPSs for a general model of negative feedback oscillators. The analytical solutions are validated by comparing them with numerical solutions. The analytical solutions explicitly show how the dynamic properties depend on the kinetic parameters. The ratio of a threshold to the amplitude has a strong impact on the period MPS. As the ratio approaches to one, the MPS increases, indicating that the period becomes more sensitive to changes in kinetic parameters. We present the first mathematical proof that the distributed time-delay mechanism contributes to making the oscillation period robust to parameter fluctuations. The MPS decreases with an increase in the feedback loop length (i.e., the number of molecular species constituting the feedback loop).

Conclusions

Since a general model of negative feedback oscillators was employed, the results shown in this paper are expected to be true for many of biological oscillators. This study strongly supports that the hypothesis that phosphorylations of clock proteins contribute to the robustness of circadian rhythms. The analytical solutions give synthetic biologists some clues to design gene oscillators with robust and desired period.

Strategy revealing phenotypic differences among synthetic oscillator designs

Considerable progress has been made in identifying and characterizing the component parts of genetic oscillators, which play central roles in all organisms. Nonlinear interaction among components is sufficiently complex that mathematical models are required to elucidate their elusive integrated behavior. Although natural and synthetic oscillators exhibit common architectures, there are numerous differences that are poorly understood. Utilizing synthetic biology to uncover basic principles of simpler circuits is a way to advance understanding of natural circadian clocks and rhythms. Following this strategy, we address the following questions: What are the implications of different architectures and molecular modes of transcriptional control for the phenotypic repertoire of genetic oscillators? Are there designs that are more realizable or robust? We compare synthetic oscillators involving one of three architectures and various combinations of the two modes of transcriptional control using a methodology that provides three innovations: a rigorous definition of phenotype, a procedure for deconstructing complex systems into qualitatively distinct phenotypes, and a graphical representation for illuminating the relationship between genotype, environment, and the qualitatively distinct phenotypes of a system. These methods provide a global perspective on the behavioral repertoire, facilitate comparisons of alternatives, and assist the rational design of synthetic gene circuitry. In particular, the results of their application here reveal distinctive phenotypes for several designs that have been studied experimentally as well as a best design among the alternatives that has yet to be constructed and tested.

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.

Forced desynchronization of activity rhythms in a model of chronic jet lag in mice

We studied locomotor activity rhythms of C57/Bl6 mice under a chronic jet lag (CJL) protocol (ChrA(6/2) ), which consisted of 6-hour phase advances of the light-dark schedule (LD) every 2 days. Through periodogram analysis, we found 2 components of the activity rhythm: a short-period component (21.01 ± 0.04 h) that was entrained by the LD schedule and a long-period component (24.68 ± 0.26 h). We developed a mathematical model comprising 2 coupled circadian oscillators that was tested experimentally with different CJL schedules. Our simulations suggested that under CJL, the system behaves as if it were under a zeitgeber with a period determined by (24 - [phase shift size/days between shifts]). Desynchronization within the system arises according to whether this effective zeitgeber is inside or outside the range of entrainment of the oscillators. In this sense, ChrA(6/2) is interpreted as a (24 - 6/2 = 21 h) zeitgeber, and simulations predicted the behavior of mice under other CJL schedules with an effective 21-hour zeitgeber. Animals studied under an asymmetric T = 21 h zeitgeber (carried out by a 3-hour shortening of every dark phase) showed 2 activity components as observed under ChrA(6/2): an entrained short-period (21.01 ± 0.03 h) and a long-period component (23.93 ± 0.31 h). Internal desynchronization was lost when mice were subjected to 9-hour advances every 3 days, a possibility also contemplated by the simulations. Simulations also predicted that desynchronization should be less prevalent under delaying than under advancing CJL. Indeed, most mice subjected to 6-hour delay shifts every 2 days (an effective 27-hour zeitgeber) displayed a single entrained activity component (26.92 ± 0.11 h). Our results demonstrate that the disruption provoked by CJL schedules is not dependent on the phase-shift magnitude or the frequency of the shifts separately but on the combination of both, through its ratio and additionally on their absolute values. In this study, we present a novel model of forced desynchronization in mice under a specific CJL schedule; in addition, our model provides theoretical tools for the evaluation of circadian disruption under CJL conditions that are currently used in circadian research.

Tuning the mammalian circadian clock: robust synergy of two loops

The circadian clock is accountable for the regulation of internal rhythms in most living organisms. It allows the anticipation of environmental changes during the day and a better adaptation of physiological processes. In mammals the main clock is located in the suprachiasmatic nucleus (SCN) and synchronizes secondary clocks throughout the body. Its molecular constituents form an intracellular network which dictates circadian time and regulates clock-controlled genes. These clock-controlled genes are involved in crucial biological processes including metabolism and cell cycle regulation. Its malfunction can lead to disruption of biological rhythms and cause severe damage to the organism. The detailed mechanisms that govern the circadian system are not yet completely understood. Mathematical models can be of great help to exploit the mechanism of the circadian circuitry. We built a mathematical model for the core clock system using available data on phases and amplitudes of clock components obtained from an extensive literature search. This model was used to answer complex questions for example: how does the degradation rate of Per affect the period of the system and what is the role of the ROR/Bmal/REV-ERB (RBR) loop? Our findings indicate that an increase in the RNA degradation rate of the clock gene Period (Per) can contribute to increase or decrease of the period--a consequence of a non-monotonic effect of Per transcript stability on the circadian period identified by our model. Furthermore, we provide theoretical evidence for a potential role of the RBR loop as an independent oscillator. We carried out overexpression experiments on members of the RBR loop which lead to loss of oscillations consistent with our predictions. These findings challenge the role of the RBR loop as a merely auxiliary loop and might change our view of the clock molecular circuitry and of the function of the nuclear receptors (REV-ERB and ROR) as a putative driving force of molecular oscillations.

Circadian regulation of sleep-wake behaviour in nocturnal rats requires multiple signals from suprachiasmatic nucleus

The dynamics of sleep and wake are strongly linked to the circadian clock. Many models have accurately predicted behaviour resulting from dynamic interactions between these two systems without specifying physiological substrates for these interactions. By contrast, recent experimental work has identified much of the relevant physiology for circadian and sleep-wake regulation, but interaction dynamics are difficult to study experimentally. To bridge these approaches, we developed a neuronal population model for the dynamic, bidirectional, neurotransmitter-mediated interactions of the sleep-wake and circadian regulatory systems in nocturnal rats. This model proposes that the central circadian pacemaker, located within the suprachiasmatic nucleus (SCN) of the hypothalamus, promotes sleep through single neurotransmitter-mediated signalling to sleep-wake regulatory populations. Feedback projections from these populations to the SCN alter SCN firing patterns and fine-tune this modulation. Although this model reproduced circadian variation in sleep-wake dynamics in nocturnal rats, it failed to describe the sleep-wake dynamics observed in SCN-lesioned rats. We thus propose two alternative, physiologically based models in which neurotransmitter- and neuropeptide-mediated signalling from the SCN to sleep-wake populations introduces mechanisms to account for the behaviour of both the intact and SCN-lesioned rat. These models generate testable predictions and offer a new framework for modelling sleep-wake and circadian interactions.

Modeling circadian clocks: Roles, advantages, and limitations

Circadian rhythms are endogenous oscillations characterized by a period of about 24h. They constitute the biological rhythms with the longest period known to be generated at the molecular level. The abundance of genetic information and the complexity of the molecular circuitry make circadian clocks a system of choice for theoretical studies. Many mathematical models have been proposed to understand the molecular regulatory mechanisms that underly these circadian oscillations and to account for their dynamic properties (temperature compensation, entrainment by light dark cycles, phase shifts by light pulses, rhythm splitting, robustness to molecular noise, intercellular synchronization). The roles and advantages of modeling are discussed and illustrated using a variety of selected examples. This survey will lead to the proposal of an integrated view of the circadian system in which various aspects (interlocked feedback loops, inter-cellular coupling, and stochasticity) should be considered together to understand the design and the dynamics of circadian clocks. Some limitations of these models are commented and challenges for the future identified.

Modeling circadian clocks: From equations to oscillations

Circadian rhythms are generated at the cellular level by a small but tightly regulated genetic network. In higher eukaryotes, interlocked transcriptional-translational feedback loops form the core of this network, which ensures the activation of the right genes (proteins) at the right time of the day. Understanding how such a complex molecular network can generate robust, self-sustained oscillations and accurately responds to signals from the environment (such as light and temperature) is greatly helped by mathematical modeling. In the present paper we review some mathematical models for circadian clocks, ranging from abstract, phenomenological models to the most detailed molecular models. We explain how the equations are derived, highlighting the challenges for the modelers, and how the models are analyzed. We show how to compute bifurcation diagrams, entrainment, and phase response curves. In the subsequent paper, we discuss, through a selection of examples, how modeling efforts have contributed to a better understanding of the dynamics of the circadian regulatory network.

Calculating activation energies for temperature compensation in circadian rhythms

Many biological species possess a circadian clock, which helps them anticipate daily variations in the environment. In the absence of external stimuli, the rhythm persists autonomously with a period of approximately 24 h. However, single pulses of light, nutrients, chemicals or temperature can shift the clock phase. In the case of light- and temperature-cycles, this allows entrainment of the clock to cycles of exactly 24 h. Circadian clocks have the remarkable property of temperature compensation, that is, the period of the circadian rhythm remains relatively constant within a physiological range of temperatures. For several organisms, temperature-regulated processes within the circadian clock have been identified in recent years. However, how these processes contribute to temperature compensation is not fully understood. Here, we theoretically investigate temperature compensation in general oscillatory systems. It is known that every oscillator can be locally temperature compensated around a reference temperature, if reactions are appropriately balanced. A balancing is always possible if the control coefficient with respect to the oscillation period of at least one reaction in the oscillator network is positive. However, for global temperature compensation, the whole physiological temperature range is relevant. Here, we use an approach which leads to an optimization problem subject to the local balancing principle. We use this approach to analyse different circadian clock models proposed in the literature and calculate activation energies that lead to temperature compensation.

Modeling light adaptation in circadian clock: prediction of the response that stabilizes entrainment

Periods of biological clocks are close to but often different from the rotation period of the earth. Thus, the clocks of organisms must be adjusted to synchronize with day-night cycles. The primary signal that adjusts the clocks is light. In Neurospora, light transiently up-regulates the expression of specific clock genes. This molecular response to light is called light adaptation. Does light adaptation occur in other organisms? Using published experimental data, we first estimated the time course of the up-regulation rate of gene expression by light. Intriguingly, the estimated up-regulation rate was transient during light period in mice as well as Neurospora. Next, we constructed a computational model to consider how light adaptation had an effect on the entrainment of circadian oscillation to 24-h light-dark cycles. We found that cellular oscillations are more likely to be destabilized without light adaption especially when light intensity is very high. From the present results, we predict that the instability of circadian oscillations under 24-h light-dark cycles can be experimentally observed if light adaptation is altered. We conclude that the functional consequence of light adaptation is to increase the adjustability to 24-h light-dark cycles and then adapt to fluctuating environments in nature.

The contributions of interlocking loops and extensive nonlinearity to the properties of circadian clock models

BACKGROUND

Sensitivity and robustness are essential properties of circadian clock systems, enabling them to respond to the environment but resist noisy variations. These properties should be recapitulated in computational models of the circadian clock. Highly nonlinear kinetics and multiple loops are often incorporated into models to match experimental time-series data, but these also impact on model properties for clock models.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we study the consequences of complicated structure and nonlinearity using simple Goodwin-type oscillators and the complex Arabidopsis circadian clock models. Sensitivity analysis of the simple oscillators implies that an interlocked multi-loop structure reinforces sensitivity/robustness properties, enhancing the response to external and internal variations. Furthermore, we found that reducing the degree of nonlinearity could sometimes enhance the robustness of models, implying that ad hoc incorporation of nonlinearity could be detrimental to a model's perceived credibility.

CONCLUSION

The correct multi-loop structure and degree of nonlinearity are therefore critical in contributing to the desired properties of a model as well as its capacity to match experimental data.

Design of regulation and dynamics in simple biochemical pathways

Complex regulation of biochemical pathways in a cell is brought about by the interaction of simpler regulatory structures. Among the basic regulatory designs, feedback inhibition of gene expression is the most common motif in gene regulation and a ubiquitous control structure found in nature. In this work, we have studied a common structural feature (delayed feedback) in gene organisation and shown, both theoretically and experimentally, its subtle but important functional role in gene expression kinetics in a negatively auto-regulated system. Using simple deterministic and stochastic models with varying levels of realism, we present detailed theoretical representations of negatively auto-regulated transcriptional circuits with increasing delays in the establishment of feedback of repression. The models of the circuits with and without delay are studied analytically as well as numerically for variation of parameters and delay lengths. The positive invariance, boundedness of the solutions, local and global asymptotic stability of both the systems around the unique positive steady state are studied analytically. Existence of transient temporal dynamics is shown mathematically. Comparison of the two types of model circuits shows that even though the long-term dynamics is stable and not affected by delays in repression, there is interesting variation in the transient dynamical features with increasing delays. Theoretical predictions are validated through experimentally constructed gene circuits of similar designs. This combined theoretical and experimental study helps delineate the opposing effects of delay-induced instability, and the stability-enhancing property of negative feedback in the pathway behaviour, and gives rationale for the abundance of similar designs in real biochemical pathways.

Entrainment Concepts Revisited

The traditional approaches to predict entrainment of circadian clocks by light are based on 2 concepts that were never successfully unified: the non-parametric approach assumes that entrainment occurs via discrete daily phase shifts while the parametric approach assumes that entrainment involves changes of the clock’s velocity. Here the authors suggest a new approach to predict and model entrainment. Unlike the traditional approaches, it does not assume a priori the mechanism of how the internal and external cycle lengths are matched (via phase shifts, velocity changes, or even other mechanisms). It is based on a circadian integrated response characteristic (CIRC) that describes how the circadian system integrates light signals at different circadian phases, without specifying exactly when and how fast its progression is affected. Light around subjective dawn compresses the internal cycle; light around subjective dusk expands it. While the phase response curve (PRC) describes the results of experiments using light stimuli (of specified duration and intensity), the CIRC reflects how the system integrates any given light profile, be it single pulses or any form of light-dark cycle (from skeleton photoperiods to natural light profiles). CIRCs are characterized by their shape (determining the extent of their dead zone) and their asymmetry (the ratio of its compressing and expanding portions). They are dimensionless (time/time), and their maximum is by definition 1. To make predictions about entrainment, the CIRC is multiplied with the light intensity/sensitivity at any given time point. Unlike the PRC and the velocity response curve, the CIRC can be assessed on the basis of entrained steady states, by modeling experimental results. The CIRC approach makes several predictions that can be tested experimentally.

Quantitative analysis of cellular metabolic dissipative, self-organized structures

One of the most important goals of the postgenomic era is understanding the metabolic dynamic processes and the functional structures generated by them. Extensive studies during the last three decades have shown that the dissipative self-organization of the functional enzymatic associations, the catalytic reactions produced during the metabolite channeling, the microcompartmentalization of these metabolic processes and the emergence of dissipative networks are the fundamental elements of the dynamical organization of cell metabolism. Here we present an overview of how mathematical models can be used to address the properties of dissipative metabolic structures at different organizational levels, both for individual enzymatic associations and for enzymatic networks. Recent analyses performed with dissipative metabolic networks have shown that unicellular organisms display a singular global enzymatic structure common to all living cellular organisms, which seems to be an intrinsic property of the functional metabolism as a whole. Mathematical models firmly based on experiments and their corresponding computational approaches are needed to fully grasp the molecular mechanisms of metabolic dynamical processes. They are necessary to enable the quantitative and qualitative analysis of the cellular catalytic reactions and also to help comprehend the conditions under which the structural dynamical phenomena and biological rhythms arise. Understanding the molecular mechanisms responsible for the metabolic dissipative structures is crucial for unraveling the dynamics of cellular life.

[Circadian rhythms and systems biology]

Cellular rhythms represent a field of choice for studies in system biology. The examples of circadian rhythms and of the cell cycle show how the experimental and modeling approaches contribute to clarify the conditions in which periodic behavior spontaneously arises in regulatory networks at the cellular level. Circadian rhythms originate from intertwined positive and negative feedback loops controlling the expression of several clock genes. Models can be used to address the dynamical bases of physiological disorders related to dysfunctions of the mammalian circadian clock. The cell cycle is driven by a network of cyclin-dependent kinases (Cdks). Modeled in the form of four modules coupled through multiple regulatory interactions, the Cdk network operates in an oscillatory manner in the presence of sufficient amounts of growth factor. For circadian rhythms and the cell cycle, as for other recently observed cellular rhythms, periodic behavior represents an emergent property of biological systems related to their regulatory structure.

[Synthetic and integrative biology]

Integrative biology currently undergoes a deep renewal as we witness the increasing influence of systems biology, which explores life's logic, and of synthetic biology, which exploits it.

Features of anti-inflammatory effects of modulated extremely high-frequency electromagnetic radiation

Using a model of acute zymosan-induced paw edema in NMRI mice, we test the hypothesis that anti-inflammatory effects of extremely high-frequency electromagnetic radiation (EHF EMR) can be essentially modified by application of pulse modulation with certain frequencies. It has been revealed that a single exposure of animals to continuous EHF EMR for 20 min reduced the exudative edema of inflamed paw on average by 19% at intensities of 0.1-0.7 mW/cm(2) and frequencies from the range of 42.2-42.6 GHz. At fixed effective carrier frequency of 42.2 GHz, the anti-inflammatory effect of EHF EMR did not depend on modulation frequencies, that is, application of different modulation frequencies from the range of 0.03-100 Hz did not lead to considerable changes in the effect level. On the contrary, at "ineffective" carrier frequencies of 43.0 and 61.22 GHz, the use of modulation frequencies of 0.07-0.1 and 20-30 Hz has allowed us to restore the effect up to a maximal level. The results obtained show the critical dependence of anti-inflammatory action of low-intensity EHF EMR on carrier and modulation frequencies. Within the framework of this study, the possibility of changing the level of expected biological effect of modulated EMR by a special selection of combination of carrier and modulation frequencies is confirmed.

Circadian clocks and phosphorylation: Insights from computational modeling

Circadian clocks are based on a molecular mechanism regulated at the transcriptional, translational and post-translational levels. Recent experimental data unravel a complex role of the phosphorylations in these clocks. In mammals, several kinases play differential roles in the regulation of circadian rhythmicity. A dysfunction in the phosphorylation of one clock protein could lead to sleep disorders such as the Familial Advanced Sleep Phase Disorder, FASPS. Moreover, several drugs are targeting kinases of the circadian clocks and can be used in cancer chronotherapy or to treat mood disorders. In Drosophila, recent experimental observations also revealed a complex role of the phosphorylations. Because of its high degree of homology with mammals, the Drosophila system is of particular interest. In the circadian clock of cyanobacteria, an atypical regulatory mechanism is based only on three clock proteins (KaiA, KaiB, KaiC) and ATP and is sufficient to produce robust temperature-compensated circadian oscillations of KaiC phosphorylation. This review will show how computational modeling has become a powerful and useful tool in investigating the regulatory mechanism of circadian clocks, but also how models can give rise to testable predictions or reveal unexpected results.

Minimum criteria for DNA damage-induced phase advances in circadian rhythms

Robust oscillatory behaviors are common features of circadian and cell cycle rhythms. These cyclic processes, however, behave distinctively in terms of their periods and phases in response to external influences such as light, temperature, nutrients, etc. Nevertheless, several links have been found between these two oscillators. Cell division cycles gated by the circadian clock have been observed since the late 1950s. On the other hand, ionizing radiation (IR) treatments cause cells to undergo a DNA damage response, which leads to phase shifts (mostly advances) in circadian rhythms. Circadian gating of the cell cycle can be attributed to the cell cycle inhibitor kinase Wee1 (which is regulated by the heterodimeric circadian clock transcription factor, BMAL1/CLK), and possibly in conjunction with other cell cycle components that are known to be regulated by the circadian clock (i.e., c-Myc and cyclin D1). It has also been shown that DNA damage-induced activation of the cell cycle regulator, Chk2, leads to phosphorylation and destruction of a circadian clock component (i.e., PER1 in Mus or FRQ in Neurospora crassa). However, the molecular mechanism underlying how DNA damage causes predominantly phase advances in the circadian clock remains unknown. In order to address this question, we employ mathematical modeling to simulate different phase response curves (PRCs) from either dexamethasone (Dex) or IR treatment experiments. Dex is known to synchronize circadian rhythms in cell culture and may generate both phase advances and delays. We observe unique phase responses with minimum delays of the circadian clock upon DNA damage when two criteria are met: (1) existence of an autocatalytic positive feedback mechanism in addition to the time-delayed negative feedback loop in the clock system and (2) Chk2-dependent phosphorylation and degradation of PERs that are not bound to BMAL1/CLK.

Molecular components and mechanisms that connect cell cycle and circadian rhythms are important for the well-being of an organism. Cell cycle machinery regulates the progress of cell growth and division while the circadian rhythm network generates an ∼24 h time-keeping mechanism that regulates the daily processes of an organism (i.e. metabolism, bowel movements, body temperature, etc.). It is observed that cell divisions usually occur during a certain time window of a day, which indicated that there are circadian-gated cell divisions. Moreover, it's been shown that mice are more prone to develop cancer when certain clock genes are mutated resulting in an arrhythmic clock. Recently, a cell cycle checkpoint regulator, Chk2, was identified as a component that influences a core clock component and creates mostly phase advances (i.e., jet lags due to traveling east) in circadian rhythms upon DNA damage. This phase response with minimum delays is an unexpected result, and the molecular mechanism behind this phenomenon remains unknown. Our computational analyses of a mathematical model reveal two molecular criteria that account for the experimentally observed phase responses of the circadian clock upon DNA damage. These results demonstrate how circadian clock regulation by cell cycle checkpoint controllers provides another layer of complexity for efficient DNA damage responses.

Design principles of biochemical oscillators

Cellular rhythms are generated by complex interactions among genes, proteins and metabolites. They are used to control every aspect of cell physiology, from signalling, motility and development to growth, division and death. We consider specific examples of oscillatory processes and discuss four general requirements for biochemical oscillations: negative feedback, time delay, sufficient 'nonlinearity' of the reaction kinetics and proper balancing of the timescales of opposing chemical reactions. Positive feedback is one mechanism to delay the negative-feedback signal. Biological oscillators can be classified according to the topology of the positive- and negative-feedback loops in the underlying regulatory mechanism.