Identification and analysis of novel tandem repeats in the cell surface proteins of archaeal and bacterial genomes using computational tools

cazy_webscraper: local compilation and interrogation of comprehensive CAZyme datasets

Carbohydrate active enzymes (CAZymes) are pivotal in biological processes including energy metabolism, cell structure maintenance, signalling, and pathogen recognition. Bioinformatic prediction and mining of CAZymes improves our understanding of these activities and enables discovery of candidates of interest for industrial biotechnology, particularly the processing of organic waste for biofuel production. CAZy (www.cazy.org) is a high-quality, manually curated, and authoritative database of CAZymes that is often the starting point for these analyses. Automated querying and integration of CAZy data with other public datasets would constitute a powerful resource for mining and exploring CAZyme diversity. However, CAZy does not itself provide methods to automate queries, or integrate annotation data from other sources (except by following hyperlinks) to support further analysis. To overcome these limitations we developed cazy_webscraper, a command-line tool that retrieves data from CAZy and other online resources to build a local, shareable and reproducible database that augments and extends the authoritative CAZy database. cazy_webscraper's integration of curated CAZyme annotations with their corresponding protein sequences, up-to-date taxonomy assignments, and protein structure data facilitates automated large-scale and targeted bioinformatic CAZyme family analysis and candidate screening. This tool has found widespread uptake in the community, with over 35 000 downloads (from April 2021 to June 2023). We demonstrate the use and application of cazy_webscraper to: (i) augment, update and correct CAZy database accessions; (ii) explore the taxonomic distribution of CAZymes recorded in CAZy, identifying under-represented taxa and unusual CAZy class distributions; and (iii) investigate three CAZymes having potential biotechnological application for degradation of biomass, but lacking a representative structure in the PDB database. We describe in general how cazy_webscraper facilitates functional, structural and evolutionary studies to aid identification of candidate enzymes for further characterization, and specifically note that CAZy provides supporting evidence for recent expansion of the Auxiliary Activities (AA) CAZy family in eukaryotes, consistent with functions potentially specific to eukaryotic lifestyles.

Identification and structural analysis of a carbohydrate-binding module specific to alginate, a representative of a new family, CBM96

Carbohydrate-binding modules (CBMs) are the noncatalytic modules that assist functions of the catalytic modules in carbohydrate-active enzymes, and they are usually discrete structural domains in larger multimodular enzymes. CBMs often occur in tandem in different alginate lyases belonging to the CBM families 13, 16, and 32. However, none of the currently known CBMs in alginate lyases specifically bind to an internal alginate chain. In our investigation of the multidomain alginate lyase Dp0100 carrying several ancillary domains, we identified an alginate-binding domain denoted TM6-N4 using protein truncation analysis. The structure of this CBM domain was determined at 1.35 Å resolution. TM6-N4 exhibited an overall β-sandwich fold architecture with two antiparallel β-sheets. We identified an extended binding groove in the CBM using site-directed mutagenesis, docking, and surface electrostatic potential analysis. Affinity analysis revealed that residues of Lys10, Lys22, Lys25, Lys27, Lys31, Arg36, and Tyr159 located on the bottom or the wall of the shallow groove are responsible for alginate binding, and isothermal titration calorimetry analyses indicated that the binding cleft consists of six subsites for sugar recognition. This substrate binding pattern is typical for type B CBM, and it represents the first CBM domain that specifically binds internal alginate chain. Phylogenetic analysis supports that TM6-N4 constitutes the founding member of a new CBM family denoted as CBM96. Our reported structure not only facilitates the investigation of the CBM-alginate ligand recognition mechanism but also inspires the utilization of the CBM domain in biotechnical applications.

Genomic Insights Into the Lifestyles of Thaumarchaeota Inside Sponges

Sponges are among the oldest metazoans and their success is partly due to their abundant and diverse microbial symbionts. They are one of the few animals that have Thaumarchaeota symbionts. Here we compare genomes of 11 Thaumarchaeota sponge symbionts, including three new genomes, to free-living ones. Like their free-living counterparts, sponge-associated Thaumarchaeota can oxidize ammonia, fix carbon, and produce several vitamins. Adaptions to life inside the sponge host include enrichment in transposases, toxin-antitoxin systems and restriction modifications systems, enrichments previously reported also from bacterial sponge symbionts. Most thaumarchaeal sponge symbionts lost the ability to synthesize rhamnose, which likely alters their cell surface and allows them to evade digestion by the host. All but one archaeal sponge symbiont encoded a high-affinity, branched-chain amino acid transporter system that was absent from the analyzed free-living thaumarchaeota suggesting a mixotrophic lifestyle for the sponge symbionts. Most of the other unique features found in sponge-associated Thaumarchaeota, were limited to only a few specific symbionts. These features included the presence of exopolyphosphatases and a glycine cleavage system found in the novel genomes. Thaumarchaeota have thus likely highly specific interactions with their sponge host, which is supported by the limited number of host sponge species to which each of these symbionts is restricted.

