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]     

Role of Partitioning-defective 1/Microtubule Affinity-regulating Kinases in the Morphogenetic Activity of Helicobacter pylori CagA

Helicobacter pylori CagA plays a key role in gastric carcinogenesis. Upon delivery into gastric epithelial cells, CagA binds and deregulates SHP-2 phosphatase, a bona fide oncoprotein, thereby causing sustained ERK activation and impaired focal adhesions. CagA also binds and inhibits PAR1b/MARK2, one of the four members of the PAR1 family of kinases, to elicit epithelial polarity defect. In nonpolarized gastric epithelial cells, CagA induces the hummingbird phenotype, an extremely elongated cell shape characterized by a rear retraction defect. This morphological change is dependent on CagA-deregulated SHP-2 and is thus thought to reflect the oncogenic potential of CagA. In this study, we investigated the role of the PAR1 family of kinases in the hummingbird phenotype. We found that CagA binds not only PAR1b but also other PAR1 isoforms, with order of strength as follows: PAR1b > PAR1d ≥ PAR1a > PAR1c. Binding of CagA with PAR1 isoforms inhibits the kinase activity. This abolishes the ability of PAR1 to destabilize microtubules and thereby promotes disassembly of focal adhesions, which contributes to the hummingbird phenotype. Consistently, PAR1 knockdown potentiates induction of the hummingbird phenotype by CagA. The morphogenetic activity of CagA was also found to be augmented through inhibition of non-muscle myosin II. Because myosin II is functionally associated with PAR1, perturbation of PAR1-regulated myosin II by CagA may underlie the defect of rear retraction in the hummingbird phenotype. Our findings reveal that CagA systemically inhibits PAR1 family kinases and indicate that malfunctioning of microtubules and myosin II by CagA-mediated PAR1 inhibition cooperates with deregulated SHP-2 in the morphogenetic activity of CagA.Infection with Helicobacter pylori strains bearing cagA (cytotoxin-associated gene A)-positive strains is the strongest risk factor for the development of gastric carcinoma, the second leading cause of cancer-related death worldwide (1–3). The cagA gene is located within a 40-kb DNA fragment, termed the cag pathogenicity island, which is specifically present in the genome of cagA-positive H. pylori strains (4–6). In addition to cagA, there are ∼30 genes in the cag pathogenicity island, many of which encode a bacterial type IV secretion system that delivers the cagA-encoded CagA protein into gastric epithelial cells (7–10). Upon delivery into gastric epithelial cells, CagA is localized to the plasma membrane, where it undergoes tyrosine phosphorylation at the C-terminal Glu-Pro-Ile-Tyr-Ala motifs by Src family kinases or c-Abl kinase (11–14). The C-terminal Glu-Pro-Ile-Tyr-Ala-containing region of CagA is noted for the structural diversity among distinct H. pylori isolates. Oncogenic potential of CagA has recently been confirmed by a study showing that systemic expression of CagA in mice induces gastrointestinal and hematological malignancies (15).When expressed in gastric epithelial cells, CagA induces morphological transformation termed the hummingbird phenotype, which is characterized by the development of one or two long and thin protrusions resembling the beak of the hummingbird. It has been thought that the hummingbird phenotype is related to the oncogenic action of CagA (7, 16–19). Pathophysiological relevance for the hummingbird phenotype in gastric carcinogenesis has recently been provided by the observation that infection with H. pylori carrying CagA with greater ability to induce the hummingbird phenotype is more closely associated with gastric carcinoma (20–23). Elevated motility of hummingbird cells (cells showing the hummingbird phenotype) may also contribute to invasion and metastasis of gastric carcinoma.In host cells, CagA interacts with the SHP-2 phosphatase, C-terminal Src kinase, and Crk adaptor in a tyrosine phosphorylation-dependent manner (16, 24, 25) and also associates with Grb2 adaptor and c-Met in a phosphorylation-independent manner (26, 27). Among these CagA targets, much attention has been focused on SHP-2 because the phosphatase has been recognized as a bona fide oncoprotein, gain-of-function mutations of which are found in various human malignancies (17, 18, 28). Stable interaction of CagA with SHP-2 requires CagA dimerization, which is mediated by a 16-amino acid CagA-multimerization (CM)² sequence present in the C-terminal region of CagA (29). Upon complex formation, CagA aberrantly activates SHP-2 and thereby elicits sustained ERK MAP kinase activation that promotes mitogenesis (30). Also, CagA-activated SHP-2 dephosphorylates and inhibits focal adhesion kinase (FAK), causing impaired focal adhesions. It has been shown previously that both aberrant ERK activation and FAK inhibition by CagA-deregulated SHP-2 are involved in induction of the hummingbird phenotype (31).Partitioning-defective 1 (PAR1)/microtubule affinity-regulating kinase (MARK) is an evolutionally conserved serine/threonine kinase originally isolated in C. elegans (32–34). Mammalian cells possess four structurally related PAR1 isoforms, PAR1a/MARK3, PAR1b/MARK2, PAR1c/MARK1, and PAR1d/MARK4 (35–37). Among these, PAR1a, PAR1b, and PAR1c are expressed in a variety of cells, whereas PAR1d is predominantly expressed in neural cells (35, 37). These PAR1 isoforms phosphorylate microtubule-associated proteins (MAPs) and thereby destabilize microtubules (35, 38), allowing asymmetric distribution of molecules that are involved in the establishment and maintenance of cell polarity.In polarized epithelial cells, CagA disrupts the tight junctions and causes loss of apical-basolateral polarity (39, 40). This CagA activity involves the interaction of CagA with PAR1b/MARK2 (19, 41). CagA directly binds to the kinase domain of PAR1b in a tyrosine phosphorylation-independent manner and inhibits the kinase activity. Notably, CagA binds to PAR1b via the CM sequence (19). Because PAR1b is present as a dimer in cells (42), CagA may passively homodimerize upon complex formation with the PAR1 dimer via the CM sequence, and this PAR1-directed CagA dimer would form a stable complex with SHP-2 through its two SH2 domains.Because of the critical role of CagA in gastric carcinogenesis (7, 16–19), it is important to elucidate the molecular basis underlying the morphogenetic activity of CagA. In this study, we investigated the role of PAR1 isoforms in induction of the hummingbird phenotype by CagA, and we obtained evidence that CagA-mediated inhibition of PAR1 kinases contributes to the development of the morphological change by perturbing microtubules and non-muscle myosin II.     

