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]     

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.     

Arachidonic Acid Stimulates Cell Adhesion through a Novel p38 MAPK-RhoA Signaling Pathway That Involves Heat Shock Protein 27

Rho GTPases are critical components of cellular signal transduction pathways. Both hyperactivity and overexpression of these proteins have been observed in human cancers and have been implicated as important factors in metastasis. We previously showed that dietary n-6 fatty acids increase cancer cell adhesion to extracellular matrix proteins, such as type IV collagen. Here we report that in MDA-MB-435 human melanoma cells, arachidonic acid activates RhoA, and inhibition of RhoA signaling with either C3 exoenzyme or dominant negative Rho blocked arachidonic acid-induced cell adhesion. Inhibition of the Rho kinase (ROCK) with either small molecule inhibitors or ROCK II-specific small interfering RNA (siRNA) blocked the fatty acid-induced adhesion. However, unlike other systems, inhibition of ROCK did not block the activation of p38 mitogen-activated protein kinase (MAPK); instead, Rho activation depended on p38 MAPK activity and the presence of heat shock protein 27 (HSP27), which is phosphorylated downstream of p38 after arachidonic acid treatment. HSP27 associated with p115RhoGEF in fatty acid-treated cells, and this association was blocked when p38 was inhibited. Furthermore, siRNA knockdown of HSP27 blocked the fatty acid-stimulated Rho activity. Expression of dominant negative p115-RhoGEF or p115RhoGEF-specific siRNA inhibited both RhoA activation and adhesion on type IV collagen, whereas a constitutively active p115RhoGEF restored the arachidonic acid stimulation in cells in which the p38 MAPK had been inhibited. These data suggest that n-6 dietary fatty acids stimulate a set of interactions that regulates cell adhesion through RhoA and ROCK II via a p38 MAPK-dependent association of HSP27 and p115RhoGEF.The ability of tumor cells to metastasize to secondary sites is a hallmark of neoplastic disease. Unfortunately, this propensity to spread is the primary cause of morbidity and death in cancer patients (1). Metastasis is clearly a highly regulated, multistep process that occurs in a spatiotemporal manner (2–4). To escape the restrictive compartment boundaries characteristic of adult tissue, separate intravasation and extravasation steps requiring alterations in co-adhesion, adhesion, invasion, and migration must occur. Execution of these biological processes, involving multiple proteins and cellular organelles, require highly coordinated cell signaling mechanisms.The Rho family of small GTPases regulates many facets of cytoskeletal rearrangements that facilitate cell attachment and migration (5–7). Rho GTPases act as molecular switches by changing from an inactive GDP-bound conformation to an active GTP-bound conformation, thereby regulating a signaling pathway. These proteins are directly regulated by Rho guanine nucleotide exchange factors (GEFs),² Rho GTPase activating proteins, and Rho GDP-dissociation inhibitors (8–12). RhoGEFs bind to the GTPase to catalyze the dissociation of GDP, allowing the binding of GTP and thereby promoting Rho activation (8). The RGS (regulators of G protein signaling) domain-containing RhoGEFs are a recently described family of GEFs. Currently, there are three members of this family, PDZ-RhoGEF, LARG, and p115RhoGEF (13–15), in which the RGS domains function as a heterotrimeric GTPase-activating domain (13, 15, 16). The RGS family of RhoGEFs has been shown to regulate Rho during several processes including cytoskeletal rearrangements, cell adhesion, and cancer progression (17–21).There is significant interplay between the activity of small GTPases and signaling derived from fatty acid metabolism (22–28). Linoleic acid, which is metabolized to arachidonic acid, is an n-6 polyunsaturated fatty acid that is present at high levels in most western diets (29). In animal models, diets high in n-6 polyunsaturated fatty acids have been shown to enhance tumor progression and metastasis (30, 31). Additionally, arachidonic acid is stored in cell membranes and is made available by phospholipases under conditions of increased inflammatory response (32). Arachidonic acid is further metabolized by cyclooxygenases (COX), lipoxygenases (LOX), and cytochrome P450 monooxygenases to yield bioactive products that have myriad effects on cells, and altered metabolism of arachidonic acid by COX, LOX, and P450 has been implicated in cancer progression (31, 33–36).We have studied mechanisms of cell adhesion using the MDA-MB-435 cells as a model of a highly metastatic human cancer cell line (37). These cells have been extensively studied for their ability to recapitulate the metastatic cascade in vivo and in vitro, although recent work indicates that the cells currently in use are most likely a human melanoma line (38). We initially observed that arachidonic acid (AA) enhanced adhesion of MDA-MB-435 cells to type IV collagen through specific integrin-mediated pathways (37). Exogenous AA led to the activation of mitogen-activated protein kinase (MAPK)-activated protein kinase 2 and the phosphorylation of heat shock protein 27 (HSP27) via a p38 MAPK-dependent process (39). Inhibition of p38 MAPK activation blocked cell adhesion as did function-blocking antibodies specific for subunits of the collagen receptor (40). More recently, we identified the key metabolite of AA (15-(S)- hydroxyeicosatetraenoic acid) and the upstream kinases (TAK1 and MKK6) that are responsible for activation of p38 MAPK in this system (41).In this study we investigated the role of Rho activation in the MDA-MB-435 cells after exposure to arachidonic acid. Several aspects of the regulation of Rho signaling in these cells provide insights into the cross-talk between important signaling pathways.     

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]     

