Modeling the VPAC2-activated cAMP/PKA signaling pathway: from receptor to circadian clock gene induction

Vasoactive intestinal peptide excites GnRH neurons via KCa3.1, a potential player in the slow afterhyperpolarization current

Vasoactive intestinal peptide (VIP) is an important component of the suprachiasmatic nucleus (SCN) which relays circadian information to neuronal populations, including GnRH neurons. Human and animal studies have shown an impact of disrupted daily rhythms (chronic shift work, temporal food restriction, clock gene disruption) on both male and female reproduction and fertility. To date, how VIP modulates GnRH neurons remains unknown. Calcium imaging and electrophysiology on primary GnRH neurons in explants and adult mouse brain slice, respectively, were used to address this question. We found VIP excites GnRH neurons via the VIP receptor, VPAC2. The downstream signaling pathway uses both Gs protein/adenylyl cyclase/protein kinase A (PKA) and phospholipase C/phosphatidylinositol 4,5-bisphosphate (PIP2) depletion. Furthermore, we identified a UCL2077-sensitive target, likely contributing to the slow afterhyperpolarization current (IAHP), as the PKA and PIP2 depletion target, and the KCa3.1 channel as a specific target. Thus, VIP/VPAC2 provides an example of Gs protein-coupled receptor-triggered excitation in GnRH neurons, modulating GnRH neurons likely via the slow IAHP. The possible identification of KCa3.1 in the GnRH neuron slow IAHP may provide a new therapeutical target for fertility treatments.

Targeting VIP and PACAP Receptor Signaling: New Insights into Designing Drugs for the PACAP Subfamily of Receptors

Pituitary Adenylate Cyclase-Activating Peptide (PACAP) and Vasoactive Intestinal Peptide (VIP) are neuropeptides involved in a diverse array of physiological and pathological processes through activating the PACAP subfamily of class B1 G protein-coupled receptors (GPCRs): VIP receptor 1 (VPAC1R), VIP receptor 2 (VPAC2R), and PACAP type I receptor (PAC1R). VIP and PACAP share nearly 70% amino acid sequence identity, while their receptors PAC1R, VPAC1R, and VPAC2R share 60% homology in the transmembrane regions of the receptor. PACAP binds with high affinity to all three receptors, while VIP binds with high affinity to VPAC1R and VPAC2R, and has a thousand-fold lower affinity for PAC1R compared to PACAP. Due to the wide distribution of VIP and PACAP receptors in the body, potential therapeutic applications of drugs targeting these receptors, as well as expected undesired side effects, are numerous. Designing selective therapeutics targeting these receptors remains challenging due to their structural similarities. This review discusses recent discoveries on the molecular mechanisms involved in the selectivity and signaling of the PACAP subfamily of receptors, and future considerations for therapeutic targeting.

Deciphering clock cell network morphology within the biological master clock, suprachiasmatic nucleus: From the perspective of circadian wave dynamics

The biological master clock, suprachiasmatic nucleus (of rat and mouse), is composed of ~10,000 clock cells which are heterogeneous with respect to their circadian periods. Despite this inhomogeneity, an intact SCN maintains a very good degree of circadian phase (time) coherence which is vital for sustaining various circadian rhythmic activities, and it is supposedly achieved by not just one but a few different cell-to-cell coupling mechanisms, among which action potential (AP)-mediated connectivity is known to be essential. But, due to technical difficulties and limitations in experiments, so far very little information is available about the morphology of the connectivity at a cellular scale. Building upon this limited amount of information, here we exhaustively and systematically explore a large pool (~25,000) of various network morphologies to come up with some plausible network features of SCN networks. All candidates under consideration reflect an experimentally obtained 'indegree distribution' as well as a 'physical range distribution of afferent clock cells.' Then, importantly, with a set of multitude criteria based on the properties of SCN circadian phase waves in extrinsically perturbed as well as in their natural states, we select out appropriate model networks: Some important measures are, 1) level of phase dispersal and direction of wave propagation, 2) phase-resetting ability of the model networks subject to external circadian forcing, and 3) decay rate of perturbation induced "phase-singularities." The successful, realistic networks have several common features: 1) "indegree" and "outdegree" should have a positive correlation; 2) the cells in the SCN ventrolateral region (core) have a much larger total degree than that of the dorsal medial region (shell); 3) The number of intra-core edges is about 7.5 times that of intra-shell edges; and 4) the distance probability density function for the afferent connections fits well to a beta function. We believe that these newly identified network features would be a useful guide for future explorations on the very much unknown AP-mediated clock cell connectome within the SCN.

