Pathway databases and tools for their exploitation: benefits, current limitations and challenges

In past years, comprehensive representations of cell signalling pathways have been developed by manual curation from literature, which requires huge effort and would benefit from information stored in databases and from automatic retrieval and integration methods. Once a reconstruction of the network of interactions is achieved, analysis of its structural features and its dynamic behaviour can take place. Mathematical modelling techniques are used to simulate the complex behaviour of cell signalling networks, which ultimately sheds light on the mechanisms leading to complex diseases or helps in the identification of drug targets. A variety of databases containing information on cell signalling pathways have been developed in conjunction with methodologies to access and analyse the data. In principle, the scenario is prepared to make the most of this information for the analysis of the dynamics of signalling pathways. However, are the knowledge repositories of signalling pathways ready to realize the systems biology promise? In this article we aim to initiate this discussion and to provide some insights on this issue.

Neutralizing antibodies in Multiple Sclerosis a model-based analysis of Interferon beta signaling pathway in macrophages

Multiple Sclerosis is the most common non-traumatic cause of neurologicaldisability in young people. There is no cure yet, and until recently, few long-termtherapies existed. Interferon beta (IFNβ) was the first treatment, and remains the mostcommonly prescribed. One of the most significant problems of IFNβ therapy is theproduction of drug specific antibodies. Up to 45% of patients develop neutralizingantibodies (NAbs) to IFNβ products. The neutralizing antibody binds to the biologicalagent preventing its interaction with its receptor, inhibiting the biological action of theprotein, which abrogates the clinical efficacy of IFNβ treatment. Interferon-betamediates its response by binding to its high affinity cell surface receptor and initiatingthe JAK/STAT signalling cascade. In this project we have analyzed the IFNβ signalingpathway in macrophages when neutralizing antibodies are present. The response tothis pathway after IFNβ stimulation shows a transient oscillatory rhythm of STAT1phosphorylation, which varies as NAbs concentration increases. To improve ourunderstanding of that behavior, we extended an existing mathematical model based onnonlinear ordinary differential equations of JAK/STAT pathway by including IFN-NAbassociation and IFN-activation receptor. Combining our theoretical model withexperimental data we could study the role of neutralizing antibodies on the molecularresponse and determine its lifetime after cytokine stimulation.

Muscle protein waste in tumor-bearing rats is effectively antagonized by a beta 2-adrenergic agonist (clenbuterol). Role of the ATP-ubiquitin-dependent proteolytic pathway.

Tissue protein hypercatabolism (TPH) is a most important feature in cancer cachexia, particularly with regard to the skeletal muscle. The rat ascites hepatoma Yoshida AH-130 is a very suitable model system for studying the mechanisms involved in the processes that lead to tissue depletion, since it induces in the host a rapid and progressive muscle waste mainly due to TPH (Tessitore, L., G. Bonelli, and F. M. Baccino. 1987. Biochem. J. 241:153-159). Detectable plasma levels of tumor necrosis factor-alpha associated with marked perturbations in the hormonal homeostasis have been shown to concur in forcing metabolism into a catabolic setting (Tessitore, L., P. Costelli, and F. M. Baccino. 1993. Br. J. Cancer. 67:15-23). The present study was directed to investigate if beta 2-adrenergic agonists, which are known to favor skeletal muscle hypertrophy, could effectively antagonize the enhanced muscle protein breakdown in this cancer cachexia model. One such agent, i.e., clenbuterol, indeed largely prevented skeletal muscle waste in AH-130-bearing rats by restoring protein degradative rates close to control values. This normalization of protein breakdown rates was achieved through a decrease of the hyperactivation of the ATP-ubiquitin-dependent proteolytic pathway, as previously demonstrated in our laboratory (Llovera, M., C. García-Martínez, N. Agell, M. Marzábal, F. J. López-Soriano, and J. M. Argilés. 1994. FEBS (Fed. Eur. Biochem. Soc.) Lett. 338:311-318). By contrast, the drug did not exert any measurable effect on various parenchymal organs, nor did it modify the plasma level of corticosterone and insulin, which were increased and decreased, respectively, in the tumor hosts. The present data give new insights into the mechanisms by which clenbuterol exerts its preventive effect on muscle protein waste and seem to warrant the implementation of experimental protocols involving the use of clenbuterol or alike drugs in the treatment of pathological states involving TPH, particularly in skeletal muscle and heart, such as in the present model of cancer cachexia.

