GlyGly-CTERM and Rhombosortase: A C-Terminal Protein Processing Signal in a Many-to-One Pairing with a Rhomboid Family Intramembrane Serine Protease

The rhomboid family of serine proteases occurs in all domains of life. Its members contain at least six hydrophobic membrane-spanning helices, with an active site serine located deep within the hydrophobic interior of the plasma membrane. The model member GlpG from Escherichia coli is heavily studied through engineered mutant forms, varied model substrates, and multiple X-ray crystal studies, yet its relationship to endogenous substrates is not well understood. Here we describe an apparent membrane anchoring C-terminal homology domain that appears in numerous genera including Shewanella, Vibrio, Acinetobacter, and Ralstonia, but excluding Escherichia and Haemophilus. Individual genomes encode up to thirteen members, usually homologous to each other only in this C-terminal region. The domain's tripartite architecture consists of motif, transmembrane helix, and cluster of basic residues at the protein C-terminus, as also seen with the LPXTG recognition sequence for sortase A and the PEP-CTERM recognition sequence for exosortase. Partial Phylogenetic Profiling identifies a distinctive rhomboid-like protease subfamily almost perfectly co-distributed with this recognition sequence. This protease subfamily and its putative target domain are hereby renamed rhombosortase and GlyGly-CTERM, respectively. The protease and target are encoded by consecutive genes in most genomes with just a single target, but far apart otherwise. The signature motif of the Rhombo-CTERM domain, often SGGS, only partially resembles known cleavage sites of rhomboid protease family model substrates. Some protein families that have several members with C-terminal GlyGly-CTERM domains also have additional members with LPXTG or PEP-CTERM domains instead, suggesting there may be common themes to the post-translational processing of these proteins by three different membrane protein superfamilies

Bioinformatic evidence for a widely distributed, ribosomally produced electron carrier precursor, its maturation proteins, and its nicotinoprotein redox partners

<p>Abstract</p> <p>Background</p> <p>Enzymes in the radical SAM (rSAM) domain family serve in a wide variety of biological processes, including RNA modification, enzyme activation, bacteriocin core peptide maturation, and cofactor biosynthesis. Evolutionary pressures and relationships to other cellular constituents impose recognizable grammars on each class of rSAM-containing system, shaping patterns in results obtained through various comparative genomics analyses.</p> <p>Results</p> <p>An uncharacterized gene cluster found in many Actinobacteria and sporadically in Firmicutes, Chloroflexi, Deltaproteobacteria, and one Archaeal plasmid contains a PqqE-like rSAM protein family that includes Rv0693 from <it>Mycobacterium tuberculosis</it>. Members occur clustered with a strikingly well-conserved small polypeptide we designate "mycofactocin," similar in size to bacteriocins and PqqA, precursor of pyrroloquinoline quinone (PQQ). Partial Phylogenetic Profiling (PPP) based on the distribution of these markers identifies the mycofactocin cluster, but also a second tier of high-scoring proteins. This tier, strikingly, is filled with up to thirty-one members per genome from three variant subfamilies that occur, one each, in three unrelated classes of nicotinoproteins. The pattern suggests these variant enzymes require not only NAD(P), but also the novel gene cluster. Further study was conducted using SIMBAL, a PPP-like tool, to search these nicotinoproteins for subsequences best correlated across multiple genomes to the presence of mycofactocin. For both the short chain dehydrogenase/reductase (SDR) and iron-containing dehydrogenase families, aligning SIMBAL's top-scoring sequences to homologous solved crystal structures shows signals centered over NAD(P)-binding sites rather than over substrate-binding or active site residues. Previous studies on some of these proteins have revealed a non-exchangeable NAD cofactor, such that enzymatic activity <it>in vitro </it>requires an artificial electron acceptor such as N,N-dimethyl-4-nitrosoaniline (NDMA) for the enzyme to cycle.</p> <p>Conclusions</p> <p>Taken together, these findings suggest that the mycofactocin precursor is modified by the Rv0693 family rSAM protein and other enzymes in its cluster. It becomes an electron carrier molecule that serves <it>in vivo </it>as NDMA and other artificial electron acceptors do <it>in vitro</it>. Subclasses from three different nicotinoprotein families show "only-if" relationships to mycofactocin because they require its presence. This framework suggests a segregated redox pool in which mycofactocin mediates communication among enzymes with non-exchangeable cofactors.</p

The B6 database: a tool for the description and classification of vitamin B6-dependent enzymatic activities and of the corresponding protein families