The VIP/VPAC1R Pathway Regulates Energy and Glucose Homeostasis by Modulating GLP-1, Glucagon, Leptin and PYY Levels in Mice

Simple Summary

The current study is the first complete characterization of the phenotypic, metabolic, calorimetric, and homeostatic effects of VPAC1R in a null murine model. To evaluate the role of VPAC1R on body phenotype, feeding behavior, glucose/energy homeostasis, metabolic rate and plasma hormones, a long-term study was conducted in VPAC1R^−/− and WT mice. The outcome data document that VPAC1R^−/− mice have altered metabolism and insulin intolerance, with significant increase of feeding bouts, reduction of total energy expenditure and respiratory gases during both the dark and light cycle, together with elevated fasting levels of GLP-1 and PYY, and higher postprandial levels of GLP-1, glucagon leptin and PYY. These findings suggests that VPAC1R controls glucose homeostasis and energy balance by regulating plasma metabolic hormones.

Abstract

Vasoactive Intestinal Peptide binds with high affinity to VPAC1R and VPAC2R, thus regulating key physiologic functions. Previously, we documented in VIP^−/− mice a leaner body phenotype and altered metabolic hormones. Past reports described in VPAC2^−/− mice impaired circadian rhythm, reduced food intake, and altered metabolism. To better define the effects of VPAC1R on body phenotype, energy/glucose homeostasis, and metabolism, we conducted a 12-week study in a VPAC1R null model. Our results reveal that VPAC1^−/− mice experienced significant metabolic alterations during the dark cycle with greater numbers of feeding bouts (p = 0.009), lower Total Energy Expenditure (p = 0.025), VO₂ (p = 0.029), and VCO₂ (p = 0.016); as well as during the light cycle with lower Total Energy Expenditure (p = 0.04), VO₂ (p = 0.044), and VCO₂ (p = 0.029). Furthermore, VPAC1^−/− mice had significantly higher levels of GLP-1 and PYY during fasting, and higher levels of GLP-1, glucagon leptin and PYY during postprandial conditions. In addition, VPAC1^−/− mice had lower levels of glucose at 60′ and 120′, as assessed by insulin tolerance test. In conclusion, this study supports a key role for VPAC1R in the regulation of body glucose/energy homeostasis and metabolism.

Competitive Tuning Among Ca2+/Calmodulin-Dependent Proteins: Analysis of in silico Model Robustness and Parameter Variability

Introduction

Calcium/Calmodulin-dependent (Ca²⁺/CaM-dependent) regulation of protein signaling has long been recognized for its importance in a number of physiological contexts. Found in almost all eukaryotic cells, Ca²⁺/CaM-dependent signaling participates in muscle development, immune responses, cardiac myocyte function and regulation of neuronal connectivity. In excitatory neurons, dynamic changes in the strength of synaptic connections, known as synaptic plasticity, occur when calcium ions (Ca²⁺) flux through NMDA receptors and bind the Ca²⁺-sensor calmodulin (CaM). Ca²⁺/CaM, in turn, regulates downstream protein signaling in actin polymerization, receptor trafficking, and transcription factor activation.The activation of downstream Ca²⁺/CaM-dependent binding proteins (CBPs) is a function of the frequency of Ca²⁺ flux, such that each CBP is preferentially "tuned" to different Ca²⁺ input signals. We have recently reported that competition among CBPs for CaM binding is alone sufficient to recreate in silico the observed in vivo frequency-dependence of several CBPs. However, CBP activation may strongly depend on the identity and concentration of proteins that constitute the competitive pool; with important implications in the regulation of CBPs in both normal and disease states.

Methods

