New decoding algorithms for Hidden Markov Models using distance measures on labellings

<p>Abstract</p> <p>Background</p> <p>Existing hidden Markov model decoding algorithms do not focus on approximately identifying the sequence feature boundaries.</p> <p>Results</p> <p>We give a set of algorithms to compute the conditional probability of all labellings "near" a reference labelling <it>λ </it>for a sequence <it>y </it>for a variety of definitions of "near". In addition, we give optimization algorithms to find the best labelling for a sequence in the robust sense of having all of its feature boundaries nearly correct. Natural problems in this domain are <it>NP</it>-hard to optimize. For membrane proteins, our algorithms find the approximate topology of such proteins with comparable success to existing programs, while being substantially more accurate in estimating the positions of transmembrane helix boundaries.</p> <p>Conclusion</p> <p>More robust HMM decoding may allow for better analysis of sequence features, in reasonable runtimes.</p

Identification of Giardia lamblia DHHC Proteins and the Role of Protein S-palmitoylation in the Encystation Process

Protein S-palmitoylation, a hydrophobic post-translational modification, is performed by protein acyltransferases that have a common DHHC Cys-rich domain (DHHC proteins), and provides a regulatory switch for protein membrane association. In this work, we analyzed the presence of DHHC proteins in the protozoa parasite Giardia lamblia and the function of the reversible S-palmitoylation of proteins during parasite differentiation into cyst. Two specific events were observed: encysting cells displayed a larger amount of palmitoylated proteins, and parasites treated with palmitoylation inhibitors produced a reduced number of mature cysts. With bioinformatics tools, we found nine DHHC proteins, potential protein acyltransferases, in the Giardia proteome. These proteins displayed a conserved structure when compared to different organisms and are distributed in different monophyletic clades. Although all Giardia DHHC proteins were found to be present in trophozoites and encysting cells, these proteins showed a different intracellular localization in trophozoites and seemed to be differently involved in the encystation process when they were overexpressed. dhhc transgenic parasites showed a different pattern of cyst wall protein expression and yielded different amounts of mature cysts when they were induced to encyst. Our findings disclosed some important issues regarding the role of DHHC proteins and palmitoylation during Giardia encystation.Fil: Merino, Maria Cecilia. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra. Universidad Nacional de Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra; ArgentinaFil: Zamponi, Nahuel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra. Universidad Nacional de Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra; ArgentinaFil: Vranych, Cecilia Verónica. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra. Universidad Nacional de Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra; ArgentinaFil: Touz, Maria Carolina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra. Universidad Nacional de Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra; ArgentinaFil: Ropolo, Andrea Silvana. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra. Universidad Nacional de Córdoba. Instituto de Investigación Médica Mercedes y Martín Ferreyra; Argentin

Measurement of the running of the QED coupling in small-angle Bhabha scattering at LEP

Using the OPAL detector at LEP, the running of the effective QED coupling alpha(t) is measured for space-like momentum transfer from the angular distribution of small-angle Bhabha scattering. In an almost ideal QED framework, with very favourable experimental conditions, we obtain: Delta alpha(-6.07GeV^2) - Delta alpha(-1.81GeV^2) = (440 pm 58 pm 43 pm 30) X 10^-5, where the first error is statistical, the second is the experimental systematic and the third is the theoretical uncertainty. This agrees with current evaluations of alpha(t).The null hypothesis that alpha remains constant within the above interval of -t is excluded with a significance above 5sigma. Similarly, our results are inconsistent at the level of 3sigma with the hypothesis that only leptonic loops contribute to the running. This is currently the most significant direct measurment where the running alpha(t) is probed differentially within the measured t range.Comment: 43 pages, 12 figures, Submitted to Euro. Phys. J.

A Measurement of Rb using a Double Tagging Method

The fraction of Z to bbbar events in hadronic Z decays has been measured by the OPAL experiment using the data collected at LEP between 1992 and 1995. The Z to bbbar decays were tagged using displaced secondary vertices, and high momentum electrons and muons. Systematic uncertainties were reduced by measuring the b-tagging efficiency using a double tagging technique. Efficiency correlations between opposite hemispheres of an event are small, and are well understood through comparisons between real and simulated data samples. A value of Rb = 0.2178 +- 0.0011 +- 0.0013 was obtained, where the first error is statistical and the second systematic. The uncertainty on Rc, the fraction of Z to ccbar events in hadronic Z decays, is not included in the errors. The dependence on Rc is Delta(Rb)/Rb = -0.056*Delta(Rc)/Rc where Delta(Rc) is the deviation of Rc from the value 0.172 predicted by the Standard Model. The result for Rb agrees with the value of 0.2155 +- 0.0003 predicted by the Standard Model.Comment: 42 pages, LaTeX, 14 eps figures included, submitted to European Physical Journal

Oxygen Reperfusion Damage in an Insect

