A new fork for clinical application: Targeting forkhead transcription factors in cancer therapy

RAS-ERK-MAPK (Mitogen-activated protein kinase) pathway plays an essential role in proliferation, differentiation, and tumor progression. In this study, we showed that ERK downregulated FOXO3a through directly interacting with and phosphorylating FOXO3a at Serine 294, Serine 344, and Serine 425. ERK-phosphorylated FOXO3a was degraded by MDM2-mediated ubiquitin-proteosome pathway. FOXO3a phosphorylation and degradation consequently promoted cell proliferation and tumorigenesis. However, the non-phosphorylated FOXO3a mutant, which was resistant to the interaction and degradation by MDM2, resulted in inhibition of tumor formation. Forkhead O transcription factors (FOXOs) are important in the regulation of cellular functions including cell cycle arrest and cell death. Perturbation of FOXOs function leads to deregulated cell proliferation and cancer. Inactivation of FOXO proteins by activation of cell survival pathways, such as PI3K/AKT/IKK, is associated with tumorigenesis. Our study will further highlight FOXOs as new therapeutic targets in a broad spectrum of cancers. ^ Chemotherapeutic drug resistance is the most concerned problem in cancer therapy as resistance ultimately leads to treatment failure of cancer patients. In another study, we showed that blocking ERK activity with AZD6244, an established MEK1/2 inhibitor currently under human cancer clinical trials, enhances FOXO3a expression in various human cancer cell lines in vitro, and also in human colon cancer cell xenografts in vivo. Knocking down FOXO3a and its downstream gene Bim impaired AZD6244-induced growth suppression, whereas restoring activation of FOXO3a sensitized human cancer cell to AZD6244-induced growth arrest and apoptosis. More importantly, AZD6244-resistant cancer cells showed impaired endogenous FOXO3a nuclear translocation, reduced FOXO3a-Bim promoter association and significantly decreased Bim expression in response to AZD6244. AZD6244-resistant cancer cells can be sensitized to API-2 (an AKT inhibitor) and LY294002 (a PI3K inhibitor) in suppressing cell growth and colony formation, these inhibitors were known to enhance FOXO3a activity/nuclear translocation through inhibiting PI3K-AKT pathway. This study reveals novel molecular mechanism contributing to AZD6244-resistance through regulation of FOXO3a activity, further provides significant clinical implication of combining AZD6244 with PI3K/AKT inhibitors for sensitizing AZD6244-resistant cancer cells by activating FOXO3a. FOXO3a activation can be an essential pharmacological target and indicator to mediate and predict AZD6244 efficacy in clinical use. ^

L -Sox5 and Sox6: Architectural transcription factors in chondrogenesis

Transcription factors often determine cell fate and tissue development. Chondrogenesis is the developmental process by which cartilages form. Recently, gene targeting studies have shown that two transcription factors, L-Sox5 and Sox6, play essential and redundant roles in chondrogenesis in vivo by converting precartilaginous cell condensations into cartilages. Both are highly similar High-Mobility-Group (HMG)-domain proteins that bind and subsequently bend DNA containing the 7bp HMG site (A/T)(A/T)CAA(A/T)G. They have no transactivation domain, but homo- and hetero-dimerize and preferentially bind DNA containing two HMG sites. They are thought to play an architectural role in transactivation by facilitating long-range DNA and protein interactions. To understand their molecular mechanism of action, we investigated how phasing, orientation, and spacing between HMG sites affect L-Sox5 and Sox6 DNA-binding. We determined that L-Sox5 and Sox6 dimers bind with high affinity to paired HMG sites in DNA rather than a single HMG site. Binding of paired sites is independent of DNA helical phasing, orientation of paired HMG sites and independent of distance up to 255 base pairs between sites. Mutational analysis demonstrated that binding of L-Sox5 and Sox6, independent of orientation of the sites, is critically dependent on the presence of paired HMG sites rather than one HMG site alone. Our data support a unique and novel model whereby L-Sox5 and Sox6 dimerize and bind DNA with pronounced spatial flexibility, possibly by a flexible hinge, and act as architectural transcription factors that bring distant DNA sites and proteins together to form higher order transcriptional complexes that are essential for the activation of their target genes in chondrogenesis. ^