Here, we extend our previous deterministic model of competition among CBPs to include phosphodiesterases, AMPAR receptors that are important in synaptic plasticity, and enzymatic function of CBPs: cAMP regulation, kinase activity, and phosphatase activity. After rigorous parameterization and validation by global sensitivity analysis using Latin Hypercube Sampling (LHS) and Partial Rank Correlation Coefficients (PRCC), we explore how perturbing the competitive pool of CBPs influences downstream signaling events. In particular, we hypothesize that although perturbations may decrease activation of one CBP, increased activation of a separate, but enzymatically-related CBP could compensate for this loss, providing a homeostatic effect.

Results and Conclusions

First we compare dynamic model output of two models: a two-state model of Ca²⁺/CaM binding and a four-state model of Ca²⁺/CaM binding. We find that a four-state model of Ca²⁺/CaM binding best captures the dynamic nature of the rapid response of CaM and CBPs to Ca²⁺ flux in the system. Using global sensitivity analysis, we find that model output is robust to parameter variability. Indeed, although variations in the expression of the CaM buffer neurogranin (Ng) may cause a decrease in Ca²⁺/CaM-dependent kinase II (CaMKII) activation, overall AMPA receptor phosphorylation is preserved; ostensibly by a concomitant increase in adenylyl cyclase 8 (AC8)-mediated activation of protein kinase A (PKA). Indeed phosphorylation of AMPAR receptors by CaMKII and PKA is robust across a wide range of Ng concentrations, though increases in AMPAR phosphorylation is seen at low Ng levels approaching zero. Our results may explain recent counter-intuitive results in neurogranin knockout mice and provide further evidence that competitive tuning is an important mechanism in synaptic plasticity. These results may be readily translated to other Ca²⁺/CaM-dependent signaling systems in other cell types and can be used to suggest targeted experimental investigation to explain counter-intuitive or unexpected downstream signaling outcomes.Tamara Kinzer-Ursem is an Assistant Professor in the Weldon School of Biomedical Engineering. She received her B.S. in Bioengineering from the University of Toledo and her M.S. and Ph.D. degrees in Chemical Engineering from the University of Michigan, and her post-doctoral training in Molecular Neuroscience at the California Institute of Technology. Prior to joining Purdue she was the Head of R&D in Biochemistry at Maven Biotechnologies and Visiting Associate in Chemical Engineering at the California Institute of Technology.Research in the Kinzer-Ursem lab focuses on developing tools to advance quantitative descriptions of cellular processes and disease within three areas of expertise: 1) Using particle diffusivity measurements to quantify biomolecular processes. Particle diffusometry is being used as a sensitive biosensor to detect the presence of pathogens in environmental and patient samples. 2) Development of novel protein tagging technologies that are used to label proteins in vivo to enable quantitative description of protein function and elucidate disease mechanisms. 3) Computational modeling of signal transduction mechanisms to understand cellular processes. Using computational techniques, we have recently described "competitive tuning" as a mechanism that might be used to regulate information transfer through protein networks, with implications in cell behavior and drug target analysis.

Circadian gating of neuronal functionality: a basis for iterative metaplasticity

Brain plasticity, the ability of the nervous system to encode experience, is a modulatory process leading to long-lasting structural and functional changes. Salient experiences induce plastic changes in neurons of the hippocampus, the basis of memory formation and recall. In the suprachiasmatic nucleus (SCN), the central circadian (~24-h) clock, experience with light at night induces changes in neuronal state, leading to circadian plasticity. The SCN's endogenous ~24-h time-generator comprises a dynamic series of functional states, which gate plastic responses. This restricts light-induced alteration in SCN state-dynamics and outputs to the nighttime. Endogenously generated circadian oscillators coordinate the cyclic states of excitability and intracellular signaling molecules that prime SCN receptivity to plasticity signals, generating nightly windows of susceptibility. We propose that this constitutes a paradigm of ~24-h iterative metaplasticity, the repeated, patterned occurrence of susceptibility to induction of neuronal plasticity. We detail effectors permissive for the cyclic susceptibility to plasticity. We consider similarities of intracellular and membrane mechanisms underlying plasticity in SCN circadian plasticity and in hippocampal long-term potentiation (LTP). The emerging prominence of the hippocampal circadian clock points to iterative metaplasticity in that tissue as well. Exploring these links holds great promise for understanding circadian shaping of synaptic plasticity, learning, and memory.

A model for the fast synchronous oscillations of firing rate in rat suprachiasmatic nucleus neurons cultured in a multielectrode array dish

