A mathematical model for the role of dopamine-D2 self-regulation in the production of ultradian rhythms.

Many self-motivated and goal-directed behaviours display highly flexible, approximately 4 hour ultradian (shorter than a day) oscillations. Despite lacking direct correspondence to physical cycles in the environment, these ultradian rhythms may be involved in optimizing functional interactions with the environment and reflect intrinsic neural dynamics. Current evidence supports a role of mesostriatal dopamine (DA) in the expression and propagation of ultradian rhythmicity, however, the biochemical processes underpinning these oscillations remain to be identified. Here, we use a mathematical model to investigate D2 autoreceptor-dependent DA self-regulation as the source of ultradian behavioural rhythms. DA concentration at the midbrain-striatal synapses is governed through a dual-negative feedback-loop structure, which naturally gives rise to rhythmicity. This model shows the propensity of striatal DA to produce an ultradian oscillation characterized by a flexible period that is highly sensitive to parameter variations. Circadian (approximately 24 hour) regulation consolidates the ultradian oscillations and alters their response to the phase-dependent, rapid-resetting effect of a transient excitatory stimulus. Within a circadian framework, the ultradian rhythm orchestrates behavioural activity and enhances responsiveness to an external stimulus. This suggests a role for the circadian-ultradian timekeeping hierarchy in governing organized behaviour and shaping daily experience through coordinating the motivation to engage in recurring, albeit not highly predictable events, such as social interactions.

Kinetic modelling of ß-cell metabolism reveals control points in the insulin-regulating pyruvate cycling pathways.

Insulin, a key hormone in the regulation of glucose homoeostasis, is secreted by pancreatic ß-cells in response to elevated glucose levels. Insulin is released in a biphasic manner in response to glucose metabolism in ß-cells. The first phase of insulin secretion is triggered by an increase in the ATP:ADP ratio; the second phase occurs in response to both a rise in ATP:ADP and other key metabolic signals, including a rise in the NADPH:NADP+ ratio. Experimental evidence indicates that pyruvate-cycling pathways play an important role in the elevation of the NADPH:NADP+ ratio in response to glucose. The authors developed a kinetic model for the tricarboxylic acid cycle and pyruvate cycling pathways. The authors successfully validated the model against experimental observations and performed a sensitivity analysis to identify key regulatory interactions in the system. The model predicts that the dicarboxylate carrier and the pyruvate transporter are the most important regulators of pyruvate cycling and NADPH production. In contrast, the analysis showed that variation in the pyruvate carboxylase flux was compensated by a response in the activity of mitochondrial isocitrate dehydrogenase (ICDm ) resulting in minimal effect on overall pyruvate cycling flux. The model predictions suggest starting points for further experimental investigation, as well as potential drug targets for the treatment of type 2 diabetes.

A multi-scale simulation of retinal physiology.

We present a detailed physiological model of the (human) retina that includes the biochemistry and electrophysiology of phototransduction, neuronal electrical coupling, and the spherical geometry of the eye. The model is a parabolic-elliptic system of partial differential equations based on the mathematical framework of the bi-domain equations, which we have generalized to account for multiple cell-types. We discretize in space with non-uniform finite differences and step through time with a custom adaptive time-stepper that employs a backward differentiation formula and an inexact Newton method. A refinement study confirms the accuracy and efficiency of our numerical method. Numerical simulations using the model compare favorably with experimental findings, such as desensitization to light stimuli and calcium buffering in photoreceptors. Other numerical simulations suggest an interplay between photoreceptor gap junctions and inner segment, but not outer segment, calcium concentration. Applications of this model and simulation include analysis of retinal calcium imaging experiments, the design of electroretinograms, the design of visual prosthetics, and studies of ephaptic coupling within the retina.

Multiple entrained oscillator model of food anticipatory circadian rhythms.

For many animal species, knowing when to look for food may be as important as knowing where to look. Rats and other species use a feeding-responsive circadian timing mechanism to anticipate, behaviorally and physiologically, a predictable daily feeding opportunity. How this mechanism for anticipating a daily meal accommodates more than one predictable mealtime is unclear. Rats were trained to press a lever for food, and then limited to one or more daily meals at fixed or systematically varying times of day. The rats were able to anticipate up to 4 of 4 daily meals at fixed times of day and two 'daily' meals recurring at 24 h and 26 h intervals. When deprived of food, in constant dark, lever pressing recurred for multiple cycles at expected mealtimes, consistent with the periodicity of the prior feeding schedule. Anticipation did not require the suprachiasmatic nucleus circadian pacemaker. The anticipation rhythms could be simulated using a Kuramoto model in which clusters of coupled oscillators entrain to specific mealtimes based on initial phase and intrinsic circadian periodicity. A flexibly coupled system of food-entrainable circadian oscillators endows rats with adaptive plasticity in daily programming of foraging activity.