Insights into transcription enhancer factor 1 (TEF-1) activity from the solution structure of the TEA domain.

Transcription enhancer factor 1 is essential for cardiac, skeletal, and smooth muscle development and uses its N-terminal TEA domain (TEAD) to bind M-CAT elements. Here, we present the first structure of TEAD and show that it is a three-helix bundle with a homeodomain fold. Structural data reveal how TEAD binds DNA. Using structure-function correlations, we find that the L1 loop is essential for cooperative loading of TEAD molecules on to tandemly duplicated M-CAT sites. Furthermore, using a microarray chip-based assay, we establish that known binding sites of the full-length protein are only a subset of DNA elements recognized by TEAD. Our results provide a model for understanding the regulation of genome-wide gene expression during development by TEA/ATTS family of transcription factors.

Human TOB, an antiproliferative transcription factor, is a poly(A)-binding protein-dependent positive regulator of cytoplasmic mRNA deadenylation.

In mammalian cells, mRNA decay begins with deadenylation, which involves two consecutive phases mediated by the PAN2-PAN3 and the CCR4-CAF1 complexes, respectively. The regulation of the critical deadenylation step and its relationship with RNA-processing bodies (P-bodies), which are thought to be a site where poly(A)-shortened mRNAs get degraded, are poorly understood. Using the Tet-Off transcriptional pulsing approach to investigate mRNA decay in mouse NIH 3T3 fibroblasts, we found that TOB, an antiproliferative transcription factor, enhances mRNA deadenylation in vivo. Results from glutathione S-transferase pull-down and coimmunoprecipitation experiments indicate that TOB can simultaneously interact with the poly(A) nuclease complex CCR4-CAF1 and the cytoplasmic poly(A)-binding protein, PABPC1. Combining these findings with those from mutagenesis studies, we further identified the protein motifs on TOB and PABPC1 that are necessary for their interaction and found that interaction with PABPC1 is necessary for TOB's deadenylation-enhancing effect. Moreover, our immunofluorescence microscopy results revealed that TOB colocalizes with P-bodies, suggesting a role of TOB in linking deadenylation to the P-bodies. Our findings reveal a new mechanism by which the fate of mammalian mRNA is modulated at the deadenylation step by a protein that recruits poly(A) nuclease(s) to the 3' poly(A) tail-PABP complex.

The molecular and genetic characterization of the MADS box transcription factor {\it Drosophila\/} MEF2

The myocyte enhancer factor (MEF)-2 family of transcription factors has been implicated in the regulation of muscle transcription in vertebrates, but the precise position of these regulators within the genetic hierarchy leading to myogenesis is unclear. The MEF2 proteins bind to a conserved A/T-rich DNA sequence present in numerous muscle-specific genes, and they are expressed in the cells of the developing somites and in the embryonic heart at the onset of muscle formation in mammals. The MEF2 genes belong to the MADS box family of transcription factors, which control specific programs of gene expression in species ranging from yeast to humans. Each MEF2 family member contains two highly conserved protein motifs, the MADS domain and the MEF2-specific domain, which together provide the MEF2 factors with their unique DNA binding and dimerization properties. In an effort to further define the function of the MEF2 proteins, and to evaluate the degree of conservation shared among these factors and the phylogenetic pathways that they regulate, we sought to identify MEF2 family members in other species. In Drosophila, a homolog of the vertebrate MEF2 genes was identified and termed D-mef2. The D-MEF2 protein binds to the consensus MEF2 element and can activate transcription through tandem copies of that site. During Drosophila embryogenesis, D-MEF2 is specific to the mesoderm germ layer of the developing embryo and becomes expressed in all muscle cell types within the embryo. The role of D-mef2 in Drosophila embryogenesis was examined by generating a loss-of-function mutation in the D-mef2 gene. In embryos homozygous for this mutant allele, somatic, cardiac, and visceral muscles fail to differentiate, but precursors of these myogenic lineages are normally specified and positioned. These results demonstrate that different muscle cell types share a common myogenic differentiation program controlled by MEF2 and suggest that this program has been conserved from Drosophila to mammals. ^