Identification of Uncharacterized Components of Prokaryotic Immune Systems and Their Diverse Eukaryotic Reformulations

Nucleotide-activated effector deployment, prototyped by interferon-dependent immunity, is a common mechanistic theme shared by immune systems of several animals and prokaryotes. Prokaryotic versions include CRISPR-Cas with the CRISPR polymerase domain, their minimal variants, and systems with second messenger oligonucleotide or dinucleotide synthetase (SMODS). Cyclic or linear oligonucleotide signals in these systems help set a threshold for the activation of potentially deleterious downstream effectors in response to invader detection. We establish such a regulatory mechanism to be a more general principle of immune systems, which can also operate independently of such messengers. Using sensitive sequence analysis and comparative genomics, we identify 12 new prokaryotic immune systems, which we unify by this principle of threshold-dependent effector activation. These display regulatory mechanisms paralleling physiological signaling based on 3'-5' cyclic mononucleotides, NAD⁺-derived messengers, two- and one-component signaling that includes histidine kinase-based signaling, and proteolytic activation. Furthermore, these systems allowed the identification of multiple new sensory signal sensory components, such as a tetratricopeptide repeat (TPR) scaffold predicted to recognize NAD⁺-derived signals, unreported versions of the STING domain, prokaryotic YEATS domains, and a predicted nucleotide sensor related to receiver domains. We also identify previously unrecognized invader detection components and effector components, such as prokaryotic versions of the Wnt domain. Finally, we show that there have been multiple acquisitions of unidentified STING domains in eukaryotes, while the TPR scaffold was incorporated into the animal immunity/apoptosis signal-regulating kinase (ASK) signalosome.IMPORTANCE Both prokaryotic and eukaryotic immune systems face the dangers of premature activation of effectors and degradation of self-molecules in the absence of an invader. To mitigate this, they have evolved threshold-setting regulatory mechanisms for the triggering of effectors only upon the detection of a sufficiently strong invader signal. This work defines general templates for such regulation in effector-based immune systems. Using this, we identify several previously uncharacterized prokaryotic immune mechanisms that accomplish the regulation of downstream effector deployment by using nucleotide, NAD⁺-derived, two-component, and one-component signals paralleling physiological homeostasis. This study has also helped identify several previously unknown sensor and effector modules in these systems. Our findings also augment the growing evidence for the emergence of key animal immunity and chromatin regulatory components from prokaryotic progenitors.

In silico functional prediction of hypothetical proteins from the core genome of Corynebacterium pseudotuberculosis biovar ovis

Corynebacterium pseudotuberculosis is a pathogen of veterinary relevance diseases, being divided into two biovars: equi and ovis; causing ulcerative lymphangitis and caseous lymphadenitis, respectively. The isolation and sequencing of C. pseudotuberculosis biovar ovis strains in the Northern and Northeastern regions of Brazil exhibited the emergence of this pathogen, which causes economic losses to small ruminant producers, and condemnation of carcasses and skins of animals. Through the pan-genomic approach, it is possible to determine and analyze genes that are shared by all strains of a species—the core genome. However, many of these genes do not have any predicted function, being characterized as hypothetical proteins (HP). In this study, we considered 32 C. pseudotuberculosis biovar ovis genomes for the pan-genomic analysis, where were identified 172 HP present in a core genome composed by 1255 genes. We are able to functionally annotate 80 sequences previously characterized as HP through the identification of structural features as conserved domains and families. Furthermore, we analyzed the physicochemical properties, subcellular localization and molecular function. Additionally, through RNA-seq data, we investigated the differential gene expression of the annotated HP. Genes inserted in pathogenicity islands had their virulence potential evaluated. Also, we have analyzed the existence of functional associations for their products based on protein–protein interaction networks, and perform the structural prediction of three targets. Due to the integration of different strategies, this study can underlie deeper in vitro researches in the characterization of these HP and the search for new solutions for combat this pathogen.

The C-terminal domain of Corynebacterium glutamicum mycoloyltransferase A is composed of five repeated motifs involved in cell wall binding and stability