When dispersed and cultured in a multielectrode dish (MED), suprachiasmatic nucleus (SCN) neurons express fast oscillations of firing rate (FOFR; fast relative to the circadian cycle), with burst duration ∼10 min, and interburst interval varying from 20 to 60 min in different cells but remaining nevertheless rather regular in individual cells. In many cases, separate neurons in distant parts of the 1 mm recording area of a MED exhibited correlated FOFR. Neither the mechanism of FOFR nor the mechanism of their synchronization among neurons is known. Based on recent data implicating vasoactive intestinal polypeptide (VIP) as a key intercellular synchronizing agent, we built a model in which VIP acts as both a feedback regulator to generate FOFR in individual neurons, and a diffusible synchronizing agent to produce coherent electrical output of a neuronal network. In our model, VIP binding to its (VPAC₂) receptors acts through G_s G-proteins to activate adenylyl cyclase (AC), increase intracellular cAMP, and open cyclic-nucleotide-gated (CNG) cation channels, thus depolarizing the cell and generating neuronal firing to release VIP. In parallel, slowly developing homologous desensitization and internalization of VPAC₂ receptors terminates elevation of cAMP and thereby provides an interpulse silent interval. Through mathematical modeling, we show that this VIP/VPAC₂/AC/cAMP/CNG-channel mechanism is sufficient for generating reliable FOFR in single neurons. When our model for FOFR is combined with a published model of synchronization of circadian rhythms based on VIP/VPAC₂ and Per gene regulation synchronization of circadian rhythms is significantly accelerated. These results suggest that (a) auto/paracrine regulation by VIP/VPAC₂ and intracellular AC/cAMP/CNG-channels are sufficient to provide robust FOFR and synchrony among neurons in a heterogeneous network, and (b) this system may also participate in synchronization of circadian rhythms.

Homologous and heterologous desensitization of guanylyl cyclase-B signaling in GH3 somatolactotropes

The guanylyl cyclases, GC-A and GC-B, are selective receptors for atrial and C-type natriuretic peptides (ANP and CNP, respectively). In the anterior pituitary, CNP and GC-B are major regulators of cGMP production in gonadotropes and yet mouse models of disrupted CNP and GC-B indicate a potential role in growth hormone secretion. In the current study, we investigate the molecular and pharmacological properties of the CNP/GC-B system in somatotrope lineage cells. Primary rat pituitary and GH3 somatolactotropes expressed functional GC-A and GC-B receptors that had similar EC50 properties in terms of cGMP production. Interestingly, GC-B signaling underwent rapid homologous desensitization in a protein phosphatase 2A (PP2A)-dependent manner. Chronic exposure to either CNP or ANP caused a significant down-regulation of both GC-A- and GC-B-dependent cGMP accumulation in a ligand-specific manner. However, this down-regulation was not accompanied by alterations in the sub-cellular localization of these receptors. Heterologous desensitization of GC-B signaling occurred in GH3 cells following exposure to either sphingosine-1-phosphate or thyrotrophin-releasing hormone (TRH). This heterologous desensitization was protein kinase C (PKC)-dependent, as pre-treatment with GF109203X prevented the effect of TRH on CNP/GC-B signaling. Collectively, these data indicate common and distinct properties of particulate guanylyl cyclase receptors in somatotropes and reveal that independent mechanisms of homologous and heterologous desensitization occur involving either PP2A or PKC. Guanylyl cyclase receptors thus represent potential novel therapeutic targets for treating growth-hormone-associated disorders.

Fast synchronous oscillations of firing rate in cultured rat suprachiasmatic nucleus neurons: possible role in circadian synchronization in the intact nucleus

The coherent circadian rhythm of the brain's master circadian pacemaker, the suprachiasmatic nucleus (SCN), is a result of synchronization of electrical activity of many SCN neurons possessing their own circadian oscillators. However, how the activity of these neurons is synchronized is not precisely known. By plotting the electrical firing rates of dispersed rat SCN neurons in multi-electrode array dishes with 20-s averaging of action-potential activity, we have investigated a novel phenomenon: fast (relative to the circadian cycle) oscillations of firing rate (FOFR) with duration of bursts ∼10min and interburst interval varying in a range from 20 to 60min in different cells, remaining nevertheless rather regular in individual cells. In many cases, separate neurons in distant parts of the 1mm recording area of an array exhibited correlated FOFR. FOFR of individual cells were positively or negatively correlated with those of other cells in a functioning neural network. Intriguingly, in occasional neuron pairs, transformation of their irregular firing to circadian peaks was accompanied by appearance of FOFR and an increase in the magnitude of firing correlation. We hypothesize that this FOFR observed in cultured SCN neurons contribute to synchronization of the circadian rhythm in the intact SCN.

