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]     

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.     

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]     

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]     

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]     

Inflammasome-dependent Caspase-1 Activation in Cervical Epithelial Cells Stimulates Growth of the Intracellular Pathogen Chlamydia trachomatis

Inflammasomes have been extensively characterized in monocytes and macrophages, but not in epithelial cells, which are the preferred host cells for many pathogens. Here we show that cervical epithelial cells express a functional inflammasome. Infection of the cells by Chlamydia trachomatis leads to activation of caspase-1, through a process requiring the NOD-like receptor family member NLRP3 and the inflammasome adaptor protein ASC. Secretion of newly synthesized virulence proteins from the chlamydial vacuole through a type III secretion apparatus results in efflux of K⁺ through glibenclamide-sensitive K⁺ channels, which in turn stimulates production of reactive oxygen species. Elevated levels of reactive oxygen species are responsible for NLRP3-dependent caspase-1 activation in the infected cells. In monocytes and macrophages, caspase-1 is involved in processing and secretion of pro-inflammatory cytokines such as interleukin-1β. However, in epithelial cells, which are not known to secrete large quantities of interleukin-1β, caspase-1 has been shown previously to enhance lipid metabolism. Here we show that, in cervical epithelial cells, caspase-1 activation is required for optimal growth of the intracellular chlamydiae.Chlamydia trachomatis is the most common cause of bacterial sexually transmitted disease in the United States, and it is the leading cause of preventable blindness in the world (1–5). Untreated, C. trachomatis infection in women can cause pelvic inflammatory disease, which can lead to infertility and ectopic pregnancy because of scarring of the ovaries and the Fallopian tubes (6). Infection by the lymphogranuloma venereum (LGV)² strain of C. trachomatis, which has become more common in North America and Europe (7, 8), is characterized by swelling and inflammation of the lymph nodes in the groin (9).Chlamydiae are intracellular pathogens that preferentially infect epithelial mucosa and have a biphasic infection cycle (10). A metabolically inactive form, the elementary body, infects the epithelial host cells through entry vesicles that avoid fusion with host cell lysosomes and develop into a membrane-bound inclusion (11–13). Despite their intravacuolar localization, chlamydiae are still able to acquire nutrients from the host cell and interact with host-cell signaling pathways (13–23). Within a few hours, the elementary bodies differentiate into larger, metabolically active reticulate bodies, which proliferate but are noninfectious. Depending on the strain of C. trachomatis, the reticulate bodies transform back into elementary bodies after 1–3 days and are released into the extracellular medium to infect other cells (11, 24, 25). Chlamydial species possess a type III secretion (T3S) system that secretes bacterial virulence factors into host cell cytosol and may control interactions between the inclusion and host-cell compartments (26).Long before the adaptive immune response is activated, infected epithelial cells produce proinflammatory cytokines and chemokines, including interleukin (IL)-6, IL-8, and granulocyte-macrophage colony-stimulating factor (27), which recruit neutrophils to the site of infection and activate other immune effector cells. However, in many cases the immune system fails to clear the infection, and the chronic release of cytokines becomes a major contributor to the scarring and damage associated with the infection (28–30).The innate immune response during C. trachomatis infection is initiated by chlamydial pathogen-associated molecular patterns, including lipopolysaccharides, which bind to pattern recognition receptors such as Toll-like receptors and cytosolic NOD-like receptors (NLRs), ultimately promoting pro-inflammatory cytokine gene expression and secretion of the cytokine proteins (31–37). However, secretion of the key pro-inflammatory cytokine IL-1β is tightly regulated (38). First, pro-IL-1β is produced following activation of pattern recognition receptor, and the precursor is then cleaved into the mature form by the pro-inflammatory cysteine protease, caspase-1 (also known as interleukin-1 converting enzyme or ICE). The mechanism by which caspase-1 is activated in response to infection or tissue damage was found to be modulated by a macromolecular protein complex termed the “inflammasome,” which consists of an NLR family member, an adaptor protein (apoptosis-associated speck-like protein containing a caspase activation recruitment domain or ASC), and an inactive caspase-1 precursor (pro-caspase-1) (39, 40). Previous studies demonstrated that IL-1β is produced in response to chlamydial infection in dendritic cells, macrophages, and monocytes (41–44). Moreover, C. trachomatis or Chlamydia caviae infection activates caspase-1 in epithelial cells or monocytes (43, 45, 46). However, whether caspase-1 activation during chlamydial infection requires the formation of an inflammasome remains unclear.Previous studies have shown that different pathogens can cause inflammasome-mediated caspase-1 activation in macrophages and monocytes (47). However, epithelial cells lining mucosal surfaces are not only the preferred target for chlamydial infection and other intracellular pathogens but also play an important role in early host immune response to infection by secreting proinflammatory cytokines and chemokines (27). Although epithelial cells are not known to secrete large amounts of IL-1β, inflammasome-dependent caspase-1 activation in epithelial cells is known to contribute to lipid metabolism and membrane regeneration in epithelial cells damaged by the membrane-disrupting toxin, aerolysin (48). As lipids are sorted from the Golgi apparatus to the chlamydial inclusion (13, 15, 49), we therefore investigated whether C. trachomatis induces caspase-1 activation in epithelial cells via the assembly of an inflammasome. We demonstrated that C. trachomatis-induced caspase-1 activation is mediated by an inflammasome containing the NLR member, NLRP3. Several studies have demonstrated the involvement of T3S apparatus in inflammasome-mediated caspase-1 activation by different pathogens in macrophages and monocytes (50–56). Therefore, we further investigated the mechanism by which C. trachomatis triggers the formation of the NLRP3 inflammasome. Our results showed that metabolically active chlamydiae, relying on their T3S apparatus, cause K⁺ efflux, which in turn leads to formation of reactive oxygen species (ROS) and ultimately NLRP3-dependent caspase-1 activation. Epithelial cells do not typically secrete large amounts of IL-1β; instead, caspase-1 activation in cervical epithelial cells contributes to development of the chlamydial inclusion.     

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.