Functional characterization of transcription factor USF

USF, Upstream Stimulatory Factor, is a family of ubiquitous transcription factors that contain highly conserved basic helix-loop-helix leucine zipper DNA binding domains and recognize the core DNA sequence CACGTG. In human and mouse, two members of the USF family, USF1 and USF2, encoded by two different genes, contribute to the USF activity. In order to gain insights into the mechanisms by which USFs function as transcriptional activators, different approaches were used to map the domains of USF2 responsible for nuclear localization and transcriptional activation. Two stretches of amino acids, one in the basic region of the DNA binding domain, the other in a highly conserved N-terminal region, were found to direct nuclear localization independently of one another. Two distinct activation domains were also identified. The first one, located in the conserved N-terminal region that overlaps the C-terminal nuclear localization signal, functioned only in the presence of an initiator element in the promoter of the reporter. The second, in a nonconserved region, activated transcription in the absence of an initiator element or when fused to a heterologous DNA binding domain. These results suggest that USF2 functions in different promoter contexts by selectively utilizing different activation domains.^ The deletion analysis of USF2 also identified two dominant negative mutants of USF, one lacking the activation domain, the other lacking the basic domain. The latter proved useful for testing the direct involvement of USFs in the transcriptional activation mediated by the viral protein IE62.^ To investigate the biological function of USFs, foci and colony formation assays were used to study the growth regulation by USFs. It was found that USFs had a strong antagonistic effect on cellular transformation mediated by the bHLH/LZ protein Myc. This effect required the DNA binding activity of either USF 1 or USF2. Moreover, USF2, but not USF1 or other mutants of USFs, was also found to have strong inhibitory effect on the cellular transformation by E1a and on the growth of HeLa cells. These results demonstrate that USFs could potentially regulate growth through two mechanisms, one by antagonizing the function of Myc in cellular transformation, the other by mediating a more general growth inhibitory effect. ^

Transcription factor SEF: Purification and functional characterization

Cytokine-induced transcription of the serum amyloid A3 (SAA3) gene promoter requires a transcriptional enhancer that contains three functional elements: two C/EBP-binding sites and a third site that interacts with a constitutively expressed transcription factor, SAA3 enhancer factor (SEF). Deletion or site-specific mutations in the SEF-binding site drastically reduced SAA3 promoter activity, strongly suggesting that SEF is important in SAA3 promoter function. To further elucidate its role in the regulation of the SAA3 gene, we purified SEF from HeLa cell nuclear extracts to near homogeneity by using conventional liquid chromatography and DNA-affinity chromatography. Ultraviolet cross-linking and Southwestern experiments indicated that SEF consisted of a single polypeptide with an apparent molecular mass of 65 kDa. Protein sequencing, oligonucleotide competition and antibody supershift experiments identified SEF as transcription factor LBP-1c/CP2/LSF. Cotransfection of SEF expression plasmid with SAA3-luciferase reporter resulted in 3- to 5-fold activation of SAA3 promoter. Interestingly, when SEF-transfected cells were treated with either conditioned medium (CM) or interleukin (IL) 1, the SAA3 promoter was synergistically activated in a dose-dependent manner. Furthermore, when SEF-binding site was mutated, the response of SAA3 promoter to IL-1 or CM stimulation was abolished or drastically decreased, suggesting that SEF may functionally cooperate with an IL-1-inducible transcription factor. Indeed, our functional studies showed that NFκB is a key transcription factor that mediates the IL-1-induced expression of SAA3 gene, and that SEF can synergize with NFκBp65 to activate SAA3 promoter. By coimmunoprecipitation experiments, we found that SEF could specifically interact with NFκBp65, and that the association of these two factors was enhanced upon IL-1 and CM stimulation. This suggests that the molecular basis for the functional synergy between SEF and NFκB may be due to the ability of SEF to physically interact with NPκB. In addition to its interaction with SEF, NFκB-dependent activation also requires the weak κB site in the C element and its interaction with C/EBP. Besides its role in regulating SAA3 gene expression, we provide evidence that SEF could also bind in a sequence-specific manner to the promoters of α2-macroglobulin, Aα fibrinogen, and 6–16 genes and to an intronic enhancer of the human Wilm's tumor 1 gene, suggesting a functional role in the regulation of these genes. By coimmunoprecipitation experiments, we determined that SEF could specifically associate with both Stat3 and Stat2 upon cytokine stimulation. To examine the functional roles of such interactions, we evaluated the effects of SEF on the transcriptional regulation of two reporter genes: Aα fibrinogen and 6–16, which are IL-6- and interferon-α-responsive, respectively. Our results showed that cotransfection of SEF expression plasmid can activate the expression of Aα fibrinogen gene and 6–16 gene. Moreover, SEF can dramatically enhance the interferon-α-induced expression of 6–16 gene and IL-6-induced expression of Aα fibrinogen gene, suggesting that SEF may functionally cooperate with ISGF3 and Stat3 to mediate interferon-α and IL-6 signaling. ^ Our findings that SEF can interact with multiple cytokine-inducible transcription factors to mediate the expression of target genes open a new avenue of investigation of cooperative transcriptional regulation of gene expression, and should further our understanding of differential gene expression in response to a specific stimulus. In summary, our data provide evidence that SEF can mediate the signaling of different cytokines by interacting with various cytokine-inducible transcription factors. ^

