A systems biology approach for the study of cumulative oncogenes with applications to the MAPK signal transduction pathway

The Extracellular signal Regulated Kinase (ERK) pathway is one of the most well-studied signaling pathways in cell cycle regulation. Disruption in the normal functioning of this pathway is linked to many forms of cancer. In a previous study [D.K. Pant, A. Ghosh, Automated oncogene detection in complex protein networks, with applications to the MAPK signal transduction pathway, Biophys. Chem. 113 (2005) 275-288.], we developed a novel approach to predict single point mutations that are likely to cause cellular transformation in signaling transduction networks. We have extended this method to study disparate pair mutation in enzyme/protein interactions and in expression levels in signal transduction pathway and have applied it to the MAPK signaling pathway to study how synergistic or cooperative mutation within signaling networks acts in unison to cause malignant transformation. The method provides a quantitative ranking of the modifier pair of ERK activation. It is seen that the highest ranking single point mutations comprise the highest ranking pair mutations. We validate some of our results with experimental literature on multiple mutations. A second order sensitivity analysis scheme is additionally used to determine the effect of correlations among mutations at different sites in the pathways.     

Imaging individual mRNA molecules using multiple singly labeled probes

We describe a method for imaging individual mRNA molecules in fixed cells by probing each mRNA species with 48 or more short, singly labeled oligonucleotide probes. This makes each mRNA molecule visible as a computationally identifiable fluorescent spot by fluorescence microscopy. We demonstrate simultaneous detection of three mRNA species in single cells and mRNA detection in yeast, nematodes, fruit fly wing discs, and mammalian cell lines and neurons.     

Assigning Quantitative Function to Post-Translational Modifications Reveals Multiple Sites of Phosphorylation That Tune Yeast Pheromone Signaling Output

Cell signaling systems transmit information by post-translationally modifying signaling proteins, often via phosphorylation. While thousands of sites of phosphorylation have been identified in proteomic studies, the vast majority of sites have no known function. Assigning functional roles to the catalog of uncharacterized phosphorylation sites is a key research challenge. Here we present a general approach to address this challenge and apply it to a prototypical signaling pathway, the pheromone response pathway in Saccharomyces cerevisiae. The pheromone pathway includes a mitogen activated protein kinase (MAPK) cascade activated by a G-protein coupled receptor (GPCR). We used published mass spectrometry-based proteomics data to identify putative sites of phosphorylation on pheromone pathway components, and we used evolutionary conservation to assign priority to a list of candidate MAPK regulatory sites. We made targeted alterations in those sites, and measured the effects of the mutations on pheromone pathway output in single cells. Our work identified six new sites that quantitatively tuned system output. We developed simple computational models to find system architectures that recapitulated the quantitative phenotypes of the mutants. Our results identify a number of putative phosphorylation events that contribute to adjust the input-output relationship of this model eukaryotic signaling system. We believe this combined approach constitutes a general means not only to reveal modification sites required to turn a pathway on and off, but also those required for more subtle quantitative effects that tune pathway output. Our results suggest that relatively small quantitative influences from individual phosphorylation events endow signaling systems with plasticity that evolution may exploit to quantitatively tailor signaling outcomes.     

Platelet-derived growth factor receptor-mediated signal transduction from endosomes

Although accumulated evidence supports the concept of endosomal signaling of receptor tyrosine kinases, most results are generated from studies of epidermal growth factor receptor (EGFR). It is not clear whether the concept of endosomal signaling could be generally applied to the other receptor tyrosine kinases. For example, platelet-derived growth factor receptor (PDGFR) is very similar to EGFR in terms of both signaling and trafficking; however, little is known about the endosomal signaling of PDGFR. In this research, we applied the same approaches from our recent studies regarding EGFR endosomal signaling to investigate the endosomal signaling of PDGFR. We showed in this communication that we are able to establish a system that allows the specific activation of endosome-associated PDGFR without the activation of the plasma membrane-associated PDGFR and without disrupting the overall endocytosis pathway. By using this system, we showed that endosomal activation of PDGFR recruits various signaling proteins including Grb2, SHC, phospholipase C-gamma1, and the p85alpha subunit of phosphatidylinositol 3-kinase into endosomes and forms signaling complexes with PDGFR. We also showed that endosomal PDGFR signaling is sufficient to activate the major signaling pathways implicated in cell proliferation and survival. Moreover, we demonstrate that endosomal PDGFR signaling is sufficient to generate physiological output including cell proliferation and cell survival.     

