ASAP1, a Phospholipid-Dependent Arf GTPase-Activating Protein That Associates with and Is Phosphorylated by Src

Membrane trafficking is regulated in part by small GTP-binding proteins of the ADP-ribosylation factor (Arf) family. Arf function depends on the controlled exchange and hydrolysis of GTP. We have purified and cloned two variants of a 130-kDa phosphatidylinositol 4,5-biphosphate (PIP₂)-dependent Arf1 GTPase-activating protein (GAP), which we call ASAP1a and ASAP1b. Both contain a pleckstrin homology (PH) domain, a zinc finger similar to that found in another Arf GAP, three ankyrin (ANK) repeats, a proline-rich region with alternative splicing and SH3 binding motifs, eight repeats of the sequence E/DLPPKP, and an SH3 domain. Together, the PH, zinc finger, and ANK repeat regions possess PIP₂-dependent GAP activity on Arf1 and Arf5, less activity on Arf6, and no detectable activity on Arl2 in vitro. The cDNA for ASAP1 was independently identified in a screen for proteins that interact with the SH3 domain of the tyrosine kinase Src. ASAP1 associates in vitro with the SH3 domains of Src family members and with the Crk adapter protein. ASAP1 coprecipitates with Src from cell lysates and is phosphorylated on tyrosine residues in cells expressing activated Src. Both coimmunoprecipitation and tyrosine phosphorylation depend on the same proline-rich class II Src SH3 binding site required for in vitro association. By directly interacting with both Arfs and tyrosine kinases involved in regulating cell growth and cytoskeletal organization, ASAP1 could coordinate membrane remodeling events with these processes.Membrane traffic, the transfer of material between membrane-bound compartments, is needed for such diverse cellular processes as secretion, endocytosis, and changes in cell shape that accompany cell growth, division, and migration (reviewed in references 84, 85, and 87). It is mediated by transport vesicles that are formed by budding from a donor membrane. The process of budding is driven by the assembly of a proteinaceous coat. Once the vesicle is formed, the coat must dissociate to permit fusion with an acceptor membrane and the consequent delivery of the vesicle’s contents. These steps are regulated in part by the Arf family of small GTP-binding proteins (reviewed in references 8, 23, 61, and 63). Arfs are highly conserved and are found in eukaryotes ranging from yeast to humans. The mammalian Arf family is divided into several classes based largely on sequence similarity: class I (Arfs 1 through 3), class II (Arfs 4 and 5), class III (Arf6), and the more distantly related Arf-like (Arl) class. By linking GTP binding and hydrolysis to coat assembly and disassembly, Arfs regulate membrane trafficking at a number of sites. Arf1 has been implicated in endoplasmic reticulum-to-Golgi and intra-Golgi transport, endosome-to-endosome fusion, and synaptic vesicle formation (8, 23, 28, 61, 63, 66). Arf6 has been implicated in regulation of membrane traffic between the plasma membrane and a specialized endocytic compartment, and its function has been linked to cytoskeletal reorganization (25, 26, 71, 73, 74). The specific sites of action of the other Arf family members are not known.The hydrolysis of GTP on Arf requires a GTPase-activating protein (GAP) (19, 61). With multiple Arfs and multiple sites of action, the existence of several unique Arf GAPs had been anticipated. A number of activities have been purified or partially purified from mammalian sources, including rat liver (19, 57, 77), rat spleen (21), and bovine brain (79), and two Arf GAP activities from rat liver have been resolved (77). They have similar Arf specificities but differ in their lipid dependencies. One of the Arf GAPs (ArfGAP/ArfGAP1, hereafter referred to as ArfGAP1) which functions in the Golgi is activated by dioleoglycerols (3, 4, 19, 40). ArfGAP1, in common with a yeast Arf GAP, GCS1 (72), contains a zinc finger domain which is required for activity (19). The second Arf GAP (ArfGAP2) is specifically activated by phosphatidylinositol 4,5-bisphosphate (PIP₂) and phosphatidic acid (PA). Based on lipid requirements, ArfGAP2 was speculated to function at the plasma membrane and be regulated independently of ArfGAP1 (77). ArfGAP1 and ArfGAP2 were antigenically distinct and, therefore, likely to be distinct gene products; however, prior to this study, only ArfGAP1 had been cloned (19).Src, a cytoplasmic tyrosine kinase with N-terminal Src homology 3 (SH3) and SH2 domains, transduces signals important for cell growth and cytoskeletal organization (12, 68, 91). A number of studies suggest that Src is also involved in regulating membrane traffic. Src associates primarily with endosomal membranes and in several cell types has been localized to specialized secretory vesicles, including synaptic vesicles (5, 20, 34, 46, 54, 69, 81). Overexpression of Src accelerates endocytosis (95). In addition, Src associates with or phosphorylates several proteins involved in membrane trafficking (5, 31, 43, 65).Here, we report the purification and cloning of a PIP₂-dependent Arf GAP, ASAP1. ASAP1 contains a zinc finger domain similar to that required for GAP activity in ArfGAP1 and GCS1. ASAP1 also contains a number of domains that are likely to be involved in regulation and/or localization: a pleckstrin homology (PH) domain, three ankyrin (ANK) repeats, a proline-rich region with SH3 binding motifs, and an SH3 domain. In addition, ASAP1 was identified independently as a binding protein for Src and was found to be phosphorylated on tyrosine in cells that express activated Src. ASAP1 also associated with the adapter protein c-Crk in vitro. ASAP1 was localized to the cytoplasm and the cell edge likely associated with the plasma membrane. We propose that ASAP1, by binding both Src and PIP₂, could coordinate membrane trafficking with cell growth or actin cytoskeleton remodeling.     