Deletion of the Chloride Transporter Slc26a7 Causes Distal Renal Tubular Acidosis and Impairs Gastric Acid Secretion

SLC26A7 (human)/Slc26a7 (mouse) is a recently identified chloride-base exchanger and/or chloride transporter that is expressed on the basolateral membrane of acid-secreting cells in the renal outer medullary collecting duct (OMCD) and in gastric parietal cells. Here, we show that mice with genetic deletion of Slc26a7 expression develop distal renal tubular acidosis, as manifested by metabolic acidosis and alkaline urine pH. In the kidney, basolateral Cl⁻/HCO₃⁻ exchange activity in acid-secreting intercalated cells in the OMCD was significantly decreased in hypertonic medium (a normal milieu for the medulla) but was reduced only mildly in isotonic medium. Changing from a hypertonic to isotonic medium (relative hypotonicity) decreased the membrane abundance of Slc26a7 in kidney cells in vivo and in vitro. In the stomach, stimulated acid secretion was significantly impaired in isolated gastric mucosa and in the intact organ. We propose that SLC26A7 dysfunction should be investigated as a potential cause of unexplained distal renal tubular acidosis or decreased gastric acid secretion in humans.The collecting duct segment of the distal kidney nephron plays a major role in systemic acid base homeostasis by acid secretion and bicarbonate absorption. The acid secretion occurs via H⁺-ATPase and H-K-ATPase into the lumen and bicarbonate is absorbed via basolateral Cl⁻/HCO₃⁻ exchangers (1–4). The tubules, which are located within the outer medullary region of the kidney collecting duct (OMCD),² have the highest rate of acid secretion among the distal tubule segments and are therefore essential to the maintenance of acid base balance (2).The gastric parietal cell is the site of generation of acid and bicarbonate through the action of cytosolic carbonic anhydrase II (5, 6). The intracellular acid is secreted into the lumen via gastric H-K-ATPase, which works in conjunction with a chloride channel and a K⁺ recycling pathway (7–10). The intracellular bicarbonate is transported to the blood via basolateral Cl⁻/HCO₃⁻ exchangers (11–14).SLC26 (human)/Slc26 (mouse) isoforms are members of a conserved family of anion transporters that display tissue-specific patterns of expression in epithelial cells (15–24). Several SLC26 members can function as chloride/bicarbonate exchangers. These include SLC26A3 (DRA), SLC26A4 (pendrin), SLC26A6 (PAT1 or CFEX), SLC26A7, and SLC26A9 (25–31). SLC26A7 and SLC26A9 can also function as chloride channels (32–34).SLC26A7/Slc26a7 is predominantly expressed in the kidney and stomach (28, 29). In the kidney, Slc26a7 co-localizes with AE1, a well-known Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of (acid-secreting) A-intercalated cells in OMCD cells (29, 35, 36) (supplemental Fig. 1). In the stomach, Slc26a7 co-localizes with AE2, a major Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of acid secreting parietal cells (28). To address the physiological function of Slc26a7 in the intact mouse, we have generated Slc26a7 ko mice. We report here that Slc26a7 ko mice exhibit distal renal tubular acidosis and impaired gastric acidification in the absence of morphological abnormalities in kidney or stomach.     

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.     

Protein Analysis of Purified Respiratory Syncytial Virus Particles Reveals an Important Role for Heat Shock Protein 90 in Virus Particle Assembly