Simultaneous visualization of neuronal protein and receptor mRNA

We describe a combined immunocytochemistry/in situ hybridization technique which allows for the simultaneous localization of protein and mRNA in a single cell. We have carried out these studies either on non-innervated skeletal myotubes or on myotubes which we have innervated with spinal cord explants or ciliary neurons. Our methods allowed us to detect acetylcholine receptor gene mRNA sequences which are expressed in low abundance within the cells and to determine the intracellular and intranuclear domains where these sequences are concentrated, as well as to identify neurons and their processes. The isotopic detection of RNA in combination with fluorescence microscopy produces high-resolution double-label images, with little background and good preservation of morphology, providing a powerful tool for detection of gene expression and protein content at the single-cell level.     

Parallel monitoring of RNA abundance,localization and compactness with correlative single molecule FISH on LR White embedded samples

Single mRNA molecules are frequently detected by single molecule fluorescence in situ hybridization (smFISH) using branched DNA technology. While providing strong and background-reduced signals, the method is inefficient in detecting mRNAs within dense structures, in monitoring mRNA compactness and in quantifying abundant mRNAs. To overcome these limitations, we have hybridized slices of high pressure frozen, freeze-substituted and LR White embedded cells (LR White smFISH). mRNA detection is physically restricted to the surface of the resin. This enables single molecule detection of RNAs with accuracy comparable to RNA sequencing, irrespective of their abundance, while at the same time providing spatial information on RNA localization that can be complemented with immunofluorescence and electron microscopy, as well as array tomography. Moreover, LR White embedding restricts the number of available probe pair recognition sites for each mRNA to a small subset. As a consequence, differences in signal intensities between RNA populations reflect differences in RNA structures, and we show that the method can be employed to determine mRNA compactness. We apply the method to answer some outstanding questions related to trans-splicing, RNA granules and mitochondrial RNA editing in single-cellular trypanosomes and we show an example of differential gene expression in the metazoan Caenorhabditis elegans.     

Three-dimensional imaging of the intracellular localization of growth hormone and prolactin and their mRNA using nanocrystal (Quantum dot) and confocal laser scanning microscopy techniques.

Semiconductor nanocrystals (Quantum dots, Qdots) have recently been used in biological research, because they do not fade on exposure to light, and they enable us to obtain multicolor imaging because of a narrow emission peak that can be excited via a single wavelength of light. There have been no reports of simultaneous localization of mRNA and protein using Qdots. We successfully applied these advantages of Qdot and confocal laser scanning microscopy (CLSM) to three-dimensional images of the intracellular localization of growth hormone and prolactin and to their mRNA. In situ hybridization and immunohistochemistry using Qdots combined with CLSM can optimally illustrate the relationship between protein and mRNA simultaneously in three dimensions. Such an approach enables us to visualize functional images of proteins in relation with mRNA synthesis and localization.     

The mitochondrial-apoptotic pathway is triggered in Xenopus mesoderm cells deprived of PDGF receptor signaling during gastrulation

Platelet-derived growth factor receptor (PDGFR) signaling is required for normal gastrulation in Xenopus laevis. Embryos deprived of PDGFR signaling develop with a range of gastrulation-specific defects including spina bifida, shortened anteroposterior axis, and reduced anterior structures. These defects arise because the involuting mesoderm fails to move appropriately. In this study, we determine that inhibition of PDGFR signaling causes prospective head mesoderm cells to appear in the blastocoel cavity at the onset of gastrulation, stage 10. These aberrant cells undergo apoptosis via the caspase 3 pathway at an embryonic checkpoint called the early gastrula transition (EGT). They are TUNEL-positive and have increased levels of caspase 3 activity compared to control embryos. Apoptotic death of these mesoderm cells can be prevented by co-injection of mRNA encoding Bcl-2 or by injection of either a general caspase inhibitor or a caspase 3-specific inhibitor. Prevention of cell death, however, is not sufficient to rescue gastrulation defects in these embryos. Based on these data, we propose that PDGFR signaling is necessary for survival of prospective head mesoderm cells, and also plays an essential role in the control of their cell movement during gastrulation.     

