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]     

Semaphorin3A Signaling Mediated by Fyn-dependent Tyrosine Phosphorylation of Collapsin Response Mediator Protein 2 at Tyrosine 32

Collapsin response mediator protein 2 (CRMP2) is an intracellular protein that mediates signaling of Semaphorin3A (Sema3A), a repulsive axon guidance molecule. Fyn, a Src-type tyrosine kinase, is involved in the Sema3A signaling. However, the relationship between CRMP2 and Fyn in this signaling pathway is still unknown. In our research, we demonstrated that Fyn phosphorylated CRMP2 at Tyr³² residues in HEK293T cells. Immunohistochemical analysis using a phospho-specific antibody at Tyr³² of CRMP showed that Tyr³²-phosphorylated CRMP was abundant in the nervous system, including dorsal root ganglion neurons, the molecular and Purkinje cell layer of adult cerebellum, and hippocampal fimbria. Overexpression of a nonphosphorylated mutant (Tyr³² to Phe³²) of CRMP2 in dorsal root ganglion neurons interfered with Sema3A-induced growth cone collapse response. These results suggest that Fyn-dependent phosphorylation of CRMP2 at Tyr³² is involved in Sema3A signaling.Collapsin response mediator proteins (CRMPs)⁴ have been identified as intracellular proteins that mediate Semaphorin3A (Sema3A) signaling in the nervous system (1). CRMP2 is one of the five members of the CRMP family. CRMPs also mediate signal transduction of NT3, Ephrin, and Reelin (2–4). CRMPs interact with several intracellular molecules, including tubulin, Numb, kinesin1, and Sra1 (5–8). CRMPs are involved in axon guidance, axonal elongation, cell migration, synapse maturation, and the generation of neuronal polarity (1, 2, 4, 5).CRMP family proteins are known to be the major phosphoproteins in the developing brain (1, 9). CRMP2 is phosphorylated by several Ser/Thr kinases, such as Rho kinase, cyclin-dependent kinase 5 (Cdk5), and glycogen synthase kinase 3β (GSK3β) (2, 10–13). The phosphorylation sites of CRMP2 by these kinases are clustered in the C terminus and have already been identified. Rho kinase phosphorylates CRMP2 at Thr⁵⁵⁵ (10). Cdk5 phosphorylates CRMP2 at Ser⁵²², and this phosphorylation is essential for sequential phosphorylations by GSK3β at Ser⁵¹⁸, Thr⁵¹⁴, and Thr⁵⁰⁹ (2, 11–13). These phosphorylations disrupt the interaction of CRMP2 with tubulin or Numb (2, 3, 13). The sequential phosphorylation of CRMP2 by Cdk5 and GSK3β is an essential step in Sema3A signaling (11, 13). Furthermore, the neurofibrillary tangles in the brains of people with Alzheimer disease contain hyperphosphorylated CRMP2 at Thr⁵⁰⁹, Ser⁵¹⁸, and Ser⁵²² (14, 15).CRMPs are also substrates of several tyrosine kinases. The phosphorylation of CRMP2 by Fes/Fps and Fer has been shown to be involved in Sema3A signaling (16, 17). Phosphorylation of CRMP2 at Tyr⁴⁷⁹ by a Src family tyrosine kinase Yes regulates CXCL12-induced T lymphocyte migration (18). We reported previously that Fyn is involved in Sema3A signaling (19). Fyn associates with PlexinA2, one of the components of the Sema3A receptor complex. Fyn also activates Cdk5 through the phosphorylation at Tyr¹⁵ of Cdk5 (19). In dorsal root ganglion (DRG) neurons from fyn-deficient mice, Sema3A-induced growth cone collapse response is attenuated compared with control mice (19). Furthermore, we recently found that Fyn phosphorylates CRMP1 and that this phosphorylation is involved in Reelin signaling (4). Although it has been shown that CRMP2 is involved in Sema3A signaling (1, 11, 13), the relationship between Fyn and CRMP2 in Sema3A signaling and the tyrosine phosphorylation site(s) of CRMPs remain unknown.Here, we show that Fyn phosphorylates CRMP2 at Tyr³². Using a phospho-specific antibody against Tyr³², we determined that the residue is phosphorylated in vivo. A nonphosphorylated mutant CRMP2Y32F inhibits Sema3A-induced growth cone collapse. These results indicate that tyrosine phosphorylation by Fyn at Tyr³² is involved in Sema3A signaling.     

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.     