BACKGROUND: Enzymes that depend on vitamin B6 (and in particular on its metabolically active form, pyridoxal 5'-phosphate, PLP) are of great relevance to biology and medicine, as they catalyze a wide variety of biochemical reactions mainly involving amino acid substrates. Although PLP-dependent enzymes belong to a small number of independent evolutionary lineages, they encompass more than 160 distinct catalytic functions, thus representing a striking example of divergent evolution. The importance and remarkable versatility of these enzymes, as well as the difficulties in their functional classification, create a need for an integrated source of information about them. DESCRIPTION: The B6 database http://bioinformatics.unipr.it/B6db contains documented B6-dependent activities and the relevant protein families, defined as monophyletic groups of sequences possessing the same enzymatic function. One or more families were associated to each of 121 PLP-dependent activities with known sequences. Hidden Markov models (HMMs) were built from family alignments and incorporated in the database. These HMMs can be used for the functional classification of PLP-dependent enzymes in genomic sets of predicted protein sequences. An example of such analyses (a census of human genes coding for PLP-dependent enzymes) is provided here, whereas many more are accessible through the database itself. CONCLUSION: The B6 database is a curated repository of biochemical and molecular information about an important group of enzymes. This information is logically organized and available for computational analyses, providing a key resource for the identification, classification and comparative analysis of B6-dependent enzymes

FastBLAST: Homology Relationships for Millions of Proteins

BackgroundAll-versus-all BLAST, which searches for homologous pairs of sequences in a database of proteins, is used to identify potential orthologs, to find new protein families, and to provide rapid access to these homology relationships. As DNA sequencing accelerates and data sets grow, all-versus-all BLAST has become computationally demanding.Methodology/principal findingsWe present FastBLAST, a heuristic replacement for all-versus-all BLAST that relies on alignments of proteins to known families, obtained from tools such as PSI-BLAST and HMMer. FastBLAST avoids most of the work of all-versus-all BLAST by taking advantage of these alignments and by clustering similar sequences. FastBLAST runs in two stages: the first stage identifies additional families and aligns them, and the second stage quickly identifies the homologs of a query sequence, based on the alignments of the families, before generating pairwise alignments. On 6.53 million proteins from the non-redundant Genbank database ("NR"), FastBLAST identifies new families 25 times faster than all-versus-all BLAST. Once the first stage is completed, FastBLAST identifies homologs for the average query in less than 5 seconds (8.6 times faster than BLAST) and gives nearly identical results. For hits above 70 bits, FastBLAST identifies 98% of the top 3,250 hits per query.Conclusions/significanceFastBLAST enables research groups that do not have supercomputers to analyze large protein sequence data sets. FastBLAST is open source software and is available at http://microbesonline.org/fastblast

Exopolysaccharide-associated protein sorting in environmental organisms: the PEP-CTERM/EpsH system. Application of a novel phylogenetic profiling heuristic

BACKGROUND: Protein translocation to the proper cellular destination may be guided by various classes of sorting signals recognizable in the primary sequence. Detection in some genomes, but not others, may reveal sorting system components by comparison of the phylogenetic profile of the class of sorting signal to that of various protein families. RESULTS: We describe a short C-terminal homology domain, sporadically distributed in bacteria, with several key characteristics of protein sorting signals. The domain includes a near-invariant motif Pro-Glu-Pro (PEP). This possible recognition or processing site is followed by a predicted transmembrane helix and a cluster rich in basic amino acids. We designate this domain PEP-CTERM. It tends to occur multiple times in a genome if it occurs at all, with a median count of eight instances; Verrucomicrobium spinosum has sixty-five. PEP-CTERM-containing proteins generally contain an N-terminal signal peptide and exhibit high diversity and little homology to known proteins. All bacteria with PEP-CTERM have both an outer membrane and exopolysaccharide (EPS) production genes. By a simple heuristic for screening phylogenetic profiles in the absence of pre-formed protein families, we discovered that a homolog of the membrane protein EpsH (exopolysaccharide locus protein H) occurs in a species when PEP-CTERM domains are found. The EpsH family contains invariant residues consistent with a transpeptidase function. Most PEP-CTERM proteins are encoded by single-gene operons preceded by large intergenic regions. In the Proteobacteria, most of these upstream regions share a DNA sequence, a probable cis-regulatory site that contains a sigma-54 binding motif. The phylogenetic profile for this DNA sequence exactly matches that of three proteins: a sigma-54-interacting response regulator (PrsR), a transmembrane histidine kinase (PrsK), and a TPR protein (PrsT). CONCLUSION: These findings are consistent with the hypothesis that PEP-CTERM and EpsH form a protein export sorting system, analogous to the LPXTG/sortase system of Gram-positive bacteria, and correlated to EPS expression. It occurs preferentially in bacteria from sediments, soils, and biofilms. The novel method that led to these findings, partial phylogenetic profiling, requires neither global sequence clustering nor arbitrary similarity cutoffs and appears to be a rapid, effective alternative to other profiling methods

HAD hydrolase function unveiled by substrate screening: enzymatic characterization of Arabidopsis thaliana subclass I phosphosugar phosphatise AtSgpp

[EN] This work presents the isolation and the biochemical characterization of the Arabidopsis thaliana gene AtSgpp. This gene shows homology with the Arabidopsis low molecular weight phosphatases AtGpp1 and AtGpp2 and the yeast counterpart GPP1 and GPP2, which have a high specificity for dl-glycerol-3-phosphate. In addition, it exhibits homology with DOG1 and DOG2 that dephosphorylate 2-deoxy-d-glucose-6-phosphate. Using a comparative genomic approach, we identified the AtSgpp gene as a conceptual translated haloacid dehalogenase-like hydrolase HAD protein. AtSgpp (locus tag At2g38740), encodes a protein with a predicted Mw of 26.7 kDa and a pI of 4.6. Its sequence motifs and expected structure revealed that AtSgpp belongs to the HAD hydrolases subfamily I, with the C1-type cap domain. In the presence of Mg2+ ions, the enzyme has a phosphatase activity over a wide range of phosphosugars substrates (pH optima at 7.0 and K (m) in the range of 3.6-7.7 mM). AtSgpp promiscuity is preferentially detectable on d-ribose-5-phosphate, 2-deoxy-d-ribose-5-phosphate, 2-deoxy-d-glucose-6-phosphate, d-mannose-6-phosphate, d-fructose-1-phosphate, d-glucose-6-phosphate, dl-glycerol-3-phosphate, and d-fructose-6-phosphate, as substrates. AtSgpp is ubiquitously expressed throughout development in most plant organs, mainly in sepal and guard cell. Interestingly, expression is affected by abiotic and biotic stresses, being the greatest under Pi starvation and cyclopentenone oxylipins induction. Based on both, substrate lax specificity and gene expression, the physiological function of AtSgpp in housekeeping detoxification, modulation of sugar-phosphate balance and Pi homeostasis, is provisionally assigned.We acknowledge Professors Montserrat Pages (CSIC Barcelona, Spain), Thomas Kupke (University of Heidelberg, Germany) and Manuel Hernandez (University Polytechnic of Valencia, Spain) for their warm support. We also thank the advice and provision of plasmid pSBETa by Dr. Florence Vignols and Yves Meyer (University of Perpignan, France); the computer software helps by Ramon Nogales-Rangel and Alexis Gonzalez-Policarpo; Eugenio Grau-Ferrando for kind advice and help for sequencing. This work was funded by the 10 month research contract MEC-FEDER to J.A.C.-M., 10 month research contract JAE-DOC to I.M.-S. and by the research project BIO2006-10138 from the MEC-FEDER of Spain to F.A.C.-M. In memoriam of Dr. Mari Cruz Cutanda-Perez.Caparrós Martín, JA.; Mccarthy Suarez, I.; Culiañez Macia, FA. (2013). HAD hydrolase function unveiled by substrate screening: enzymatic characterization of Arabidopsis thaliana subclass I phosphosugar phosphatise AtSgpp. Planta. 237(4):943-954. https://doi.org/10.1007/s00425-012-1809-5S9439542374Allen KN, Dunaway-Mariano D (2004) Phosphoryl group transfer: evolution of a catalytic scaffold. Trends Biochem Sci 29:495–503Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403–410Ames BN (1966) Assay of inorganic phosphate, total phosphate, and phosphatases. Methods Enzymol 8:115–118Böhmer M, Schroeder JI (2011) Quantitative transcriptomic analysis of abscisic acid-induced and reactive oxygen species-dependent expression changes and proteomic profiling in Arabidopsis suspension cells. Plant J 67:105–118Bradford MM (1976) A rapid and sensitive method for the quantization of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 72:248–254Burroughs AM, Allen KN, Dunaway-Mariano D, Aravind L (2006) Evolutionary genomics of the HAD superfamily: understanding the structural adaptations and catalytic diversity in a superfamily of phosphoesterases and allied enzymes. J Mol Biol 361:1003–1034Caparrós-Martín JA, Reiland S, Köchert K, Cutanda MC, Culiáñez-Macia FA (2007) Arabidopsis thaliana AtGpp 1 and AtGpp2: two novel low molecular weight phosphatases involved in plant glycerol metabolism. Plant Mol Biol 63:505–517Collet JF, Stroobant V, Pirard M, Delpierre G, Van Schaftingen E (1998) A new class of phosphotransferases phosphorylated on an aspartate residue in an amino-terminal DXDX(T/V) motif. J Biol Chem 273:14107–14112Corpet F, Servantm F, Gouzy J, Kahn D (2000) ProDom and ProDom-CG: tools for protein domain analysis and whole genome comparisons. Nucleic Acids Res 28:267–269Cutanda MC (2003) Effect of altering levels of hexoses phosphate in carbohydrate metabolism and glucose signalling in yeast and plants. PhD thesis, Polytechnic University of Valencia, Valencia, SpainHiggins D, Thompson J, Gibson T, Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673–4680Koonin EV, Tatusov RL (1994) Computer analysis of bacterial haloacid dehalogenases defines a large superfamily of hydrolases with diverse specificity. Application of an iterative approach to database search. J Mol Biol 244:125–132Kupke T, Caparrós-Martín JA, Malquichagua Salazar KJ, Culiàñez-Macià FA (2009) Biochemical and physiological characterization of Arabidopsis thaliana AtCoAse: a Nudix CoA hydrolyzing protein that improves plant development. Physiol Plant 135:365–378Kuznetsova E, Proudfoot M, Sanders SA, Reinking J, Savchenko A, Arrowsmith CH, Edwards AM, Yakunin AF (2005) Enzyme genomics: application of general enzymatic screens to discover new enzymes. FEMS Microbiol Rev 29:263–279Kuznetsova E, Proudfoo M, Gonzalez CF, Brown G, Omelchenko MV, Borozan I, Carmel L, Wolf YI, Mori H, Savchenko AV, Arrowsmith CH, Koonin EV, Edwards AM, Yakunin AF (2006) Genome-wide analysis of substrate specificities of the Escherichia coli haloacid dehalogenase-like phosphatase family. J Biol Chem 281:36149–36161Lahiri SD, Zhang G, Dai J, Dunaway-Mariano D, Allen KN (2004) Analysis of the substrate specificity loop of the HAD superfamily cap domain. Biochemistry 43:2812–2820Lahiri SD, Zhang G, Dunaway-Mariano D, Allen KN (2006) Diversification of function in the haloacid dehalogenase enzyme superfamily: the role of the cap domain in hydrolytic phosphorus—carbon bond cleavage. Bioorganic Chem 34:394–409Lambert C, Leonard N, De Bolle X, Depiereux E (2002) ESyPred3D: prediction of proteins 3D structures. Bioinformatics 18:1250–1256Lu Z, Dunaway-Mariano D, Allen KN (2005) HAD superfamily phosphotransferase substrate diversification: structure and function analysis of HAD subclass IIB sugar phosphatase BT4131. Biochemistry 44:8684–8696Lu Z, Dunaway-Mariano D, Allen KN (2008) The catalytic scaffold of the haloalkanoic acid dehalogenase enzyme superfamily acts as a mold for the trigonal bipyramidal transition state. Proc Natl Acad Sci USA 105:5687–5692Maniatis T, Fritsch EF, Sambrook J (1982) Molecular cloning: a laboratory manual. Cold Spring Harbor Laboratory Press, Cold Spring HarborMorais MC, Zhang W, Baker AS, Zhang G, Dunaway-Mariano D, Allen KN (2000) The crystal structure of Bacillus cereus phosphonoacetaldehyde hydrolase: insight into catalysis of phosphorus bond cleavage and catalytic diversification within the HAD enzyme superfamily. Biochemistry 39:10385–10396Morais MC, Zhang G, Zhang W, Olsen DB, Dunaway-Mariano D, Allen KN (2004) X-ray crystallographic and site-directed mutagenesis analysis of the mechanism of Schiff-base formation in phosphonoacetaldehyde hydrolase catalysis. J Biol Chem 279:9353–9361Mueller WS, Hilbert B, Dueckershoff K, Roitsch T, Krischke M, Mueller MJ, Berger S (2008) General detoxification and stress responses are mediated by oxidized lipids through TGA transcription factors in Arabidopsis. Plant Cell 20:768–785Murashige T, Skoog F (1962) A revised medium for rapid growth and bioassays with tobacco cultures. Physiol Plant 15:473–497Norbeck J, Pahlman AK, Akhtar N, Blomberg A, Adler L (1996) Purification and characterization of two isoenzymes of dl-glycerol-3-phosphatase from Saccharomyces cerevisiae. Identification of the corresponding GPP1 and GPP2 genes and evidence for osmotic regulation of Gpp 2p expression by the osmosensing mitogen-activated protein kinase signal transduction pathway. J Biol Chem 271:13875–13881Rández-Gil F, Blasco A, Prieto JA, Sanz P (1995) DOGR1 and DOGR2: two genes from Saccharomyces cerevisiae that confer 2-deoxyglucose resistance when overexpressed. Yeast 11:1233–1240Rao KN, Kumaran D, Seetharaman J, Bonanno JB, Burley SK, Swaminathan S (2006) Crystal structure of trehalose-6-phosphate phosphatase-related protein: biochemical and biological implications. Protein Sci 15:1735–1744Schagger H, von Jagow G (1987) Tricine-sodium dodecyl sulfatepolyacrylamide gel electrophoresis for the separation of proteins in the range from 1 to 100 kDa. Anal Biochem 166:368–379Schenk PM, Baumann S, Mattes R, Steinbiss HH (1995) Improved high-level expression system for eukaryotic genes in Escherichia coli using T7 RNA polymerase and rare ArgtRNAs. Biotechniques 19:196–200Selengut JD (2001) MDP-1 is a new and distinct member of the haloacid dehalogenase family of aspartate-dependent phosphohydrolases. Biochemistry 40:12704–12711Selengut JD, Levine RL (2000) MDP-1: a novel eukaryotic magnesium-dependent phosphatase. Biochemistry 39:8315–8324Shin DH, Roberts A, Jancarik J, Yocota H, Kim R, Wemmer DE, Kim S-H (2003) Crystal structure of a phosphatase with a unique substrate binding domain from Thermotoga maritime. Protein Sci 12:1464–1472Sussman I, Avron M (1981) Characterization and partial puri-fication of dl-glycerol-1-phosphatase from Dunaliella salina. Biochim Biophys Acta 661:199–204The Arabidopsis Genome Initiative (2000) Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408:796–815Tremblay LW, Dunaway-Mariano D, Allen KN (2006) Structure and activity analyses of Escherichia coli K-12 NagD provide insight into the evolution of biochemical function in the haloalkanoic acid dehalogenase superfamily. Biochemistry 45:1183–1193Vicient CM, Delseny M (1999) Isolation of total RNA from Arabidopsis thaliana seeds. Anal Biochem 268:412–413Wang W, Cho HS, Kim R, Jancarik J, Yokota H, Nguyen HH, Grigoriev IV, Wemmer DE, Kim S-H (2002) Structural characterization of the reaction pathway in phosphoserine phosphatase: crystallographic “snapshots” of intermediate states. J Mol Biol 319:421–431Zhang G, Mazurkie AS, Dunaway-Mariano D, Allen KN (2002) Kinetic evidence for a substrate-induced fit in phosphonoacetaldehyde hydrolase catalysis. Biochemistry 41:13370–13377Zhang G, Morais MC, Dai J, Zhang W, Dunaway-Mariano D, Allen KN (2004) Investigation of metal Ion binding in phosphonoacetaldehyde hydrolase identifies sequence markers for metal-activated enzymes of the HAD enzyme superfamily. Biochemistry 43:4990–4997Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W (2004) GENEVESTIGATOR: Arabidopsis microarray database and analysis toolbox. Plant Physiol 136:2621–263

Predicted Relative Metabolomic Turnover (PRMT): determining metabolic turnover from a coastal marine metagenomic dataset

We present an approach in which the semantics of an XML language is defined by means of a transformation from an XML document model (an XML schema) to an application specific model. The application specific model implements the intended behavior of documents written in the language. A transformation is specified in a model transformation language used in the Model Driven Architecture (MDA) approach for software development. Our approach provides a better separation of three concerns found in XML applications: syntax, syntax processing logic and intended meaning of the syntax. It frees the developer of low-level syntactical details and improves the adaptability and reusability of XML applications. Declarative transformation rules and the explicit application model provide a finer control over the application parts affected by adaptations. Transformation rules and the application model for an XML language may be composed with the corresponding rules and application models defined for other XML languages. In that way we achieve reuse and composition of XML applications