A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

COP 9 SIGNALOSOME SUBUNIT 6 STABILIZES COP1, A NOVEL E3 UBIQUITIN LIGASE FOR 14-3-3σ

14-3-3σ, a gene upregulated by p53 in response to DNA damage, exists as part of a positive-feedback loop which activates p53 and is a human cancer epithelial marker downregulated in various cancer types. 14-3-3σ levels are critical for maintaining p53 activity in response to DNA damage and regulating signal mediator such as Akt. Here, we identify Mammalian Constitutive Photomorphogenic 1 (COP1) as a novel E3 ubiquitin ligase for targeting 14-3-3σ through proteasome degradation. We show for the first time that COP9 signalosome subunit 6 (CSN6) associates with COP1 and is involved in 14-3-3σ ubiquitin-mediated degradation. Mechanistic studies show that CSN6 expression leads to stabilization of COP1 through reducing COP1 self-ubiquitination and decelerating COP1’s turnover rate. We also show that CSN6-mediated 14-3-3σ ubiquitination is compromised when COP1 is knocked down. Thus, CSN6 mediates 14-3-3σ ubiquitination through enhancing COP1 stability. Subsequently, we show that CSN6 causes 14-3-3σ downregulation, thereby activating Akt and promoting cell survival by suppressing FOXO, an Akt target, transcriptional activity. Also, CSN6 overexpression leads to increased cell growth, transformation and promotes tumorigenicity. Significantly, 14-3-3σ expression can correct the abnormalities mediated by CSN6 expression. These data suggest that the CSN6-COP1 axis is involved in 14-3-3σ degradation, and that deregulation of this axis will promote cell growth and tumorigenicity.

Constitutive NF-kappaB and NFAT activation leads to stimulation of the BLyS survival pathway in aggressive B-cell lymphomas.

B-lymphocyte stimulator (BLyS), a relatively recently recognized member of the tumor necrosis factor ligand family (TNF), is a potent cell-survival factor expressed in many hematopoietic cells. BLyS binds to 3 TNF-R receptors, TACI, BCMA, BAFF-R, to regulate B-cell survival, differentiation, and proliferation. The mechanisms involved in BLYS gene expression and regulation are still incompletely understood. In this study, we examined BLYS gene expression, function, and regulation in B-cell non-Hodgkin lymphoma (NHL-B) cells. Our studies indicate that BLyS is constitutively expressed in aggressive NHL-B cells, including large B-cell lymphoma (LBCL) and mantle cell lymphoma (MCL), playing an important role in the survival and proliferation of malignant B cells. We found that 2 important transcription factors, NF-kappaB and NFAT, are involved in regulating BLyS expression through at least one NF-kappaB and 2 NFAT binding sites in the BLYS promoter. We also provide evidence suggesting that the constitutive activation of NF-kappaB and BLyS in NHL-B cells forms a positive feedback loop associated with lymphoma cell survival and proliferation. Our findings indicate that constitutive NF-kappaB and NFAT activations are crucial transcriptional regulators of the BLyS survival pathway in malignant B cells that could be therapeutic targets in aggressive NHL-B.

BLyS/BAFF-R signaling pathway in human aggressive non Hodgkin's lymphoma B cell and normal peripheral blood B lymphocytes