The genomes of Corynebacteriales contain several genes encoding mycoloyltransferases (Myt) that are specific cell envelope enzymes essential for the biogenesis of the outer membrane. MytA is a major mycoloyltransferase of Corynebacterium glutamicum, displaying an N-terminal domain with esterase activity and a C-terminal extension containing a conserved repeated Leu-Gly-Phe-Pro (LGFP) sequence motif of unknown function. This motif is highly conserved in Corynebacteriales and found associated with cell wall hydrolases and with proteins of unknown function. In this study, we determined the crystal structure of MytA and found that its C-terminal domain is composed of five LGFP motifs and forms a long stalk perpendicular to the N-terminal catalytic α/β-hydrolase domain. The LGFP motifs are composed of a 4-stranded β-fold and occupy alternating orientations along the axis of the stalk. Multiple acetate binding pockets were identified in the stalk, which could correspond to putative ligand-binding sites. By using various MytA mutants and complementary in vitro and in vivo approaches, we provide evidence that the C-terminal LGFP domain interacts with the cell wall peptidoglycan-arabinogalactan polymer. We also show that the C-terminal LGFP domain is not required for the activity of MytA but rather contributes to the overall integrity of the cell envelope.

Satellites in the prokaryote world

Background

Satellites or tandem repeats are very abundant in many eukaryotic genomes. Occasionally they have been reported to be present in some prokaryotes, but to our knowledge there is no general comparative study on their occurrence. For this reason we present here an overview of the distribution and properties of satellites in a set of representative species. Our results provide novel insights into the evolutionary relationship between eukaryotes, Archaea and Bacteria.

Results

We have searched all possible satellites present in the NCBI reference group of genomes in Archaea (142 species) and in Bacteria (119 species), detecting 2735 satellites in Archaea and 1067 in Bacteria. We have found that the distribution of satellites is very variable in different organisms. The archaeal Methanosarcina class stands out for the large amount of satellites in their genomes. Satellites from a few species have similar characteristics to those in eukaryotes, but most species have very few satellites: only 21 species in Archaea and 18 in Bacteria have more than 4 satellites/Mb. The distribution of satellites in these species is reminiscent of what is found in eukaryotes, but we find two significant differences: most satellites have a short length and many of them correspond to segments of genes coding for amino acid repeats. Transposition of non-coding satellites throughout the genome occurs rarely: only in the bacteria Leptospira interrogans and the archaea Methanocella conradii we have detected satellite families of transposed satellites with long repeats.

Conclusions

Our results demonstrate that the presence of satellites in the genome is not an exclusive feature of eukaryotes. We have described a few prokaryotes which do contain satellites. We present a discussion on their eventual evolutionary significance.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1504-2) contains supplementary material, which is available to authorized users.

Sequence and structural analysis of fibronectin-binding protein reveals importance of multiple intrinsic disordered tandem repeats

The location of certain amino acid sequences like repeats along the polypeptide chain is very important in the context of forming the overall shape of the protein molecule which in fact determines its function. In gram-positive bacteria, fibronectin-binding protein (FnBP) is one such repeat containing protein, and it is a cell wall-attached protein responsible for various acute infections in human. Several studies on sequence, structure, and function of fibronectin-binding regions of FnBPs were reported; however, no detailed study was carried out on the full-length protein sequence. In the present study, we have made a thorough sequence and structure analysis on FnBP_A of Staphylococcus aureus and explored the presence of dual ligand-binding ability of fibrinogen (fg)-binding region and its molecular recognition processes. Multiple sequence alignment and protein-protein docking analysis reveal the regions which are likely involved in dual ligand binding. Further analysis of docking of FnBP_A fg-binding region and fn N-terminal modules suggests that if the latter binds to the fg-binding region of FnBP_A, it would inhibit the subsequent binding of fg because of steric hindrance. The sequence analysis further suggests that the abundance of disorder promoting residue glutamic acid and dual personality (both order/disorder promoting) residue threonine in tandem repeats of FnBP_A and B proteins possibly would help the molecule to undergo a conformational change while binding with fn by β-zipper mechanism. The segment-based power spectral analysis was carried out which helps to understand the distribution of hydrophobic residues along the sequence particularly in intrinsic disordered tandem repeats. The results presented here will help to understand the role of internal repeats and intrinsic disorder in the molecular recognition process of a pathogenic cell surface protein.

The Rv3799-Rv3807 Gene Cluster in Mycobacterium Tuberculosis Genome Corresponds to the ‘Ancient Conserved Region’ in CMN Mycolyltransferases

We have identified based on gene cluster analysis that the genes between Rv3799–Rv3807 in M. tuberculosis have orthologs in Corynebacteria, Mycobacteria and Nocardia (CMN) genomes. Therefore, this gene cluster possibly corresponds to the ‘Ancient Conserved Region’ of CMN mycolyltransferases. The evolutionary trace analysis suggests that twelve amino acid residues; Leu39, Trp51, Pro71, Trp82, Trp97, Phe100, Gly124, Ser126, Asp192, Glu230, Gly260 and Trp264 are ‘absolutely conserved’. These amino acid residues constitute the active site and conserved hydrophobic tunnel in CMN mycolyltransferases.

