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.

Conversion of myoblasts to physiologically active neuronal phenotype.

Repressor element 1 (RE1)-silencing transcription factor (REST)/neuron-restrictive silencer factor (NRSF) can repress several terminal neuronal differentiation genes by binding to a specific DNA sequence (RE1/neuron-restrictive silencer element [NRSE]) present in their regulatory regions. REST-VP16 binds to the same RE1/NRSE, but activates these REST/NRSF target genes. However, it is unclear whether REST-VP16 expression is sufficient to cause formation of functional neurons either from neural stem cells or from heterologous stem cells. Here we show that the expression of REST-VP16 in myoblasts grown under muscle differentiation conditions blocked entry into the muscle differentiation pathway, countered endogenous REST/NRSF-dependent repression, activated the REST/NRSF target genes, and, surprisingly, activated other neuronal differentiation genes and converted the myoblasts to a physiologically active neuronal phenotype. Furthermore, in vitro differentiated neurons produced by REST-VP16-expressing myoblasts, when injected into mouse brain, survived, incorporated into the normal brain, and did not form tumors. This is the first instance in which myoblasts were converted to a neuronal phenotype. Our results suggest that direct activation of REST/NRSF target genes with a single transgene, REST-VP16, is sufficient to activate other terminal neuronal differentiation genes and to override the muscle differentiation pathways, and they suggest that this approach provides an efficient way of triggering neuronal differentiation in myoblasts and possibly other stem cells.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

ELUCIDATING FUNCTIONAL ROLES FOR MYOGENIN IN ADULT SKELETAL MUSCLE METABOLISM, EXERCISE CAPACITY, AND REGENERATION

The four basic helix-loop-helix myogenic transcription factors, myogenin, Myf5, MRF4, and MyoD are critical for embryonic skeletal muscle development. Myogenin is necessary for the terminal differentiation of myoblasts into myofibers during embryogenesis, but little is known about the roles played by myogenin in adult skeletal muscle function and metabolism. Furthermore, while metabolism is a well-studied physiological process, how it is regulated at the transcriptional level remains poorly understood. In this study, my aim was to determine the function of myogenin in adult skeletal muscle metabolism, exercise capacity, and regeneration. To investigate this, I utilized a mouse strain harboring the Myogflox allele and a Cre recombinase transgene, enabling the efficient deletion of myogenin in the adult mouse. Myogflox/flox mice were stressed physically through involuntary treadmill running and by breeding them with a strain harboring the Duchenne’s muscular dystrophy (DMDmdx) allele. Surprisingly, Myog-deleted animals exhibited an enhanced capacity for exercise, running farther and faster than their wild-type counterparts. Increased lactate production and utilization of glucose as a fuel source indicated that Myog-deleted animals exhibited an increased glycolytic flux. Hypoglycemic Myog-deleted mice no longer possessed the ability to outrun their wild-type counterparts, implying the ability of these animals to further deplete their glucose reserves confers their enhanced exercise capacity. Moreover, Myog-deleted mice exhibited an enhanced response to long-term exercise training. The mice developed a greater proportion of type 1 oxidative muscle fibers, and displayed increased levels of succinate dehydrogenase activity, indicative of increased oxidative metabolism. Mdx:Myog-deleted mice exhibited a similar phenotype, outperforming their mdx counterparts, although lagging behind wild-type animals. The morphology of muscle tissue from mdx:Myog-deleted mice appears to mimic that of mdx animals, indicating that myogenin is dispensable for adult skeletal muscle regeneration. Through global gene expression profiling and quantitative (q)RT-PCR, I identified a unique set of putative myogenin-dependent genes involved in regulating metabolic processes. These data suggest myogenin’s functions during adulthood are distinctly different than those during embryogenesis, and myogenin acts as a high-level transcription factor regulating metabolic activity in adult skeletal muscle.

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.

NF-kappaB and AP-1 connection: mechanism of NF-kappaB-dependent regulation of AP-1 activity.