Identification and Characterization of RBEL1 Subfamily of GTPases in the Ras Superfamily Involved in Cell Growth Regulation

Recently, we reported the identification of a novel gene named RBEL1 (Rab-like protein 1) and characterized its two encoded isoforms, RBEL1A and RBEL1B, that function as novel GTPases of Ras superfamily. Here we report the identification of two additional splice variants of RBEL1 that we have named RBEL1C and -D. All four RBEL1 isoforms (A, B, C, and D) have identical N termini harboring the Rab-like GTPase domains but contain variable C termini. Although all isoforms can be detected in both cytoplasm and nucleus, RBEL1A is predominantly cytoplasmic, whereas RBEL1B is mostly nuclear. RBEL1C and -D, by contrast, are evenly distributed between the cytoplasm and nucleus. Furthermore, all four RBEL1 proteins are also capable of associating with cellular membrane. The RBEL1 proteins also exhibit a unique nucleotide-binding potential and, whereas the larger A and B isoforms are mainly GTP-bound, the smaller C and D variants bind to both GTP and GDP. Furthermore, a regulatory region at amino acid position 236–302 immediately adjacent to the GTP-binding domain is important for GTP-binding potential of RBEL1A, because deletion of this region converts RBEL1A from predominantly GTP-bound to GDP-bound. RBEL1 knockdown via RNA interference results in marked cell growth suppression, which is associated with morphological and biochemical features of apoptosis as well as inhibition of extracellular signal-regulated kinase phosphorylation. Taken together, our results indicate that RBEL1 proteins are linked to cell growth and survival and possess unique biochemical, cellular, and functional characteristics and, therefore, appear to form a novel subfamily of GTPases within the Ras superfamily.The Ras superfamily is known to comprise five structurally distinct subfamilies of small GTPases, including Ras, Rho, Rab, Sar1/Arf, and Ran, and each subfamily of these GTPases possess distinct functions in the regulation of a variety of cellular processes such as cell proliferation, cell differentiation, cytoskeletal organization, protein transport, and trafficking (1–4). The Ras subfamily of GTPases (N-, H-, and K-Ras) function predominantly in relaying signals from receptors at the plasma membrane and modulating cell signaling pathways that regulate cell proliferation, differentiation, and survival (5). Ran GTPase, on other hand, is a key regulator of nucleocytoplasmic transport that regulates protein transport across the nuclear pore complex (6, 7). The Rab subfamily is the largest subfamily among the Ras superfamily and contains more than 60 members. The key functions of the Rab GTPases are to regulate protein exocytic and endocytic pathways and modulate intracellular protein transport/trafficking (8–13).In general, the Ras superfamily GTPases cycle between an active GTP-bound state and an inactive GDP-bound state. There are five N-terminal motifs involved in the binding and hydrolysis of GTP that are highly conserved among all GTPases: G1 (GXXXXGK(S/T)), G2 (T), G3 (DXXG), G4 ((N/T)(K/Q)XD), and G5 (EXSAX). Each sequence has particular functions involved in binding nucleotides (GTP or GDP) and facilitating hydrolysis (4, 14, 15). In general, the intrinsic GTPase activity (converting GTP to GDP) and exchange of GDP for GTP are slow processes for these GTPases and thus require regulatory proteins such as GTPase-activating proteins and GDP/GTP exchange factors to facilitate these processes (16–18).For the last two decades, the Ras superfamily has been a major focus in the cancer field as many of the members are either mutated or dysregulated in cancer. The founding members of the Ras superfamily, H-Ras and K-Ras, were first identified as viral oncogenes (1, 4). Later studies demonstrated that mutations of the Ras proteins (H-, N-, and K-Ras) occur frequently in human cancers, and the mutations identified are mostly clustered within the GTP-binding domains of the proteins thus locking Ras proteins in a GTP-bound configuration. GTP-bound Ras is constitutively active; it constantly activates its effector proteins to transduce cell proliferative signals (1, 4). Unlike Ras subfamily genes, mutations occurring in Rab and Rab-like genes are less common, yet alterations in gene expression of a number of Rab genes have been reported in multiple human malignancies. For example, Rab25 overexpression has been linked to prostate cancer progression (19). Rab2 overexpression has been found in lung adenomas and adenocarcinomas (20). In addition, alterations in Rab gene expression have also been linked to cancer drug resistance. For instance, resistance to the anticancer drug doxorubicin in MCF-7 cells has been linked with reduced expression of Rab6C, and introduction of exogenous Rab6C restores drug sensitivity (21).We have recently reported the identification two novel Ras superfamily GTPases, RBEL1A and RBEL1B (22). RBEL1A and RBEL1B are two splice variants of the RBEL1 gene and are highly homologous to the Rab and Ran GTPases within their N-terminal GTP-binding domains (22). Our studies show that both RBEL1A and -B predominantly bind to GTP. A single point mutation (T57N) in the GTP-binding domain of RBEL1A and -B abolishes their ability to bind to both GTP and GDP. Both RBEL1A and RBEL1B localize in the nucleus as well as in the cytosol. Whereas RBEL1A is predominantly cytosolic, RBEL1B is primarily nuclear. Interestingly, our studies also suggested that nucleotide (GTP or GDP)-binding could be important for the nuclear distribution of RBEL1B, because the nucleotide binding-deficient mutant form (T57N) of RBEL1B did not reside in the nucleus but rather became largely cytosolic (22).In our continuous efforts to fully elucidate the function of RBEL1, we have identified two additional splice variants that we have named RBEL1C and RBEL1D. Here we report further characterization of all four RBEL1 splice variants in terms of their GTPase activities, subcellular localizations, regulations, and potential functions. Our results indicate that RBEL1 GTPases, although sharing some common features with other Ras superfamily members, also harbor unique characteristics that are significantly different from other Ras superfamily GTPases. Based on our findings, we suggest that RBEL1 proteins appear to form a novel subfamily of GTPases within the Ras superfamily.     