Application of an integrated physical and functional screening approach to identify inhibitors of the Wnt pathway

Large‐scale proteomic approaches have been used to study signaling pathways. However, identification of biologically relevant hits from a single screen remains challenging due to limitations inherent in each individual approach. To overcome these limitations, we implemented an integrated, multi‐dimensional approach and used it to identify Wnt pathway modulators. The LUMIER protein–protein interaction mapping method was used in conjunction with two functional screens that examined the effect of overexpression and siRNA‐mediated gene knockdown on Wnt signaling. Meta‐analysis of the three data sets yielded a combined pathway score (CPS) for each tested component, a value reflecting the likelihood that an individual protein is a Wnt pathway regulator. We characterized the role of two proteins with high CPSs, Ube2m and Nkd1. We show that Ube2m interacts with and modulates β‐catenin stability, and that the antagonistic effect of Nkd1 on Wnt signaling requires interaction with Axin, itself a negative pathway regulator. Thus, integrated physical and functional mapping in mammalian cells can identify signaling components with high confidence and provides unanticipated insights into pathway regulators.     

Sequence-Based Antigenic Change Prediction by a Sparse Learning Method Incorporating Co-Evolutionary Information

Rapid identification of influenza antigenic variants will be critical in selecting optimal vaccine candidates and thus a key to developing an effective vaccination program. Recent studies suggest that multiple simultaneous mutations at antigenic sites accumulatively enhance antigenic drift of influenza A viruses. However, pre-existing methods on antigenic variant identification are based on analyses from individual sites. Because the impacts of these co-evolved sites on influenza antigenicity may not be additive, it will be critical to quantify the impact of not only those single mutations but also multiple simultaneous mutations or co-evolved sites. Here, we developed and applied a computational method, AntigenCO, to identify and quantify both single and co-evolutionary sites driving the historical antigenic drifts. AntigenCO achieved an accuracy of up to 90.05% for antigenic variant prediction, significantly outperforming methods based on single sites. AntigenCO can be useful in antigenic variant identification in influenza surveillance.     

Tyrosine phosphorylation mapping of the epidermal growth factor receptor signaling pathway.

Phosphorylation is one of the most common forms of protein modification. The most frequent targets for protein phosphorylation in eukaryotes are serine and threonine residues, although tyrosine residues also undergo phosphorylation. Many of the currently applied methods for the detection and localization of protein phosphorylation sites are mass spectrometry-based and are biased against the analysis of tyrosine-phosphorylated residues because of the stability and low reactivity of phosphotyrosines. To overcome this lack of sensitive methods for the detection of phosphotyrosine-containing peptides, we have recently developed a method that is not affected by the more predominant threonine or serine phosphorylation within cells. It is based on the specific detection of immonium ion of phosphotyrosine at 216.043 Da and does not require prior knowledge of the protein sequence. In this report, we describe the first application of this new method in a proteomic strategy. Using anti-phosphotyrosine antibodies for immunoprecipitation and one-dimensional gel electrophoresis, we have identified 10 proteins in the epidermal growth factor receptor signaling pathway, of which 8 have been shown previously to be involved in epidermal growth factor signaling. Most importantly, in addition to several known tyrosine phosphorylation sites, we have identified five novel sites on SHIP-2, Hrs, Cbl, STAM, and STAM2, most of which were not predicted to be phosphorylated. Because of its sensitivity and selectivity, this approach will be useful in proteomic approaches to study tyrosine phosphorylation in a number of signal transduction pathways.     

TLR-4 and sustained calcium agonists synergistically produce eicosanoids independent of protein synthesis in RAW264.7 cells