Alteration of CEACAM1 gene transcription in prostate cancer cells: The role of AR, Sp2, and HDAC

Cell to cell adhesion molecule (CEACAM1), a type II tumor suppressor, has been found to be down-regulated in prostate cancer cells. The mechanism that causes CEACAM1's down-regulation in tumorigenesis is unknown. Here we show that the transcriptional activity of CEACAM1 is decreased in prostate cancer cells. This decrease is not due to methylation of the CEACAM1's promoter, but rather to the alteration of transcription factors regulating CEACAM1 expression. ^ Since androgen/androgen receptors (AR) are potent regulators of prostate growth and differentiation, their role on CEACAM1 gene transcription was examined. The androgen receptor could directly increase CEACAM1 transcriptional activity in a ligand dependent manner by interacting with an AR consensus element that resides in the CEACAM1 promoter. However, AR binding to the CEACAM1 promoter is not related to the loss of CEACAM1 during prostate cancer progression. ^ Further analysis enabled us to determine the particular region in the CEACAM1 promoter that mediates a decrease in CEACAM1 transcriptional activity in prostate cancer cells. Upon further examination, we found that this CEACAM1 promoter region interacts with the Sp1, Sp2, and Sp3 transcription factors. However, only Sp2 expression was found to increase in prostate cancer cells. Inhibiting Sp2 from binding to the CEACAM1 promoter caused an increase in CEACAM1 transcriptional activity in prostate cancer cells. In addition, over-expressing Sp2 in normal prostate cells resulted in a decrease in CEACAM1 transcriptional activity and endogenous protein expression. These observations suggest that Sp2 is a transcription repressor of CEACAM1. Furthermore, prostate cancer cells treated with trichostatin A (TSA), a specific histone deacetylase (HDAC) inhibitor, activated CEACAM1 transcriptional activity. This implies that HDACs are involved in CEACAM1 transcriptional activity. Mutation of the Sp2 DNA binding region on the CEACAM1 promoter inhibited TSA activation of CEACAM1 transcriptional activity. This indicates that HDACs inhibit CEACAM1 transcriptional activity through Sp2. Base on these results, we propose that Sp2 is critical for down-regulating CEACAM1 expression, and one mechanism by which Sp2 represses CEACAM1 expression is by recruiting HDAC to the CEACAM1 promoter in prostate cancer cells. Collectively, these findings provide novel insights into mechanisms that cause the down-regulation of CEACAM1 expression in prostate cancer cells. ^

Regulation of chromatin structure, Smad binding and AFP expression by the Forkhead box transcription factor: Foxa1

