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.     

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.     

Human Disease-causing Mutations Disrupt an N-C-terminal Interaction and Channel Function of Bestrophin 1

Mutations in the human bestrophin 1 (hBest1) chloride channel cause Best vitelliform macular dystrophy. Although mutations in its transmembrane domains were found to alter biophysical properties of the channel, the mechanism for disease-causing mutations in its N and C termini remains elusive. We hypothesized that these mutations lead to channel dysfunction through disruption of an N-C-terminal interaction. Here, we present data demonstrating that hBest1 N and C termini indeed interact both in vivo and in vitro. In addition, using a spectrum-based fluorescence resonance energy transfer method, we showed that functional hBest1 channels in the plasma membrane were multimers. Disease-causing mutations in the N terminus (R19C, R25C, and K30C) and the C terminus (G299E, D301N, and D312N) caused channel dysfunction and disruption of the N-C interaction. Consistent with the functional and biochemical results, mutants D301N and D312N clearly reduced fluorescence resonance energy transfer signal, indicating that the N-C interaction was indeed perturbed. These results suggest that hBest1 functions as a multimer in the plasma membrane, and disruption of the N-C interaction by mutations leads to hBest1 channel dysfunction.Extensive evidence has shown that bestrophins are anion channels when expressed heterologously in cell lines (1–9) and represent a type of endogenous calcium-activated chlorine channel in several cell types (10–12). Mutations disrupting hBest1 channel function lead to Best vitelliform macular dystrophy (Best disease) (13). Bestrophin 1 channels have been reported to exist as oligomers (2, 14). Results from immunoprecipitation experiments suggested that hBest1 exists as tetramers or pentamers when expressed heterologously (2), whereas a hydrodynamic study based on the calculations of the molecular mass of the Best1-detergent complex concluded that native porcine Best1 channels are dimers (14). Besides the uncertainty about the subunit stoichiometry, it is still unclear how bestrophin subunits assemble and interact in functional channels.Understanding how disease-causing mutations might adversely affect the assembly and interaction of bestrophin subunits is important for elucidating the molecular mechanism underlying the Best disease. The hBest1 protein contains six transmembrane domains (TMD1–6) with intracellularly located N and C termini. There are three “hot spots” of disease-causing mutations in hBest1: the N-terminal region, TMD2, and the C terminus proximal to TMD6 (13, 15). Multiple lines of evidence support that mutations in TMD2 alter the biophysical properties of the channel (3, 4, 7). How mutations distributed around the N- and C-terminal regions cause Best disease is less clear, although strong evidence recently suggested that certain mutations in the C terminus disrupt hBest1 channel gating by Ca²⁺ (16). N-C-terminal interaction of multimeric channels has been demonstrated to be involved in the activation of many channel types including inward rectifier K⁺ channels (17) and cyclic nucleotide-gated channels (18–20). In the present study we tested the hypothesis that an interaction between the N and C termini plays an important role in normal hBest1 channel function, and weakening or disruption of this interaction by mutations leads to channel dysfunction.To test our hypothesis, we introduced disease-causing mutations into the N- or C-terminal regions of hBest1 expressed in HEK293 cells and tested the effects of mutations on subunit interaction and channel function with a combination of electrophysiological, biochemical, and optical methods. We found that these mutations not only disrupted channel function but also caused N- and C-terminal dissociation both in vitro and in vivo. The findings suggest a novel molecular mechanism for the mutations in hBest1 to cause human Best disease.     

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.     

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]     

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.     

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]