In this study, we used imaging and proteomics to identify the presence of virus-associated cellular proteins that may play a role in respiratory syncytial virus (RSV) maturation. Fluorescence microscopy of virus-infected cells revealed the presence of virus-induced cytoplasmic inclusion bodies and mature virus particles, the latter appearing as virus filaments. In situ electron tomography suggested that the virus filaments were complex structures that were able to package multiple copies of the virus genome. The virus particles were purified, and the protein content was analyzed by one-dimensional nano-LC MS/MS. In addition to all the major virus structural proteins, 25 cellular proteins were also detected, including proteins associated with the cortical actin network, energy pathways, and heat shock proteins (HSP70, HSC70, and HSP90). Representative actin-associated proteins, HSC70, and HSP90 were selected for further biological validation. The presence of β-actin, filamin-1, cofilin-1, HSC70, and HSP90 in the virus preparation was confirmed by immunoblotting using relevant antibodies. Immunofluorescence microscopy of infected cells stained with antibodies against relevant virus and cellular proteins confirmed the presence of these cellular proteins in the virus filaments and inclusion bodies. The relevance of HSP90 to virus infection was examined using the specific inhibitors 17-N-Allylamino-17-demethoxygeldanamycin. Although virus protein expression was largely unaffected by these drugs, we noted that the formation of virus particles was inhibited, and virus transmission was impaired, suggesting an important role for HSP90 in virus maturation. This study highlights the utility of proteomics in facilitating both our understanding of the role that cellular proteins play during RSV maturation and, by extrapolation, the identification of new potential targets for antiviral therapy.Respiratory syncytial virus (RSV)¹ belongs to the paramyxovirus group of viruses, and it is the most important respiratory virus causing lower respiratory tract infection in young children and neonates. The mature RSV particle comprises a ribonucleoparticle (RNP) core formed by the interaction between the viral genomic RNA (vRNA), the nucleocapsid (N) protein (42 kDa), the phospho (P) protein (35 kDa), and the large (L) protein (250 kDa). The RNP core is visualized by electron microscopy as a strand of repeating N protein subunits that form a herringbone-like structure of ∼10–20 nm in diameter (1). Although the minimal functional polymerase activity requires an association between the N, P, and L proteins and the virus genome vRNA (2–4), additional viral proteins called the M2-1 protein (22 kDa), M2-2 protein, and M protein (28 kDa) regulate the activity of the polymerase (5–8). The virus is surrounded by a lipid envelope that is formed from the host cell during the budding process in which the three virus membrane proteins are inserted. The G protein (90 kDa) mediates attachment of the virus to the cell during virus entry (9), and the fusion (F) protein (10) mediates the fusion of the virus and host cell membranes during virus entry, whereas the role of the SH protein is currently unknown. In addition, two non-structural proteins called NS1 and NS2, which are thought not to be present in the virus particle but play a role in countering the host innate immune response (11), are expressed.During virus infection two distinct virus structures are formed, virus filaments and inclusion bodies. The virus filaments are membrane-bound structures that are ∼150–200 nm thick and can be up to 6 μm in length (1, 12–16); they form at the sites of virus assembly and are the progeny viruses. The inclusion bodies form in the cytoplasm and can be several μm in diameter, consisting of accumulations of RNP cores (17–19). Inclusion bodies are found in all RSV-infected tissue culture cells, and they have also been observed in biopsy material isolated from RSV-infected patients (20) suggesting a clinical relevance. Although the cellular processes that lead to assembly of the mature virus filaments are still poorly understood, the involvement of lipid raft microdomains and the cortical cytoskeleton network appear to play an important role in this process (16, 21–25). For example, rhoA kinase is a raft-associated signaling molecule that is involved in regulating actin structure (26), and it has been implicated in virus filament formation (27, 28). Virus filament formation also requires phosphoinositide 3-kinase (PI3K) activity (25, 29, 30); PI3K is a raft-associated kinase activated by rhoA kinase (31). The identification of cellular proteins that interact with the virus particles should further facilitate the identification of the cellular pathways that are involved in RSV maturation. In this study, we examined virus-host cell interactions during RSV assembly using a combination of advanced imaging techniques and analyzed the protein content of purified virus particles by proteomics technology. Our analysis provides evidence that cellular proteins that regulate actin structures in the cell may also play an important role in formation of infectious RSV particles, and that the HSP90 protein plays an important role in the virus assembly process.