Sip1, a Novel RS Domain-Containing Protein Essential for Pre-mRNA Splicing

Previous studies have shown that protein-protein interactions among splicing factors may play an important role in pre-mRNA splicing. We report here identification and functional characterization of a new splicing factor, Sip1 (SC35-interacting protein 1). Sip1 was initially identified by virtue of its interaction with SC35, a splicing factor of the SR family. Sip1 interacts with not only several SR proteins but also with U1-70K and U2AF65, proteins associated with 5′ and 3′ splice sites, respectively. The predicted Sip1 sequence contains an arginine-serine-rich (RS) domain but does not have any known RNA-binding motifs, indicating that it is not a member of the SR family. Sip1 also contains a region with weak sequence similarity to the Drosophila splicing regulator suppressor of white apricot (SWAP). An essential role for Sip1 in pre-mRNA splicing was suggested by the observation that anti-Sip1 antibodies depleted splicing activity from HeLa nuclear extract. Purified recombinant Sip1 protein, but not other RS domain-containing proteins such as SC35, ASF/SF2, and U2AF65, restored the splicing activity of the Sip1-immunodepleted extract. Addition of U2AF65 protein further enhanced the splicing reconstitution by the Sip1 protein. Deficiency in the formation of both A and B splicing complexes in the Sip1-depleted nuclear extract indicates an important role of Sip1 in spliceosome assembly. Together, these results demonstrate that Sip1 is a novel RS domain-containing protein required for pre-mRNA splicing and that the functional role of Sip1 in splicing is distinct from those of known RS domain-containing splicing factors.Pre-mRNA splicing takes place in spliceosomes, the large RNA-protein complexes containing pre-mRNA, U1, U2, U4/6, and U5 small nuclear ribonucleoprotein particles (snRNPs), and a large number of accessory protein factors (for reviews, see references 21, 22, 37, 44, and 48). It is increasingly clear that the protein factors are important for pre-mRNA splicing and that studies of these factors are essential for further understanding of molecular mechanisms of pre-mRNA splicing.Most mammalian splicing factors have been identified by biochemical fractionation and purification (3, 15, 19, 31–36, 45, 69–71, 73), by using antibodies recognizing splicing factors (8, 9, 16, 17, 61, 66, 67, 74), and by sequence homology (25, 52, 74).Splicing factors containing arginine-serine-rich (RS) domains have emerged as important players in pre-mRNA splicing. These include members of the SR family, both subunits of U2 auxiliary factor (U2AF), and the U1 snRNP protein U1-70K (for reviews, see references 18, 41, and 59). Drosophila alternative splicing regulators transformer (Tra), transformer 2 (Tra2), and suppressor of white apricot (SWAP) also contain RS domains (20, 40, 42). RS domains in these proteins play important roles in pre-mRNA splicing (7, 71, 75), in nuclear localization of these splicing proteins (23, 40), and in protein-RNA interactions (56, 60, 64). Previous studies by us and others have demonstrated that one mechanism whereby SR proteins function in splicing is to mediate specific protein-protein interactions among spliceosomal components and between general splicing factors and alternative splicing regulators (1, 1a, 6, 10, 27, 63, 74, 77). Such protein-protein interactions may play critical roles in splice site recognition and association (for reviews, see references 4, 18, 37, 41, 47 and 59). Specific interactions among the splicing factors also suggest that it is possible to identify new splicing factors by their interactions with known splicing factors.Here we report identification of a new splicing factor, Sip1, by its interaction with the essential splicing factor SC35. The predicted Sip1 protein sequence contains an RS domain and a region with sequence similarity to the Drosophila splicing regulator, SWAP. We have expressed and purified recombinant Sip1 protein and raised polyclonal antibodies against the recombinant Sip1 protein. The anti-Sip1 antibodies specifically recognize a protein migrating at a molecular mass of approximately 210 kDa in HeLa nuclear extract. The anti-Sip1 antibodies sufficiently deplete Sip1 protein from the nuclear extract, and the Sip1-depleted extract is inactive in pre-mRNA splicing. Addition of recombinant Sip1 protein can partially restore splicing activity to the Sip1-depleted nuclear extract, indicating an essential role of Sip1 in pre-mRNA splicing. Other RS domain-containing proteins, including SC35, ASF/SF2, and U2AF65, cannot substitute for Sip1 in reconstituting splicing activity of the Sip1-depleted nuclear extract. However, addition of U2AF65 further increases splicing activity of Sip1-reconstituted nuclear extract, suggesting that there may be a functional interaction between Sip1 and U2AF65 in nuclear extract.     