The deleterious effects of anoxia followed by reperfusion with oxygen in higher animals including mammals are well known. A convenient and genetically well characterized small-animal model that exhibits reproducible, quantifiable oxygen reperfusion damage is currently lacking. Here we describe the dynamics of whole-organism metabolic recovery from anoxia in an insect, Drosophila melanogaster, and report that damage caused by oxygen reperfusion can be quantified in a novel but straightforward way. We monitored CO2 emission (an index of mitochondrial activity) and water vapor output (an index of neuromuscular control of the spiracles, which are valves between the outside air and the insect's tracheal system) during entry into, and recovery from, rapid-onset anoxia exposure with durations ranging from 7.5 to 120 minutes. Anoxia caused a brief peak of CO2 output followed by knock-out. Mitochondrial respiration ceased and the spiracle constrictor muscles relaxed, but then re-contracted, presumably powered by anaerobic processes. Reperfusion to sustained normoxia caused a bimodal re-activation of mitochondrial respiration, and in the case of the spiracle constrictor muscles, slow inactivation followed by re-activation. After long anoxia durations, both the bimodality of mitochondrial reactivation and the recovery of spiracular control were impaired. Repeated reperfusion followed by episodes of anoxia depressed mitochondrial respiratory flux rates and damaged the integrity of the spiracular control system in a dose-dependent fashion. This is the first time that physiological evidence of oxygen reperfusion damage has been described in an insect or any invertebrate. We suggest that some of the traditional approaches of insect respiratory biology, such as quantifying respiratory water loss, may facilitate using D. melanogaster as a convenient, well-characterized experimental model for studying the underlying biology and mechanisms of ischemia and reperfusion damage and its possible mitigation

Measurement of the B+ and B-0 lifetimes and search for CP(T) violation using reconstructed secondary vertices

The lifetimes of the B+ and B-0 mesons, and their ratio, have been measured in the OPAL experiment using 2.4 million hadronic Z(0) decays recorded at LEP. Z(0) --> b (b) over bar decays were tagged using displaced secondary vertices and high momentum electrons and muons. The lifetimes were then measured using well-reconstructed charged and neutral secondary vertices selected in this tagged data sample. The results aretau(B+) = 1.643 +/- 0.037 +/- 0.025 pstau(Bo) = 1.523 +/- 0.057 +/- 0.053 pstau(B+)/tau(Bo) = 1.079 +/- 0.064 +/- 0.041,where in each case the first error is statistical and the second systematic.A larger data sample of 3.1 million hadronic Z(o) decays has been used to search for CP and CPT violating effects by comparison of inclusive b and (b) over bar hadron decays, No evidence fur such effects is seen. The CP violation parameter Re(epsilon(B)) is measured to be Re(epsilon(B)) = 0.001 +/- 0.014 +/- 0.003and the fractional difference between b and (b) over bar hadron lifetimes is measured to(Delta tau/tau)(b) = tau(b hadron) - tau((b) over bar hadron)/tau(average) = -0.001 +/- 0.012 +/- 0.008

Bacterial Transmembrane Proteins that Lack N-Terminal Signal Sequences

Tail-anchored membrane proteins (TAMPs), a class of proteins characterized by their lack of N-terminal signal sequence and Sec-independent membrane targeting, play critical roles in apoptosis, vesicle trafficking and other vital processes in eukaryotic organisms. Until recently, this class of membrane proteins has been unknown in bacteria. Here we present the results of bioinformatic analysis revealing proteins that are superficially similar to eukaryotic TAMPs in the bacterium Streptomyces coelicolor. We demonstrate that at least four of these proteins are bona fide membrane-spanning proteins capable of targeting to the membrane in the absence of their N-terminus and the C-terminal membrane-spanning domain is sufficient for membrane targeting. Several of these proteins, including a serine/threonine kinase and the SecE component of the Sec translocon, are widely conserved in bacteria

Efficiency of Purine Utilization by Helicobacter pylori: Roles for Adenosine Deaminase and a NupC Homolog

The ability to synthesize and salvage purines is crucial for colonization by a variety of human bacterial pathogens. Helicobacter pylori colonizes the gastric epithelium of humans, yet its specific purine requirements are poorly understood, and the transport mechanisms underlying purine uptake remain unknown. Using a fully defined synthetic growth medium, we determined that H. pylori 26695 possesses a complete salvage pathway that allows for growth on any biological purine nucleobase or nucleoside with the exception of xanthosine. Doubling times in this medium varied between 7 and 14 hours depending on the purine source, with hypoxanthine, inosine and adenosine representing the purines utilized most efficiently for growth. The ability to grow on adenine or adenosine was studied using enzyme assays, revealing deamination of adenosine but not adenine by H. pylori 26695 cell lysates. Using mutant analysis we show that a strain lacking the gene encoding a NupC homolog (HP1180) was growth-retarded in a defined medium supplemented with certain purines. This strain was attenuated for uptake of radiolabeled adenosine, guanosine, and inosine, showing a role for this transporter in uptake of purine nucleosides. Deletion of the GMP biosynthesis gene guaA had no discernible effect on mouse stomach colonization, in contrast to findings in numerous bacterial pathogens. In this study we define a more comprehensive model for purine acquisition and salvage in H. pylori that includes purine uptake by a NupC homolog and catabolism of adenosine via adenosine deaminase

Positively Selected Codons in Immune-Exposed Loops of the Vaccine Candidate OMP-P1 of Haemophilus influenzae

The high levels of variation in surface epitopes can be considered as an evolutionary hallmark of immune selection. New computational tools enable analysis of this variation by identifying codons that exhibit high rates of amino acid changes relative to the synonymous substitution rate. In the outer membrane protein P1 of Haemophilus influenzae, a vaccine candidate for nontypeable strains, we identified four codons with this attribute in domains that did not correspond to known or assumed B- and T-cell epitopes of OMP-P1. These codons flank hypervariable domains and do not appear to be false positives as judged from parsimony and maximum likelihood analyses. Some closely spaced positively selected codons have been previously considered part of a transmembrane domain, which would render this region unsuited for inclusion in a vaccine. Secondary structure analysis, three-dimensional structural database searches, and homology modeling using FadL of E. coli as a structural homologue, however, revealed that all positively selected codons are located in or near extracellular looping domains. The spacing and level of diversity of these positively selected and exposed codons in OMP-P1 suggest that vaccine targets based on these and conserved flanking residues may provide broad coverage in H. influenzae