Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Semaphorin3A Signaling Mediated by Fyn-dependent Tyrosine Phosphorylation of Collapsin Response Mediator Protein 2 at Tyrosine 32

Collapsin response mediator protein 2 (CRMP2) is an intracellular protein that mediates signaling of Semaphorin3A (Sema3A), a repulsive axon guidance molecule. Fyn, a Src-type tyrosine kinase, is involved in the Sema3A signaling. However, the relationship between CRMP2 and Fyn in this signaling pathway is still unknown. In our research, we demonstrated that Fyn phosphorylated CRMP2 at Tyr³² residues in HEK293T cells. Immunohistochemical analysis using a phospho-specific antibody at Tyr³² of CRMP showed that Tyr³²-phosphorylated CRMP was abundant in the nervous system, including dorsal root ganglion neurons, the molecular and Purkinje cell layer of adult cerebellum, and hippocampal fimbria. Overexpression of a nonphosphorylated mutant (Tyr³² to Phe³²) of CRMP2 in dorsal root ganglion neurons interfered with Sema3A-induced growth cone collapse response. These results suggest that Fyn-dependent phosphorylation of CRMP2 at Tyr³² is involved in Sema3A signaling.Collapsin response mediator proteins (CRMPs)⁴ have been identified as intracellular proteins that mediate Semaphorin3A (Sema3A) signaling in the nervous system (1). CRMP2 is one of the five members of the CRMP family. CRMPs also mediate signal transduction of NT3, Ephrin, and Reelin (2–4). CRMPs interact with several intracellular molecules, including tubulin, Numb, kinesin1, and Sra1 (5–8). CRMPs are involved in axon guidance, axonal elongation, cell migration, synapse maturation, and the generation of neuronal polarity (1, 2, 4, 5).CRMP family proteins are known to be the major phosphoproteins in the developing brain (1, 9). CRMP2 is phosphorylated by several Ser/Thr kinases, such as Rho kinase, cyclin-dependent kinase 5 (Cdk5), and glycogen synthase kinase 3β (GSK3β) (2, 10–13). The phosphorylation sites of CRMP2 by these kinases are clustered in the C terminus and have already been identified. Rho kinase phosphorylates CRMP2 at Thr⁵⁵⁵ (10). Cdk5 phosphorylates CRMP2 at Ser⁵²², and this phosphorylation is essential for sequential phosphorylations by GSK3β at Ser⁵¹⁸, Thr⁵¹⁴, and Thr⁵⁰⁹ (2, 11–13). These phosphorylations disrupt the interaction of CRMP2 with tubulin or Numb (2, 3, 13). The sequential phosphorylation of CRMP2 by Cdk5 and GSK3β is an essential step in Sema3A signaling (11, 13). Furthermore, the neurofibrillary tangles in the brains of people with Alzheimer disease contain hyperphosphorylated CRMP2 at Thr⁵⁰⁹, Ser⁵¹⁸, and Ser⁵²² (14, 15).CRMPs are also substrates of several tyrosine kinases. The phosphorylation of CRMP2 by Fes/Fps and Fer has been shown to be involved in Sema3A signaling (16, 17). Phosphorylation of CRMP2 at Tyr⁴⁷⁹ by a Src family tyrosine kinase Yes regulates CXCL12-induced T lymphocyte migration (18). We reported previously that Fyn is involved in Sema3A signaling (19). Fyn associates with PlexinA2, one of the components of the Sema3A receptor complex. Fyn also activates Cdk5 through the phosphorylation at Tyr¹⁵ of Cdk5 (19). In dorsal root ganglion (DRG) neurons from fyn-deficient mice, Sema3A-induced growth cone collapse response is attenuated compared with control mice (19). Furthermore, we recently found that Fyn phosphorylates CRMP1 and that this phosphorylation is involved in Reelin signaling (4). Although it has been shown that CRMP2 is involved in Sema3A signaling (1, 11, 13), the relationship between Fyn and CRMP2 in Sema3A signaling and the tyrosine phosphorylation site(s) of CRMPs remain unknown.Here, we show that Fyn phosphorylates CRMP2 at Tyr³². Using a phospho-specific antibody against Tyr³², we determined that the residue is phosphorylated in vivo. A nonphosphorylated mutant CRMP2Y32F inhibits Sema3A-induced growth cone collapse. These results indicate that tyrosine phosphorylation by Fyn at Tyr³² is involved in Sema3A signaling.     