Nuclear factor kappaB (NF-kappaB) and activator protein 1 (AP-1) transcription factors regulate many important biological and pathological processes. Activation of NF-kappaB is regulated by the inducible phosphorylation of NF-kappaB inhibitor IkappaB by IkappaB kinase. In contrast, Fos, a key component of AP-1, is primarily transcriptionally regulated by serum responsive factors (SRFs) and ternary complex factors (TCFs). Despite these different regulatory mechanisms, there is an intriguing possibility that NF-kappaB and AP-1 may modulate each other, thus expanding the scope of these two rapidly inducible transcription factors. To determine whether NF-kappaB activity is involved in the regulation of fos expression in response to various stimuli, we analyzed activity of AP-1 and expression of fos, fosB, fra-1, fra-2, jun, junB, and junD, as well as AP-1 downstream target gene VEGF, using MDAPanc-28 and MDAPanc-28/IkappaBalphaM pancreatic tumor cells and wild-type, IKK1-/-, and IKK2-/- murine embryonic fibroblast cells. Our results show that elk-1, a member of TCFs, is one of the NF-kappaB downstream target genes. Inhibition of NF-kappaB activity greatly decreased expression of elk-1. Consequently, the reduced level of activated Elk-1 protein by extracellular signal-regulated kinase impeded constitutive, serum-, and superoxide-inducible c-fos expression. Thus, our study revealed a distinct and essential role of NF-kappaB in participating in the regulation of elk-1, c-fos, and VEGF expression.

Abnormal expression of REST/NRSF and Myc in neural stem/progenitor cells causes cerebellar tumors by blocking neuronal differentiation.

Medulloblastoma, one of the most malignant brain tumors in children, is thought to arise from undifferentiated neural stem/progenitor cells (NSCs) present in the external granule layer of the cerebellum. However, the mechanism of tumorigenesis remains unknown for the majority of medulloblastomas. In this study, we found that many human medulloblastomas express significantly elevated levels of both myc oncogenes, regulators of neural progenitor proliferation, and REST/NRSF, a transcriptional repressor of neuronal differentiation genes. Previous studies have shown that neither c-Myc nor REST/NRSF alone could cause tumor formation. To determine whether c-Myc and REST/NRSF act together to cause medulloblastomas, we used a previously established cell line derived from external granule layer stem cells transduced with activated c-myc (NSC-M). These immortalized NSCs were able to differentiate into neurons in vitro. In contrast, when the cells were engineered to express a doxycycline-regulated REST/NRSF transgene (NSC-M-R), they no longer underwent terminal neuronal differentiation in vitro. When injected into intracranial locations in mice, the NSC-M cells did not form tumors either in the cerebellum or in the cerebral cortex. In contrast, the NSC-M-R cells did produce tumors in the cerebellum, the site of human medulloblastoma formation, but not when injected into the cerebral cortex. Furthermore, the NSC-M-R tumors were blocked from terminal neuronal differentiation. In addition, countering REST/NRSF function blocked the tumorigenic potential of NSC-M-R cells. To our knowledge, this is the first study in which abnormal expression of a sequence-specific DNA-binding transcriptional repressor has been shown to contribute directly to brain tumor formation. Our findings indicate that abnormal expression of REST/NRSF and Myc in NSCs causes cerebellum-specific tumors by blocking neuronal differentiation and thus maintaining the "stemness" of these cells. Furthermore, these results suggest that such a mechanism plays a role in the formation of human medulloblastoma.

FoxM1B transcriptionally regulates vascular endothelial growth factor expression and promotes the angiogenesis and growth of glioma cells.

We previously found that FoxM1B is overexpressed in human glioblastomas and that forced FoxM1B expression in anaplastic astrocytoma cells leads to the formation of highly angiogenic glioblastoma in nude mice. However, the molecular mechanisms by which FoxM1B enhances glioma angiogenesis are currently unknown. In this study, we found that vascular endothelial growth factor (VEGF) is a direct transcriptional target of FoxM1B. FoxM1B overexpression increased VEGF expression, whereas blockade of FoxM1 expression suppressed VEGF expression in glioma cells. Transfection of FoxM1 into glioma cells directly activated the VEGF promoter, and inhibition of FoxM1 expression by FoxM1 siRNA suppressed VEGF promoter activation. We identified two FoxM1-binding sites in the VEGF promoter that specifically bound to the FoxM1 protein. Mutation of these FoxM1-binding sites significantly attenuated VEGF promoter activity. Furthermore, FoxM1 overexpression increased and inhibition of FoxM1 expression suppressed the angiogenic ability of glioma cells. Finally, an immunohistochemical analysis of 59 human glioblastoma specimens also showed a significant correlation between FoxM1 overexpression and elevated VEGF expression. Our findings provide both clinical and mechanistic evidence that FoxM1 contributes to glioma progression by enhancing VEGF gene transcription and thus tumor angiogenesis.

Eomesodermin, a target gene of Pou4f2, is required for retinal ganglion cell and optic nerve development in the mouse.