B-lymphocyte stimulator (BLyS also called BAFF), is a potent cell survival factor expressed in many hematopoietic cells. BLyS levels are elevated in the serum of non-Hodgkin lymphoma (NHL) patients, and have been reported to be associated with disease progression, and prognosis. To understand the mechanisms involved in BLyS gene expression and regulation, we examined expression, function, and regulation of the BLyS gene in B cell non-Hodgkin's lymphoma (NHL-B) cells. BLyS is constitutively expressed in aggressive NHL-B cells including large B cell lymphoma (LBCL) and mantle cell lymphoma (MCL) contributing to survival and proliferation of malignant B cells. Two important transcription factors, NF-κB and NFAT, were found to be involved in regulating BLyS expression through at least one NF-κB and two NFAT binding sites in the BLyS promoter. Further study indicates that the constitutive activation of NF-κB and BLyS in NHL-B cells forms a positive feedback loop contributing to cell survival and proliferation. In order to further investigate BLyS signaling pathway, we studied the function of BAFF-R, a major BLyS receptor, on B cells survival and proliferation. Initial study revealed that BAFF-R was also found in the nucleus, in addition to its presence on plasma membrane of B cells. Nuclear presentation of BAFF-R can be increased by anti-IgM and soluble BLyS treatment in normal peripheral B lymphocytes. Inhibition of BLyS expression decreases nuclear BAFF-R level in LBCL cells. Furthermore, we showed that BAFF-R translocated to nucleus through the classic karyopherin pathway. A candidate nuclear localization sequence (NLS) was identified in the BAFF-R protein sequence and mutation of this putative NLS can block BAFF-R entering nucleus and LBCL cell proliferation. Further study showed that BAFF-R co-localized with NF-κB family member, c-rel in the nucleus. We also found BAFF-R mediated transcriptional activity, which could be increased by c-rel. We also found that nuclear BAFF-R could bind to the NF-κB binding site on the promoters of NF-κB target genes such as BLyS, CD154, Bcl-xL, Bfl-1/A1 and IL-8. These findings indicate that BAFF-R may also promote survival and proliferation of normal B cells and NHL-B cells by directly functioning as a transcriptional co-factor with NF-κB family member. ^

TRANSCRIPTIONAL AND POST-TRANSLATIONAL MECHANISMS CONTRIBUTE TO MAINTENANCE OF REST IN NEURAL TUMORS

The RE-1 silencing transcription factor (REST) is an important regulator of normal nervous system development. It negatively regulates neuronal lineage specification in neural progenitors by binding to its consensus RE-1 element(s) located in the regulatory region of its target neuronal differentiation genes. The developmentally coordinated down-regulation of REST mRNA and protein in neural progenitors triggers terminal neurogenesis. REST is overexpressed in pediatric neural tumors such as medulloblastoma and neuroblastoma and is associated with poor neuronal differentiation. High REST protein correlate with poor prognosis for patients with medulloblastoma, however similar studies have not been done with neuroblastoma patients. Mechanism(s) underlying elevated REST levels medulloblastoma and neuroblastoma are unclear, and is the focus of this thesis project. We discovered that transcriptional and post-translational mechanisms govern REST mis-regulation in medulloblastoma and neuroblastoma. In medulloblastoma, REST transcript is aberrantly elevated in a subset of patient samples. Using loss of function and gain of function experiments, we provide evidence that the Hairy Enhancer of Split (HES1) protein represses REST transcription in medulloblastoma cell lines, modulates the expression of neuronal differentiation genes, and alters the survival potential of these cells in vitro. We also show that REST directly represses its own expression in an auto-regulatory feedback loop. Interestingly, our studies identified a novel interaction between REST and HES1. We also observed their co-occupancy at the RE-1 sites, thereby suggesting potential for co-regulation of REST expression. Our pharmacological studies in neuroblastoma using retinoic acid revealed that REST levels are controlled by transcriptional and post-transcriptional mechanisms. Post-transcriptional mechanisms are mediated by modulation of E3 ligase or REST, SCFβ-TRCP, and contribute to resistance of some cells to retinoic acid treatment.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

Interlinked dual-time feedback loops can enhance robustness to stochasticity and persistence of memory.

Multiple interlinked positive feedback loops shape the stimulus responses of various biochemical systems, such as the cell cycle or intracellular Ca2+ release. Recent studies with simplified models have identified two advantages of coupling fast and slow feedback loops. This dual-time structure enables a fast response while enhancing resistances of responses and bistability to stimulus noise. We now find that (1) the dual-time structure similarly confers resistance to internal noise due to molecule number fluctuations, and (2) model variants with altered coupling, which better represent some specific biochemical systems, share all the above advantages. We also develop a similar bistable model with coupling of a fast autoactivation loop to a slow loop. This model's topology was suggested by positive feedback proposed to play a role in long-term synaptic potentiation (LTP). The advantages of fast response and noise resistance are also present in this autoactivation model. Empirically, LTP develops resistance to reversal over approximately 1h . The model suggests this resistance may result from increased amounts of synaptic kinases involved in positive feedback.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.