A Phenomenological Mouse Circadian Pacemaker Model.

Mathematical models have been used extensively in chronobiology to explore characteristics of biological clocks. In particular, for human circadian studies, the Kronauer model has been modified multiple times to describe rhythm production and responses to sensory input. This phenomenological model comprises a single set of parameters which can simulate circadian responses in humans under a variety of environmental conditions. However, corresponding models for nocturnal rodents commonly used in circadian rhythm studies are not available and may require new parameter values for different species and even strains. Moreover, due to a considerable variation in experimental data collected from mice of the same strain, within and across laboratories, a range of valid parameters is essential. This study develops a Kronauer-like model for mice by re-fitting relevant parameters to published phase response curve and period data using total least squares. Local parameter sensitivity analysis and parameter distributions determine the parameter ranges that give a near-identical model and data distribution of periods. However, the model required further parameter adjustments to match characteristics of other mouse strains, implying that the model itself detects changes in the core processes of rhythm generation and control. The model is a useful tool to understand and interpret future mouse circadian clock experiments.

Promoter methylation in a mixed feedback loop circadian clock model.

We investigate how epigenetic modifications to clock gene promoters affect transcriptomic activity in the circadian clock. Motivated by experimental observations that link DNA methylation with the behavior of the clock, we introduce and analyze an extension of the mixed feedback loop (MFL) model of François and Hakim. We extend the original model to include an additional methylated promoter state and allow for reversible protein sequestration, an important feature for circadian applications. First, working with the general form of the MFL model, we find that the qualitative behavior of the model is dictated by the promoter state with the highest transcription rate. We then build on the work of Kim and Forger, who analyzed the stability of the mammalian circadian clock by using a reduced form of the MFL model. We present a rigorous procedure for translating between the MFL model and the reduction of Kim and Forger. We then propose a model reduction more appropriate for the study of oscillatory promoter states, making use of a fully coupled quasi-steady-state approximation rather than the standard partially uncoupled quasi-steady-state approach. Working with the novel reduced form of the model, we find substantial differences in the transcription function and show that, although methylation contributes to period control, excessive methylation can abolish rhythmicity. Altogether our results show that even in a minimal clock model, DNA methylation has a nontrivial influence on the system's ability to oscillate.

M1-Type, but Not M4-Type, Melanopsin Ganglion Cells Are Physiologically Tuned to the Central Circadian Clock.

Proper circadian photoentrainment is crucial for the survival of many organisms. In mammals, intrinsically photosensitive retinal ganglion cells (ipRGCs) can use the photopigment melanopsin to sense light independently from rod and cone photoreceptors and send this information to many brain nuclei such as the suprachiasmatic nucleus (SCN), the site of the central circadian pacemaker. Here, we measure ionic currents and develop mathematical models of the electrical activity of two types of ipRGCs: M1, which projects to the SCN, and M4, which does not. We illustrate how their ionic properties differ, mainly how ionic currents generate lower spike rates and depolarization block in M1 ipRGCs. Both M1 and M4 cells have large geometries and project to higher visual centers of the brain via the optic nerve. Using a partial differential equation model, we show how axons of M1 and M4 cells faithfully convey information from the soma to the synapse even when the signal at the soma is attenuated due to depolarization block. Finally, we consider an ionic model of circadian photoentrainment from ipRGCs synapsing on SCN neurons and show how the properties of M1 ipRGCs are tuned to create accurate transmission of visual signals from the retina to the central pacemaker, whereas M4 ipRGCs would not evoke nearly as efficient a postsynaptic response. This work shows how ipRGCs and SCN neurons' electrical activities are tuned to allow for accurate circadian photoentrainment.

Targeted modification of the Per2 clock gene alters circadian function in mPer2luciferase (mPer2Luc) mice.

Modification of the Per2 clock gene in mPer2Luc reporter mice significantly alters circadian function. Behavioral period in constant dark is lengthened, and dissociates into two distinct components in constant light. Rhythms exhibit increased bimodality, enhanced phase resetting to light pulses, and altered entrainment to scheduled feeding. Mechanistic mathematical modelling predicts that enhanced protein interactions with the modified mPER2 C-terminus, combined with differential clock regulation among SCN subregions, can account for effects on circadian behavior via increased Per2 transcript and protein stability. PER2::LUC produces greater suppression of CLOCK:BMAL1 E-box activity than PER2. mPer2Luc carries a 72 bp deletion in exon 23 of Per2, and retains a neomycin resistance cassette that affects rhythm amplitude but not period. The results show that mPer2Luc acts as a circadian clock mutation illustrating a need for detailed assessment of potential impacts of c-terminal tags in genetically modified animal models.

Multiplexing Visual Signals in the Suprachiasmatic Nuclei.

