Unlimited multistability in multisite phosphorylation systems

Reversible phosphorylation on serine, threonine and tyrosine is the most widely studied posttranslational modification of proteins (1, 2). The number of phosphorylated sites on a protein (n) shows a significant increase from prokaryotes, with n less than or equal to 7 sites, to eukaryotes, with examples having n greater than or equal to 150 sites (3). Multisite phosphorylation has many roles (4, 5) and site conservation indicates that increasing numbers of sites cannot be due merely to promiscuous phosphorylation. A substrate with n sites has an exponential number (2^n) of phospho-forms and individual phospho-forms may have distinct biological effects (6, 7). The distribution of these phospho-forms and how this distribution is regulated have remained unknown. Here we show that, when kinase and phosphatase act in opposition on a multisite substrate, the system can exhibit distinct stable phospho-form distributions at steady state and that the maximum number of such distributions increases with n. Whereas some stable distributions are focused on a single phospho-form, others are more diffuse, giving the phospho-proteome the potential to behave as a fluid regulatory network able to encode information and flexibly respond to varying demands. Such plasticity may underlie complex information processing in eukaryotic cells (8) and suggests a functional advantage in having many sites. Our results follow from the unusual geometry of the steady-state phospho-form concentrations, which we show to constitute a rational algebraic curve, irrespective of n. We thereby reduce the complexity of calculating steady states from simulating 3 times 2^n differential equations to solving two algebraic equations, while treating parameters symbolically. We anticipate that these methods can be extended to systems with multiple substrates and multiple enzymes catalysing different modifications, as found in posttranslational modification 'codes' (9) such as the histone code (10, 11). Whereas simulations struggle with exponentially increasing molecular complexity, mathematical methods of the kind developed here can provide a new language in which to articulate the principles of cellular information processing (12)

Climbing Fiber Burst Size and Olivary Sub-threshold Oscillations in a Network Setting

The inferior olivary nucleus provides one of the two main inputs to the cerebellum: the so-called climbing fibers. Activation of climbing fibers is generally believed to be related to timing of motor commands and/or motor learning. Climbing fiber spikes lead to large all-or-none action potentials in cerebellar Purkinje cells, overriding any other ongoing activity and silencing these cells for a brief period of time afterwards. Empirical evidence shows that the climbing fiber can transmit a short burst of spikes as a result of an olivary cell somatic spike, potentially increasing the information being transferred to the cerebellum per climbing fiber activation. Previously reported results from in vitro studies suggested that the information encoded in the climbing fiber burst is related to the occurrence of the spike relative to the ongoing sub-threshold membrane potential oscillation of the olivary cell, i.e. that the phase of the oscillation is reflected in the size of the climbing fiber burst. We used a detailed three-compartmental model of an inferior olivary cell to further investigate the possible factors determining the size of the climbing fiber burst. Our findings suggest that the phase-dependency of the burst size is present but limited and that charge flow between soma and dendrite is a major determinant of the climbing fiber burst. From our findings it follows that phenomena such as cell ensemble synchrony can have a big effect on the climbing fiber burst size through dendrodendritic gap-junctional coupling between olivary cells

The disruption of proteostasis in neurodegenerative diseases

Cells count on surveillance systems to monitor and protect the cellular proteome which, besides being highly heterogeneous, is constantly being challenged by intrinsic and environmental factors. In this context, the proteostasis network (PN) is essential to achieve a stable and functional proteome. Disruption of the PN is associated with aging and can lead to and/or potentiate the occurrence of many neurodegenerative diseases (ND). This not only emphasizes the importance of the PN in health span and aging but also how its modulation can be a potential target for intervention and treatment of human diseases.info:eu-repo/semantics/publishedVersio

Protein Kinase A Binds and Activates Heat Shock Factor 1