Sgt1 Dimerization Is Negatively Regulated by Protein Kinase CK2-mediated Phosphorylation at Ser361

The kinetochore, which consists of centromere DNA and structural proteins, is essential for proper chromosome segregation in eukaryotes. In budding yeast, Sgt1 and Hsp90 are required for the binding of Skp1 to Ctf13 (a component of the core kinetochore complex CBF3) and therefore for the assembly of CBF3. We have previously shown that Sgt1 dimerization is important for this kinetochore assembly mechanism. In this study, we report that protein kinase CK2 phosphorylates Ser³⁶¹ on Sgt1, and this phosphorylation inhibits Sgt1 dimerization.The kinetochore is a structural protein complex located in the centromeric region of the chromosome coupled to spindle microtubules (1, 2). The kinetochore generates a signal to arrest cells during mitosis when it is not properly attached to microtubules, thereby preventing chromosome missegregation, which can lead to aneuploidy (3, 4). The molecular structure of the kinetochore complex of the budding yeast Saccharomyces cerevisiae has been well characterized; it is composed of more than 70 proteins, many of which are conserved in mammals (2).The centromere DNA in the budding yeast is a 125-bp region that contains three conserved regions, CDEI, CDEII, and CDEIII (5, 6). CDEIII (25 bp) is essential for centromere function (7) and is bound to a key component of the centromere, the CBF3 complex. The CBF3 complex contains four proteins, Ndc10, Cep3, Ctf13 (8–15), and Skp1 (14, 15), all essential for viability. Mutations in any of the CBF3 proteins abolish the ability of CDEIII to bind to CBF3 (16, 17). All of the kinetochore proteins, except the CDEI-binding Cbf1 (18–20), localize to the kinetochores in a CBF3-dependent manner (2). Thus, CBF3 is a fundamental kinetochore complex, and its mechanism of assembly is of great interest.We have previously found that Sgt1 and Skp1 activate Ctf13; thus, they are required for assembly of the CBF3 complex (21). The molecular chaperone Hsp90 is also required to form the active Ctf13-Skp1 complex (22). Sgt1 has two highly conserved motifs that are required for protein-protein interaction: the tetratricopeptide repeat (21) and the CHORD protein and Sgt1-specific motif. We and others have found that both domains are important for the interaction of Sgt1 with Hsp90 (23–26), which is required for assembly of the core kinetochore complex. This interaction is an initial step in kinetochore activation (24, 26, 27), which is conserved between yeast and humans (28, 29).We have recently shown that Sgt1 dimerization is important for Sgt1-Skp1 binding and therefore for kinetochore assembly (30). In this study, we have found that protein kinase CK2 phosphorylates Sgt1 at Ser³⁶¹, and this phosphorylation inhibits Sgt1 dimerization. Therefore, CK2 appears to regulate kinetochore assembly negatively in budding yeast.     

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]     

Dbf4-Cdc7 Phosphorylation of Mcm2 Is Required for Cell Growth