New Insights into How the Rho Guanine Nucleotide Dissociation Inhibitor Regulates the Interaction of Cdc42 with Membranes

The subcellular localization of the Rho family GTPases is of fundamental importance to their proper functioning in cells. The Rho guanine nucleotide dissociation inhibitor (RhoGDI) plays a key regulatory role by influencing the cellular localization of Rho GTPases and is essential for the transforming activity of oncogenic forms of Cdc42. However, the mechanism by which RhoGDI helps Cdc42 to undergo the transition between a membrane-associated protein and a soluble (cytosolic) species has been poorly understood. Here, we examine how RhoGDI influences the binding of Cdc42 to lipid bilayers. Despite having similar affinities for the signaling-inactive (GDP-bound) and signaling-active (GTP-bound) forms of Cdc42 in solution, we show that when RhoGDI interacts with Cdc42 along the membrane surface, it has a much higher affinity for GDP-bound Cdc42 compared with its GTP-bound counterpart. Interestingly, the rate for the dissociation of Cdc42·RhoGDI complexes from membranes is unaffected by the nucleotide-bound state of Cdc42. Moreover, the membrane release of Cdc42·RhoGDI complexes occurs at a similar rate as the release of Cdc42 alone, with the major effect of RhoGDI being to impede the re-association of Cdc42 with membranes. These findings lead us to propose a new model for how RhoGDI influences the ability of Cdc42 to move between membranes and the cytosol, which highlights the role of the membrane in helping RhoGDI to distinguish between the GDP- and GTP-bound forms of Cdc42 and holds important implications for how it functions as a key regulator of the cellular localization and signaling activities of this GTPase.The Rho family GTPases are a tightly regulated class of signaling proteins that controls a number of important cellular processes. Known most prominently for their ability to remodel the actin cytoskeleton in mammalian cells (1–3), members of this GTPase family have been shown to play essential roles in cell migration, epithelial cell polarization, phagocytosis, and cell cycle progression (4–11). The Rho family member Cdc42 was discovered for its essential role in bud formation in Saccharomyces cerevisiae (12). However, after its identification in higher organisms (13), Cdc42 has been implicated in a diverse array of signaling pathways including those involved in the regulation of cell growth and in the induction of malignant transformation (14). Indeed, point mutations which enable Cdc42 to undergo the spontaneous exchange of GDP for GTP cause NIH3T3 cells to form colonies in soft agar and grow in low serum, two hallmarks of cellular transformation (15). The introduction of activated Cdc42 mutants into nude mice gives rise to tumor formation (16). Moreover, cellular transformation by oncogenic Ras, one of the most commonly mutated proteins in human cancers, requires the activation of Cdc42 (17).At the molecular level, there are a number of mechanisms that possibly contribute to the roles played by Cdc42 in cell growth control and cellular transformation. These include the ability of Cdc42 to activate the c-Jun NH₂-terminal kinase and p38/Mpk2 signaling pathways (18–20) as well as spatially regulate proteins implicated in the establishment of microtubule-dependent cell polarity including glycogen synthase kinase-3β and adenomatous polyposis coli (21), extend the lifetime of epidermal growth factor receptor-signaling activities by sequestering Cbl, a ubiquitin E3 ligase (22), and influence intracellular trafficking events (23, 24). To mediate such a wide range of cellular responses, two parameters must be properly regulated; that is, the activation state of Cdc42 and its subcellular localization. As is the case with other GTPases, the activation of Cdc42 occurs as an outcome of GDP-GTP exchange, which then enables it to undergo high affinity interactions with effector proteins (25–27). Upon the hydrolysis of GTP to GDP, Cdc42 is converted back to a signaling-inactive state. Two families of proteins work in opposing fashion to regulate the GTP-binding/GTPase cycle of Cdc42. GTPase-activating proteins recognize the GTP-bound form of Cdc42 and accelerate the hydrolysis of GTP to GDP, rendering Cdc42 inactive (28, 29). Guanine nucleotide exchange factors (GEFs)² stimulate the dissociation of GDP from Cdc42, thereby promoting the formation of its signaling-active, GTP-bound state (29, 30).Of equal importance to its activation status is the spatial regulation of Cdc42. This is highly contingent on the particular cellular membranes that serve as sites of binding and/or recruitment of Cdc42 (31–33). The vast majority of in vitro studies performed on Cdc42 have been carried out in the absence of lipids, which is an important omission considering that virtually all of the physiological functions of Cdc42 occur on a membrane surface (34). Cdc42, along with most other Rho family GTPases, undergoes a series of carboxyl-terminal modifications which result in the covalent attachment of a 20-carbon geranylgeranyl lipid anchor (35–37). Directly preceding this lipid tail is a sequence of basic residues that further stabilizes the association of Cdc42 with the membrane surface (31, 33, 38). A ubiquitously expressed 22-kDa protein called Rho guanine nucleotide dissociation inhibitor (RhoGDI) was found to form a soluble (cytosolic) complex with Cdc42 and other Rho GTPases and to apparently promote their release from membranes (39, 40). RhoGDI was originally discovered and named for its ability to block the GEF- and EDTA-stimulated nucleotide exchange activity of Rho family GTPases (39, 41, 42) and then subsequently shown to inhibit the GTP-hydrolytic activity of Cdc42 (43) and to be capable of interacting with the GDP- and GTP-bound forms of Cdc42 in solution with equal affinity (44). The x-ray crystal structure of a complex between RhoGDI and Cdc42-GDP revealed two types of binding interactions (45). An amino-terminal regulatory arm of RhoGDI was shown to form a helix-loop-helix motif that binds to both of the switch domains of Cdc42, leading to the inhibition of GTP hydrolysis and GDP dissociation (45, 46). The carboxyl-terminal two-thirds of RhoGDI assumes an immunoglobulin-like domain, forming a hydrophobic pocket that in effect provides a membrane substitute for the geranylgeranyl moiety of Cdc42. After release from membranes, the lipid anchor of Cdc42 binds in the hydrophobic pocket of RhoGDI, thereby helping to maintain Cdc42 in solution (45–47).Prior work from our laboratory has demonstrated an essential role for RhoGDI in Cdc42-mediated cellular transformation. Based on the x-ray crystal structure for the Cdc42·RhoGDI complex, Arg-66 of Cdc42 makes multiple contacts with RhoGDI. When this residue was changed to alanine, Cdc42 was unable to bind to RhoGDI but was still capable of interacting with its other regulatory and effector proteins. Interestingly, when the R66A mutant of Cdc42 was examined in the constitutively active Cdc42(F28L) background, the resulting Cdc42 double mutant was no longer able to transform cells (48). Knocking down RhoGDI by small interfering RNA also blocked transformation by Cdc42. These findings highlighted a key role for RhoGDI in the ability of Cdc42 to stimulate signaling pathways of importance to cellular transformation, presumably by influencing the membrane association of Cdc42 and ensuring its proper cellular localization.In the present study we have set out to better understand how RhoGDI regulates the signaling functions of Cdc42 and, in particular, how RhoGDI affects the association of Cdc42 with membranes. We show how the membrane plays a previously unappreciated role in allowing RhoGDI to distinguish between the signaling-inactive (GDP-bound) and signaling-active (GTP-bound) forms of Cdc42. By assaying the binding of Cdc42 to insect cell membranes and compositionally defined liposomes through different approaches including a sensitive, real-time fluorescence resonance energy transfer (FRET) readout, we have been able to establish how RhoGDI influences the ability of Cdc42 to transition between a membrane-bound and soluble species. This has led us to propose a new mechanism describing how RhoGDI performs its important regulatory function.     

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]