Systems-wide Analysis of K-Ras,Cdc42, and PAK4 Signaling by Quantitative Phosphoproteomics

Although K-Ras, Cdc42, and PAK4 signaling are commonly deregulated in cancer, only a few studies have sought to comprehensively examine the spectrum of phosphorylation-mediated signaling downstream of each of these key signaling nodes. In this study, we completed a label-free quantitative analysis of oncogenic K-Ras, activated Cdc42, and PAK4-mediated phosphorylation signaling, and report relative quantitation of 2152 phosphorylated peptides on 1062 proteins. We define the overlap in phosphopeptides regulated by K-Ras, Cdc42, and PAK4, and find that perturbation of these signaling components affects phosphoproteins associated with microtubule depolymerization, cytoskeletal organization, and the cell cycle. These findings provide a resource for future studies to characterize novel targets of oncogenic K-Ras signaling and validate biomarkers of PAK4 inhibition.The Ras oncoproteins are small monomeric GTPases that transduce mitogenic signals from cell surface receptor tyrosine kinases (RTKs) to intracellular serine/threonine kinases. Approximately thirty percent of human tumors harbor a somatic gain-of-function mutation in one of three RAS genes, resulting in the constitutive activation of Ras signaling and the aberrant hyperactivation of growth-promoting effector pathways (1). Designing therapeutic agents that directly target Ras has been challenging (2, 3), and thus clinical development efforts have focused on targeting effector pathways downstream of Ras. The Raf-MEK-ERK and PI3K-Akt effector pathways have been extensively studied and several small molecule inhibitors targeting these pathways are currently under clinical evaluation (4, 5). However, biochemical studies and mouse models indicate that several additional effector pathways are essential for Ras-driven transformation and tumorigenesis (6–11). Hence, a comprehensive characterization of these effector pathways may reveal additional druggable targets.The Rho GTPase Cdc42 lies downstream of Ras (12–14) and regulates many cellular processes that are commonly perturbed in cancer, including migration, polarization, and proliferation (15) (Fig. 1A). Importantly, Cdc42 is overexpressed in several types of human cancer (16–20) and is required for Ras-driven cellular transformation (13, 21, 22). Recent studies show that genetic ablation of Cdc42 impairs Ras-driven tumorigenesis (13), indicating the potential of Cdc42 and its effectors as drug targets in Ras mutant tumors.Open in a separate windowFig. 1.Experimental workflow. A, K-Ras is a small GTPase that regulates the activity of a variety of downstream proteins including the Rho GTPase Cdc42. The PAK4 serine/threonine kinase is a direct effector of Cdc42 and regulates actin reorganization, microtubule stability, and cell polarity. B, To measure large-scale phosphorylation changes induced by constitutive K-Ras or Cdc42 signaling or PAK4 ablation, the quantitative label-free PTMscan® approach was employed (Cell Signaling Technology). Briefly, for each condition extracted proteins were digested with trypsin and separated from non-peptide material by solid-phase extraction with Sep-Pak C18 cartridges. Three phosphorylation motif antibodies were used serially to isolate phosphorylated peptides in independent immunoaffinity purifications (CDK substrate motif [K R]-pS-P-X-[K R], CK substrate motif pT-[D E]-X-[D E], PKD substrate motif l-X-R-X-X-p[S T]). The samples were run in duplicate and tandem mass spectra were collected with an LTQ-Orbitrap hybrid mass spectrometer. pLPC is an empty vector control.In particular, the p21-activated kinases (PAKs) are Cdc42 effectors that have generated significant interest (23, 24), as they are central components of key oncogenic signaling pathways and regulate cytoskeletal organization, cell migration, and nuclear signaling (25). The PAK family is comprised of six members and is subdivided into two groups (Groups I and II) based on sequence and structural homology. Group I PAKs (PAK1–3) are relatively well characterized, however, much less is known regarding the function and regulation of Group II PAKs (PAK4–6). The kinase domains of Group I and II PAKs share only about 50% identity, suggesting the two groups may recognize distinct substrates and govern unique cellular processes (26).The Group II PAK family member PAK4 is of particular interest as it is overexpressed or genetically amplified in several lung, colon, prostate, pancreas, and breast tumor cell lines and samples (26–30). Furthermore, functional studies have implicated PAK4 in cell transformation, cell invasion, and migration (27, 31). Xenograft studies in athymic mice show an important role for PAK4 in mediating Cdc42- or K-Ras-driven tumor formation, highlighting a critical role for Pak4 downstream of these GTPases (32). Given its roles in transformation, tumorigenesis, and oncogenic signaling, there is significant interest in targeting PAK4 therapeutically (23). PAK4 binds and phosphorylates several proteins involved in cytoskeletal organization and apoptosis, including Lim domain kinase 1 (LIMK1) (33), guanine nucleotide exchange factor-H1 (GEF-H1) (34), Raf-1 (35), and Bad (36). However, the Group I PAK family member PAK1 also phosphorylates several of these PAK4 targets (37). Thus, there remains a need to identify robust and selective pharmacodynamic biomarkers for PAK4 inhibition.Despite the importance of PAK4 and its upstream regulators in cancer development, few studies have sought to comprehensively characterize the spectrum of K-Ras, Cdc42, or PAK4 mediated phosphorylation signaling (37–39). Recent developments in mass spectrometry allow the in-depth identification and quantitation of thousands of phosphorylation sites (40–43). The majority of large-scale efforts have aimed to identify the basal phosphoproteomes of different species (44, 45) or tissues (46) to characterize global steady-state phosphorylation. However, this methodology can also be applied to quantify perturbed phosphorylation regulation in cancer signaling pathways (40, 47–49), and has the potential to reveal novel biomarkers of oncogenic signaling.In this study, we completed a label-free quantitative analysis of K-Ras, Cdc42, and PAK4 phosphorylation signaling using the PTMScan® method, which has proven as robust and reproducible quantitation technology (50, 51). We quantified phosphorylation levels in wild-type and PAK4 knockout NIH3T3 cells expressing oncogenic K-Ras, activated Cdc42, or an empty vector control to elucidate the molecular pathways and functions modulated by these key signaling proteins. We report relative quantitation of 2152 phosphorylated peptides on 1062 proteins among the different conditions, and find that many of the regulated phosphoproteins are associated with microtubule depolymerization, cytoskeletal organization, and the cell cycle. To our knowledge, our study is the first to examine the overlap among signaling networks regulated by K-Ras, Cdc42, and PAK4, and provides a resource for future studies to further interrogate the perturbation of this signaling pathway.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