Muscle protein waste in tumor-bearing rats is effectively antagonized by a beta 2-adrenergic agonist (clenbuterol). Role of the ATP-ubiquitin-dependent proteolytic pathway.

Tissue protein hypercatabolism (TPH) is a most important feature in cancer cachexia, particularly with regard to the skeletal muscle. The rat ascites hepatoma Yoshida AH-130 is a very suitable model system for studying the mechanisms involved in the processes that lead to tissue depletion, since it induces in the host a rapid and progressive muscle waste mainly due to TPH (Tessitore, L., G. Bonelli, and F. M. Baccino. 1987. Biochem. J. 241:153-159). Detectable plasma levels of tumor necrosis factor-alpha associated with marked perturbations in the hormonal homeostasis have been shown to concur in forcing metabolism into a catabolic setting (Tessitore, L., P. Costelli, and F. M. Baccino. 1993. Br. J. Cancer. 67:15-23). The present study was directed to investigate if beta 2-adrenergic agonists, which are known to favor skeletal muscle hypertrophy, could effectively antagonize the enhanced muscle protein breakdown in this cancer cachexia model. One such agent, i.e., clenbuterol, indeed largely prevented skeletal muscle waste in AH-130-bearing rats by restoring protein degradative rates close to control values. This normalization of protein breakdown rates was achieved through a decrease of the hyperactivation of the ATP-ubiquitin-dependent proteolytic pathway, as previously demonstrated in our laboratory (Llovera, M., C. García-Martínez, N. Agell, M. Marzábal, F. J. López-Soriano, and J. M. Argilés. 1994. FEBS (Fed. Eur. Biochem. Soc.) Lett. 338:311-318). By contrast, the drug did not exert any measurable effect on various parenchymal organs, nor did it modify the plasma level of corticosterone and insulin, which were increased and decreased, respectively, in the tumor hosts. The present data give new insights into the mechanisms by which clenbuterol exerts its preventive effect on muscle protein waste and seem to warrant the implementation of experimental protocols involving the use of clenbuterol or alike drugs in the treatment of pathological states involving TPH, particularly in skeletal muscle and heart, such as in the present model of cancer cachexia.

Muscle protein waste in tumor-bearing rats is effectively antagonized by a beta 2-adrenergic agonist (clenbuterol). Role of the ATP-ubiquitin-dependent proteolytic pathway.

Tissue protein hypercatabolism (TPH) is a most important feature in cancer cachexia, particularly with regard to the skeletal muscle. The rat ascites hepatoma Yoshida AH-130 is a very suitable model system for studying the mechanisms involved in the processes that lead to tissue depletion, since it induces in the host a rapid and progressive muscle waste mainly due to TPH (Tessitore, L., G. Bonelli, and F. M. Baccino. 1987. Biochem. J. 241:153-159). Detectable plasma levels of tumor necrosis factor-alpha associated with marked perturbations in the hormonal homeostasis have been shown to concur in forcing metabolism into a catabolic setting (Tessitore, L., P. Costelli, and F. M. Baccino. 1993. Br. J. Cancer. 67:15-23). The present study was directed to investigate if beta 2-adrenergic agonists, which are known to favor skeletal muscle hypertrophy, could effectively antagonize the enhanced muscle protein breakdown in this cancer cachexia model. One such agent, i.e., clenbuterol, indeed largely prevented skeletal muscle waste in AH-130-bearing rats by restoring protein degradative rates close to control values. This normalization of protein breakdown rates was achieved through a decrease of the hyperactivation of the ATP-ubiquitin-dependent proteolytic pathway, as previously demonstrated in our laboratory (Llovera, M., C. García-Martínez, N. Agell, M. Marzábal, F. J. López-Soriano, and J. M. Argilés. 1994. FEBS (Fed. Eur. Biochem. Soc.) Lett. 338:311-318). By contrast, the drug did not exert any measurable effect on various parenchymal organs, nor did it modify the plasma level of corticosterone and insulin, which were increased and decreased, respectively, in the tumor hosts. The present data give new insights into the mechanisms by which clenbuterol exerts its preventive effect on muscle protein waste and seem to warrant the implementation of experimental protocols involving the use of clenbuterol or alike drugs in the treatment of pathological states involving TPH, particularly in skeletal muscle and heart, such as in the present model of cancer cachexia.