The mechanisms regulating retinal ganglion cell (RGC) development are crucial for retinogenesis and for the establishment of normal vision. However, these mechanisms are only vaguely understood. RGCs are the first neuronal lineage to segregate from pluripotent progenitors in the developing retina. As output neurons, RGCs display developmental features very distinct from those of the other retinal cell types. To better understand RGC development, we have previously constructed a gene regulatory network featuring a hierarchical cascade of transcription factors that ultimately controls the expression of downstream effector genes. This has revealed the existence of a Pou domain transcription factor, Pou4f2, that occupies a key node in the RGC gene regulatory network and that is essential for RGC differentiation. However, little is known about the genes that connect upstream regulatory genes, such as Pou4f2 with downstream effector genes responsible for RGC differentiation. The purpose of this study was to characterize the retinal function of eomesodermin (Eomes), a T-box transcription factor with previously unsuspected roles in retinogenesis. We show that Eomes is expressed in developing RGCs and is a mediator of Pou4f2 function. Pou4f2 directly regulates Eomes expression through a cis-regulatory element within a conserved retinal enhancer. Deleting Eomes in the developing retina causes defects reminiscent of those in Pou4f2(-/-) retinas. Moreover, myelin ensheathment in the optic nerves of Eomes(-/-) embryos is severely impaired, suggesting that Eomes regulates this process. We conclude that Eomes is a crucial regulator positioned immediately downstream of Pou4f2 and is required for RGC differentiation and optic nerve development.

REST maintains self-renewal and pluripotency of embryonic stem cells.

The neuronal repressor REST (RE1-silencing transcription factor; also called NRSF) is expressed at high levels in mouse embryonic stem (ES) cells, but its role in these cells is unclear. Here we show that REST maintains self-renewal and pluripotency in mouse ES cells through suppression of the microRNA miR-21. We found that, as with known self-renewal markers, the level of REST expression is much higher in self-renewing mouse ES cells than in differentiating mouse ES (embryoid body, EB) cells. Heterozygous deletion of Rest (Rest+/-) and its short-interfering-RNA-mediated knockdown in mouse ES cells cause a loss of self-renewal-even when these cells are grown under self-renewal conditions-and lead to the expression of markers specific for multiple lineages. Conversely, exogenously added REST maintains self-renewal in mouse EB cells. Furthermore, Rest+/- mouse ES cells cultured under self-renewal conditions express substantially reduced levels of several self-renewal regulators, including Oct4 (also called Pou5f1), Nanog, Sox2 and c-Myc, and exogenously added REST in mouse EB cells maintains the self-renewal phenotypes and expression of these self-renewal regulators. We also show that in mouse ES cells, REST is bound to the gene chromatin of a set of miRNAs that potentially target self-renewal genes. Whereas mouse ES cells and mouse EB cells containing exogenously added REST express lower levels of these miRNAs, EB cells, Rest+/- ES cells and ES cells treated with short interfering RNA targeting Rest express higher levels of these miRNAs. At least one of these REST-regulated miRNAs, miR-21, specifically suppresses the self-renewal of mouse ES cells, corresponding to the decreased expression of Oct4, Nanog, Sox2 and c-Myc. Thus, REST is a newly discovered element of the interconnected regulatory network that maintains the self-renewal and pluripotency of mouse ES cells.

FoxM1B regulates NEDD4-1 expression, leading to cellular transformation and full malignant phenotype in immortalized human astrocytes.

Our recent studies have shown that the FoxM1B transcription factor is overexpressed in human glioma tissues and that the level of its expression correlates directly with glioma grade. However, whether FoxM1B plays a role in the early development of glioma (i.e., in transformation) is unknown. In this study, we found that the FoxM1B molecule causes cellular transformation and tumor formation in normal human astrocytes (NHA) immortalized by p53 and pRB inhibition. Moreover, brain tumors that arose from intracranial injection of FoxM1B-expressing immortalized NHAs displayed glioblastoma multiforme (GBM) phenotypes, suggesting that FoxM1B overexpression in immortalized NHAs not only transforms the cells but also leads to GBM formation. Mechanistically, our results showed that overexpression of FoxM1B upregulated NEDD4-1, an E3 ligase that mediates the degradation and downregulation of phosphatase and tensin homologue (PTEN) in multiple cell lines. Decreased PTEN in turn resulted in the hyperactivation of Akt, which led to phosphorylation and cytoplasmic retention of FoxO3a. Blocking Akt activation with phosphoinositide 3-kinase/Akt inhibitors inhibited the FoxM1B-induced transformation of immortalized NHAs. Furthermore, overexpression of FoxM1B in immortalized NHAs increased the expression of survivin, cyclin D1, and cyclin E, which are important molecules for tumor growth. Collectively, these results indicate that overexpression of FoxM1B, in cooperation with p53 and pRB inhibition in NHA cells, promotes astrocyte transformation and GBM formation through multiple mechanisms.