Mycoloyltransferases: A large and major family of enzymes shaping the cell envelope of Corynebacteriales

Mycobacterium and Corynebacterium are important genera of the Corynebacteriales order, the members of which are characterized by an atypical diderm cell envelope. Indeed the cytoplasmic membrane of these bacteria is surrounded by a thick mycolic acid-arabinogalactan-peptidoglycan (mAGP) covalent polymer. The mycolic acid-containing part of this complex associates with other lipids (mainly trehalose monomycolate (TMM) and trehalose dimycolate (TDM)) to form an outer membrane. The metabolism of mycolates in the cell envelope is governed by esterases called mycoloyltransferases that catalyze the transfer of mycoloyl chains from TMM to another TMM molecule or to other acceptors such as the terminal arabinoses of arabinogalactan or specific polypeptides. In this review we present an overview of this family of Corynebacteriales enzymes, starting with their expression, localization, structure and activity to finally discuss their putative functions in the cell. In addition, we show that Corynebacteriales possess multiple mycoloyltransferases encoding genes in their genome. The reason for this multiplicity is not known, as their function in mycolates biogenesis appear to be only partially redundant. It is thus possible that, in some species living in specific environments, some mycoloyltransferases have evolved to gain some new functions. In any case, the few characterized mycoloyltransferases are very important for the bacterial physiology and are also involved in adaptation in the host where they constitute major secreted antigens. Although not discussed in this review, all these functions make them interesting targets for the discovery of new antibiotics and promising vaccines candidates. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Features of two new proteins with OmpA-like domains identified in the genome sequences of Leptospira interrogans

Leptospirosis is an acute febrile disease caused by pathogenic spirochetes of the genus Leptospira. It is considered an important re-emerging infectious disease that affects humans worldwide. The knowledge about the mechanisms by which pathogenic leptospires invade and colonize the host remains limited since very few virulence factors contributing to the pathogenesis of the disease have been identified. Here, we report the identification and characterization of two new leptospiral proteins with OmpA-like domains. The recombinant proteins, which exhibit extracellular matrix-binding properties, are called Lsa46 - LIC13479 and Lsa77 - LIC10050 (Leptospiral surface adhesins of 46 and 77 kDa, respectively). Attachment of Lsa46 and Lsa77 to laminin was specific, dose dependent and saturable, with KD values of 24.3 ± 17.0 and 53.0 ± 17.5 nM, respectively. Lsa46 and Lsa77 also bind plasma fibronectin, and both adhesins are plasminogen (PLG)-interacting proteins, capable of generating plasmin (PLA) and as such, increase the proteolytic ability of leptospires. The proteins corresponding to Lsa46 and Lsa77 are present in virulent L. interrogans L1-130 and in saprophyte L. biflexa Patoc 1 strains, as detected by immunofluorescence. The adhesins are recognized by human leptospirosis serum samples at the onset and convalescent phases of the disease, suggesting that they are expressed during infection. Taken together, our data could offer valuable information to the understanding of leptospiral pathogenesis.

Systematic analysis of compositional order of proteins reveals new characteristics of biological functions and a universal correlate of macroevolution

We present a novel analysis of compositional order (CO) based on the occurrence of Frequent amino-acid Triplets (FTs) that appear much more than random in protein sequences. The method captures all types of proteomic compositional order including single amino-acid runs, tandem repeats, periodic structure of motifs and otherwise low complexity amino-acid regions. We introduce new order measures, distinguishing between ‘regularity’, ‘periodicity’ and ‘vocabulary’, to quantify these phenomena and to facilitate the identification of evolutionary effects. Detailed analysis of representative species across the tree-of-life demonstrates that CO proteins exhibit numerous functional enrichments, including a wide repertoire of particular patterns of dependencies on regularity and periodicity. Comparison between human and mouse proteomes further reveals the interplay of CO with evolutionary trends, such as faster substitution rate in mouse leading to decrease of periodicity, while innovation along the human lineage leads to larger regularity. Large-scale analysis of 94 proteomes leads to systematic ordering of all major taxonomic groups according to FT-vocabulary size. This is measured by the count of Different Frequent Triplets (DFT) in proteomes. The latter provides a clear hierarchical delineation of vertebrates, invertebrates, plants, fungi and prokaryotes, with thermophiles showing the lowest level of FT-vocabulary. Among eukaryotes, this ordering correlates with phylogenetic proximity. Interestingly, in all kingdoms CO accumulation in the proteome has universal characteristics. We suggest that CO is a genomic-information correlate of both macroevolution and various protein functions. The results indicate a mechanism of genomic ‘innovation’ at the peptide level, involved in protein elongation, shaped in a universal manner by mutational and selective forces.