Comment on 'The Molecular Evolutionary Patterns of the Insulin/FOXO Signaling Pathway'

Letter to the Editor on Wang M, Wang Q, Wang Z, Zhang X, Pan Y. The molecular evolutionary patterns of the insulin/FOXO signaling pathway

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment and innate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints. Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection are frequent on genes that are known to be at bifurcation points, and that are identified as being in key position by a network-level analysis such as MGAT3 and GCS1. These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

Pint lincRNA connects the p53 pathway with epigenetic silencing by the Polycomb repressive complex 2

BACKGROUND: The p53 transcription factor is located at the core of a complex wiring of signaling pathways that are critical for the preservation of cellular homeostasis. Only recently it has become clear that p53 regulates the expression of several long intergenic noncoding RNAs (lincRNAs). However, relatively little is known about the role that lincRNAs play in this pathway. RESULTS: Here we characterize a lincRNA named Pint (p53 induced noncoding transcript). We show that Pint is a ubiquitously expressed lincRNA that is finely regulated by p53. In mouse cells, Pint promotes cell proliferation and survival by regulating the expression of genes of the TGF-β, MAPK and p53 pathways. Pint is a nuclear lincRNA that directly interacts with the Polycomb repressive complex 2 (PRC2), and is required for PRC2 targeting of specific genes for H3K27 tri-methylation and repression. Furthermore, Pint functional activity is highly dependent on PRC2 expression. We have also identified Pint human ortholog (PINT), which presents suggestive analogies with the murine lincRNA. PINT is similarly regulated by p53, and its expression significantly correlates with the same cellular pathways as the mouse ortholog, including the p53 pathway. Interestingly, PINT is downregulated in colon primary tumors, while its overexpression inhibits the proliferation of tumor cells, suggesting a possible role as tumor suppressor. CONCLUSIONS: Our results reveal a p53 autoregulatory negative mechanism where a lincRNA connects p53 activation with epigenetic silencing by PRC2. Additionally, we show analogies and differences between the murine and human orthologs, identifying a novel tumor suppressor candidate lincRNA.

In silico pathway reconstruction: Iron-sulfur cluster biogenesis in Saccharomyces cerevisiae

Background: Current advances in genomics, proteomics and other areas of molecular biology make the identification and reconstruction of novel pathways an emerging area of great interest. One such class of pathways is involved in the biogenesis of Iron-Sulfur Clusters (ISC). Results: Our goal is the development of a new approach based on the use and combination of mathematical, theoretical and computational methods to identify the topology of a target network. In this approach, mathematical models play a central role for the evaluation of the alternative network structures that arise from literature data-mining, phylogenetic profiling, structural methods, and human curation. As a test case, we reconstruct the topology of the reaction and regulatory network for the mitochondrial ISC biogenesis pathway in S. cerevisiae. Predictions regarding how proteins act in ISC biogenesis are validated by comparison with published experimental results. For example, the predicted role of Arh1 and Yah1 and some of the interactions we predict for Grx5 both matches experimental evidence. A putative role for frataxin in directly regulating mitochondrial iron import is discarded from our analysis, which agrees with also published experimental results. Additionally, we propose a number of experiments for testing other predictions and further improve the identification of the network structure. Conclusion: We propose and apply an iterative in silico procedure for predictive reconstruction of the network topology of metabolic pathways. The procedure combines structural bioinformatics tools and mathematical modeling techniques that allow the reconstruction of biochemical networks. Using the Iron Sulfur cluster biogenesis in S. cerevisiae as a test case we indicate how this procedure can be used to analyze and validate the network model against experimental results. Critical evaluation of the obtained results through this procedure allows devising new wet lab experiments to confirm its predictions or provide alternative explanations for further improving the models.

Saccharomyces cerevisiae Grx6 and Grx7 are monothiol glutaredoxins associated with the early secretory pathway