Transcription factors must be able to access their DNA binding sites to either activate or repress transcription. However, DNA wrapping and compaction into chromatin occludes most binding sites from ready access by proteins. Pioneer transcription factors are capable of binding their DNA elements within a condensed chromatin context and then reducing the level of nucleosome occupancy so that the chromatin structure is more accessible. This altered accessibility increases the probability of other transcription factors binding to their own DNA binding elements. My hypothesis is that Foxa1, a ‘pioneer’ transcription factor, activates alpha-fetoprotein (AFP) expression by binding DNA in a chromatinized environment, reducing the nucleosome occupancy and facilitating binding of additional transcription factors.^ Using retinoic-acid differentiated mouse embryonic stem cells, we illustrate a mechanism for activation of the tumor marker AFP by the pioneer transcription factor Foxa1 and TGF-β downstream effector transcription factors Smad2 and Smad4. In differentiating embryonic stem cells, binding of the Foxa1 forkhead box transcription factor to chromatin reduces nucleosome occupancy and levels of linker histone H1 at the AFP distal promoter. The more accessible DNA is subsequently bound by the Smad2 and Smad4 transcription factors, concurrent with activation of transcription. Chromatin immunoprecipitation analyses combined with siRNA-mediated knockdown indicate that Smad protein binding and the reduction of nucleosome occupancy at the AFP distal promoter is dependent on Foxa1. In addition to facilitating transcription factor binding, Foxa1 is also associated with histone modifications related to active gene expression. Acetylation of lysine 9 on histone H3, a mark that is associated active transcription, is dependent on Foxa1, while methylation of H3K4, also associated with active transcription, is independent of Foxa1. I propose that Foxa1 potentiates a region of chromatin to respond to Smad proteins, leading to active expression of AFP.^ These studies demonstrate one mechanism whereby a transcription factor can alter the accessibility of additional transcription factors to chromatin, by altering nucleosome positions. Specifically, Foxa1 exposes DNA so that Smad4 can bind to its regulatory element and activate transcription of the tumor-marker gene AFP.^

Regulation of myogenesis by growth signals -growth factors, oncogenes and protein kinases

Establishment of a myogenic phenotype involves antagonism between cell proliferation and differentiation. The recent identification of the MyoD family of muscle-specific transcription factors provides opportunities to dissect at the molecular level the mechanisms through which defined cell type-specific transcription factors respond to environmental cues and regulate differentiation programs. This project is aimed at elucidation of the molecular mechanism whereby growth factors repress myogenesis. Initial studies demonstrated that nuclear oncogenes such as c-fos, junB and c-jun are immediate early genes that respond to serum and TGF-$\beta$. Using the muscle creatine kinase (MCK) enhancer linked to the reporter gene CAT as a marker for differentiation, we showed that transcriptional function of myogenin can be disrupted in the presence of c-Fos, JunB and cjun. In contrast, JunD, which shares DNA-binding specificity with JunB and c-Jun but is expressed constitutively in muscle cells, failed to show the inhibition. The repression by Fos and Jun is targeted at KE-2 motif, the same sequence that mediates myogenin-dependent activation and muscle-specific transactivation. Deletion analysis indicated that the transactivation domain of c-Jun at the N-terminus is responsible for the repression. Considering that myogenin is a phosphoprotein and cAMP and TPA are able to regulate myogenesis, we examined whether constitutively active protein kinase C (PKC) and protein kinase A (PKA) could substitute for exogenous growth factors and prevent transcription activation by myogenin. Indeed, the basic region of myogenin is phosphorylated by PKC at a threonine that is conserved in all members of the MyoD family. Phosphorylation at this site attenuates DNA binding activity of myogenin. Protein kinase A can also phosphorylate myogenin in a region adjacent to the DNA binding domain. However, phosphorylation at this site is insufficient to abrogate myogenin's DNA binding capacity, suggesting that PKA and PKC may affect myogenin transcriptional activity through different mechanisms. These findings provide insight into the mechanisms through which growth factor signals negatively regulate the muscle differentiation program and contribute to an understanding of signal transducing pathways between the cell membrane and nucleus. ^

Functional studies of a LIM -homeodomain transcription factor lmx1b in murine mid -hindbrain and limb development