Arachidonic acid is released by phospholipase A(2) and converted into hundreds of distinct bioactive mediators by a variety of cyclooxygenases (COX), lipoxygenases (LO), and cytochrome P450s. Because of the size and diversity of the eicosanoid class of signaling molecules produced, a thorough and systematic investigation of these biological processes requires the simultaneous quantitation of a large number of eicosanoids in a single analysis. We have developed a robust liquid chromatography/tandem mass spectrometry method that can identify and quantitate over 60 different eicosanoids in a single analysis, and we applied it to agonist-stimulated RAW264.7 murine macrophages. Fifteen different eicosanoids produced through COX and 5-LO were detected either intracellularly or in the media following stimulation with 16 different agonists of Toll-like receptors (TLR), G protein-coupled receptors, and purinergic receptors. No significant differences in the COX metabolite profiles were detected using the different agonists; however, we determined that only agonists creating a sustained Ca(2+) influx were capable of activating the 5-LO pathway in these cells. Synergy between Ca(2+) and TLR pathways was detected and discovered to be independent of NF-kappaB-induced protein synthesis. This demonstrates that TLR induction of protein synthesis and priming for enhanced phospholipase A(2)-mediated eicosanoid production work through two distinct pathways.     

LPA-induced mutually exclusive subcellular localization of active RhoA and Arp2 mRNA revealed by sequential FRET and FISH

We previously demonstrated that mRNAs for the subunits of the Arp2/3 complex localize to protrusions in fibroblasts (Mingle et al. in J Cell Sci 118:2425–2433, 2005). However, the signaling pathway that regulates Arp2/3 complex mRNA localization remains unknown. In this study we have identified lysophosphatidic acid (LPA) as a potent inducer of Arp2 mRNA localization to protrusions in fibroblasts via the RhoA-ROCK pathway. As RhoA is known to be activated locally in the cells, we sought to understand how spatial activation of Rho affects Arp2 mRNA localization. By sequentially performing fluorescence resonance energy transfer (FRET) and fluorescence in situ hybridization (FISH), we have visualized active RhoA and Arp2 mRNA in the same cells. Upon LPA stimulation, approximately two times more cells than those in the serum-free medium showed mutually exclusive localization of active RhoA and Arp2 mRNA. These results demonstrate the importance of localized activation of Rho in Arp2 mRNA localization and provide new insights as to how Rho regulates Arp2/3 complex mRNA localization. To our best knowledge, this is the first report in which FRET and FISH are combined to detect localized protein activity and mRNA in the same cells. This method should be easily adopted for the detection of other fluorescence protein based biosensors and DNA/RNA in the same cells.     

Double and triple in situ hybridization for coexpression studies: combined fluorescent and chromogenic detection of neuropeptide Y (NPY) and serotonin receptor subtype mRNAs expressed at different abundance levels

Multiple fluorescence in situ hybridization is the method of choice for studies aimed at determining simultaneous production of signal transduction molecules and neuromodulators in neurons. In our analyses of the monoamine receptor mRNA expression of peptidergic neurons in the rat telencephalon, double tyramide-signal-amplified fluorescence in situ hybridization delivered satisfactory results for coexpression analysis of neuropeptide Y (NPY) and serotonin receptor 2C (5-HT2C) mRNA, a receptor subtype expressed at high-to-moderate abundance in the regions analyzed. However, expression of 5-HT1A mRNA, which is expressed at comparatively low abundance in many telencephalic areas, could not be unequivocally identified in NPY mRNA-reactive neurons due to high background and poor signal-to-noise ratio in fluorescent receptor mRNA detections. Parallel chromogenic in situ hybridization provided clear labeling for 5-HT1A mRNA and additionally offered the possibility to monitor the chromogen deposition at regular time intervals to determine the optimal signal-to-noise ratio. We first developed a double labeling protocol combining fluorescence and chromogenic in situ hybridization and subsequently expanded this variation to combine double fluorescence and chromogenic in situ hybridization for triple labelings. With this method, we documented expression of 5-HT2C and/or 5-HT1A in subpopulations of telencephalic NPY-producing neurons. The method developed in the present study appears suitable for conventional light and fluorescence microscopy, combines advantages of fluorescence and chromogenic in situ hybridization protocols and thus provides a reliable non-radioactive alternative to previously published multiple labeling methods for coexpression analyses in which one mRNA species requires highly sensitive detection.