Variations in compositionally ordered (CO) sections of proteins, such as amino acid runs, tandem repeats and low complexity regions, are often considered as a third type of genomic variation along with SNP and CNV. At the microevolutionary scale, they are involved in the rapid evolution of numerous biological functions and the development of novel phenotypic complex traits, including disease in human, in particular neurodegeneration and cancer. At the macroevolutionary scale, the best discriminating proteomic factor between super-kingdoms is the prevalence of CO proteins in eukaryotes. The analysis of CO structures has so far been quite eclectic. Here we introduce a novel unifying methodology, accounting for all types of low-complexity regions and repetitive phenomena, including the existence of large periodic structures in protein sequences. We define new CO measures providing insights into the correlation of CO with protein function and with evolution. In particular, a large-scale analysis of 94 proteomes shows that the CO vocabulary of frequently appearing amino acid triplets serves as a measure of taxonomic ordering separating major clades from each other. It unravels a missing genomic correlate of macroevolution and serves as a novel phylogenetic tool. This suggests that major CO generation occurs during the creation of a completely new species, i.e. during macroevolutionary events.

Genome sequence of the Nocardia bacteriophage NBR1

We here characterize a novel bacteriophage (NBR1) that is lytic for Nocardia otitidiscaviarum and N. brasiliensis. NBR1 is a member of the family Siphoviridae and appears to have a structurally more complex tail than previously reported Siphoviridae phages. NBR1 has a linear genome of 46,140 bp and a sequence that appears novel when compared to those of other phage sequences in GenBank. Annotation of the genome reveals 68 putative open reading frames. The phage genome organization appears to be similar to other Siphoviridae phage genomes in that it has a modular arrangement.

Genome and proteome analysis of phage E3 infecting the soil-borne actinomycete Rhodococcus equi

We report on the characterization and genomic analysis of bacteriophage E3 isolated from soil and propagating in Rhodococcus equi strains. Phage E3 has a circular genome of 142 563 bp and is the first Myoviridae reported for the genus Rhodococcus and for a non-mycobacterial actinomycete. Phylogenetic analyses placed E3 in a distinct Myoviridae clade together with Mycobacterium phages Bxz1 and Myrna. The highly syntenic genomes of this myoviridal group comprise vertically evolving core phage modules flanked by hyperplastic regions specific to each phage and rich in horizontally acquired DNA. The hyperplastic regions contain numerous tRNA genes in the mycobacteriophages which are absent in E3, possibly reflecting bacterial host-specific translation-related phage fitness constraints associated with rate-limiting tRNAs. A structural proteome analysis identified 28 E3 polypeptides, including 15 not previously known to be virion-associated proteins. The E3 genome and comparative analysis provide insight into short-term genome evolution and adaptive plasticity in tailed phages from the environmental microbiome.

Mycobacteriophage endolysins: diverse and modular enzymes with multiple catalytic activities

The mycobacterial cell wall presents significant challenges to mycobacteriophages – viruses that infect mycobacterial hosts – because of its unusual structure containing a mycolic acid-rich mycobacterial outer membrane attached to an arabinogalactan layer that is in turn linked to the peptidoglycan. Although little is known about how mycobacteriophages circumvent these barriers during the process of infection, destroying it for lysis at the end of their lytic cycles requires an unusual set of functions. These include Lysin B proteins that cleave the linkage of mycolic acids to the arabinogalactan layer, chaperones required for endolysin delivery to peptidoglycan, holins that regulate lysis timing, and the endolysins (Lysin As) that hydrolyze peptidoglycan. Because mycobacterial peptidoglycan contains atypical features including 3→3 interpeptide linkages, it is not surprising that the mycobacteriophage endolysins also have non-canonical features. We present here a bioinformatic dissection of these lysins and show that they are highly diverse and extensively modular, with an impressive number of domain organizations. Most contain three domains with a novel N-terminal predicted peptidase, a centrally located amidase, muramidase, or transglycosylase, and a C-terminal putative cell wall binding domain.

Genome sequence and characterization of the Tsukamurella bacteriophage TPA2

