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.     

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]     

Mtv-1 Superantigen Trafficks Independently of Major Histocompatibility Complex Class II Directly to the B-Cell Surface by the Exocytic Pathway

Presentation of the Mtv-1 superantigen (vSag1) to specific V_β-bearing T cells requires association with major histocompatibility complex class II molecules. The intracellular route by which vSag1 trafficks to the cell surface and the site of vSag1-class II complex assembly in antigen-presenting B lymphocytes have not been determined. Here, we show that vSag1 trafficks independently of class II to the plasma membrane by the exocytic secretory pathway. At the surface of B cells, vSag1 associates primarily with mature peptide-bound class II αβ dimers, which are stable in sodium dodecyl sulfate. vSag1 is unstable on the cell surface in the absence of class II, and reagents that alter the surface expression of vSag1 and the conformation of class II molecules affect vSag1 stimulation of superantigen reactive T cells.

T lymphocytes respond to peptide antigens presented by either major histocompatibility complex (MHC) class I or class II molecules. Many viruses have evolved sophisticated strategies that interfere with antigen presentation by infected cells in order to escape recognition by T lymphocytes. Most strategies studied rely on disrupting MHC class I presentation, either by affecting components of the processing machinery that generate and transport viral peptides into the endoplasmic reticulum (ER) or by retarding transport or targeting class I molecules into the degradation pathway (for a review, see reference 73).In contrast, mouse mammary tumor virus (MMTV) utilizes T-cell stimulation to promote its life cycle. MMTVs encode within their 3′ long terminal repeat a viral superantigen (vSag), and coexpression of the Sag glycoprotein with MHC class II molecules on the surface of virally infected B cells induces V_β-specific T-cell stimulation, generating an immune response which is critical for amplification of MMTV and ensures vertical transmission of virus to the next generation (13, 29, 30). In the absence of B cells, MHC class II, or Sag-reactive T cells, the infection is short-lived (5, 6, 24, 28). The assembly and functional expression of vSag-class II complexes are therefore essential to the viral life cycle. When inherited as germ line elements, Mtv proviruses expressing vSags during ontogeny trigger V_β-specific clonal elimination of immature T cells and profoundly shape the T-cell repertoire (for a review, see reference 1).vSags are type II integral membrane glycoproteins (14, 36). They possess up to six potential N-linked glycosylation sites, and carbohydrate addition is essential for vSag stability and activity (45). Their protein sequence is highly conserved among all MMTV strains except at the C-terminal 29 to 32 residues, which vary and confer T-cell V_β specificity (77). Biochemical analyses of vSag7 (minor lymphocyte stimulating locus 1, Mls-1^a) molecular forms after transfection into a murine B-cell line have identified a predominant 45-kDa endo-β-N-acetylglucosaminidase H (endo H)-sensitive ER-resident glycoprotein, as well as multiple highly glycosylated forms (74). It is thought that an 18-kDa C-terminal fragment binds MHC class II products (75). It has also been suggested that vSags associate weakly with class II in the ER and that proteolytic processing is required for the efficient assembly of vSag-class II complexes for presentation to T cells (46, 49, 75). As yet, the intracellular route that vSags take to the cell surface, the compartment in which they bind class II, and whether they associate with peptide-loaded class II dimers have been enigmatic.Newly synthesized MHC class II αβ heterodimers assemble with invariant chain (Ii), a type II integral membrane protein, to form an oligomeric complex in the ER (37). Ii prevents class II heterodimers from binding peptides in the ER and Golgi complex (55), and signals in its cytoplasmic tail sort the complex into the endocytic pathway (4, 42). In this acidic, protease-rich compartment, Ii is degraded and class II binds antigenic peptides. After the formation of peptide-class II dimers, the complexes are exported to the plasma membrane (8, 48). In the absence of Ii, class II αβ heterodimers exhibit defective post-ER transport, and their conversion into functionally mature, sodium dodecyl sulfate (SDS)-stable compact dimers by peptide antigens is affected (7, 16, 22, 70).A specialized endosomal compartment where class II peptide loading occurs, termed the MHC class II-enriched compartment (MIIC or CIIV), has been found recently in antigen-presenting cells (2, 50, 53, 58, 68, 71). Whether nascent Ii-class II complexes traffic directly to the MIIC from the trans-Golgi network (TGN) or transit first to early endosomes, either directly or via the cell surface, before entering late endocytic vesicles and MIIC is still under debate (26, 56, 57). Transport by all these routes most probably occurs to ensure the capture and loading of antigenic peptides throughout the endocytic pathway (12). MIIC vesicles are positive for lysosome-associated membrane proteins (LAMPs) and cathepsin D and are enriched for HLA-DM or H-2M (18, 32, 59), proteins that facilitate the catalytic exchange of class II-associated invariant peptide chain (CLIP) for antigenic peptides (19, 61, 62). The ultrastructural colocalization of DM with intracellular peptide-class II complexes suggests that the MIIC is a main site where class II dimers bind exogenous and endogenous peptide antigens (47, 58).Determining the route by which vSag protein(s) trafficks to the cell surface and the cellular location where vSag1 processing and assembly with class II molecules occurs is central to understanding the mechanism whereby vSags activate T cells to maintain the viral life cycle. It has been unclear whether vSags traffic independently by the constitutive exocytic pathway or with class II and Ii to the MIIC before reaching the cell surface. Reagents that alter class II expression have been shown to affect vSag presentation (43, 46). Furthermore, mice lacking Ii show reduced intrathymic V_β-specific T-cell deletion (70), suggesting that Ii may play a role, either by ensuring proper maturation of class II dimers or by targeting vSag-class II complexes to the MIIC, in promoting efficient vSag-induced immune responses.To investigate these issues, we used immunochemical detection of vSag1 protein in combination with subcellular fractionation and surface reexpression assays. We show that class II is required for stable vSag1 surface expression. vSag1 trafficks directly to the cell surface independently of class II, and reagents that alter the conversion of newly synthesized class II into peptide-loaded SDS-stable dimers affect functional vSag1 surface expression.     

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]     

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]     

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.