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.     

PTIP Regulates 53BP1 and SMC1 at the DNA Damage Sites

Although PTIP is implicated in the DNA damage response, through interactions with 53BP1, the function of PTIP in the DNA damage response remain elusive. Here, we show that RNF8 controls DNA damage-induced nuclear foci formation of PTIP, which in turn regulates 53BP1 localization to the DNA damage sites. In addition, SMC1, a substrate of ATM, could not be phosphorylated at the DNA damage sites in the absence of PTIP. The PTIP-dependent pathway is important for DNA double strand breaks repair and DNA damage-induced intra-S phase checkpoint activation. Taken together, these results suggest that the role of PTIP in the DNA damage response is downstream of RNF8 and upstream of 53BP1. Thus, PTIP regulates 53BP1-dependent signaling pathway following DNA damage.The DNA damage response pathways are signal transduction pathways with DNA damage sensors, mediators, and effectors, which are essential for maintaining genomic stability (1–3). Following DNA double strand breaks, histone H2AX at the DNA damage sites is rapidly phosphorylated by ATM/ATR/DNAPK (4–10), a family homologous to phosphoinositide 3-kinases (11, 12). Subsequently, phospho-H2AX (γH2AX) provides the platform for accumulation of a larger group of DNA damage response factors, such as MDC1, BRCA1, 53BP1, and the MRE11·RAD50·NBS1 complex (13, 14), at the DNA damage sites. Translocalization of these proteins to the DNA double strand breaks (DSBs)³ facilitates DNA damage checkpoint activation and enhances the efficiency of DNA damage repair (14, 15).Recently, PTIP (Pax2 transactivation domain-interacting protein, or Paxip) has been identified as a DNA damage response protein and is required for cell survival when exposed to ionizing radiation (IR) (1, 16–18). PTIP is a 1069-amino acid nuclear protein and has been originally identified in a yeast two-hybrid screening as a partner of Pax2 (19). Genetic deletion of the PTIP gene in mice leads to early embryonic lethality at embryonic day 8.5, suggesting that PTIP is essential for early embryonic development (20). Structurally, PTIP contains six tandem BRCT (BRCA1 carboxyl-terminal) domains (16–18, 21). The BRCT domain is a phospho-group binding domain that mediates protein-protein interactions (17, 22, 23). Interestingly, the BRCT domain has been found in a large number of proteins involved in the cellular response to DNA damages, such as BRCA1, MDC1, and 53BP1 (7, 24–29). Like other BRCT domain-containing proteins, upon exposure to IR, PTIP forms nuclear foci at the DSBs, which is dependent on its BRCT domains (16–18). By protein affinity purification, PTIP has been found in two large complexes. One includes the histone H3K4 methyltransferase ALR and its associated cofactors, the other contains DNA damage response proteins, including 53BP1 and SMC1 (30, 31). Further experiments have revealed that DNA damage enhances the interaction between PTIP and 53BP1 (18, 31).To elucidate the DNA damage response pathways, we have examined the upstream and downstream partners of PTIP. Here, we report that PTIP is downstream of RNF8 and upstream of 53BP1 in response to DNA damage. Moreover, PTIP and 53BP1 are required for the phospho-ATM association with the chromatin, which phosphorylates SMC1 at the DSBs. This PTIP-dependent pathway is involved in DSBs repair.     

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]     

Algorithms for Finding Small Attractors in Boolean Networks

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

Proinsulin and the Genetics of Diabetes Mellitus

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

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

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

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]     

In Vivo Identification of Nuclear Factors Interacting with the Conserved Elements of Box C/D Small Nucleolar RNAs