Multicellular model for intercellular synchronization in circadian neural networks

We developed a multicellular model characterized by a high degree of heterogeneity to investigate possible mechanisms that underlie circadian network synchronization and rhythmicity in the suprachiasmatic nucleus (SCN). We populated a two-dimensional grid with 400 model neurons coupled via γ-aminobutyric acid (GABA) and vasoactive intestinal polypeptide (VIP) neurotransmitters through a putative Ca(2+) mediated signaling cascade to investigate their roles in gene expression and electrical firing activity of cell populations. As observed experimentally, our model predicted that GABA would affect the amplitude of circadian oscillations but not synchrony among individual oscillators. Our model recapitulated experimental findings of decreased synchrony and average periods, loss of rhythmicity, and reduced circadian amplitudes as VIP signaling was eliminated. In addition, simulated increases of VIP reduced periodicity and synchrony. We therefore postulated a physiological range of VIP within which the system is able to produce sustained and synchronized oscillations. Our model recapitulated experimental findings of diminished amplitudes and periodicity with decreasing intracellular Ca(2+) concentrations, suggesting that such behavior could be due to simultaneous decrease of individual oscillation amplitudes and population synchrony. Simulated increases in Cl(-) levels resulted in increased Cl(-) influx into the cytosol, a decrease of inhibitory postsynaptic currents, and ultimately a shift of GABA-elicited responses from inhibitory to excitatory. The simultaneous reduction of IPSCs and increase in membrane resting potential produced GABA dose-dependent increases in firing rates across the population, as has been observed experimentally. By integrating circadian gene regulation and electrophysiology with intracellular and intercellular signaling, we were able to develop the first (to our knowledge) multicellular model that allows the effects of clock genes, electrical firing, Ca(2+), GABA, and VIP on circadian system behavior to be predicted.

Multiscale complexity in the mammalian circadian clock

The field of systems biology studies how the interactions among individual components (e.g. genes and proteins) yield interesting and complex behavior. The circadian (daily) timekeeping system in mammals is an ideal system to study complexity because of its many biological scales (from genes to animal behavior). A wealth of data at each of these scales has recently been discovered. Within each scale, modeling can advance our understanding of challenging problems that arise in studying mammalian timekeeping. However, future work must focus on bridging the multiple spatial and temporal scales in the modeling of SCN network. Here we review recent advances, and then delve into a few areas that are promising research directions. We also discuss the flavor of modeling needed (simple or detailed) as well as new techniques that are needed to meet the challenges in modeling data across scales.

Robust dynamic balance of AP-1 transcription factors in a neuronal gene regulatory network

BACKGROUND

The octapeptide Angiotensin II is a key hormone that acts via its receptor AT1R in the brainstem to modulate the blood pressure control circuits and thus plays a central role in the cardiac and respiratory homeostasis. This modulation occurs via activation of a complex network of signaling proteins and transcription factors, leading to changes in levels of key genes and proteins. AT1R initiated activity in the nucleus tractus solitarius (NTS), which regulates blood pressure, has been the subject of extensive molecular analysis. But the adaptive network interactions in the NTS response to AT1R, plausibly related to the development of hypertension, are not understood.

RESULTS

We developed and analyzed a mathematical model of AT1R-activated signaling kinases and a downstream gene regulatory network, with structural basis in our transcriptomic data analysis and literature. To our knowledge, our report presents the first computational model of this key regulatory network. Our simulations and analysis reveal a dynamic balance among distinct dimers of the AP-1 family of transcription factors. We investigated the robustness of this behavior to simultaneous perturbations in the network parameters using a novel multivariate approach that integrates global sensitivity analysis with decision-tree methods. Our analysis implicates a subset of Fos and Jun dependent mechanisms, with dynamic sensitivities shifting from Fos-regulating kinase (FRK)-mediated processes to those downstream of c-Jun N-terminal kinase (JNK). Decision-tree analysis indicated that while there may be a large combinatorial functional space feasible for neuronal states and parameters, the network behavior is constrained to a small set of AP-1 response profiles. Many of the paths through the combinatorial parameter space lead to a dynamic balance of AP-1 dimer forms, yielding a robust AP-1 response counteracting the biological variability.