Localized Diacylglycerol-dependent Stimulation of Ras and Rap1 during Phagocytosis

We describe a role for diacylglycerol in the activation of Ras and Rap1 at the phagosomal membrane. During phagocytosis, Ras density was similar on the surface and invaginating areas of the membrane, but activation was detectable only in the latter and in sealed phagosomes. Ras activation was associated with the recruitment of RasGRP3, a diacylglycerol-dependent Ras/Rap1 exchange factor. Recruitment to phagosomes of RasGRP3, which contains a C1 domain, parallels and appears to be due to the formation of diacylglycerol. Accordingly, Ras and Rap1 activation was precluded by antagonists of phospholipase C and of diacylglycerol binding. Ras is dispensable for phagocytosis but controls activation of extracellular signal-regulated kinase, which is partially impeded by diacylglycerol inhibitors. By contrast, cross-activation of complement receptors by stimulation of Fcγ receptors requires Rap1 and involves diacylglycerol. We suggest a role for diacylglycerol-dependent exchange factors in the activation of Ras and Rap1, which govern distinct processes induced by Fcγ receptor-mediated phagocytosis to enhance the innate immune response.Receptors that interact with the constant region of IgG (FcγR)⁴ mediate the recognition and elimination of soluble immune complexes and particles coated (opsonized) with immunoglobulins. Clustering of FcγR on the surface of leukocytes upon attachment to multivalent ligands induces their activation and subsequent internalization. Soluble immune complexes are internalized by endocytosis, a clathrin- and ubiquitylation-dependent process (1). In contrast, large, particulate complexes like IgG-coated pathogens are ingested by phagocytosis, a process that is contingent on extensive actin polymerization that drives the extension of pseudopods (2). In parallel with the internalization of the opsonized targets, cross-linking of phagocytic receptors triggers a variety of other responses that are essential components of the innate immune response. These include degranulation, activation of the respiratory burst, and the synthesis and release of multiple inflammatory agents (3, 4).Like T and B cell receptors, FcγR possesses an immunoreceptor tyrosine-based activation motif that is critical for signal transduction (3, 4). Upon receptor clustering, tyrosyl residues of the immunoreceptor tyrosine-based activation motif are phosphorylated by Src family kinases, thereby generating a docking site for Syk, a tyrosine kinase of the ZAP70 family (3, 4). The recruitment and activation of Syk in turn initiates a cascade of events that include activation of Tec family kinases, Rho- and ARF-family GTPases, phosphatidylinositol 3-kinase, phospholipase Cγ (PLCγ), and a multitude of additional effectors that together remodel the underlying cytoskeleton, culminating in internalization of the bound particle (5, 6).Phosphoinositide metabolism is thought to be critical for FcγR-induced phagocytosis (7, 8). Highly localized and very dynamic phosphoinositide changes have been observed at sites of phagocytosis: phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P₂) undergoes a transient accumulation at the phagocytic cup, which is rapidly superseded by its complete elimination from the nascent phagosome (7). The secondary disappearance of PtdIns(4,5)P₂ is attributable in part to the localized generation of phosphatidylinositol 3,4,5-trisphosphate, which has been reported to accumulate at sites of phagocytosis (9). Activation of PLCγ is also believed to contribute to the acute disappearance of PtdIns(4,5)P₂ in nascent phagosomes. Indeed, the generation of diacylglycerol (DAG) and inositol 1,4,5-trisphosphate has been detected by chemical means during FcγR-evoked particle ingestion (10, 11). Moreover, imaging experiments revealed that DAG appears at the time and at the precise site where PtdIns(4,5)P₂ is consumed (7).Two lines of evidence suggest that the DAG generated upon engagement of phagocytic receptors modulates particle engulfment. First, antagonists of PLC severely impair phagocytosis by macrophages (7, 12). This inhibition is not mimicked by preventing the associated [Ca²⁺] transient, suggesting that DAG, and not inositol 1,4,5-trisphosphate, is the crucial product of the PLC (13). Second, the addition of exogenous DAG or phorbol esters, which mimic the actions of endogenous DAG, augment phagocytosis (14, 15).Selective recognition of DAG by cellular ligands is generally mediated by specific regions of its target proteins, called C1 domains (16). Proteins bearing C1 domains include, most notably, members of the classical and novel families of protein kinase C (PKC), making them suitable candidates to account for the DAG dependence of phagocytosis. Indeed, PKCα, a classical isoform, and PKCϵ and PKCδ, both novel isoforms, are recruited to phagosomes (12, 15, 17, 18). Although the role of the various PKC isoforms in particle engulfment has been equivocal over the years, Cheeseman et al. (12) convincingly demonstrated that PKCϵ contributes to particle uptake in a PLC- and DAG-dependent manner.PKCs are not the sole proteins bearing DAG-binding C1 domains. Similar domains are also found in several other proteins, including members of the RasGRP family, chimaerins, and Munc-13 (19–21). One or more of these could contribute to the complex set of responses elicited by FcγR-induced DAG production. The RasGRP proteins are a class of exchange factors for the Ras/Rap family of GTPases (22). There are four RasGRP proteins (RasGRP1 to -4), and emerging evidence has implicated RasGRP1 and RasGRP3 in T and B cell receptor signaling (23–27).The possible role of DAG-mediated signaling pathways other than PKC in phagocytosis and the subsequent inflammatory response has not been explored. Here, we provide evidence that DAG stimulates Ras and Rap1 at sites of phagocytosis, probably through RasGRPs. Last, the functional consequences of Ras and Rap1 activation were analyzed.     

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]     

