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]     

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]     

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.     

Aberrant Splice Variants of HAS1 (Hyaluronan Synthase 1) Multimerize with and Modulate Normally Spliced HAS1 Protein: A POTENTIAL MECHANISM PROMOTING HUMAN CANCER*

Most human genes undergo alternative splicing, but aberrant splice forms are hallmarks of many cancers, usually resulting from mutations initiating abnormal exon skipping, intron retention, or the introduction of a new splice sites. We have identified a family of aberrant splice variants of HAS1 (the hyaluronan synthase 1 gene) in some B lineage cancers, characterized by exon skipping and/or partial intron retention events that occur either together or independently in different variants, apparently due to accumulation of inherited and acquired mutations. Cellular, biochemical, and oncogenic properties of full-length HAS1 (HAS1-FL) and HAS1 splice variants Va, Vb, and Vc (HAS1-Vs) are compared and characterized. When co-expressed, the properties of HAS1-Vs are dominant over those of HAS1-FL. HAS1-FL appears to be diffusely expressed in the cell, but HAS1-Vs are concentrated in the cytoplasm and/or Golgi apparatus. HAS1-Vs synthesize detectable de novo HA intracellularly. Each of the HAS1-Vs is able to relocalize HAS1-FL protein from diffuse cytoskeleton-anchored locations to deeper cytoplasmic spaces. This HAS1-Vs-mediated relocalization occurs through strong molecular interactions, which also serve to protect HAS1-FL from its otherwise high turnover kinetics. In co-transfected cells, HAS1-FL and HAS1-Vs interact with themselves and with each other to form heteromeric multiprotein assemblies. HAS1-Vc was found to be transforming in vitro and tumorigenic in vivo when introduced as a single oncogene to untransformed cells. The altered distribution and half-life of HAS1-FL, coupled with the characteristics of the HAS1-Vs suggest possible mechanisms whereby the aberrant splicing observed in human cancer may contribute to oncogenesis and disease progression.About 70–80% of human genes undergo alternative splicing, contributing to proteomic diversity and regulatory complexities in normal development (1). About 10% of mutations listed so far in the Human Gene Mutation Database (HGMD) of “gene lesions responsible for human inherited disease” were found to be located within splice sites. Furthermore, it is becoming increasingly apparent that aberrant splice variants, generated mostly due to splicing defects, play a key role in cancer. Germ line or acquired genomic changes (mutations) in/around splicing elements (2–4) promote aberrant splicing and aberrant protein isoforms.Hyaluronan (HA)³ is synthesized by three different plasma membrane-bound hyaluronan synthases (1, 2, and 3). HAS1 undergoes alternative and aberrant intronic splicing in multiple myeloma, producing truncated variants termed Va, Vb, and Vc (5, 6), which predicted for poor survival in a cohort of multiple myeloma patients (5). Our work suggests that this aberrant splicing arises due to inherited predispositions and acquired mutations in the HAS1 gene (7). Cancer-related, defective mRNA splicing caused by polymorphisms and/or mutations in splicing elements often results in inactivation of tumor suppressor activity (e.g. HRPT2 (8, 9), PTEN (10), MLHI (11–14), and ATR (15)) or generation of dominant negative inhibitors (e.g. CHEK2 (16) and VWOX (17)). In breast cancer, aberrantly spliced forms of progesterone and estrogen receptors are found (reviewed in Ref. 3). Intronic mutations inactivate p53 through aberrant splicing and intron retention (18). Somatic mutations with the potential to alter splicing are frequent in some cancers (19–25). Single nucleotide polymorphisms in the cyclin D1 proto-oncogene predispose to aberrant splicing and the cyclin D1b intronic splice variant (26–29). Cyclin D1b confers anchorage independence, is tumorogenic in vivo, and is detectable in human tumors (30), but as yet no clinical studies have confirmed an impact on outcome. On the other hand, aberrant splicing of HAS1 shows an association between aberrant splice variants and malignancy, suggesting that such variants may be potential therapeutic targets and diagnostic indicators (19, 31–33). Increased HA expression has been associated with malignant progression of multiple tumor types, including breast, prostate, colon, glioma, mesothelioma, and multiple myeloma (34). The three mammalian HA synthase (HAS) isoenzymes synthesize HA and are integral transmembrane proteins with a probable porelike structural assembly (35–39). Although in humans, the three HAS genes are located on different chromosomes (hCh19, hCh8, and hCh16, respectively) (40), they share a high degree of sequence homology (41, 42). HAS isoenzymes synthesize a different size range of HA molecules, which exhibit different functions (43, 44). HASs contribute to a variety of cancers (45–55). Overexpression of HASs promotes growth and/or metastatic development in fibrosarcoma, prostate, and mammary carcinoma, and the removal of the HA matrix from a migratory cell membrane inhibits cell movement (45, 53). HAS2 confers anchorage independence (56). Our work has shown aberrant HAS1 splicing in multiple myeloma (5) and Waldenstrom''s macroglobulinemia (6). HAS1 is overexpressed in colon (57), ovarian (58), endometrial (59), mesothelioma (60), and bladder cancers (61). A HAS1 splice variant is detected in bladder cancer (61).Here, we characterize molecular and biochemical characteristics of HAS1 variants (HAS1-Vs) (5), generated by aberrant splicing. Using transient transfectants and tagged HAS1 family constructs, we show that HAS1-Vs differ in cellular localization, de novo HA localization, and turnover kinetics, as compared with HAS1-FL, and dominantly influence HAS1-FL when co-expressed. HAS1-Vs proteins form intra- and intermolecular associations among themselves and with HAS1-FL, including covalent interactions and multimer formation. HAS1-Vc supports vigorous cellular transformation of NIH3T3 cells in vitro, and HAS1-Vc-transformed NIH3T3 cells are tumorogenic in vivo.     

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]     