CONCLUSIONS

Based on the simulation and analysis results, we demonstrate that a dynamic balance among distinct dimers of the AP-1 family of transcription factors underlies the robust activation of neuronal gene expression in the NTS response to AT1R activation. Such a differential sensitivity to limited set of mechanisms is likely to underlie the stable homeostatic physiological response.

Emergence of noise-induced oscillations in the central circadian pacemaker

Bmal1 is an essential transcriptional activator within the mammalian circadian clock. We report here that the suprachiasmatic nucleus (SCN) of Bmal1-null mutant mice, unexpectedly, generates stochastic oscillations with periods that overlap the circadian range. Dissociated SCN neurons expressed fluctuating levels of PER2 detected by bioluminescence imaging but could not generate circadian oscillations intrinsically. Inhibition of intercellular communication or cyclic-AMP signaling in SCN slices, which provide a positive feed-forward signal to drive the intracellular negative feedback loop, abolished the stochastic oscillations. Propagation of this feed-forward signal between SCN neurons then promotes quasi-circadian oscillations that arise as an emergent property of the SCN network. Experimental analysis and mathematical modeling argue that both intercellular coupling and molecular noise are required for the stochastic rhythms, providing a novel biological example of noise-induced oscillations. The emergence of stochastic circadian oscillations from the SCN network in the absence of cell-autonomous circadian oscillatory function highlights a previously unrecognized level of circadian organization.

Fathead minnow steroidogenesis: in silico analyses reveals tradeoffs between nominal target efficacy and robustness to cross-talk

Background

Interpreting proteomic and genomic data is a major challenge in predictive ecotoxicology that can be addressed by a systems biology approach. Mathematical modeling provides an organizational platform to consolidate protein dynamics with possible genomic regulation. Here, a model of ovarian steroidogenesis in the fathead minnow, Pimephales promelas, (FHM) is developed to evaluate possible transcriptional regulation of steroid production observed in microarray studies.

Results

The model was developed from literature sources, integrating key signaling components (G-protein and PKA activation) with their ensuing effect on steroid production. The model properly predicted trajectory behavior of estradiol and testosterone when fish were exposed to fadrozole, a specific aromatase inhibitor, but failed to predict the steroid hormone behavior occurring one week post-exposure as well as the increase in steroid levels when the stressor was removed. In vivo microarray data implicated three modes of regulation which may account for over-production of steroids during a depuration phase (when the stressor is removed): P450 enzyme up-regulation, inhibin down-regulation, and luteinizing hormone receptor up-regulation. Simulation studies and sensitivity analysis were used to evaluate each case as possible source of compensation to endocrine stress.

Conclusions

Simulation studies of the testosterone and estradiol response to regulation observed in microarray data supported the hypothesis that the FHM steroidogenesis network compensated for endocrine stress by modulating the sensitivity of the ovarian network to global cues coming from the hypothalamus and pituitary. Model predictions of luteinizing hormone receptor regulation were consistent with depuration and in vitro data. These results challenge the traditional approach to network elucidation in systems biology. Generally, the most sensitive interactions in a network are targeted for further elucidation but microarray evidence shows that homeostatic regulation of the steroidogenic network is likely maintained by a mildly sensitive interaction. We hypothesize that effective network elucidation must consider both the sensitivity of the target as well as the target's robustness to biological noise (in this case, to cross-talk) when identifying possible points of regulation.

Light-dark cycle memory in the mammalian suprachiasmatic nucleus