A PH Domain in the Arf GTPase-activating Protein (GAP) ARAP1 Binds Phosphatidylinositol 3,4,5-Trisphosphate and Regulates Arf GAP Activity Independently of Recruitment to the Plasma Membranes

ARAP1 is a phosphatidylinositol 3,4,5-trisphosphate (PtdIns(3,4,5)P₃)-dependent Arf GTPase-activating protein (GAP) with five PH domains that regulates endocytic trafficking of the epidermal growth factor receptor (EGFR). Two tandem PH domains are immediately N-terminal of the Arf GAP domain, and one of these fits the consensus sequence for PtdIns(3,4,5)P₃ binding. Here, we tested the hypothesis that PtdIns(3,4,5)P₃-dependent recruitment mediated by the first PH domain of ARAP1 regulates the in vivo and in vitro function of ARAP1. We found that PH1 of ARAP1 specifically bound to PtdIns(3,4,5)P₃, but with relatively low affinity (≈1.6 μm), and the PH domains did not mediate PtdIns(3,4,5)P₃-dependent recruitment to membranes in cells. However, PtdIns(3,4,5)P₃ binding to the PH domain stimulated GAP activity and was required for in vivo function of ARAP1 as a regulator of endocytic trafficking of the EGFR. Based on these results, we propose a variation on the model for the function of phosphoinositide-binding PH domains. In our model, ARAP1 is recruited to membranes independently of PtdIns(3,4,5)P₃, the subsequent production of which triggers enzymatic activity.Pleckstrin homology (PH)² domains are a common structural motif encoded by the human genome (1, 2). Approximately 10% of PH domains bind to phosphoinositides. These PH domains are thought to mediate phosphoinositide-dependent recruitment to membranes (1–3). Most PH domains likely have functions other than or in addition to phosphoinositide binding. For example, PH domains have been found to bind to protein and DNA (4–12). In addition, some PH domains have been found to be structurally and functionally integrated with adjacent domains (13, 14). A small fraction of PH domain-containing proteins (about 9% of the human proteins) have multiple PH domains arranged in tandem, which have been proposed to function as adaptors but have only been examined in one protein (15, 16). Arf GTPase-activating proteins (GAPs) of the ARAP family are phosphatidylinositol 3,4,5-trisphosphate (PtdIns(3,4,5)P₃)-dependent Arf GAPs with tandem PH domains (17, 18). The function of specific PH domains in regulating Arf GAP activity and for biologic activity has not been described.Arf GAPs are proteins that induce the hydrolysis of GTP bound to Arfs (19–23). The Arf proteins are members of the Ras superfamily of GTP-binding proteins (24–27). The six Arf proteins in mammals (five in humans) are divided into three classes based on primary sequence: Arf1, -2, and -3 are class 1, Arf4 and -5 are class 2, and Arf6 is class 3 (23, 24, 27–29). Class 1 and class 3 Arf proteins have been studied more extensively than class 2. They have been found to regulate membrane traffic and the actin cytoskeleton.The Arf GAPs are a family of proteins with diverse domain structures (20, 21, 23, 30). ARAPs, the most structurally complex of the Arf GAPs, contain, in addition to an Arf GAP domain, the sterile α motif (SAM), five PH, Rho GAP, and Ras association domains (17, 18, 31, 32). The first and second and the third and fourth PH domains are tandem (Fig. 1). The first and third PH domains of the ARAPs fit the consensus for PtdIns(3,4,5)P₃ binding (33–35). ARAPs have been found to affect actin and membrane traffic (21, 23). ARAP3 regulates growth factor-induced ruffling of porcine aortic endothelial cells (31, 36, 37). The function is dependent on the Arf GAP and Rho GAP domains. ARAP2 regulates focal adhesions, an actin cytoskeletal structure (17). ARAP2 function requires Arf GAP activity and a Rho GAP domain capable of binding RhoA·GTP. ARAP1 has been found to have a role in membrane traffic (18). The protein associates with pre-early endosomes involved in the attenuation of EGFR signals. The function of the tandem PH domains in the ARAPs has not been examined.Open in a separate windowFIGURE 1.ARAP1 binding to phospholipids. A, schematic of the recombinant proteins used in this study. Domain abbreviations: Ank, ankyrin repeat; PLCδ-PH, PH domain of phospholipase C δ; RA, Ras association motif; RhoGAP, Rho GTPase-activating domain. B, ARAP1 phosphoinositide binding specificity. 500 nm PH1-Ank recombinant protein was incubated with sucrose-loaded LUVs formed by extrusion through a 1-μm pore filter. LUVs contained PtdIns alone or PtdIns with 2.5 μm PtdIns(3,4,5)P₃, 2.5 μm PtdIns(3)P, 2.5 μm PtdIns(4)P, 2.5 μm PtdIns(5)P, 2.5 μm PtdIns(3,4)P₂, 2.5 μm PtdIns(3,5)P₂, or 2.5 μm PtdIns(4,5)P₂ with a total phosphoinositide concentration of 50 μm and a total phospholipid concentration of 500 μm. Vesicles were precipitated by ultracentrifugation, and associated proteins were separated by SDS-PAGE. The amount of precipitated protein was determined by densitometry of the Coomassie Blue-stained gels with standards on each gel. C, PtdIns(3,4,5)P₃-dependent binding of ARAP1 to LUVs. 1 μm PH1-Ank or ArfGAP-Ank recombinant protein was incubated with 1 mm sucrose-loaded LUVs formed by extrusion through a 1-μm pore size filter containing varying concentration of PtdIns(3,4,5)P₃. Precipitation of LUVs and analysis of associated proteins were performed as described in B. The average ± S.E. of three independent experiments is presented.Here we investigated the role of the first two PH domains of ARAP1 for catalysis and in vivo function. The first PH domain fits the consensus sequence for PtdIns(3,4,5)P₃ binding (33–35). The second does not fit a phosphoinositide binding consensus but is immediately N-terminal to the GAP domain. We have previously reported that the PH domain that occurs immediately N-terminal of the Arf GAP domain of ASAP1 is critical for the catalytic function of the protein (38, 39). We tested the hypothesis that the two PH domains of ARAP1 function independently; one recruits ARAP1 to PtdIns(3,4,5)P₃-rich membranes, and the other functions with the catalytic domain. As predicted, PH1 interacted specifically with PtdIns(3,4,5)P₃, and PH2 did not. However, both PH domains contributed to catalysis independently of recruitment to membranes. None of the PH domains in ARAP1 mediated PtdIns(3,4,5)P₃-dependent targeting to plasma membranes (PM). PtdIns(3,4,5)P₃ stimulated GAP activity, and the ability to bind PtdIns(3,4,5)P₃ was required for ARAP1 to regulate membrane traffic. We propose that ARAP1 is recruited independently of PtdIns(3,4,5)P₃ to the PM where PtdIns(3,4,5)P₃ subsequently regulates its GAP activity to control endocytic events.     

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]     

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]     