Biochemical Characterization and Function of Complexes Formed by Hyaluronan and the Heavy Chains of Inter-��-inhibitor (HC��HA) Purified from Extracts of Human Amniotic Membrane

Clinically, amniotic membrane (AM) suppresses inflammation, scarring, and angiogenesis. AM contains abundant hyaluronan (HA) but its function in exerting these therapeutic actions remains unclear. Herein, AM was extracted sequentially with buffers A, B, and C, or separately by phosphate-buffered saline (PBS) alone. Agarose gel electrophoresis showed that high molecular weight (HMW) HA (an average of ∼3000 kDa) was predominantly extracted in isotonic Extract A (70.1 ± 6.0%) and PBS (37.7 ± 3.2%). Western blot analysis of these extracts with hyaluronidase digestion or NaOH treatment revealed that HMW HA was covalently linked with the heavy chains (HCs) of inter-α-inhibitor (IαI) via a NaOH-sensitive bond, likely transferred by the tumor necrosis factor-α stimulated gene-6 protein (TSG-6). This HC·HA complex (nHC·HA) could be purified from Extract PBS by two rounds of CsCl/guanidine HCl ultracentrifugation as well as in vitro reconstituted (rcHC·HA) by mixing HMW HA, serum IαI, and recombinant TSG-6. Consistent with previous reports, Extract PBS suppressed transforming growth factor-β1 promoter activation in corneal fibroblasts and induced mac ro phage apo pto sis. However, these effects were abolished by hyaluronidase digestion or heat treatment. More importantly, the effects were retained in the nHC·HA or rcHC·HA. These data collectively suggest that the HC·HA complex is the active component in AM responsible in part for clinically observed anti-inflammatory and anti-scarring actions.Hyaluronan (HA)⁴ is widely distributed in extracellular matrices, tissues, body fluids, and even in intracellular compartments (reviewed in Refs. 1 and 2). The molecular weight of HA ranges from 200 to 10,000 kDa depending on the source (3), but can also exist as smaller fragments and oligosaccharides under certain physiological or pathological conditions (1). Investigations over the last 15 years have suggested that low M_r HA can induce the gene expression of proinflammatory mediators and proangiogenesis, whereas high molecular weight (HMW) HA inhibits these processes (4–7).Several proteins have been shown to bind to HA (8) such as aggrecan (9), cartilage link protein (10), versican (11), CD44 (12, 13), inter-α-inhibitor (IαI) (14, 15), and tumor necrosis factor-α stimulated gene-6 protein (TSG-6) (16, 17). IαI consists of two heavy chains (HCs) (HC1 and HC2), both of which are linked through ester bonds to a chondroitin sulfate chain that is attached to the light chain, i.e. bikunin. Among all HA-binding proteins, only the HCs of IαI have been clearly demonstrated to be covalently coupled to HA (14, 18). However, TSG-6 has also been reported to form stable, possibly covalent, complexes with HA, either alone (19, 20) or when associated with HC (21).The formation of covalent bonds between HCs and HA is mediated by TSG-6 (22–24) where its expression is often induced by inflammatory mediators such as tumor necrosis factor-α and interleukin-1 (25, 26). TSG-6 is also expressed in inflammatory-like processes, such as ovulation (21, 27, 28) and cervical ripening (29). TSG-6 interacts with both HA (17) and IαI (21, 24, 30–33), and is essential for covalently transferring HCs on to HA (22–24). The TSG-6-mediated formation of the HC·HA complex has been demonstrated to play a crucial role in female fertility in mice. The HC·HA complex is an integral part of an expanded extracellular “cumulus” matrix around the oocyte, which plays a critical role in successful ovulation and fertilization in vivo (22, 34). HC·HA complexes have also been found at sites of inflammation (35–38) where its pro- or anti-inflammatory role remain arguable (39, 40).Immunostaining reveals abundant HA in the avascular stromal matrix of the AM (41, 42).⁵ In ophthalmology, cryopreserved AM has been widely used as a surgical graft for ocular surface reconstruction and exerts clinically observable actions to promote epithelial wound healing and to suppress inflammation, scarring, and angiogenesis (for reviews see Refs. 43–45). However, it is not clear whether HA in AM forms HC·HA complex, and if so whether such an HC·HA complex exerts any of the above therapeutic actions. To address these questions, we extracted AM with buffers of increasing salt concentration. Because HMW HA was found to form the HC·HA complex and was mainly extractable by isotonic solutions, we further purified it from the isotonic AM extract and reconstituted it in vitro from three defined components, i.e. HMW HA, serum IαI, and recombinant TSG-6. Our results showed that the HC·HA complex is an active component in AM responsible for the suppression of TGF-β1 promoter activity, linkable to the scarring process noted before by AM (46–48) and by the AM soluble extract (49), as well as for the promotion of macrophage death, linkable to the inflammatory process noted by AM (50) and the AM soluble extract (51).