This thesis is centered on applying molecular genetics to study pattern formation during animal development. More specifically, this thesis describes the functional studies of a LIM-homeodomain gene called lmx1b during murine embryogenesis. Lmx1b expression is restricted to the mid-hindbrain junction as well as to the dorsal mesenchyme of the limb, suggesting important functions during mid-hindbrain and limb development. To test these possibilities, lmx1b homozygous mutant mice were generated and their limb and CNS phenotypes examined. Lmx1b homozygous mutant mice exhibit a large reduction of mid-hindbrain structures, and that their limbs are symmetrical along the dorsal-ventral axis as the result of a dorsal to ventral transformation. Taken together, these studies define essential functions for lmx1b in mid-hindbrain patteming and in dorsal limb cell fate determination. However, the molecular mechanisms which accounts for these phenotypes are unknown, and whether lmx1b has same or distinctive functions during the mid-hindbrain and limb development is also unclear. ^ Recently, insight into molecular mechanisms of mid-hindbrain patterning and limb development has resulted from the identification of several factors with restricted expression patterns within these regions. These include the secreted factors wnt-1, fgf-8, wnt-7a and the transcription factors pax-2, and en-1. Targeted disruption of any of these genes in mice suggests that these genes might be involved in similar regulatory pathways. Analysis of the expression of these genes in lmx1b mutants demonstrates that lmxlb is not required for the initiation, but is required to maintain their expression at the mid-hindbrain junction. Thus, lmxlb is not required for specifying mid-hindbrain cell fates, rather, it functions to ensure the establishment or maintenance of a proper organizing center at the mid-hindbrain junction. Interestingly, lmxlb functions cell non-autonomously in chimera analysis, which indicates that lmx1b might regulate the expression of secreted factors such as wnt-1 and/or fgf-8 in the organizing center. In contrast, lmx1b functions cell autonomously in the dorsal limb to govern dorsal ventral limb development and its expression is regulated by with wnt-7a and en-1. However, single and double mutant analysis suggest that all three genes have partially overlapping functions as well as independent functions. The results point toward a complicated network of cross-talks among all three limb axes. ^

NEW TARGET GENES FOR TUMOR SUPPRESSORS p53 AND p73 IN REGENERATING LIVER

The p53-family of proteins regulates expression of target genes during tissue development and differentiation. Within the p53-family, p53 and p73 have hepatic-specific functions in development and tumor suppression. Despite a growing list of p53/p73 target genes, very few of these have been studied in vivo, and the knowledge regarding functions of p53 and p73 in normal tissues remains limited. p53+/-p73+/- mice develop hepatocellular carcinoma (HCC), whereas overexpression of p53 in human HCC leads to tumor regression. However, the mechanism of p53/p73 function in liver remains poorly characterized. Here, the model of mouse liver regeneration is used to identify new target genes for p53/p73 in normal quiescent vs. proliferating cells. In response to surgical removal of ~2/3 of liver mass (partial hepatectomy, PH), the remaining hepatocytes exit G0 of cell cycle and undergo proliferation to reestablish liver mass. The hypothesis tested in this work is that p53/p73 functions in cell cycle arrest, apoptosis and senescence are repressed during liver regeneration, and reactivated at the end of the regenerative response. Chromatin immunoprecipitation (ChIP), with a p73-antibody, was used to probe arrayed genomic sequences (ChIP-chip) and uncover 158 potential targets of p73-regulation in normal liver. Global microarray analysis of mRNA levels, at T=0-48h following PH, revealed sets of genes that change expression during regeneration. Eighteen p73-bound genes changed expression after PH. Four of these genes, Foxo3, Jak1, Pea15, and Tuba1 have p53 response elements (p53REs), identified in silico within the upstream regulatory region. Forkhead transcription factor Foxo3 is the most responsive gene among transcription factors with altered expression during regenerative, cellular proliferation. p53 and p73 bind a Foxo3 p53RE and maintain active expression in quiescent liver. During liver regeneration, binding of p53 and p73, recruitment of acetyltransferase p300, and an active chromatin structure of Foxo3 are disrupted, alongside loss of Foxo3 expression. These parameters of Foxo3 regulation are reestablished at completion of liver growth and regeneration, supporting a temporary suspension of p53 and p73 regulatory functions in normal cells during tissue regeneration.

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

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.