Arachidonic Acid Stimulates Cell Adhesion through a Novel p38 MAPK-RhoA Signaling Pathway That Involves Heat Shock Protein 27

Rho GTPases are critical components of cellular signal transduction pathways. Both hyperactivity and overexpression of these proteins have been observed in human cancers and have been implicated as important factors in metastasis. We previously showed that dietary n-6 fatty acids increase cancer cell adhesion to extracellular matrix proteins, such as type IV collagen. Here we report that in MDA-MB-435 human melanoma cells, arachidonic acid activates RhoA, and inhibition of RhoA signaling with either C3 exoenzyme or dominant negative Rho blocked arachidonic acid-induced cell adhesion. Inhibition of the Rho kinase (ROCK) with either small molecule inhibitors or ROCK II-specific small interfering RNA (siRNA) blocked the fatty acid-induced adhesion. However, unlike other systems, inhibition of ROCK did not block the activation of p38 mitogen-activated protein kinase (MAPK); instead, Rho activation depended on p38 MAPK activity and the presence of heat shock protein 27 (HSP27), which is phosphorylated downstream of p38 after arachidonic acid treatment. HSP27 associated with p115RhoGEF in fatty acid-treated cells, and this association was blocked when p38 was inhibited. Furthermore, siRNA knockdown of HSP27 blocked the fatty acid-stimulated Rho activity. Expression of dominant negative p115-RhoGEF or p115RhoGEF-specific siRNA inhibited both RhoA activation and adhesion on type IV collagen, whereas a constitutively active p115RhoGEF restored the arachidonic acid stimulation in cells in which the p38 MAPK had been inhibited. These data suggest that n-6 dietary fatty acids stimulate a set of interactions that regulates cell adhesion through RhoA and ROCK II via a p38 MAPK-dependent association of HSP27 and p115RhoGEF.The ability of tumor cells to metastasize to secondary sites is a hallmark of neoplastic disease. Unfortunately, this propensity to spread is the primary cause of morbidity and death in cancer patients (1). Metastasis is clearly a highly regulated, multistep process that occurs in a spatiotemporal manner (2–4). To escape the restrictive compartment boundaries characteristic of adult tissue, separate intravasation and extravasation steps requiring alterations in co-adhesion, adhesion, invasion, and migration must occur. Execution of these biological processes, involving multiple proteins and cellular organelles, require highly coordinated cell signaling mechanisms.The Rho family of small GTPases regulates many facets of cytoskeletal rearrangements that facilitate cell attachment and migration (5–7). Rho GTPases act as molecular switches by changing from an inactive GDP-bound conformation to an active GTP-bound conformation, thereby regulating a signaling pathway. These proteins are directly regulated by Rho guanine nucleotide exchange factors (GEFs),² Rho GTPase activating proteins, and Rho GDP-dissociation inhibitors (8–12). RhoGEFs bind to the GTPase to catalyze the dissociation of GDP, allowing the binding of GTP and thereby promoting Rho activation (8). The RGS (regulators of G protein signaling) domain-containing RhoGEFs are a recently described family of GEFs. Currently, there are three members of this family, PDZ-RhoGEF, LARG, and p115RhoGEF (13–15), in which the RGS domains function as a heterotrimeric GTPase-activating domain (13, 15, 16). The RGS family of RhoGEFs has been shown to regulate Rho during several processes including cytoskeletal rearrangements, cell adhesion, and cancer progression (17–21).There is significant interplay between the activity of small GTPases and signaling derived from fatty acid metabolism (22–28). Linoleic acid, which is metabolized to arachidonic acid, is an n-6 polyunsaturated fatty acid that is present at high levels in most western diets (29). In animal models, diets high in n-6 polyunsaturated fatty acids have been shown to enhance tumor progression and metastasis (30, 31). Additionally, arachidonic acid is stored in cell membranes and is made available by phospholipases under conditions of increased inflammatory response (32). Arachidonic acid is further metabolized by cyclooxygenases (COX), lipoxygenases (LOX), and cytochrome P450 monooxygenases to yield bioactive products that have myriad effects on cells, and altered metabolism of arachidonic acid by COX, LOX, and P450 has been implicated in cancer progression (31, 33–36).We have studied mechanisms of cell adhesion using the MDA-MB-435 cells as a model of a highly metastatic human cancer cell line (37). These cells have been extensively studied for their ability to recapitulate the metastatic cascade in vivo and in vitro, although recent work indicates that the cells currently in use are most likely a human melanoma line (38). We initially observed that arachidonic acid (AA) enhanced adhesion of MDA-MB-435 cells to type IV collagen through specific integrin-mediated pathways (37). Exogenous AA led to the activation of mitogen-activated protein kinase (MAPK)-activated protein kinase 2 and the phosphorylation of heat shock protein 27 (HSP27) via a p38 MAPK-dependent process (39). Inhibition of p38 MAPK activation blocked cell adhesion as did function-blocking antibodies specific for subunits of the collagen receptor (40). More recently, we identified the key metabolite of AA (15-(S)- hydroxyeicosatetraenoic acid) and the upstream kinases (TAK1 and MKK6) that are responsible for activation of p38 MAPK in this system (41).In this study we investigated the role of Rho activation in the MDA-MB-435 cells after exposure to arachidonic acid. Several aspects of the regulation of Rho signaling in these cells provide insights into the cross-talk between important signaling pathways.     

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]     

