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.

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.

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.

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.

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.

A model of late long-term potentiation simulates aspects of memory maintenance.

Late long-term potentiation (L-LTP) denotes long-lasting strengthening of synapses between neurons. L-LTP appears essential for the formation of long-term memory, with memories at least partly encoded by patterns of strengthened synapses. How memories are preserved for months or years, despite molecular turnover, is not well understood. Ongoing recurrent neuronal activity, during memory recall or during sleep, has been hypothesized to preferentially potentiate strong synapses, preserving memories. This hypothesis has not been evaluated in the context of a mathematical model representing ongoing activity and biochemical pathways important for L-LTP. In this study, ongoing activity was incorporated into two such models - a reduced model that represents some of the essential biochemical processes, and a more detailed published model. The reduced model represents synaptic tagging and gene induction simply and intuitively, and the detailed model adds activation of essential kinases by Ca(2+). Ongoing activity was modeled as continual brief elevations of Ca(2+). In each model, two stable states of synaptic strength/weight resulted. Positive feedback between synaptic weight and the amplitude of ongoing Ca(2+) transients underlies this bistability. A tetanic or theta-burst stimulus switches a model synapse from a low basal weight to a high weight that is stabilized by ongoing activity. Bistability was robust to parameter variations in both models. Simulations illustrated that prolonged periods of decreased activity reset synaptic strengths to low values, suggesting a plausible forgetting mechanism. However, episodic activity with shorter inactive intervals maintained strong synapses. Both models support experimental predictions. Tests of these predictions are expected to further understanding of how neuronal activity is coupled to maintenance of synaptic strength. Further investigations that examine the dynamics of activity and synaptic maintenance can be expected to help in understanding how memories are preserved for up to a lifetime in animals including humans.

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. ^

DETERMINING THE GENOTYPE-PHENOTYPE CONNECTION IN SYNTHETIC INDUCIBLE GENE EXPRESSION SYSTEMS

Introduction Gene expression is an important process whereby the genotype controls an individual cell’s phenotype. However, even genetically identical cells display a variety of phenotypes, which may be attributed to differences in their environment. Yet, even after controlling for these two factors, individual phenotypes still diverge due to noisy gene expression. Synthetic gene expression systems allow investigators to isolate, control, and measure the effects of noise on cell phenotypes. I used mathematical and computational methods to design, study, and predict the behavior of synthetic gene expression systems in S. cerevisiae, which were affected by noise. Methods I created probabilistic biochemical reaction models from known behaviors of the tetR and rtTA genes, gene products, and their gene architectures. I then simplified these models to account for essential behaviors of gene expression systems. Finally, I used these models to predict behaviors of modified gene expression systems, which were experimentally verified. Results Cell growth, which is often ignored when formulating chemical kinetics models, was essential for understanding gene expression behavior. Models incorporating growth effects were used to explain unexpected reductions in gene expression noise, design a set of gene expression systems with “linear” dose-responses, and quantify the speed with which cells explored their fitness landscapes due to noisy gene expression. Conclusions Models incorporating noisy gene expression and cell division were necessary to design, understand, and predict the behaviors of synthetic gene expression systems. The methods and models developed here will allow investigators to more efficiently design new gene expression systems, and infer gene expression properties of TetR based systems.

Dysregulated CD40-NF-kappaB signaling in the pathophysiology of aggressive non -Hodgkin's lymphoma B cells

Non-Hodgkin's Lymphomas (NHL) are a group (>30) of important human lymphoid cancers that unlike other tumors today, are showing a marked increase in incidence. The lack of insight to the pathogenesis of B-cell NHL poses a significant problem in the early detection and effective treatment of these malignancies. This study shows that large B-cell lymphoma (LBCL) cells, the most common type of B-cell NHL (account for more than 30% of cases), have developed a novel mechanism for autonomous neoplastic B cell growth. We have identified that the key transcription factor NF-κB, is constitutively activated in LBCL cell lines and primary biopsy-derived LBCL cells, suggesting that they are autonomously activated, and do not require accessory T-cell signaling for cell growth and survival. Further studies have indicated that LBCL cells ectopically express an important T-cell associated co-mitogenic factor, CD154 (CD40 ligand), that is able to internally activate the CD401NF-κB pathway, through constitutive binding to its cognate receptor, CD40, on the lymphoma cell surface. CD40 activation triggers the formation of a “Signalosome” comprising virtually the entire canonical CD40/NF-κB signaling pathway that is anchored by CD40 in plasma membrane lipid rafts. The CD40 Signalosome is vulnerable to interdiction by antibody against CD40 that disrupts the Signalosome and induces cell death in the malignant cells. In addition to constitutive NF-κB activation, we have found that the nuclear factor of activated T cells (NFAT) transcription factor is also constitutively activated in LBCL cells. We have demonstrated that the constitutively active NFATc1 and c-rel members of the NFAT and NF-κB families of transcription factors, respectively, interact with each other, bind to the CD154 promoter, and synergistically activate CD154 gene transcription. Down-regulation of NFATc1 and c-rel with small interfering RNA inhibits CD154 gene transcription and lymphoma cell growth. Our findings suggest that continuous CD40 activation not only provides dysregulated proliferative stimuli for lymphoma cell growth and extended tumor cell survival, but also allows continuous regeneration of the CD40 ligand in the lymphoma cell and thereby recharges the system through a positive feedback mechanism. Targeting the CD40/NF-κB signaling pathway could provide potential therapeutic modalities for LBCL cells in the future. ^