Sp1 up-regulates cAMP-response-element-binding protein expression during retinoic acid-induced mucous differentiation of normal human bronchial epithelial cells.

CREB [CRE (cAMP-response element)-binding protein] is an important transcription factor that is differentially regulated in cells of various types. We recently reported that RA (retinoic acid) rapidly activates CREB without using RARs (RA receptors) or RXRs (retinoid X receptors) in NHTBE cells (normal human tracheobronchial epithelial cells). However, little is known about the role of RA in the physiological regulation of CREB expression in the early mucous differentiation of NHTBE cells. In the present study, we report that RA up-regulates CREB gene expression and that, using 5'-serial deletion promoter analysis and mutagenesis analyses, two Sp1 (specificity protein 1)-binding sites located at nt -217 and -150, which flank the transcription initiation site, are essential for RA induction of CREB gene transcription. Furthermore, we found that CREs located at nt -119 and -98 contributed to basal promoter activity. Interestingly, RA also up-regulated Sp1 in a time- and dose-dependent manner. Knockdown of endogenous Sp1 using siRNA (small interfering RNA) decreased RA-induced CREB gene expression. However, the converse was not true: knockdown of CREB using CREB siRNA did not affect RA-induced Sp1 gene expression. We conclude that RA up-regulates CREB gene expression during the early stage of NHTBE cell differentiation and that RA-inducible Sp1 plays a major role in up-regulating human CREB gene expression. This result implies that co-operation of these two transcription factors plays a crucial role in mediating early events of normal mucous cell differentiation of bronchial epithelial cells.

Role of NonO-histone interaction in TNFalpha-suppressed prolyl-4-hydroxylase alpha1.

Inflammation is a key process in cardiovascular diseases. The extracellular matrix (ECM) of the vasculature is a major target of inflammatory cytokines, and TNFalpha regulates ECM metabolism by affecting collagen production. In this study, we have examined the pathways mediating TNFalpha-induced suppression of prolyl-4 hydroxylase alpha1 (P4Halpha1), the rate-limiting isoform of P4H responsible for procollagen hydroxylation, maturation, and organization. Using human aortic smooth muscle cells, we found that TNFalpha activated the MKK4-JNK1 pathway, which induced histone (H) 4 lysine 12 acetylation within the TNFalpha response element in the P4Halpha1 promoter. The acetylated-H4 then recruited a transcription factor, NonO, which, in turn, recruited HDACs and induced H3 lysine 9 deacetylation, thereby inhibiting transcription of the P4Halpha1 promoter. Furthermore, we found that TNFalpha oxidized DJ-1, which may be essential for the NonO-P4Halpha1 interaction because treatment with gene specific siRNA to knockout DJ-1 eliminated the TNFalpha-induced NonO-P4Halpha1 interaction and its suppression. Our findings may be relevant to aortic aneurysm and dissection and the stability of the fibrous cap of atherosclerotic plaque in which collagen metabolism is important in arterial remodeling. Defining this cytokine-mediated regulatory pathway may provide novel molecular targets for therapeutic intervention in preventing plaque rupture and acute coronary occlusion.

Rtr1 is the Saccharomyces cerevisiae homolog of a novel family of RNA polymerase II-binding proteins.

Cells must rapidly sense and respond to a wide variety of potentially cytotoxic external stressors to survive in a constantly changing environment. In a search for novel genes required for stress tolerance in Saccharomyces cerevisiae, we identified the uncharacterized open reading frame YER139C as a gene required for growth at 37 degrees C in the presence of the heat shock mimetic formamide. YER139C encodes the closest yeast homolog of the human RPAP2 protein, recently identified as a novel RNA polymerase II (RNAPII)-associated factor. Multiple lines of evidence support a role for this gene family in transcription, prompting us to rename YER139C RTR1 (regulator of transcription). The core RNAPII subunits RPB5, RPB7, and RPB9 were isolated as potent high-copy-number suppressors of the rtr1Delta temperature-sensitive growth phenotype, and deletion of the nonessential subunits RPB4 and RPB9 hypersensitized cells to RTR1 overexpression. Disruption of RTR1 resulted in mycophenolic acid sensitivity and synthetic genetic interactions with a number of genes involved in multiple phases of transcription. Consistently, rtr1Delta cells are defective in inducible transcription from the GAL1 promoter. Rtr1 constitutively shuttles between the cytoplasm and nucleus, where it physically associates with an active RNAPII transcriptional complex. Taken together, our data reveal a role for members of the RTR1/RPAP2 family as regulators of core RNAPII function.