Saccharomyces cerevisiae Grx6 and Grx7 are two monothiol glutaredoxins whose active-site sequences (CSYS and CPYS, respectively) are reminiscent of the CPYC active-site sequence of classical dithiol glutaredoxins. Both proteins contain an N-terminal transmembrane domain which is responsible for their association to membranes of the early secretory pathway vesicles, facing the luminal side. Thus, Grx6 localizes at the endoplasmic reticulum and Golgi compartments, while Grx7 is mostly at the Golgi. Expression of GRX6 is modestly upregulated by several stresses (calcium, sodium, and peroxides) in a manner dependent on the Crz1-calcineurin pathway. Some of these stresses also upregulate GRX7 expression under the control of the Msn2/4 transcription factor. The N glycosylation inhibitor tunicamycin induces the expression of both genes along with protein accumulation. Mutants lacking both glutaredoxins display reduced sensitivity to tunicamycin, although the drug is still able to manifest its inhibitory effect on a reporter glycoprotein. Grx6 and Grx7 have measurable oxidoreductase activity in vivo, which is increased in the presence of tunicamycin. Both glutaredoxins could be responsible for the regulation of the sulfhydryl oxidative state at the oxidant conditions of the early secretory pathway vesicles. However, the differences in location and expression responses against stresses suggest that their functions are not totally overlapping.

GD3 Synthase Overexpression Sensitizes Hepatocarcinoma Cells to Hypoxia and Reduces Tumor Growth by Suppressing the cSrc/NF-κB Survival Pathway

Abstract Background: Hypoxia-mediated HIF-1a stabilization and NF-kB activation play a key role in carcinogenesis by fostering cancer cell survival, angiogenesis and tumor invasion. Gangliosides are integral components of biological membranes with an increasingly recognized role as signaling intermediates. In particular, ganglioside GD3 has been characterized as a proapoptotic lipid effector by promoting cell death signaling and suppression of survival pathways. Thus, our aim was to analyze the role of GD3 in hypoxia susceptibility of hepatocarcinoma cells and in vivo tumor growth. Methodology/Principal Findings: We generated and characterized a human hepatocarcinoma cell line stably expressing GD3 synthase (Hep3B-GD3), which catalyzes the synthesis of GD3 from GM3. Despite increased GD3 levels (2-3 fold), no significant changes in cell morphology or growth were observed in Hep3B-GD3 cells compared to wild type Hep3B cells under normoxia. However, exposure of Hep3B-GD3 cells to hypoxia (2% O2) enhanced reactive oxygen species (ROS) generation, resulting in decreased cell survival, with similar findings observed in Hep3B cells exposed to increasing doses of exogenous GD3. In addition, hypoxia-induced c-Src phosphorylation at tyrosine residues, NF-kB activation and subsequent expression of Mn-SOD were observed in Hep3B cells but not in Hep3B-GD3 cells. Moreover, MnTBAP, an antioxidant with predominant SOD mimetic activity, reduced ROS generation, protecting Hep3B-GD3 cells from hypoxia-induced death. Finally, lower tumor growth, higher cell death and reduced Mn-SOD expression were observed in Hep3B-GD3 compared to Hep3B tumor xenografts. Conclusion: These findings underscore a role for GD3 in hypoxia susceptibility by disabling the c-Src/NF-kB survival pathway resulting in lower Mn-SOD expression, which may be of relevance in hepatocellular carcinoma therapy.

Hepatocyte nuclear factor-4 alpha regulates the human apolipoprotein AV gene: Identification of a novel response element and involvement in the control by peroxisome proliferator-activated receptor-gamma coactivator-1 alpha, AMP-activated protein kinase, and mitogen-activated protein kinase pathway

The recently discovered apolipoprotein AV (apoAV) gene has been reported to be a key player in modulating plasma triglyceride levels. Here we identify the hepatocyte nuclear factor-4 (HNF-4 ) as a novel regulator of human apoAV gene. Inhibition of HNF-4 expression by small interfering RNA resulted in down-regulation of apoAV. Deletion, mutagenesis, and binding assays revealed that HNF-4 directly regulates human apoAV promoter through DR1 [a direct repeat separated by one nucleotide (nt)], and via a novel element for HNF-4 consisting of an inverted repeat separated by 8 nt (IR8). In addition, we show that the coactivator peroxisome proliferator-activated receptor- coactivator-1 was capable of stimulating the HNF-4 -dependent transactivation of apoAV promoter. Furthermore, analyses in human hepatic cells demonstrated that AMP-activated protein kinase (AMPK) and the MAPK signaling pathway regulate human apoAV expression and suggested that this regulation may be mediated, at least in part, by changes in HNF-4 . Intriguingly, EMSAs and mice with a liver-specific disruption of the HNF-4 gene revealed a species-distinct regulation of apoAV by HNF-4 , which resembles that of a subset of HNF-4 target genes. Taken together, our data provide new insights into the binding properties and the modulation of HNF-4 and underscore the role of HNF-4 in regulating triglyceride metabolism.