The formation of stable foam in activated sludge plants is a global problem for which control is difficult. These foams are often stabilized by hydrophobic mycolic acid-synthesizing Actinobacteria, among which are Tsukamurella spp. This paper describes the isolation from activated sludge of the novel double-stranded DNA phage TPA2. This polyvalent Siphoviridae family phage is lytic for most Tsukamurella species. Whole-genome sequencing reveals that the TPA2 genome is circularly permuted (61,440 bp) and that 70% of its sequence is novel. We have identified 78 putative open reading frames, 95 pairs of inverted repeats, and 6 palindromes. The TPA2 genome has a modular gene structure that shares some similarity to those of Mycobacterium phages. A number of the genes display a mosaic architecture, suggesting that the TPA2 genome has evolved at least in part from genetic recombination events. The genome sequence reveals many novel genes that should inform any future discussion on Tsukamurella phage evolution.

Evolution of diverse cell division and vesicle formation systems in Archaea

Recently a novel cell division system comprised of homologues of eukaryotic ESCRT-III (endosomal sorting complex required for transport III) proteins was discovered in the hyperthermophilic crenarchaeote Sulfolobus acidocaldarius. On the basis of this discovery, we undertook a comparative genomic analysis of the machineries for cell division and vesicle formation in Archaea. Archaea possess at least three distinct membrane remodelling systems: the FtsZ-based bacterial-type system, the ESCRT-III-based eukaryote-like system and a putative novel system that uses an archaeal actin-related protein. Many archaeal genomes encode assortments of components from different systems. Evolutionary reconstruction from these findings suggests that the last common ancestor of the extant Archaea possessed a complex membrane remodelling apparatus, different components of which were lost during subsequent evolution of archaeal lineages. By contrast, eukaryotes seem to have inherited all three ancestral systems.

A framework for classification of prokaryotic protein kinases

BACKGROUND

Overwhelming majority of the Serine/Threonine protein kinases identified by gleaning archaeal and eubacterial genomes could not be classified into any of the well known Hanks and Hunter subfamilies of protein kinases. This is owing to the development of Hanks and Hunter classification scheme based on eukaryotic protein kinases which are highly divergent from their prokaryotic homologues. A large dataset of prokaryotic Serine/Threonine protein kinases recognized from genomes of prokaryotes have been used to develop a classification framework for prokaryotic Ser/Thr protein kinases.

METHODOLOGY/PRINCIPAL FINDINGS

We have used traditional sequence alignment and phylogenetic approaches and clustered the prokaryotic kinases which represent 72 subfamilies with at least 4 members in each. Such a clustering enables classification of prokaryotic Ser/Thr kinases and it can be used as a framework to classify newly identified prokaryotic Ser/Thr kinases. After series of searches in a comprehensive sequence database we recognized that 38 subfamilies of prokaryotic protein kinases are associated to a specific taxonomic level. For example 4, 6 and 3 subfamilies have been identified that are currently specific to phylum proteobacteria, cyanobacteria and actinobacteria respectively. Similarly subfamilies which are specific to an order, sub-order, class, family and genus have also been identified. In addition to these, we also identify organism-diverse subfamilies. Members of these clusters are from organisms of different taxonomic levels, such as archaea, bacteria, eukaryotes and viruses.

CONCLUSION/SIGNIFICANCE

Interestingly, occurrence of several taxonomic level specific subfamilies of prokaryotic kinases contrasts with classification of eukaryotic protein kinases in which most of the popular subfamilies of eukaryotic protein kinases occur diversely in several eukaryotes. Many prokaryotic Ser/Thr kinases exhibit a wide variety of modular organization which indicates a degree of complexity and protein-protein interactions in the signaling pathways in these microbes.

A deficiency in arabinogalactan biosynthesis affects Corynebacterium glutamicum mycolate outer membrane stability

Corynebacterineae is a specific suborder of Gram-positive bacteria that includes Mycobacterium tuberculosis and Corynebacterium glutamicum. The ultrastructure of the cell envelope is very atypical. It is composed of a heteropolymer of peptidoglycan and arabinogalactan (AG) covalently associated to an outer membrane. Five arabinosyltransferases are involved in the biosynthesis of AG in C. glutamicum. AftB catalyzes the transfer of Araf (arabinofuranosyl) onto the arabinan domain of the arabinogalactan to form terminal beta(1 --> 2)-linked Araf residues. Here we show that Delta aftB cells lack half of the arabinogalactan mycoloylation sites but are still able to assemble an outer membrane. In addition, we show that a Delta aftB mutant grown on a rich medium has a perturbed cell envelope and sheds a significant amount of membrane fragments in the external culture medium. These fragments contain mono- and dimycolate of trehalose and PorA/H, the major porin of C. glutamicum, but lack conventional phospholipids that typify the plasma membrane, suggesting that they are derived from the atypical mycolate outer membrane of the cell envelope. This is the first report of outer membrane destabilization in the Corynebacterineae, and it suggests that a strong interaction between the mycolate outer membrane and the underlying polymer is essential for cell envelope integrity. The presence of outer membrane-derived fragments (OMFs) in the external medium of the Delta aftB mutant is also a very promising tool for outer membrane characterization. Indeed, fingerprint analysis of major OMF-associated proteins has already led to the identification of 3 associated mycoloyltransferases and an unknown protein with a C-terminal hydrophobic anchoring domain reminiscent of that found for the S-layer protein PS2 of C. glutamicum.