High-definition De Novo Sequencing of Crustacean Hyperglycemic Hormone (CHH)-family Neuropeptides

A complete understanding of the biological functions of large signaling peptides (>4 kDa) requires comprehensive characterization of their amino acid sequences and post-translational modifications, which presents significant analytical challenges. In the past decade, there has been great success with mass spectrometry-based de novo sequencing of small neuropeptides. However, these approaches are less applicable to larger neuropeptides because of the inefficient fragmentation of peptides larger than 4 kDa and their lower endogenous abundance. The conventional proteomics approach focuses on large-scale determination of protein identities via database searching, lacking the ability for in-depth elucidation of individual amino acid residues. Here, we present a multifaceted MS approach for identification and characterization of large crustacean hyperglycemic hormone (CHH)-family neuropeptides, a class of peptide hormones that play central roles in the regulation of many important physiological processes of crustaceans. Six crustacean CHH-family neuropeptides (8–9.5 kDa), including two novel peptides with extensive disulfide linkages and PTMs, were fully sequenced without reference to genomic databases. High-definition de novo sequencing was achieved by a combination of bottom-up, off-line top-down, and on-line top-down tandem MS methods. Statistical evaluation indicated that these methods provided complementary information for sequence interpretation and increased the local identification confidence of each amino acid. Further investigations by MALDI imaging MS mapped the spatial distribution and colocalization patterns of various CHH-family neuropeptides in the neuroendocrine organs, revealing that two CHH-subfamilies are involved in distinct signaling pathways.Neuropeptides and hormones comprise a diverse class of signaling molecules involved in numerous essential physiological processes, including analgesia, reward, food intake, learning and memory (1). Disorders of the neurosecretory and neuroendocrine systems influence many pathological processes. For example, obesity results from failure of energy homeostasis in association with endocrine alterations (2, 3). Previous work from our lab used crustaceans as model organisms found that multiple neuropeptides were implicated in control of food intake, including RFamides, tachykinin related peptides, RYamides, and pyrokinins (4–6).Crustacean hyperglycemic hormone (CHH)¹ family neuropeptides play a central role in energy homeostasis of crustaceans (7–17). Hyperglycemic response of the CHHs was first reported after injection of crude eyestalk extract in crustaceans. Based on their preprohormone organization, the CHH family can be grouped into two sub-families: subfamily-I containing CHH, and subfamily-II containing molt-inhibiting hormone (MIH) and mandibular organ-inhibiting hormone (MOIH). The preprohormones of the subfamily-I have a CHH precursor related peptide (CPRP) that is cleaved off during processing; and preprohormones of the subfamily-II lack the CPRP (9). Uncovering their physiological functions will provide new insights into neuroendocrine regulation of energy homeostasis.Characterization of CHH-family neuropeptides is challenging. They are comprised of more than 70 amino acids and often contain multiple post-translational modifications (PTMs) and complex disulfide bridge connections (7). In addition, physiological concentrations of these peptide hormones are typically below picomolar level, and most crustacean species do not have available genome and proteome databases to assist MS-based sequencing.MS-based neuropeptidomics provides a powerful tool for rapid discovery and analysis of a large number of endogenous peptides from the brain and the central nervous system. Our group and others have greatly expanded the peptidomes of many model organisms (3, 18–33). For example, we have discovered more than 200 neuropeptides with several neuropeptide families consisting of as many as 20–40 members in a simple crustacean model system (5, 6, 25–31, 34). However, a majority of these neuropeptides are small peptides with 5–15 amino acid residues long, leaving a gap of identifying larger signaling peptides from organisms without sequenced genome. The observed lack of larger size peptide hormones can be attributed to the lack of effective de novo sequencing strategies for neuropeptides larger than 4 kDa, which are inherently more difficult to fragment using conventional techniques (34–37). Although classical proteomics studies examine larger proteins, these tools are limited to identification based on database searching with one or more peptides matching without complete amino acid sequence coverage (36, 38).Large populations of neuropeptides from 4–10 kDa exist in the nervous systems of both vertebrates and invertebrates (9, 39, 40). Understanding their functional roles requires sufficient molecular knowledge and a unique analytical approach. Therefore, developing effective and reliable methods for de novo sequencing of large neuropeptides at the individual amino acid residue level is an urgent gap to fill in neurobiology. In this study, we present a multifaceted MS strategy aimed at high-definition de novo sequencing and comprehensive characterization of the CHH-family neuropeptides in crustacean central nervous system. The high-definition de novo sequencing was achieved by a combination of three methods: (1) enzymatic digestion and LC-tandem mass spectrometry (MS/MS) bottom-up analysis to generate detailed sequences of proteolytic peptides; (2) off-line LC fractionation and subsequent top-down MS/MS to obtain high-quality fragmentation maps of intact peptides; and (3) on-line LC coupled to top-down MS/MS to allow rapid sequence analysis of low abundance peptides. Combining the three methods overcomes the limitations of each, and thus offers complementary and high-confidence determination of amino acid residues. We report the complete sequence analysis of six CHH-family neuropeptides including the discovery of two novel peptides. With the accurate molecular information, MALDI imaging and ion mobility MS were conducted for the first time to explore their anatomical distribution and biochemical properties.     