Lipid phosphate phosphatase 3 participates in transport carrier formation and protein trafficking in the early secretory pathway

The inhibition of phosphatidic acid phosphatase (PAP) activity by propanolol indicates that diacylglycerol (DAG) is required for the formation of transport carriers at the Golgi and for retrograde trafficking to the ER. Here we report that the PAP2 family member lipid phosphate phosphatase 3 (LPP3, also known as PAP2b) localizes in compartments of the secretory pathway from ER export sites to the Golgi complex. The depletion of human LPP3: (i) reduces the number of tubules generated from the ER-Golgi intermediate compartment and the Golgi, with those formed from the Golgi being longer in LPP3-silenced cells than in control cells; (ii) impairs the Rab6-dependent retrograde transport of Shiga toxin subunit B from the Golgi to the ER, but not the anterograde transport of VSV-G or ssDsRed; and (iii) induces a high accumulation of Golgi-associated membrane buds. LPP3 depletion also reduces levels of de novo synthesized DAG and the Golgi-associated DAG contents. Remarkably, overexpression of a catalytically inactive form of LPP3 mimics the effects of LPP3 knockdown on Rab6-dependent retrograde transport. We conclude that LPP3 participates in the formation of retrograde transport carriers at the ER-Golgi interface, where it transitorily cycles, and during its route to the plasma membrane.

Pint lincRNA connects the p53 pathway with epigenetic silencing by the Polycomb repressive complex 2

BACKGROUND: The p53 transcription factor is located at the core of a complex wiring of signaling pathways that are critical for the preservation of cellular homeostasis. Only recently it has become clear that p53 regulates the expression of several long intergenic noncoding RNAs (lincRNAs). However, relatively little is known about the role that lincRNAs play in this pathway. RESULTS: Here we characterize a lincRNA named Pint (p53 induced noncoding transcript). We show that Pint is a ubiquitously expressed lincRNA that is finely regulated by p53. In mouse cells, Pint promotes cell proliferation and survival by regulating the expression of genes of the TGF-β, MAPK and p53 pathways. Pint is a nuclear lincRNA that directly interacts with the Polycomb repressive complex 2 (PRC2), and is required for PRC2 targeting of specific genes for H3K27 tri-methylation and repression. Furthermore, Pint functional activity is highly dependent on PRC2 expression. We have also identified Pint human ortholog (PINT), which presents suggestive analogies with the murine lincRNA. PINT is similarly regulated by p53, and its expression significantly correlates with the same cellular pathways as the mouse ortholog, including the p53 pathway. Interestingly, PINT is downregulated in colon primary tumors, while its overexpression inhibits the proliferation of tumor cells, suggesting a possible role as tumor suppressor. CONCLUSIONS: Our results reveal a p53 autoregulatory negative mechanism where a lincRNA connects p53 activation with epigenetic silencing by PRC2. Additionally, we show analogies and differences between the murine and human orthologs, identifying a novel tumor suppressor candidate lincRNA.

Three-step pathway engineering results in more incidence rate and higher emission of nerolidol and improved attraction of Diadegma semiclausum

The concentration and ratio of terpenoids in the headspace volatile blend of plants have a fundamental role in the communication of plants and insects. The sesquiterpene (E)-nerolidol is one of the important volatiles with effect on beneficial carnivores for biologic pest management in the field. To optimize de novo biosynthesis and reliable and uniform emission of (E)-nerolidol, we engineered different steps of the (E)-nerolidol biosynthesis pathway in Arabidopsis thaliana. Introduction of a mitochondrial nerolidol synthase gene mediates de novo emission of (E)-nerolidol and linalool. Co-expression of the mitochondrial FPS1 and cytosolic HMGR1 increased the number of emitting transgenic plants (incidence rate) and the emission rate of both volatiles. No association between the emission rate of transgenic volatiles and their growth inhibitory effect could be established. (E)-Nerolidol was to a large extent metabolized to non-volatile conjugates.