BACKGROUND. Many inducible transcription factors are regulated through batteries of posttranslational modifications that couple their activity to inducing stimuli. We have studied such regulation of Heat Shock Factor 1 (HSF1), a key protein in control of the heat shock response, and a participant in carcinogenisis, neurological health and aging. As the mechanisms involved in the intracellular regulation of HSF1 in good health and its dysregulation in disease are still incomplete we are investigating the role of posttranslational modifications in such regulation. METHODOLOGY/PRINCIPAL FINDINGS. In a proteomic study of HSF1 binding partners, we have discovered its association with the pleiotropic protein kinase A (PKA). HSF1 binds avidly to the catalytic subunit of PKA, (PKAca) and becomes phosphorylated on a novel serine phosphorylation site within its central regulatory domain (serine 320 or S320), both in vitro and in vivo. Intracellular PKAca levels and phosphorylation of HSF1 at S320 were both required for HSF1 to be localized to the nucleus, bind to response elements in the promoter of an HSF1 target gene (hsp70.1) and activate hsp70.1 after stress. Reduction in PKAca levels by small hairpin RNA led to HSF1 exclusion from the nucleus, its exodus from the hsp70.1 promoter and decreased hsp70.1 transcription. Likewise, null mutation of HSF1 at S320 by alanine substitution for serine led to an HSF1 species excluded from the nucleus and deficient in hsp70.1 activation. CONCLUSIONS. These findings of PKA regulation of HSF1 through S320 phosphorylation add to our knowledge of the signaling networks converging on this factor and may contribute to elucidating its complex roles in the stress response and understanding HSF1 dysregulation in disease.National Institutes of Health (2RO1CA047407, RO1CA077465

A classification-based framework for predicting and analyzing gene regulatory response

BACKGROUND: We have recently introduced a predictive framework for studying gene transcriptional regulation in simpler organisms using a novel supervised learning algorithm called GeneClass. GeneClass is motivated by the hypothesis that in model organisms such as Saccharomyces cerevisiae, we can learn a decision rule for predicting whether a gene is up- or down-regulated in a particular microarray experiment based on the presence of binding site subsequences ("motifs") in the gene's regulatory region and the expression levels of regulators such as transcription factors in the experiment ("parents"). GeneClass formulates the learning task as a classification problem — predicting +1 and -1 labels corresponding to up- and down-regulation beyond the levels of biological and measurement noise in microarray measurements. Using the Adaboost algorithm, GeneClass learns a prediction function in the form of an alternating decision tree, a margin-based generalization of a decision tree. METHODS: In the current work, we introduce a new, robust version of the GeneClass algorithm that increases stability and computational efficiency, yielding a more scalable and reliable predictive model. The improved stability of the prediction tree enables us to introduce a detailed post-processing framework for biological interpretation, including individual and group target gene analysis to reveal condition-specific regulation programs and to suggest signaling pathways. Robust GeneClass uses a novel stabilized variant of boosting that allows a set of correlated features, rather than single features, to be included at nodes of the tree; in this way, biologically important features that are correlated with the single best feature are retained rather than decorrelated and lost in the next round of boosting. Other computational developments include fast matrix computation of the loss function for all features, allowing scalability to large datasets, and the use of abstaining weak rules, which results in a more shallow and interpretable tree. We also show how to incorporate genome-wide protein-DNA binding data from ChIP chip experiments into the GeneClass algorithm, and we use an improved noise model for gene expression data. RESULTS: Using the improved scalability of Robust GeneClass, we present larger scale experiments on a yeast environmental stress dataset, training and testing on all genes and using a comprehensive set of potential regulators. We demonstrate the improved stability of the features in the learned prediction tree, and we show the utility of the post-processing framework by analyzing two groups of genes in yeast — the protein chaperones and a set of putative targets of the Nrg1 and Nrg2 transcription factors — and suggesting novel hypotheses about their transcriptional and post-transcriptional regulation. Detailed results and Robust GeneClass source code is available for download from

IL-6 Stabilizes Twist and Enhances Tumor Cell Motility in Head and Neck Cancer Cells through Activation of Casein Kinase 2

BACKGROUND: Squamous cell carcinoma of the head and neck (SCCHN) is the seventh most common cancer worldwide. Unfortunately, the survival of patients with SCCHN has not improved in the last 40 years, and thus new targets for therapy are needed. Recently, elevations in serum level of interleukin 6 (IL-6) and expression of Twist in tumor samples were found to be associated with poor clinical outcomes in multiple types of cancer, including SCCHN. Although Twist has been proposed as a master regulator of epithelial-mesenchymal transition and metastasis in cancers, the mechanisms by which Twist levels are regulated post-translationally are not completely understood. Tumor progression is characterized by the involvement of cytokines and growth factors and Twist induction has been connected with a number of these signaling pathways including IL-6. Since many of the effects of IL-6 are mediated through activation of protein phosphorylation cascades, this implies that Twist expression must be under a tight control at the post-translational level in order to respond in a timely manner to external stimuli. METHODOLOGY/PRINCIPAL FINDINGS: Our data show that IL-6 increases Twist expression via a transcription-independent mechanism in many SCCHN cell lines. Further investigation revealed that IL-6 stabilizes Twist in SCCHN cell lines through casein kinase 2 (CK2) phosphorylation of Twist residues S18 and S20, and that this phosphorylation inhibits degradation of Twist. Twist phosphorylation not only increases its stability but also enhances cell motility. Thus, post-translational modulation of Twist contributes to its tumor-promoting properties. CONCLUSIONS/SIGNIFICANCE: Our study shows Twist expression can be regulated at the post-translational level through phosphorylation by CK2, which increases Twist stability in response to IL-6 stimulation. Our findings not only provide novel mechanistic insights into post-translational regulation of Twist but also suggest that CK2 may be a viable therapeutic target in SCCHN

The Acute Phase Protein Ceruloplasmin as a Non-Invasive Marker of Pseudopregnancy, Pregnancy, and Pregnancy Loss in the Giant Panda

After ovulation, non-pregnant female giant pandas experience pseudopregnancy. During pseudopregnancy, non-pregnant females exhibit physiological and behavioral changes similar to pregnancy. Monitoring hormonal patterns that are usually different in pregnant mammals are not effective at determining pregnancy status in many animals that undergo pseudopregnancy, including the giant panda. Therefore, a physiological test to distinguish between pregnancy and pseudopregnancy in pandas has eluded scientists for decades. We examined other potential markers of pregnancy and found that activity of the acute phase protein ceruloplasmin increases in urine of giant pandas in response to pregnancy. Results indicate that in term pregnancies, levels of active urinary ceruloplasmin were elevated the first week of pregnancy and remain elevated until 20–24 days prior to parturition, while no increase was observed during the luteal phase in known pseudopregnancies. Active ceruloplasmin also increased during ultrasound-confirmed lost pregnancies; however, the pattern was different compared to term pregnancies, particularly during the late luteal phase. In four out of the five additional reproductive cycles included in the current study where females were bred but no birth occurred, active ceruloplasmin in urine increased during the luteal phase. Similar to the known lost pregnancies, the temporal pattern of change in urinary ceruloplasmin during the luteal phase deviated from the term pregnancies suggesting that these cycles may have also been lost pregnancies. Among giant pandas in captivity, it has been presumed that there is a high rate of pregnancy loss and our results are the first to provide evidence supporting this notion

Protein Scaffolds Can Enhance the Bistability of Multisite Phosphorylation Systems

The phosphorylation of a substrate at multiple sites is a common protein modification that can give rise to important structural and electrostatic changes. Scaffold proteins can enhance protein phosphorylation by facilitating an interaction between a protein kinase enzyme and its target substrate. In this work we consider a simple mathematical model of a scaffold protein and show that under specific conditions, the presence of the scaffold can substantially raise the likelihood that the resulting system will exhibit bistable behavior. This phenomenon is especially pronounced when the enzymatic reactions have sufficiently large KM, compared to the concentration of the target substrate. We also find for a closely related model that bistable systems tend to have a specific kinetic conformation. Using deficiency theory and other methods, we provide a number of necessary conditions for bistability, such as the presence of multiple phosphorylation sites and the dependence of the scaffold binding/unbinding rates on the number of phosphorylated sites

Regulation of Cyclooxygenase-2 Expression by Heat: A Novel Aspect of Heat Shock Factor 1 Function in Human Cells

The heat-shock response, a fundamental defense mechanism against proteotoxic stress, is regulated by a family of heat-shock transcription factors (HSF). In humans HSF1 is considered the central regulator of heat-induced transcriptional responses. The main targets for HSF1 are specific promoter elements (HSE) located upstream of heat-shock genes encoding cytoprotective heat-shock proteins (HSP) with chaperone function. In addition to its cytoprotective function, HSF1 was recently hypothesized to play a more complex role, regulating the expression of non-HSP genes; however, the non-canonical role of HSF1 is still poorly understood. Herein we report that heat-stress promotes the expression of cyclooxygenase-2 (COX-2), a key regulator of inflammation controlling prostanoid and thromboxane synthesis, resulting in the production of high levels of prostaglandin-E2 in human cells. We show that heat-induced COX-2 expression is regulated at the transcriptional level via HSF1-mediated signaling and identify, by in-vitro reporter gene activity assay and deletion-mutant constructs analysis, the COX-2 heat-responsive promoter region and a new distal cis-acting HSE located at position −2495 from the transcription start site. As shown by ChIP analysis, HSF1 is recruited to the COX-2 promoter rapidly after heat treatment; by using shRNA-mediated HSF1 suppression and HSE-deletion from the COX-2 promoter, we demonstrate that HSF1 plays a central role in the transcriptional control of COX-2 by heat. Finally, COX-2 transcription is also induced at febrile temperatures in endothelial cells, suggesting that HSF1-dependent COX-2 expression could contribute to increasing blood prostaglandin levels during fever. The results identify COX-2 as a human non-classical heat-responsive gene, unveiling a new aspect of HSF1 function

Pancreatic islet-autonomous insulin and smoothened-mediated signalling modulate identity changes of glucagon+ α-cells

The mechanisms that restrict regeneration and maintain cell identity following injury are poorly characterized in higher vertebrates. Following β-cell loss, 1-2% of the glucagon-producing α-cells spontaneously engage in insulin production in mice. Here we explore the mechanisms inhibiting α-cell plasticity. We show that adaptive α-cell identity changes are constrained by intra-islet insulin- and Smoothened-mediated signalling, among others. The combination of β-cell loss or insulin-signalling inhibition, with Smoothened inactivation in α- or δ-cells, stimulates insulin production in more α-cells. These findings suggest that the removal of constitutive 'brake signals' is crucial to neutralize the refractoriness to adaptive cell-fate changes. It appears that the maintenance of cell identity is an active process mediated by repressive signals, which are released by neighbouring cells and curb an intrinsic trend of differentiated cells to change