The suprachiasmatic nuclei (SCN), the site of the mammalian circadian (daily) pacemaker, contains thousands of interconnected neurons, some of which receive direct retinal input. Here, we study the fast (<1 s) responses of SCN neurons to visual stimuli with a large-scale mathematical model tracking the ionic currents and voltage of all SCN neurons. We reconstruct the SCN network connectivity and reject 99.99% of theoretically possible SCN networks by requiring that the model reproduces experimentally determined receptive fields of SCN neurons. The model shows how the SCN neuronal network can enhance circadian entrainment by sensitizing a population of neurons in the ventral SCN to irradiance. This SCN network also increases the spatial acuity of neurons and increases the accuracy of a simulated subconscious spatial visual task. We hypothesize that much of the fast electrical activity within the SCN is related to the processing of spatial information.

Responses to Spatial Contrast in the Mouse Suprachiasmatic Nuclei.

A direct retinal projection targets the suprachiasmatic nucleus (SCN) (an important hypothalamic control center). The accepted function of this projection is to convey information about ambient light (irradiance) to synchronize the SCN's endogenous circadian clock with local time and drive the diurnal variations in physiology and behavior [1-4]. Here, we report that it also renders the SCN responsive to visual images. We map spatial receptive fields (RFs) for SCN neurons and find that only a minority are excited (or inhibited) by light from across the scene as expected for irradiance detectors. The most commonly encountered units have RFs with small excitatory centers, combined with very extensive inhibitory surrounds that reduce their sensitivity to global changes in light in favor of responses to spatial patterns. Other units have larger excitatory RF centers, but these always cover a coherent region of visual space, implying visuotopic order at the single-unit level. Approximately 75% of light-responsive SCN units modulate their firing according to simple spatial patterns (drifting or inverting gratings) without changes in irradiance. The time-averaged firing rate of the SCN is modestly increased under these conditions, but including spatial contrast did not significantly alter the circadian phase resetting efficiency of light. Our data indicate that the SCN contains information about irradiance and spatial patterns. This newly appreciated sensory capacity provides a mechanism by which behavioral and physiological systems downstream of the SCN could respond to visual images [5].

It is not the parts, but how they interact that determines the behaviour of circadian rhythms across scales and organisms.

Biological rhythms, generated by feedback loops containing interacting genes, proteins and/or cells, time physiological processes in many organisms. While many of the components of the systems that generate biological rhythms have been identified, much less is known about the details of their interactions. Using examples from the circadian (daily) clock in three organisms, Neurospora, Drosophila and mouse, we show, with mathematical models of varying complexity, how interactions among (i) promoter sites, (ii) proteins forming complexes, and (iii) cells can have a drastic effect on timekeeping. Inspired by the identification of many transcription factors, for example as involved in the Neurospora circadian clock, that can both activate and repress, we show how these multiple actions can cause complex oscillatory patterns in a transcription-translation feedback loop (TTFL). Inspired by the timekeeping complex formed by the NMO-PER-TIM-SGG complex that regulates the negative TTFL in the Drosophila circadian clock, we show how the mechanism of complex formation can determine the prevalence of oscillations in a TTFL. Finally, we note that most mathematical models of intracellular clocks model a single cell, but compare with experimental data from collections of cells. We find that refitting the most detailed model of the mammalian circadian clock, so that the coupling between cells matches experimental data, yields different dynamics and makes an interesting prediction that also matches experimental data: individual cells are bistable, and network coupling removes this bistability and causes the network to be more robust to external perturbations. Taken together, we propose that the interactions between components in biological timekeeping systems are carefully tuned towards proper function. We also show how timekeeping can be controlled by novel mechanisms at different levels of organization.

Population density approach for discrete mRNA distributions in generalized switching models for stochastic gene expression.

We present a generalization of a population density approach for modeling and analysis of stochastic gene expression. In the model, the gene of interest fluctuates stochastically between an inactive state, in which transcription cannot occur, and an active state, in which discrete transcription events occur; and the individual mRNA molecules are degraded stochastically in an independent manner. This sort of model in simplest form with exponential dwell times has been used to explain experimental estimates of the discrete distribution of random mRNA copy number. In our generalization, the random dwell times in the inactive and active states, T_{0} and T_{1}, respectively, are independent random variables drawn from any specified distributions. Consequently, the probability per unit time of switching out of a state depends on the time since entering that state. Our method exploits a connection between the fully discrete random process and a related continuous process. We present numerical methods for computing steady-state mRNA distributions and an analytical derivation of the mRNA autocovariance function. We find that empirical estimates of the steady-state mRNA probability mass function from Monte Carlo simulations of laboratory data do not allow one to distinguish between underlying models with exponential and nonexponential dwell times in some relevant parameter regimes. However, in these parameter regimes and where the autocovariance function has negative lobes, the autocovariance function disambiguates the two types of models. Our results strongly suggest that temporal data beyond the autocovariance function is required in general to characterize gene switching.