CPDadh: a new peptidase family homologous to the cysteine protease domain in bacterial MARTX toxins

A cysteine protease domain (CPD) has been recently discovered in a group of multifunctional, autoprocessing RTX toxins (MARTX) and Clostridium difficile toxins A and B. These CPDs (referred to as CPDmartx) autocleave the toxins to release domains with toxic effects inside host cells. We report identification and computational analysis of CPDadh, a new cysteine peptidase family homologous to CPDmartx. CPDadh and CPDmartx share a Rossmann-like structural core and conserved catalytic residues. In bacteria, domains of the CPDadh family are present at the N-termini of a diverse group of putative cell-cell interaction proteins and at the C-termini of some RHS (recombination hot spot) proteins. In eukaryotes, catalytically inactive members of the CPDadh family are found in cell surface protein NELF (nasal embryonic LHRH factor) and some putative signaling proteins.

S-layer, surface-accessible, and concanavalin A binding proteins of Methanosarcina acetivorans and Methanosarcina mazei

The outermost cell envelope structure of many archaea and bacteria contains a proteinaceous lattice termed the surface layer or S-layer. It is typically composed of only one or two abundant, often posttranslationally modified proteins that self-assemble to form the highly organized arrays. Surprisingly, over 100 proteins were annotated to be S-layer components in the archaeal species Methanosarcina acetivorans C2A and Methanosarcina mazei Gö1, reflecting limitations of current predictions. An in vivo biotinylation methodology was devised to affinity tag surface-exposed proteins while overcoming unique challenges in working with these fragile organisms. Cells were adapted to growth under N2 fixing conditions, thus, minimizing free amines reactive to the NHS-label, and high pH media compatible with the acylation chemistry was used. A 3-phase separation procedure was employed to isolate intact, labeled cells from lysed-cell derived proteins. Streptavidin affinity enrichment followed by stringent wash conditions removed nonspecifically bound proteins. This methodology revealed S-layer proteins in M. acetivorans C2A and M. mazei Gö1 to be MA0829 and MM1976, respectively. Each was demonstrated to exist as multiple glycosylated forms using SDS-PAGE coupled with glycoprotein-specific staining, and by interaction with the lectin, Concanavalin A. A number of additional surface-exposed proteins and glycoproteins were identified and included all three subunits of the thermosome: the latter suggests that the chaperonin complex is both surface- and cytoplasmically localized. This approach provides an alternative strategy to study surface proteins in the archaea.

Identification of the gene for disaggregatase from Methanosarcina mazei

The gene sequences encoding disaggregatase (Dag), the enzyme responsible for dispersion of cell aggregates of Methanosarcina mazei to single cells, were determined for three strains of M. mazei (S-6(T), LYC and TMA). The dag genes of the three strains were 3234 bp in length and had almost the same sequences with 97% amino acid sequence identities. Dag was predicted to comprise 1077 amino acid residues and to have a molecular mass of 120 kDa containing three repeats of the DNRLRE domain in the C terminus, which is specific to the genus Methanosarcina and may be responsible for structural organization and cell wall function. Recombinant Dag was overexpressed in Escherichia coli and preparations of the expressed protein exhibited enzymatic activity. The RT-PCR analysis showed that dag was transcribed to mRNA in M. mazei LYC and indicated that the gene was expressed in vivo. This is the first time the gene involved in the morphological change of Methanosarcina spp. from aggregate to single cells has been identified.

Cell surface proteins in archaeal and bacterial genomes comprising "LVIVD", "RIVW" and "LGxL" tandem sequence repeats are predicted to fold as beta-propeller

Proteins that share even low sequence homologies are known to adopt similar folds. The beta-propeller structural motif is one such example. Identifying sequences that adopt a beta-propeller fold is useful to annotate protein structure and function. Often, tandem sequence repeats provide the necessary signal for identifying beta-propellers in proteins. In our recent analysis to identify cell surface proteins in archaeal and bacterial genomes, we identified some proteins that contain novel tandem repeats "LVIVD", "RIVW" and "LGxL". In this work, based on protein fold predictions and three-dimensional comparative modeling methods, we predicted that these repeat types fold as beta-propeller. Further, the evolutionary trace analysis of all proteins constituting amino acid sequence repeats in beta-propellers suggest that the novel repeats have diverged from a common ancestor.

Classification of proteins based on minimal modular repeats: lessons from nature in protein design