The U16 small nucleolar RNA (snoRNA) is encoded by the third intron of the L1 (L4, according to the novel nomenclature) ribosomal protein gene of Xenopus laevis and originates from processing of the pre-mRNA in which it resides. The U16 snoRNA belongs to the box C/D snoRNA family, whose members are known to assemble in ribonucleoprotein particles (snoRNPs) containing the protein fibrillarin. We have utilized U16 snoRNA in order to characterize the factors that interact with the conserved elements common to the other members of the box C/D class. In this study, we have analyzed the in vivo assembly of U16 snoRNP particles in X. laevis oocytes and identified the proteins which interact with the RNA by label transfer after UV cross-linking. This analysis revealed two proteins, of 40- and 68-kDa apparent molecular size, which require intact boxes C and D together with the conserved 5′,3′-terminal stem for binding. Immunoprecipitation experiments showed that the p40 protein corresponds to fibrillarin, indicating that this protein is intimately associated with the RNA. We propose that fibrillarin and p68 represent the RNA-binding factors common to box C/D snoRNPs and that both proteins are essential for the assembly of snoRNP particles and the stabilization of the snoRNA.One of the most interesting recent findings related to ribosome biogenesis has been the identification of a large number of small RNAs localized in the nucleolus (snoRNAs). So far, more than 60 snoRNAs have been identified in vertebrates (17), and more than 30 have been identified in yeast (2). The total number of snoRNAs is not known, but it is likely to be close to 200 (33, 38). These snoRNAs, with the exception of the mitochondrial RNA processing (MRP) species (38), can be grouped into two major families on the basis of conserved structural and sequence elements. The first group includes molecules referred to as box C/D snoRNAs, whereas the second one comprises the species belonging to the box H/ACA family (2, 15).The two families differ in many aspects. The box C/D snoRNAs are functionally heterogeneous. Most of them function as antisense RNAs in site-specific ribose methylation of the pre-rRNA (1, 10, 17, 26); a minority have been shown to play a direct role in pre-rRNA processing in both yeast and metazoan cells (11, 21). The box C/D snoRNAs play their role by means of unusually long (up to 21 contiguous nucleotides) regions of complementarity to highly conserved sequences of 28S and 18S rRNAs (1). In contrast, several members of the H/ACA RNA family have been shown to direct site-specific isomerization of uridines into pseudouridines and to display shorter regions of complementarity to rRNA (14, 24). Mutational analysis suggests that H/ACA snoRNAs can also play a role as antisense RNAs by base pairing with complementary regions on rRNA (15, 24).Another difference between the two families can be seen by comparison of secondary structures. A Y-shaped motif, where a 5′,3′-terminal stem adjoins the C and D conserved elements, has been proposed for many box C/D snoRNAs (16, 26, 40, 42), whereas box H/ACA snoRNAs have been proposed to fold into two conserved hairpin structures connected by a single-stranded hinge region, followed by a short 3′ tail (15).Despite these differences, analogies have been found in the roles played by the conserved box elements. Mutational analysis and competition experiments indicated that C/D and H/ACA boxes are required both for processing and stable accumulation of the mature snoRNA, suggesting that they represent binding sites for specific trans-acting factors (2, 3, 8, 15, 16, 28, 36, 41).All snoRNAs are associated with proteins to form specific ribonucleoparticles (snoRNPs). The study of these particles began only recently, and so far, very few aspects of their structure and biosynthesis have been clarified. The only detailed analysis performed was on the mammalian U3 (19) and the yeast snR30 (20) snoRNPs. Of the identified components, a few appear to be more general factors: fibrillarin, which was shown to be associated with C/D snoRNPs (3, 4, 8, 13, 28, 31, 39), and the nucleolar protein GAR1, which was found associated with H/ACA snoRNAs in yeast (20). Just as the study of small nuclear RNP (snRNP) particles was crucial to the understanding of the splicing process, a detailed structural and functional analysis of snoRNP particles will be essential to elucidate the complex process of ribosome biosynthesis.In this study, we have analyzed the snoRNP assembly of wild-type and mutant U16 snoRNAs by following the kinetics of complex formation in the in vivo system of the Xenopus laevis oocyte. By a UV cross-linking technique, we have identified two proteins, of 40- and 68-kDa apparent molecular mass, which require intact boxes C and D together with the terminal stem for their binding. The 40-kDa species is specifically recognized by fibrillarin antibodies, indicating that this protein is intimately associated with the RNA.