The Dbf4-Cdc7 kinase (DDK) is required for the activation of the origins of replication, and DDK phosphorylates Mcm2 in vitro. We find that budding yeast Cdc7 alone exists in solution as a weakly active multimer. Dbf4 forms a likely heterodimer with Cdc7, and this species phosphorylates Mcm2 with substantially higher specific activity. Dbf4 alone binds tightly to Mcm2, whereas Cdc7 alone binds weakly to Mcm2, suggesting that Dbf4 recruits Cdc7 to phosphorylate Mcm2. DDK phosphorylates two serine residues of Mcm2 near the N terminus of the protein, Ser-164 and Ser-170. Expression of mcm2-S170A is lethal to yeast cells that lack endogenous MCM2 (mcm2Δ); however, this lethality is rescued in cells harboring the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Mcm2 is required for cell growth.The Cdc7 protein kinase is required throughout the yeast S phase to activate origins (1, 2). The S phase cyclin-dependent kinase also activates yeast origins of replication (3–5). It has been proposed that Dbf4 activates Cdc7 kinase in S phase, and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6). However, it is not known how Dbf4-Cdc7 (DDK)² acts during S phase to trigger the initiation of DNA replication. DDK has homologs in other eukaryotic species, and the role of Cdc7 in activation of replication origins during S phase may be conserved (7–10).The Mcm2-7 complex functions with Cdc45 and GINS to unwind DNA at a replication fork (11–15). A mutation of MCM5 (mcm5-bob1) bypasses the cellular requirements for DBF4 and CDC7 (16), suggesting a critical physiologic interaction between Dbf4-Cdc7 and Mcm proteins. DDK phosphorylates Mcm2 in vitro with proteins purified from budding yeast (17, 18) or human cells (19). Furthermore, there are mutants of MCM2 that show synthetic lethality with DBF4 mutants (6, 17), suggesting a biologically relevant interaction between DBF4 and MCM2. Nevertheless, the physiologic role of DDK phosphorylation of Mcm2 is a matter of dispute. In human cells, replacement of MCM2 DDK-phosphoacceptor residues with alanines inhibits DNA replication, suggesting that Dbf4-Cdc7 phosphorylation of Mcm2 in humans is important for DNA replication (20). In contrast, mutation of putative DDK phosphorylation sites at the N terminus of Schizosaccharomyces pombe Mcm2 results in viable cells, suggesting that phosphorylation of S. pombe Mcm2 by DDK is not critical for cell growth (10).In budding yeast, Cdc7 is present at high levels in G₁ and S phase, whereas Dbf4 levels peak in S phase (18, 21, 22). Furthermore, budding yeast DDK binds to chromatin during S phase (6), and it has been shown that Dbf4 is required for Cdc7 binding to chromatin in budding yeast (23, 24), fission yeast (25), and Xenopus (9). Human and fission yeast Cdc7 are inert on their own (7, 8), but Dbf4-Cdc7 is active in phosphorylating Mcm proteins in budding yeast (6, 26), fission yeast (7), and human (8, 10). Based on these data, it has been proposed that Dbf4 activates Cdc7 kinase in S phase and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6, 9, 18, 21–24). However, a mechanistic analysis of how Dbf4 activates Cdc7 has not yet been accomplished. For example, the multimeric state of the active Dbf4-Cdc7 complex is currently disputed. A heterodimer of fission yeast Cdc7 (Hsk1) in complex with fission yeast Dbf4 (Dfp1) can phosphorylate Mcm2 (7). However, in budding yeast, oligomers of Cdc7 exist in the cell (27), and Dbf4-Cdc7 exists as oligomers of 180 and 300 kDa (27).DDK phosphorylates the N termini of human Mcm2 (19, 20, 28), human Mcm4 (10), budding yeast Mcm4 (26), and fission yeast Mcm6 (10). Although the sequences of the Mcm N termini are poorly conserved, the DDK sites identified in each study have neighboring acidic residues. The residues of budding yeast Mcm2 that are phosphorylated by DDK have not yet been identified.In this study, we find that budding yeast Cdc7 is weakly active as a multimer in phosphorylating Mcm2. However, a low molecular weight form of Dbf4-Cdc7, likely a heterodimer, has a higher specific activity for phosphorylation of Mcm2. Dbf4 or DDK, but not Cdc7, binds tightly to Mcm2, suggesting that Dbf4 recruits Cdc7 to Mcm2. DDK phosphorylates two serine residues of Mcm2, Ser-164 and Ser-170, in an acidic region of the protein. Mutation of Ser-170 is lethal to yeast cells, but this phenotype is rescued by the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Ser-170 of Mcm2 is required for budding yeast growth.     

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]     

A Ca2+ Release-activated Ca2+ (CRAC) Modulatory Domain (CMD) within STIM1 Mediates Fast Ca2+-dependent Inactivation of ORAI1 Channels