Nucleophosmin and Nucleolin Regulate K-Ras Plasma Membrane Interactions and MAPK Signal Transduction

The spatial organization of Ras proteins into nanoclusters on the inner leaflet of the plasma membrane is essential for high fidelity signaling through the MAPK pathway. Here we identify two selective regulators of K-Ras nanoclustering from a proteomic screen for K-Ras interacting proteins. Nucleophosmin (NPM) and nucleolin are predominantly localized to the nucleolus but also have extranuclear functions. We show that a subset of NPM and nucleolin localizes to the inner leaflet of plasma membrane and forms specific complexes with K-Ras but not other Ras isoforms. Active GTP-loaded and inactive GDP-loaded K-Ras both interact with NPM, although NPM-K-Ras binding is increased by growth factor receptor activation. NPM and nucleolin both stabilize K-Ras levels on the plasma membrane, but NPM concurrently increases the clustered fraction of GTP-K-Ras. The increase in nanoclustered GTP-K-Ras in turn enhances signal gain in the MAPK pathway. In summary these results reveal novel extranucleolar functions for NPM and nucleolin as regulators of K-Ras nanocluster formation and activation of the MAPK pathway. The study also identifies a new class of K-Ras nanocluster regulator that operates independently of the structural scaffold galectin-3.Ras proteins are small GTPases that function as molecular switches on the inner leaflet of the plasma membrane, conveying extracellular signals to the cell interior. Ras proteins are critical regulators of signal transduction pathways controlling key cell fates such as cell growth, differentiation, and apoptosis. Deregulation of these pathways results in aberrant cell growth and tumor formation. Mutations that render Ras constitutively active are found in ∼15% of human cancers, making Ras one of the most clinically significant proteins in human carcinogenesis. Oncogenic mutations are most prevalent in the K-Ras gene, accounting for a large proportion of solid tumors including 90% of pancreatic cancer, 50% of colon cancer, and 30% of non-small cell lung cancer (1, 2).The three major Ras isoforms, H-, N-, and K-Ras generate distinct signal outputs in intact cells, signifying specific roles for each isoform. These functional differences stem from significant sequence divergence in the Ras C-terminal 25 amino acids of the hypervariable region (HVR)³ that directs post-translation attachment of different lipid anchors. The minimal membrane anchor of H-Ras comprises two palmitate groups and a farnesyl group, whereas K-Ras is tethered by a farnesyl group and a polybasic domain (3, 4). These minimal anchors, together with flanking protein sequences and the G-domain, interact with lipids and proteins of the plasma membrane, driving the Ras isoforms into spatially and structurally distinct nanodomains on the plasma membrane (5, 6). Ras lateral segregation is further modulated by the activation state of Ras; active GTP-loaded H-Ras is organized in cholesterol-independent nanoclusters, whereas inactive GDP-loaded H-Ras is arrayed in cholesterol-dependent nanoclusters (5, 7–9). Recent work has also shown that GTP-K-Ras clusters into nanodomains that are spatially distinct from GDP-K-Ras, although both types of nanocluster are cholesterol-independent and actin-dependent (7, 9). K-Ras-GTP nanoclustering, however, is regulated by galectin-3, which operates as a nanodomain scaffold (10, 11).Ras-GTP nanoclusters are the sites of Raf/MEK and ERK recruitment to the plasma membrane. Scaffolding all components of the MAPK module within nanoclusters rewires the biochemistry to generate a digital ERKpp output. The operation of Ras-GTP nanoclusters as highly sensitive digital switches is critical to deliver high fidelity signal transmission across the plasma membrane (12–14). A key parameter in epidermal growth factor (EGF) receptor to MAPK signal transmission is the fraction of Ras-GTP that forms nanoclusters; this clustered fraction sets the gain for cellular MAPK signaling (15, 16).NPM (also known as B23) and nucleolin are multifunctional phosphoproteins predominately localized to the nucleolus that play key roles in ribosome biogenesis (17–19). For example, NPM exhibits ribonuclease activity and preferentially cleaves pre-rRNA. NPM and nucleolin also have functions outside of the nucleolus. Both proteins shuttle between the nucleolus and the cytoplasm (20), and this shuttling may allow NPM to operate as molecular chaperone (21). In addition cytosolic NPM is involved in centrosome duplication (22). Like Ras proteins, NPM and nucleolin regulate cell proliferation and transformation and are overexpressed in multiple cancers (23). However, the physiological role of NPM in carcinogenesis remains controversial because it has been described as both an oncogene and a tumor suppressor (23).In this study we identify NPM and nucleolin as proteins that interact specifically with K-Ras but not H-Ras. Furthermore we definitively identify a subset of NPM and nucleolin on the inner leaflet of the plasma membrane where both proteins interact with K-Ras. Importantly, NPM and nucleolin stabilize K-Ras levels on the plasma membrane, leading to an increase in the K-Ras clustered fraction, which amplifies signal output from the MAPK pathway. Combined, our data indicate that NPM and nucleolin play a critical role in signal transduction via the MAPK pathway.