The mammalian circadian oscillator, or suprachiasmatic nucleus (SCN), contains several thousand clock neurons in its ventrolateral division, many of which are spontaneous oscillators with period lengths that range from 22 to 28 h. In complete darkness, this network synchronizes through the exchange of action potentials that release vasoactive intestinal polypeptide, striking a compromise, free-running period close to 24 h long. We entrained Siberian hamsters to various light-dark cycles and then tracked their activity into constant darkness to show that they retain a memory of the previous light-dark cycle before returning to their own free-running period. Employing Leloup-Goldbeter mammalian clock neurons we model the ventrolateral SCN network and show that light acting weakly upon a strongly rhythmic vasoactive intestinal polypeptide oscillation can explain the observed light-dark cycle memory. In addition, light is known to initiate a mitogen-activated protein kinase signaling cascade that induces transcription of both per and mkp1 phosphatase. We show that the ensuing phosphatase-kinase interaction can account for the dead zone in the mammalian phase response curve and hypothesize that the SCN behaves like a lock-in amplifier to entrain to the light edges of the circadian day.

Synchrony and entrainment properties of robust circadian oscillators

Systems theoretic tools (i.e. mathematical modelling, control, and feedback design) advance the understanding of robust performance in complex biological networks. We highlight phase entrainment as a key performance measure used to investigate dynamics of a single deterministic circadian oscillator for the purpose of generating insight into the behaviour of a population of (synchronized) oscillators. More specifically, the analysis of phase characteristics may facilitate the identification of appropriate coupling mechanisms for the ensemble of noisy (stochastic) circadian clocks. Phase also serves as a critical control objective to correct mismatch between the biological clock and its environment. Thus, we introduce methods of investigating synchrony and entrainment in both stochastic and deterministic frameworks, and as a property of a single oscillator or population of coupled oscillators.

A molecular model for intercellular synchronization in the mammalian circadian clock

The mechanisms and consequences of synchrony among heterogeneous oscillators are poorly understood in biological systems. We present a multicellular, molecular model of the mammalian circadian clock that incorporates recent data implicating the neurotransmitter vasoactive intestinal polypeptide (VIP) as the key synchronizing agent. The model postulates that synchrony arises among circadian neurons because they release VIP rhythmically on a daily basis and in response to ambient light. Two basic cell types, intrinsically rhythmic pacemakers and damped oscillators, are assumed to arise from a distribution of Period gene transcription rates. Postsynaptic neurons show time-of-day dependent responses to VIP binding through a signaling cascade that activates Period mRNA transcription. The heterogeneous cell ensemble model self-synchronizes, entrains to ambient light-dark cycles, and desynchronizes in constant bright light or upon removal of VIP signaling. The degree of synchronicity observed depends on cell-specific features (e.g., mean and variability of parameters within the rhythm-generating loop), in addition to the more commonly studied effect of intercellular coupling strength. These simulations closely replicate experimental data and predict that heterogeneous oscillations (e.g., sustained, damped, and arrhythmic) arise from small differences in the molecular parameters between cells, that damped oscillators participate in entrainment and synchrony of the ensemble of cells, and that constant light desynchronizes oscillators by maximizing VIP release.

GABA and Gi/o differentially control circadian rhythms and synchrony in clock neurons

Neurons in the mammalian suprachiasmatic nuclei (SCN) generate daily rhythms in physiology and behavior, but it is unclear how they maintain and synchronize these rhythms in vivo. We hypothesized that parallel signaling pathways in the SCN are required to synchronize rhythms in these neurons for coherent output. We recorded firing and clock-gene expression patterns while blocking candidate signaling pathways for at least 8 days. GABA(A) and GABA(B) antagonism increased circadian peak firing rates and rhythm precision of cultured SCN neurons, but G(i/o) did not impair synchrony or rhythmicity. In contrast, inhibiting G(i/o) with pertussis toxin abolished rhythms in most neurons and desynchronized the population, phenocopying the loss of vasoactive intestinal polypeptide (VIP). Daily VIP receptor agonist treatment restored synchrony and rhythmicity to VIP(-/-) SCN cultures during continuous GABA receptor antagonism but not during G(i/o) blockade. Pertussis toxin did not affect circadian cycling of the liver, suggesting that G(i/o) plays a specialized role in maintaining SCN rhythmicity. We conclude that endogenous GABA controls the amplitude of SCN neuronal rhythms by reducing daytime firing, whereas G(i/o) signaling suppresses nighttime firing, and it is necessary for synchrony among SCN neurons. We propose that G(i/o), not GABA activity, converges with VIP signaling to maintain and coordinate rhythms among SCN neurons.