Tyrosine Phosphorylation of Growth Factor Receptor-bound Protein-7 by Focal Adhesion Kinase in the Regulation of Cell Migration, Proliferation, and Tumorigenesis

We have previously reported that growth factor receptor-bound protein-7 (Grb7), an Src-homology 2 (SH2)-containing adaptor protein, enables interaction with focal adhesion kinase (FAK) to regulate cell migration in response to integrin activation. To further elucidate the signaling events mediated by FAK·Grb7 complexes in promoting cell migration and other cellular functions, we firstly examined the phos pho ryl a ted tyrosine site(s) of Grb7 by FAK using an in vivo mutagenesis. We found that FAK was capable of phos pho rylating at least 2 of 12 tyrosine residues within Grb7, Tyr-188 and Tyr-338. Moreover, mutations converting the identified Tyr to Phe inhibited integrin-dependent cell migration as well as impaired cell proliferation but not survival compared with the wild-type control. Interestingly, the above inhibitory effects caused by the tyrosine phos pho ryl a tion-deficient mutants are probably attributed to their down-regulation of phospho-Tyr-397 of FAK, thereby implying a mechanism by competing with wild-type Grb7 for binding to FAK. Consequently, these tyrosine phos pho ryl a tion-deficient mutants evidently altered the phospho-Tyr-118 of paxillin and phos pho ryl a tion of ERK1/2 but less on phospho-Ser-473 of AKT, implying their involvement in the FAK·Grb7-mediated cellular functions. Additionally, we also illustrated that the formation of FAK·Grb7 complexes and Grb7 phos pho ryl a tion by FAK in an integrin-dependent manner were essential for cell migration, proliferation and anchorage-independent growth in A431 epidermal carcinoma cells, indicating the importance of FAK·Grb7 complexes in tumorigenesis. Our data provide a better understanding on the signal transduction event for FAK·Grb7-mediated cellular functions as well as to shed light on a potential therapeutic in cancers.Growth factor receptor bound protein-7 (Grb7)² is initially identified as a SH2 domain-containing adaptor protein bound to the activated EGF receptor (1). Grb7 is composed of an N-terminal proline-rich region, following a putative RA (Ras-associating) domain and a central PH (pleckstrin homology) domain and a BPS motif (between PH and SH2 domains), and a C-terminal SH2 domain (2–6). Despite the lack of enzymatic activity, the presence of multiple protein-protein interaction domains allows Grb7 family adaptor proteins to participate in versatile signal transduction pathways and, therefore, to regulate many cellular functions (4–6). A number of signaling molecules has been reported to interact with these featured domains, although most of the identified Grb7 binding partners are mediated through its SH2 domain. For example, the SH2 domain of Grb7 has been demonstrated to be capable of binding to the phospho-tyrosine sites of EGF receptor (1), ErbB2 (7), ErbB3 and ErbB4 (8), Ret (9), platelet-derived growth factor receptor (10), insulin receptor (11), SHPTP2 (12), Tek/Tie2 (13), caveolin (14), c-Kit (15), EphB1 (16), G6f immunoreceptor protein (17), Rnd1 (18), Shc (7), FAK (19), and so on. The proceeding α-helix of the PH domain of Grb7 is the calmodulin-binding domain responsible for recruiting Grb7 to plasma membrane in a Ca²⁺-dependent manner (20), and the association between the PH domain of Grb7 and phosphoinositides is required for the phosphorylation by FAK (21). Two additional proteins, NIK (nuclear factor κB-inducing kinase) and FHL2 (four and half lim domains isoform 2), in association with the GM region (Grb and Mig homology region) of Grb7 are also reported, although the physiological functions for these interactions remain unknown (22, 23). Recently, other novel roles in translational controls and stress responses through the N terminus of Grb7 are implicated for the findings of Grb7 interacting with the 5′-untranslated region of capped targeted KOR (kappa opioid receptor) mRNA and the Hu antigen R of stress granules in an FAK-mediated phosphorylation manner (24, 25).Unlike its member proteins Grb10 and Grb14, the role of Grb7 in cell migration is unambiguous and well documented. This is supported by a series of studies. Firstly, Grb7 family members share a significantly conserved molecular architecture with the Caenorhabditis elegans Mig-10 protein, which is involved in neuronal cell migration during embryonic development (4, 5, 26), suggesting that Grb7 may play a role in cell migration. Moreover, Grb7 is often co-amplified with Her2/ErbB2 in certain human cancers and tumor cell lines (7, 27, 28), and its overexpression resulted in invasive and metastatic consequences of various cancers and tumor cells (23, 29–33). On the contrary, knocking down Grb7 by RNA interference conferred to an inhibitory outcome of the breast cancer motility (34). Furthermore, interaction of Grb7 with autophosphorylated FAK at Tyr-397 could promote integrin-mediated cell migration in NIH 3T3 and CHO cells, whereas overexpression of its SH2 domain, an dominant negative mutant of Grb7, inhibited cell migration (19, 35). Recruitment and phosphorylation of Grb7 by EphB1 receptors enhanced cell migration in an ephrin-dependent manner (16). Recently, G7–18NATE, a selective Grb7-SH2 domain affinity cyclic peptide, was demonstrated to efficiently block cell migration of tumor cells (32, 36). In addition to cell migration, Grb7 has been shown to play a role in a variety of physiological and pathological events, for instance, kidney development (37), tumorigenesis (7, 14, 38–41), angiogenic activity (20), proliferation (34, 42, 43), anti-apoptosis (44), gene expression regulation (24), Silver-Russell syndrome (45), rheumatoid arthritis (46), atopic dermatitis (47), and T-cell activation (17, 48). Nevertheless, it remains largely unknown regarding the downstream signaling events of Grb7-mediated various functions. In particular, given the role of Grb7 as an adaptor molecule and its SH2 domain mainly interacting with upstream regulators, it will be interesting to identify potential downstream effectors through interacting with the functional GM region or N-terminal proline-rich region.In this report, we identified two tyrosine phosphorylated sites of Grb7 by FAK and deciphered the signaling targets downstream through these phosphorylated tyrosine sites to regulate various cellular functions such as cell migration, proliferation, and survival. In addition, our study sheds light on tyrosine phosphorylation of Grb7 by FAK involved in tumorigenesis.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]     