Proteins containing internal repeats within their primary sequence have received increased attention recently, as the extent of their presence in various organisms is recognized more fully, and their role in evolution is more thoroughly studied. Presented here is a technique used to detect and classify proteins based on a modular evolutionary phenomenon that results in a series of small internal repeats. The parameters chosen are based on a minimum segment of seven residues that result in simple functional scaffolds. The genomes and corresponding proteomes of a variety of eubacteria and archaea have been analyzed using an algorithm that searches prokaryotic genomes for proteins containing small conserved repeats assembled in a modular fashion similar to a recently characterized protein from the organism Nitrosomonas europaea. This analysis has revealed additional proteins present in N. europaea with similar modular characteristics. A further survey of a variety of organisms demonstrates that this evolutionary pathway has been utilized in other organisms as well, to yield a broad assortment of small modular proteins. A thorough description of the sequential characteristics of these modular proteins follows, along with a selection and discussion of the various proteins uncovered through this expanded search and analysis. Several databases of the proteins uncovered from this work and the program used to perform the search are available.

The Rv3799-Rv3807 gene cluster in Mycobacterium tuberculosis genome corresponds to the 'Ancient Conserved Region' in CMN mycolyltransferases

We have identified based on gene cluster analysis that the genes between Rv3799-Rv3807 in M. tuberculosis have orthologs in Corynebacteria, Mycobacteria and Nocardia (CMN) genomes. Therefore, this gene cluster possibly corresponds to the 'Ancient Conserved Region' of CMN mycolyltransferases. The evolutionary trace analysis suggests that twelve amino acid residues; Leu39, Trp51, Pro71, Trp82, Trp97, Phe100, Gly124, Ser126, Asp192, Glu230, Gly260 and Trp264 are 'absolutely conserved'. These amino acid residues constitute the active site and conserved hydrophobic tunnel in CMN mycolyltransferases.

Analysis and modeling of mycolyl-transferases in the CMN group

Mycolyl-transferases are a family of proteins that are specifically present in the CMN (Corynebacterium, Mycobacterium and Nocardia) genera and are responsible for the synthesis of cell wall components. We modeled the three-dimensional structures of mycolyl-transfersases from Corynebacterium and Nocardia using homology modeling methods based on the crystal structures of mycolyl-transferases from M. tuberculosis. Comparison of the models revealed significant differences in their substrate binding site. Some mycolyl-transferases identified by the following Gene Ids: Nfa25110, Nfa45560, Nfa7210, Nfa38260, Nfa32420, Nfa23770, Nfa43800, Nfa30260, Dip0365, Ncgl0987, Ce1488, Ncgl0885, Ce0984, Ncgl2101, Ncgl0336, Ce0356 are associated with a relatively larger substrate binding site and amino acid residue mutations (D40N, R43D/G, S236N/A) are likely to affect binding to trehalose.

Three-dimensional models and structure analysis of corynemycolyltransferases in Corynebacterium glutamicum and Corynebacterium efficiens

The corynemycolyltransferase proteins were identified from Corynebacterium glutamicum and Corynebacterium efficiens genomes using computational tools available in the public domain. Three-dimensional models were constructed for corynemycolyltransferases based on the crystal structures of related mycolyltransferases in Mycobacterium tuberculosis using the comparative modeling methods. The corynemycolyltransferases share overall an alpha/beta-fold characteristic of the mycolyltransferases despite low sequence identity (<20%) shared by some of the corynemycolyltransferases. However, a significant difference is observed in the region between amino acid residues Trp82-Trp97 and Ala222-Asn223 corresponding to mycolyltransferases. The specificity pockets defined by interactions with the trehalose substrate observed in the crystal structure complex of Ag85B mycolyltransferase (PDB code: 1F0P) suggests that trehalose may not bind some corynemycolyltransferases. This is due to critical mutations in corynemycolyltransferase binding subsites that lead to loss of equivalent side-chain interactions with trehalose and unfavorable steric interactions, particularly, in the case of cmytC gene and the protein corresponding to the gene identifier CE0356 with the equivalent Ala222-Asn223 "long insertion loop". Further, the fibronectin binding region (Phe58-Val69), in mycolyltransferases associated with mediating host-pathogen interactions in M. tuberculosis comprises amino acid residue mutations in the corresponding region in the soil bacterium--Corynebacterium corynemycolyltransferases, that suggest a different epitope and therefore possible lack of binding to fibronectin. The corynemycolyltransferase cmytA responsible for the cell shape formation and for maintaining the cell surface integrity is associated with a C-terminal domain that we have recently shown to comprise tandem amino acid sequence repeats that is likely to be associated with a regular secondary structural motif.