STIM1 and ORAI1, the two limiting components in the Ca²⁺ release-activated Ca²⁺ (CRAC) signaling cascade, have been reported to interact upon store depletion, culminating in CRAC current activation. We have recently identified a modulatory domain between amino acids 474 and 485 in the cytosolic part of STIM1 that comprises 7 negatively charged residues. A STIM1 C-terminal fragment lacking this domain exhibits enhanced interaction with ORAI1 and 2–3-fold higher ORAI1/CRAC current densities. Here we focused on the role of this CRAC modulatory domain (CMD) in the fast inactivation of ORAI1/CRAC channels, utilizing the whole-cell patch clamp technique. STIM1 mutants either with C-terminal deletions including CMD or with 7 alanines replacing the negative amino acids within CMD gave rise to ORAI1 currents that displayed significantly reduced or even abolished inactivation when compared with STIM1 mutants with preserved CMD. Consistent results were obtained with cytosolic C-terminal fragments of STIM1, both in ORAI1-expressing HEK 293 cells and in RBL-2H3 mast cells containing endogenous CRAC channels. Inactivation of the latter, however, was much more pronounced than that of ORAI1. The extent of inactivation of ORAI3 channels, which is also considerably more prominent than that of ORAI1, was also substantially reduced by co-expression of STIM1 constructs missing CMD. Regarding the dependence of inactivation on Ca²⁺, a decrease in intracellular Ca²⁺ chelator concentrations promoted ORAI1 current fast inactivation, whereas Ba²⁺ substitution for extracellular Ca²⁺ completely abrogated it. In summary, CMD within the STIM1 cytosolic part provides a negative feedback signal to Ca²⁺ entry by triggering fast Ca²⁺-dependent inactivation of ORAI/CRAC channels.The Ca²⁺ release-activated Ca²⁺ (CRAC)⁵ channel is one of the best characterized store-operated entry pathways (1–7). Substantial efforts have led to identification of two key components of the CRAC channel machinery: the stromal interaction molecule 1 (STIM1), which is located in the endoplasmic reticulum and acts as a Ca²⁺ sensor (8–10), and ORAI1/CRACM1, the pore-forming subunit of the CRAC channel (11–13). Besides ORAI1, two further homologues named ORAI2 and ORAI3 belong to the ORAI channel family (12, 14).STIM1 senses endoplasmic reticulum store depletion primarily by its luminal EF-hand in its N terminus (8, 15), redistributes close to the plasma membrane, where it forms puncta-like structures, and co-clusters with ORAI1, leading to inward Ca²⁺ currents (12, 16–19). The STIM1 C terminus, located in the cytosol, contains two coiled-coil regions overlapping with an ezrin-radixin-moesin (ERM)-like domain followed by a serine/proline- and a lysine-rich region (2, 8, 20–22). Three recent studies have described the essential ORAI-activating region within the ERM domain, termed SOAR (Stim ORAI-activating region) (23), OASF (ORAI-activating small fragment) (24), and CAD (CRAC-activating domain) (25), including the second coiled coil domain and the following ∼55 amino acids. We and others have provided evidence that store depletion leads to a dynamic coupling of STIM1 to ORAI1 (26–28) that is mediated by a direct interaction of the STIM1 C terminus with ORAI1 C terminus probably involving the putative coiled-coil domain in the latter (27).Furthermore, different groups have proven that the C terminus of STIM1 is sufficient to activate CRAC as well as ORAI1 channels independent of store depletion (22–25, 27, 29). We have identified that OASF-(233–474) or shorter fragments exhibit further enhanced coupling to ORAI1 resulting in 3-fold increased constitutive Ca²⁺ currents. A STIM1 fragment containing an additional cluster of anionic amino acids C-terminal to position 474 displays weaker interaction with ORAI1 as well as reduced Ca²⁺ current comparable with that mediated by wild-type STIM1 C terminus. Hence, we have suggested that these 11 amino acids (474–485) act in a modulatory manner onto ORAI1; however, their detailed mechanistic impact within the STIM1/ORAI1 signaling machinery has remained so far unclear.In this study, we focused on the impact of this negative cluster on fast inactivation of STIM1-mediated ORAI Ca²⁺ currents. Lis et al. (30) have shown that all three ORAI homologues display distinct inactivation profiles, where ORAI2 and ORAI3 show a much more pronounced fast inactivation than ORAI1. Moreover, it has been reported (31) that different expression levels of STIM1 to ORAI1 affect the properties of CRAC current inactivation. Yamashita et al. (32) have demonstrated a linkage between the selectivity filter of ORAI1 and its Ca²⁺-dependent fast inactivation. Here we provide evidence that a cluster of acidic residues within the C terminus of STIM1 is involved in the fast inactivation of ORAI1 and further promotes that of ORAI3 and native CRAC currents.