Increased Protein Stability of FGF1 Can Compensate for Its Reduced Affinity for Heparin

Human FGF1 (fibroblast growth factor 1) is a powerful signaling molecule with a short half-life in vivo and a denaturation temperature close to physiological. Binding to heparin increases the stability of FGF1 and is believed to be important in the formation of FGF1·fibroblast growth factor receptor (FGFR) active complex. In order to reveal the function of heparin in FGF1·FGFR complex formation and signaling, we constructed several FGF1 variants with reduced affinity for heparin and with diverse stability. We determined their biophysical properties and biological activities as well as their ability to translocate across cellular membranes. Our study showed that increased thermodynamic stability of FGF1 nicely compensates for decreased binding of heparin in FGFR activation, induction of DNA synthesis, and cell proliferation. By stepwise introduction of stabilizing mutations into the K118E (K132E) FGF1 variant that shows reduced affinity for heparin and is inactive in stimulation of DNA synthesis, we were able to restore the full mitogenic activity of this mutant. Our results indicate that the main role of heparin in FGF-induced signaling is to protect this naturally unstable protein against heat and/or proteolytic degradation and that heparin is not essential for a direct FGF1-FGFR interaction and receptor activation.FGF1 (fibroblast growth factor 1) belongs to a family of polypeptide growth factors comprising in humans 22 structurally related proteins (1, 2). The signaling induced by the growth factor leads to a wide range of cellular responses during development as well as in adult life, such as growth regulation, differentiation, survival, stress response, migration, and proliferation of different cell types (3). The biological activity of FGF1 is exerted through binding to four high affinity cell surface receptors (FGFR1–4), resulting in receptor dimerization and transphosphorylation in its tyrosine kinase domain (4, 5). The activated FGFR³ induces cellular response by initiating several signaling cascades, including mitogen-activated protein kinase (MAPK), phosphoinositide 3-kinase/Akt, and phospholipase C-γ (PLC-γ) pathways (6).In addition to FGFRs, FGF1 binds to heparan sulfates (HS) associated with proteoglycans at the cell surface and in the extracellular matrix (7). Among the physiological sugars, the highest affinity for FGF1 is shown by heparin, a widely used linear, highly sulfated polysaccharide composed of 2-O-sulfated iduronic acid and 6-O-sulfated, N-sulfated glucosamine units (8).Despite many years of research, there is still controversy regarding the molecular role of heparin/HS in FGF1- and FGF2-induced signaling. Thus, the question of whether or not the linkage of two molecules of the growth factor by heparin/HS is an absolute prerequisite for induction of FGFR dimerization is still open. Numerous studies have concluded that the presence of heparin/HS is obligatory for FGF signaling. It is widely believed that heparin/HS is directly involved in receptor dimerization and is critical for mitogenic response stimulated by the growth factor (4, 6, 8–10).On the other hand, several authors working on FGF1 and FGF2 have suggested that there is no mandatory requirement for heparin for the assembly and activation of the FGF·FGFR complex. They imply that heparin only plays a role in association of two molecules of the growth factor and therefore facilitates their binding to FGFR (11). It has been reported that FGF1 and FGF2 can interact with the FGFR and trigger phosphorylation of p42/44 MAPK and activation of other signaling pathways even in the absence of HS (12–16).The accepted role of heparin/HS in FGF1 signaling is to prevent the degradation of the growth factor (17). The interaction with heparin or HS protects FGF1 against heat, acidic pH, and proteases (18, 19). HS also seems to regulate the activity of different FGFs by creating their local reservoir and generating a concentration gradient of the growth factor (6, 17).The binding of FGF1 to heparin/HS is mediated by specific residues forming a positively charged patch on the protein surface (20, 21). The major contribution is made by Lys¹¹⁸ (Lys¹³² in the full-length numbering system), which was identified by Harper and Lobb (22), and Lys¹¹² and Arg¹²² (23, 24). Additional residues of FGF1 involved in the interaction with heparin are the positively charged Lys¹¹³, Arg¹¹⁹, and Lys¹²⁸ and the polar Asn¹⁸, Asn¹¹⁴, and Gln¹²⁷ (20, 21). Site-directed mutagenesis and other studies have revealed the importance of Lys¹¹⁸ not only in heparin binding but also for the biological function of FGF1 (22, 25, 26). It was shown that the K118E (K132E) mutant is inactive in stimulation of DNA synthesis, although its affinity for FGFR and the ability to activate signaling cascades is not reduced (27, 28). Despite extensive research, the reason for the lack of mitogenic potential of K118E FGF1 is still not clear.In this paper, we verified the function of heparin in FGF1·FGFR complex formation and signaling by constructing several FGF1 mutants with reduced affinity for heparin. To recover the stability of these variants, which could no longer be stabilized by heparin, we supplemented them stepwise with stabilizing mutations (29). We analyzed thoroughly their biological activity and their ability to translocate across cellular membranes (30–34). Interestingly, the full mitogenic activity of the K118E FGF1 variant was restored by the introduced stabilizing mutations.Our results indicate that the main role of heparin in FGF-induced signaling is to protect this naturally unstable protein against heat denaturation and proteolytic degradation and that the increased stability of the growth factor can compensate for reduced heparin binding.