HHsenser: exhaustive transitive profile search using HMM-HMM comparison

CoCoNuTs are a diverse subclass of Type IV restriction systems predicted to target RNA

A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote coiled-coil nuclease tandems (CoCoNuTs) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with three distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF (CRISPR-Associated Rossmann Fold) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over.

CoCoNuTs: A diverse subclass of Type IV restriction systems predicted to target RNA

A comprehensive census of McrBC systems, among the most common forms of prokaryotic Type IV restriction systems, followed by phylogenetic analysis, reveals their enormous abundance in diverse prokaryotes and a plethora of genomic associations. We focus on a previously uncharacterized branch, which we denote CoCoNuTs (coiled-coil nuclease tandems) for their salient features: the presence of extensive coiled-coil structures and tandem nucleases. The CoCoNuTs alone show extraordinary variety, with 3 distinct types and multiple subtypes. All CoCoNuTs contain domains predicted to interact with translation system components, such as OB-folds resembling the SmpB protein that binds bacterial transfer-messenger RNA (tmRNA), YTH-like domains that might recognize methylated tmRNA, tRNA, or rRNA, and RNA-binding Hsp70 chaperone homologs, along with RNases, such as HEPN domains, all suggesting that the CoCoNuTs target RNA. Many CoCoNuTs might additionally target DNA, via McrC nuclease homologs. Additional restriction systems, such as Type I RM, BREX, and Druantia Type III, are frequently encoded in the same predicted superoperons. In many of these superoperons, CoCoNuTs are likely regulated by cyclic nucleotides, possibly, RNA fragments with cyclic termini, that bind associated CARF (CRISPR-Associated Rossmann Fold) domains. We hypothesize that the CoCoNuTs, together with the ancillary restriction factors, employ an echeloned defense strategy analogous to that of Type III CRISPR-Cas systems, in which an immune response eliminating virus DNA and/or RNA is launched first, but then, if it fails, an abortive infection response leading to PCD/dormancy via host RNA cleavage takes over.

A Structural Approach to Anti-Virulence: A Discovery Pipeline

The anti-virulence strategy is designed to prevent bacterial virulence factors produced by pathogenic bacteria from initiating and sustaining an infection. One family of bacterial virulence factors is the mono-ADP-ribosyltransferase toxins, which are produced by pathogens as tools to compromise the target host cell. These toxins are bacterial enzymes that exploit host cellular NAD+ as the donor substrate to modify an essential macromolecule acceptor target in the host cell. This biochemical reaction modifies the target macromolecule (often protein or DNA) and functions in a binary fashion to turn the target activity on or off by blocking or impairing a critical process or pathway in the host. A structural biology approach to the anti-virulence method to neutralize the cytotoxic effect of these factors requires the search and design of small molecules that bind tightly to the enzyme active site and prevent catalytic function essentially disarming the pathogen. This method requires a high-resolution structure to serve as the model for small molecule inhibitor development, which illuminates the path to drug development. This alternative strategy to antibiotic therapy represents a paradigm shift that may circumvent multi-drug resistance in the offending microbe through anti-virulence therapy. In this report, the rationale for the anti-virulence structural approach will be discussed along with recent efforts to apply this method to treat honey bee diseases using natural products.

Genomic and molecular evidence reveals novel pathways associated with cell surface polysaccharides in bacteria

Amino acid (acyl carrier protein) ligases (AALs) are a relatively new family of bacterial amino acid adenylating enzymes with unknown function(s). Here, genomic enzymology tools that employ sequence similarity networks and genome context analyses were used to hypothesize the metabolic function(s) of AALs. In over 50% of species, aal and its cognate acyl carrier protein (acp) genes, along with three more genes, formed a highly conserved AAL cassette. AAL cassettes were strongly associated with surface polysaccharide gene clusters in Proteobacteria and Actinobacteria, yet were prevalent among soil and rhizosphere-associated α- and β-Proteobacteria, including symbiotic α- and β-rhizobia and some Mycolata. Based on these associations, AAL cassettes were proposed to encode a noncanonical Acp-dependent polysaccharide modification route. Genomic-inferred predictions were substantiated by published experimental evidence, revealing a role for AAL cassettes in biosynthesis of biofilm-forming exopolysaccharide in pathogenic Burkholderia and expression of aal and acp genes in nitrogen-fixing Rhizobium bacteroids. Aal and acp genes were associated with dltBD-like homologs that modify cell wall teichoic acids with d-alanine, including in Paenibacillus and certain other bacteria. Characterization of pathways that involve AAL and Acp may lead to developing new plant and human disease-controlling agents as well as strains with improved nitrogen fixation capacity.

Several New Putative Bacterial ADP-Ribosyltransferase Toxins Are Revealed from In Silico Data Mining, Including the Novel Toxin Vorin, Encoded by the Fire Blight Pathogen Erwinia amylovora

Mono-ADP-ribosyltransferase (mART) toxins are secreted by several pathogenic bacteria that disrupt vital host cell processes in deadly diseases like cholera and whooping cough. In the last two decades, the discovery of mART toxins has helped uncover the mechanisms of disease employed by pathogens impacting agriculture, aquaculture, and human health. Due to the current abundance of mARTs in bacterial genomes, and an unprecedented availability of genomic sequence data, mART toxins are amenable to discovery using an in silico strategy involving a series of sequence pattern filters and structural predictions. In this work, a bioinformatics approach was used to discover six bacterial mART sequences, one of which was a functional mART toxin encoded by the plant pathogen, Erwinia amylovora, called Vorin. Using a yeast growth-deficiency assay, we show that wild-type Vorin inhibited yeast cell growth, while catalytic variants reversed the growth-defective phenotype. Quantitative mass spectrometry analysis revealed that Vorin may cause eukaryotic host cell death by suppressing the initiation of autophagic processes. The genomic neighbourhood of Vorin indicated that it is a Type-VI-secreted effector, and co-expression experiments showed that Vorin is neutralized by binding of a cognate immunity protein, VorinI. We demonstrate that Vorin may also act as an antibacterial effector, since bacterial expression of Vorin was not achieved in the absence of VorinI. Vorin is the newest member of the mART family; further characterization of the Vorin/VorinI complex may help refine inhibitor design for mART toxins from other deadly pathogens.

Modified base-binding EVE and DCD domains: striking diversity of genomic contexts in prokaryotes and predicted involvement in a variety of cellular processes

BACKGROUND

DNA and RNA of all cellular life forms and many viruses contain an expansive repertoire of modified bases. The modified bases play diverse biological roles that include both regulation of transcription and translation, and protection against restriction endonucleases and antibiotics. Modified bases are often recognized by dedicated protein domains. However, the elaborate networks of interactions and processes mediated by modified bases are far from being completely understood.

RESULTS

We present a comprehensive census and classification of EVE domains that belong to the PUA/ASCH domain superfamily and bind various modified bases in DNA and RNA. We employ the "guilt by association" approach to make functional inferences from comparative analysis of bacterial and archaeal genomes, based on the distribution and associations of EVE domains in (predicted) operons and functional networks of genes. Prokaryotes encode two classes of EVE domain proteins, slow-evolving and fast-evolving ones. Slow-evolving EVE domains in α-proteobacteria are embedded in conserved operons, potentially involved in coupling between translation and respiration, cytochrome c biogenesis in particular, via binding 5-methylcytosine in tRNAs. In β- and γ-proteobacteria, the conserved associations implicate the EVE domains in the coordination of cell division, biofilm formation, and global transcriptional regulation by non-coding 6S small RNAs, which are potentially modified and bound by the EVE domains. In eukaryotes, the EVE domain-containing THYN1-like proteins have been reported to inhibit PCD and regulate the cell cycle, potentially, via binding 5-methylcytosine and its derivatives in DNA and/or RNA. We hypothesize that the link between PCD and cytochrome c was inherited from the α-proteobacterial and proto-mitochondrial endosymbiont and, unexpectedly, could involve modified base recognition by EVE domains. Fast-evolving EVE domains are typically embedded in defense contexts, including toxin-antitoxin modules and type IV restriction systems, suggesting roles in the recognition of modified bases in invading DNA molecules and targeting them for restriction. We additionally identified EVE-like prokaryotic Development and Cell Death (DCD) domains that are also implicated in defense functions including PCD. This function was inherited by eukaryotes, but in animals, the DCD proteins apparently were displaced by the extended Tudor family proteins, whose partnership with Piwi-related Argonautes became the centerpiece of the Piwi-interacting RNA (piRNA) system.

CONCLUSIONS

Recognition of modified bases in DNA and RNA by EVE-like domains appears to be an important, but until now, under-appreciated, common denominator in a variety of processes including PCD, cell cycle control, antivirus immunity, stress response, and germline development in animals.

A Graph-Based Approach for Detecting Sequence Homology in Highly Diverged Repeat Protein Families

Reconstructing evolutionary relationships in repeat proteins is notoriously difficult due to the high degree of sequence divergence that typically occurs between duplicated repeats. This is complicated further by the fact that proteins with a large number of similar repeats are more likely to produce significant local sequence alignments than proteins with fewer copies of the repeat motif. Furthermore, biologically correct sequence alignments are sometimes impossible to achieve in cases where insertion or translocation events disrupt the order of repeats in one of the sequences being aligned. Combined, these attributes make traditional phylogenetic methods for studying protein families unreliable for repeat proteins, due to the dependence of such methods on accurate sequence alignment.We present here a practical solution to this problem, making use of graph clustering combined with the open-source software package HH-suite, which enables highly sensitive detection of sequence relationships. Carrying out multiple rounds of homology searches via alignment of profile hidden Markov models, large sets of related proteins are generated. By representing the relationships between proteins in these sets as graphs, subsequent clustering with the Markov cluster algorithm enables robust detection of repeat protein subfamilies.

Large-scale mapping of bioactive peptides in structural and sequence space

Health-enhancing potential bioactive peptide (BP) has driven an interest in food proteins as well as in the development of predictive methods. Research in this area has been especially active to use them as components in functional foods. Apparently, BPs do not have a given biological function in the containing proteins and they do not evolve under independent evolutionary constraints. In this work we performed a large-scale mapping of BPs in sequence and structural space. Using well curated BP deposited in BIOPEP database, we searched for exact matches in non-redundant sequences databases. Proteins containing BPs, were used in fold-recognition methods to predict the corresponding folds and BPs occurrences were mapped. We found that fold distribution of BP occurrences possibly reflects sequence relative abundance in databases. However, we also found that proteins with 5 or more than 5 BP in their sequences correspond to well populated protein folds, called superfolds. Also, we found that in well populated superfamilies, BPs tend to adopt similar locations in the protein fold, suggesting the existence of hotspots. We think that our results could contribute to the development of new bioinformatics pipeline to improve BP detection.

A Completely Reimplemented MPI Bioinformatics Toolkit with a New HHpred Server at its Core

The MPI Bioinformatics Toolkit (https://toolkit.tuebingen.mpg.de) is a free, one-stop web service for protein bioinformatic analysis. It currently offers 34 interconnected external and in-house tools, whose functionality covers sequence similarity searching, alignment construction, detection of sequence features, structure prediction, and sequence classification. This breadth has made the Toolkit an important resource for experimental biology and for teaching bioinformatic inquiry. Recently, we replaced the first version of the Toolkit, which was released in 2005 and had served around 2.5 million queries, with an entirely new version, focusing on improved features for the comprehensive analysis of proteins, as well as on promoting teaching. For instance, our popular remote homology detection server, HHpred, now allows pairwise comparison of two sequences or alignments and offers additional profile HMMs for several model organisms and domain databases. Here, we introduce the new version of our Toolkit and its application to the analysis of proteins.

Phylogenomics of Cas4 family nucleases

Background

The Cas4 family endonuclease is a component of the adaptation module in many variants of CRISPR-Cas adaptive immunity systems. Unlike most of the other Cas proteins, Cas4 is often encoded outside CRISPR-cas loci (solo-Cas4) and is also found in mobile genetic elements (MGE-Cas4).

Results

As part of our ongoing investigation of CRISPR-Cas evolution, we explored the phylogenomics of the Cas4 family. About 90% of the archaeal genomes encode Cas4 compared to only about 20% of the bacterial genomes. Many archaea encode both the CRISPR-associated form (CAS-Cas4) and solo-Cas4, whereas in bacteria, this combination is extremely rare. The solo-cas4 genes are over-represented in environmental bacteria and archaea with small genomes that typically lack CRISPR-Cas, suggesting that Cas4 could perform uncharacterized defense or repair functions in these microbes. Phylogenomic analysis indicates that both the CRISPR-associated cas4 genes are often transferred horizontally but almost exclusively, as part of the adaptation module. The evolutionary integrity of the adaptation module sharply contrasts the rampant shuffling of CRISPR-cas modules whereby a given variant of the adaptation module can combine with virtually any effector module. The solo-cas4 genes evolve primarily via vertical inheritance and are subject only to occasional horizontal transfer. The selection pressure on cas4 genes does not substantially differ between CAS-Cas4 and solo-cas4, and is close to the genomic median. Thus, cas4 genes, similarly to cas1 and cas2, evolve similarly to ‘regular’ microbial genes involved in various cellular functions, showing no evidence of direct involvement in virus-host arms races. A notable feature of the Cas4 family evolution is the frequent recruitment of cas4 genes by various mobile genetic elements (MGE), particularly, archaeal viruses. The functions of Cas4 in these elements are unknown and potentially might involve anti-defense roles.

Conclusions

Unlike most of the other Cas proteins, Cas4 family members are as often encoded by stand-alone genes as they are incorporated in CRISPR-Cas systems. In addition, cas4 genes were repeatedly recruited by MGE, perhaps, for anti-defense functions. Experimental characterization of the solo and MGE-encoded Cas4 nucleases is expected to reveal currently uncharacterized defense and anti-defense systems and their interactions with CRISPR-Cas systems.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1081-1) contains supplementary material, which is available to authorized users.

Bacterial kinesin light chain (Bklc) links the Btub cytoskeleton to membranes

Bacterial kinesin light chain is a TPR domain-containing protein encoded by the bklc gene, which co-localizes with the bacterial tubulin (btub) genes in a conserved operon in Prosthecobacter. Btub heterodimers show high structural homology with eukaryotic tubulin and assemble into head-to-tail protofilaments. Intriguingly, Bklc is homologous to the light chain of the microtubule motor kinesin and could thus represent an additional eukaryotic-like cytoskeletal element in bacteria. Using biochemical characterization as well as cryo-electron tomography we show here that Bklc interacts specifically with Btub protofilaments, as well as lipid vesicles and could thus play a role in anchoring the Btub filaments to the membrane protrusions in Prosthecobacter where they specifically localize in vivo. This work sheds new light into possible ways in which the microtubule cytoskeleton may have evolved linking precursors of microtubules to the membrane via the kinesin moiety that in today’s eukaryotic cytoskeleton links vesicle-packaged cargo to microtubules.

Tuned SMC Arms Drive Chromosomal Loading of Prokaryotic Condensin

SMC proteins support vital cellular processes in all domains of life by organizing chromosomal DNA. They are composed of ATPase “head” and “hinge“ dimerization domains and a connecting coiled-coil “arm.” Binding to a kleisin subunit creates a closed tripartite ring, whose ∼47-nm-long SMC arms act as barrier for DNA entrapment. Here, we uncover another, more active function of the bacterial Smc arm. Using high-throughput genetic engineering, we resized the arm in the range of 6–60 nm and found that it was functional only in specific length regimes following a periodic pattern. Natural SMC sequences reflect these length constraints. Mutants with improper arm length or peptide insertions in the arm efficiently target chromosomal loading sites and hydrolyze ATP but fail to use ATP hydrolysis for relocation onto flanking DNA. We propose that SMC arms implement force transmission upon nucleotide hydrolysis to mediate DNA capture or loop extrusion.

•

Short and long but not intermediate-length Smc coiled-coil arms are functional

•

Smc complexes with improper arms accumulate at chromosomal loading sites

•

Smc arms are functional units linking ATP hydrolysis to an essential DNA transaction

•

Pro- and eukaryotic SMC sequences reflect similar periodic length constraints

The MPI bioinformatics Toolkit as an integrative platform for advanced protein sequence and structure analysis

The MPI Bioinformatics Toolkit (http://toolkit.tuebingen.mpg.de) is an open, interactive web service for comprehensive and collaborative protein bioinformatic analysis. It offers a wide array of interconnected, state-of-the-art bioinformatics tools to experts and non-experts alike, developed both externally (e.g. BLAST+, HMMER3, MUSCLE) and internally (e.g. HHpred, HHblits, PCOILS). While a beta version of the Toolkit was released 10 years ago, the current production-level release has been available since 2008 and has serviced more than 1.6 million external user queries. The usage of the Toolkit has continued to increase linearly over the years, reaching more than 400 000 queries in 2015. In fact, through the breadth of its tools and their tight interconnection, the Toolkit has become an excellent platform for experimental scientists as well as a useful resource for teaching bioinformatic inquiry to students in the life sciences. In this article, we report on the evolution of the Toolkit over the last ten years, focusing on the expansion of the tool repertoire (e.g. CS-BLAST, HHblits) and on infrastructural work needed to remain operative in a changing web environment.

Properties and Phylogeny of 76 Families of Bacterial and Eukaryotic Organellar Outer Membrane Pore-Forming Proteins

We here report statistical analyses of 76 families of integral outer membrane pore-forming proteins (OMPPs) found in bacteria and eukaryotic organelles. 47 of these families fall into one superfamily (SFI) which segregate into fifteen phylogenetic clusters. Families with members of the same protein size, topology and substrate specificities often cluster together. Virtually all OMPP families include only proteins that form transmembrane pores. Nine such families, all of which cluster together in the SFI phylogenetic tree, contain both α- and β-structures, are multi domain, multi subunit systems, and transport macromolecules. Most other SFI OMPPs transport small molecules. SFII and SFV homologues derive from Actinobacteria while SFIII and SFIV proteins derive from chloroplasts. Three families of actinobacterial OMPPs and two families of eukaryotic OMPPs apparently consist primarily of α-helices (α-TMSs). Of the 71 families of (putative) β-barrel OMPPs, only twenty could not be assigned to a superfamily, and these derived primarily from Actinobacteria (1), chloroplasts (1), spirochaetes (8), and proteobacteria (10). Proteins were identified in which two or three full length OMPPs are fused together. Family characteristic are described and evidence agrees with a previous proposal suggesting that many arose by adjacent β-hairpin structural unit duplications.

3697G>A in MT-ND1 is a causative mutation in mitochondrial disease

Mitochondrial diseases are a group of clinically heterogeneous disorders that can be difficult to diagnose. We report a two and a half year old girl with clinical symptoms compatible with Leigh disease but with no definitive diagnosis. Using next generation sequencing we found that mutation 3697G>A was responsible for the patient's clinical symptoms. Corroboration was performed via segregation analysis in mother and sister and by evolutionary analysis that showed that the mutation is located in a highly conserved region across a wide range of species. Functional analyses corroborated the mutation effect and indicated that the pathophysiological alterations were partially restored by Coenzyme Q10. In addition, we proposed that the presence of the mutation at high frequencies causes the phenotype in the patient, while other family members with intermediate levels of heteroplasmy are symptoms-free.

Recent Mobility of Casposons, Self-Synthesizing Transposons at the Origin of the CRISPR-Cas Immunity

Casposons are a superfamily of putative self-synthesizing transposable elements that are predicted to employ a homolog of Cas1 protein as a recombinase and could have contributed to the origin of the CRISPR-Cas adaptive immunity systems in archaea and bacteria. Casposons remain uncharacterized experimentally, except for the recent demonstration of the integrase activity of the Cas1 homolog, and given their relative rarity in archaea and bacteria, original comparative genomic analysis has not provided direct indications of their mobility. Here, we report evidence of casposon mobility obtained by comparison of the genomes of 62 strains of the archaeon Methanosarcina mazei. In these genomes, casposons are variably inserted in three distinct sites indicative of multiple, recent gains, and losses. Some casposons are inserted into other mobile genetic elements that might provide vehicles for horizontal transfer of the casposons. Additionally, many M. mazei genomes contain previously undetected solo terminal inverted repeats that apparently are derived from casposons and could resemble intermediates in CRISPR evolution. We further demonstrate the sequence specificity of casposon insertion and note clear parallels with the adaptation mechanism of CRISPR-Cas. Finally, besides identifying additional representatives in each of the three originally defined families, we describe a new, fourth, family of casposons.

Outer membrane protein OmpQ of Bordetella bronchiseptica is required for mature biofilm formation

Bordetella bronchiseptica, an aerobic Gram-negative bacterium, is capable of colonizing the respiratory tract of diverse animals and chronically persists inside the hosts by forming biofilm. Most known virulence factors in Bordetella species are regulated by the BvgAS two-component transduction system. The Bvg-activated proteins play a critical role during host infection. OmpQ is an outer membrane porin protein which is expressed under BvgAS control. Here, we studied the contribution of OmpQ to the biofilm formation process by B. bronchiseptica. We found that the lack of expression of OmpQ did not affect the growth kinetics and final biomass of B. bronchiseptica under planktonic growth conditions. The ΔompQ mutant strain displayed no differences in attachment level and in early steps of biofilm formation. However, deletion of the ompQ gene attenuated the ability of B. bronchiseptica to form a mature biofilm. Analysis of ompQ gene expression during the biofilm formation process by B. bronchiseptica showed a dynamic expression pattern, with an increase of biofilm culture at 48 h. Moreover, we demonstrated that the addition of serum anti-OmpQ had the potential to reduce the biofilm biomass formation in a dose-dependent manner. In conclusion, we showed for the first time, to the best of our knowledge, evidence of the contribution of OmpQ to a process of importance for B. bronchiseptica pathobiology. Our results indicate that OmpQ plays a role during the biofilm development process, particularly at later stages of development, and that this porin could be a potential target for strategies of biofilm formation inhibition.

Discovery and Functional Characterization of Diverse Class 2 CRISPR-Cas Systems

Microbial CRISPR-Cas systems are divided into Class 1, with multisubunit effector complexes, and Class 2, with single protein effectors. Currently, only two Class 2 effectors, Cas9 and Cpf1, are known. We describe here three distinct Class 2 CRISPR-Cas systems. The effectors of two of the identified systems, C2c1 and C2c3, contain RuvC-like endonuclease domains distantly related to Cpf1. The third system, C2c2, contains an effector with two predicted HEPN RNase domains. Whereas production of mature CRISPR RNA (crRNA) by C2c1 depends on tracrRNA, C2c2 crRNA maturation is tracrRNA independent. We found that C2c1 systems can mediate DNA interference in a 5'-PAM-dependent fashion analogous to Cpf1. However, unlike Cpf1, which is a single-RNA-guided nuclease, C2c1 depends on both crRNA and tracrRNA for DNA cleavage. Finally, comparative analysis indicates that Class 2 CRISPR-Cas systems evolved on multiple occasions through recombination of Class 1 adaptation modules with effector proteins acquired from distinct mobile elements.

Cpf1 is a single RNA-guided endonuclease of a class 2 CRISPR-Cas system

The microbial adaptive immune system CRISPR mediates defense against foreign genetic elements through two classes of RNA-guided nuclease effectors. Class 1 effectors utilize multi-protein complexes, whereas class 2 effectors rely on single-component effector proteins such as the well-characterized Cas9. Here, we report characterization of Cpf1, a putative class 2 CRISPR effector. We demonstrate that Cpf1 mediates robust DNA interference with features distinct from Cas9. Cpf1 is a single RNA-guided endonuclease lacking tracrRNA, and it utilizes a T-rich protospacer-adjacent motif. Moreover, Cpf1 cleaves DNA via a staggered DNA double-stranded break. Out of 16 Cpf1-family proteins, we identified two candidate enzymes from Acidaminococcus and Lachnospiraceae, with efficient genome-editing activity in human cells. Identifying this mechanism of interference broadens our understanding of CRISPR-Cas systems and advances their genome editing applications.

Vaccinia virus protein A49 is an unexpected member of the B-cell Lymphoma (Bcl)-2 protein family

Vaccinia virus (VACV) encodes several proteins that inhibit activation of the proinflammatory transcription factor nuclear factor κB (NF-κB). VACV protein A49 prevents translocation of NF-κB to the nucleus by sequestering cellular β-TrCP, a protein required for the degradation of the inhibitor of κB. A49 does not share overall sequence similarity with any protein of known structure or function. We solved the crystal structure of A49 from VACV Western Reserve to 1.8 Å resolution and showed, surprisingly, that A49 has the same three-dimensional fold as Bcl-2 family proteins despite lacking identifiable sequence similarity. Whereas Bcl-2 family members characteristically modulate cellular apoptosis, A49 lacks a surface groove suitable for binding BH3 peptides and does not bind proapoptotic Bcl-2 family proteins Bax or Bak. The N-terminal 17 residues of A49 do not adopt a single well ordered conformation, consistent with their proposed role in binding β-TrCP. Whereas pairs of A49 molecules interact symmetrically via a large hydrophobic surface in crystallo, A49 does not dimerize in solution or in cells, and we propose that this hydrophobic interaction surface may mediate binding to a yet undefined cellular partner. A49 represents the eleventh VACV Bcl-2 family protein and, despite these proteins sharing very low sequence identity, structure-based phylogenetic analysis shows that all poxvirus Bcl-2 proteins are structurally more similar to each other than they are to any cellular or herpesvirus Bcl-2 proteins. This is consistent with duplication and diversification of a single BCL2 family gene acquired by an ancestral poxvirus.

The thalidomide-binding domain of cereblon defines the CULT domain family and is a new member of the β-tent fold

Despite having caused one of the greatest medical catastrophies of the last century through its teratogenic side-effects, thalidomide continues to be an important agent in the treatment of leprosy and cancer. The protein cereblon, which forms an E3 ubiquitin ligase compex together with damaged DNA-binding protein 1 (DDB1) and cullin 4A, has been recently indentified as a primary target of thalidomide and its C-terminal part as responsible for binding thalidomide within a domain carrying several invariant cysteine and tryptophan residues. This domain, which we name CULT (cereblon domain of unknown activity, binding cellular ligands and thalidomide), is also found in a family of secreted proteins from animals and in a family of bacterial proteins occurring primarily in δ-proteobacteria. Its nearest relatives are yippee, a highly conserved eukaryotic protein of unknown function, and Mis18, a protein involved in the priming of centromeres for recruitment of CENP-A. Searches for distant homologs point to an evolutionary relationship of CULT, yippee, and Mis18 to proteins sharing a common fold, which consists of two four-stranded β-meanders packing at a roughly right angle and coordinating a zinc ion at their apex. A β-hairpin inserted into the first β-meander extends across the bottom of the structure towards the C-terminal edge of the second β-meander, with which it forms a cradle-shaped binding site that is topologically conserved in all members of this fold. We name this the β-tent fold for the striking arrangement of its constituent β-sheets. The fold has internal pseudosymmetry, raising the possibility that it arose by duplication of a subdomain-sized fragment.

Diversity and variability of NOD-like receptors in fungi

Nucleotide-binding oligomerization domain (NOD)-like receptors (NLRs) are intracellular receptors that control innate immunity and other biotic interactions in animals and plants. NLRs have been characterized in plant and animal lineages, but in fungi, this gene family has not been systematically described. There is however previous indications of the involvement of NLR-like genes in nonself recognition and programmed cell death in fungi. We have analyzed 198 fungal genomes for the presence of NLRs and have annotated a total of 5,616 NLR candidates. We describe their phylogenetic distribution, domain organization, and evolution. Fungal NLRs are characterized by a great diversity of domain organizations, suggesting frequently occurring combinatorial assortments of different effector, NOD and repeat domains. The repeat domains are of the WD, ANK, and TPR type; no LRR motifs were found. As previously documented for WD-repeat domains of fungal NLRs, TPR, and ANK repeats evolve under positive selection and show highly conserved repeats and repeat length polymorphism, suggesting the possibility of concerted evolution of these repeats. We identify novel effector domains not previously found associated with NLRs, whereas others are related to effector domains of plant or animals NLRs. In particular, we show that the HET domain found in fungal NLRs may be related to Toll/interleukin-1 receptor domains found in animal and plant immune receptors. This description of fungal NLR repertoires reveals both similarities and differences with plant and animals NLR collections, highlights the importance of domain reassortment and repeat evolution and provides a novel entry point to explore the evolution of NLRs in eukaryotes.

Synthesis of bacterial polysaccharides via the Wzx/Wzy-dependent pathway

The surfaces of bacteria mediate a multitude of functions in the environment and in an infected host, including adhesion to both biotic and abiotic substrata, motility, immune system interaction and (or) activation, biofilm formation, and cell–cell communication, with many of these features directly influenced by cell-surface glycans. In both Gram-negative and Gram-positive bacteria, the majority of cell-surface polysaccharides are produced via the Wzx/Wzy-dependent assembly pathway; these glycans include heteropolymeric O-antigen, enterobacterial common antigen, exopolysaccharide, spore coat, and capsule in diverse bacteria. The key components of this assembly pathway are the integral inner membrane Wzx flippase, Wzy polymerase, and Wzz chain-length regulator proteins, which until recently have resisted detailed structural and functional characterization. In this review, we have provided a comprehensive synthesis of the latest structural and mechanistic data for each protein, as well as an examination of substrate specificity for each assembly step and complex formation between the constituent proteins. To complement the unprecedented explosion of genomic-sequencing data for bacteria, we have also highlighted both classical and state-of-the-art methods by which encoded Wzx, Wzy, and Wzz proteins can be reliably identified and annotated, using the model Gram-negative bacterium Pseudomonas aeruginosa as an example data set. Lastly, we outline future avenues of research, with the aim of stimulating researchers to take the next steps in investigating the function of, and interplay between, the constituents of this widespread assembly scheme.

morFeus: a web-based program to detect remotely conserved orthologs using symmetrical best hits and orthology network scoring

Background

Searching the orthologs of a given protein or DNA sequence is one of the most important and most commonly used Bioinformatics methods in Biology. Programs like BLAST or the orthology search engine Inparanoid can be used to find orthologs when the similarity between two sequences is sufficiently high. They however fail when the level of conservation is low. The detection of remotely conserved proteins oftentimes involves sophisticated manual intervention that is difficult to automate.

Results

Here, we introduce morFeus, a search program to find remotely conserved orthologs. Based on relaxed sequence similarity searches, morFeus selects sequences based on the similarity of their alignments to the query, tests for orthology by iterative reciprocal BLAST searches and calculates a network score for the resulting network of orthologs that is a measure of orthology independent of the E-value. Detecting remotely conserved orthologs of a protein using morFeus thus requires no manual intervention. We demonstrate the performance of morFeus by comparing it to state-of-the-art orthology resources and methods. We provide an example of remotely conserved orthologs, which were experimentally shown to be functionally equivalent in the respective organisms and therefore meet the criteria of the orthology-function conjecture.

Conclusions

Based on our results, we conclude that morFeus is a powerful and specific search method for detecting remotely conserved orthologs. morFeus is freely available at http://bio.biochem.mpg.de/morfeus/. Its source code is available from Sourceforge.net (https://sourceforge.net/p/morfeus/).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-263) contains supplementary material, which is available to authorized users.

Inhibition of Mycobacterium tuberculosis PknG by non-catalytic rubredoxin domain specific modification: reaction of an electrophilic nitro-fatty acid with the Fe-S center

PknG from Mycobacterium tuberculosis is a Ser/Thr protein kinase that regulates key metabolic processes within the bacterial cell as well as signaling pathways from the infected host cell. This multidomain protein has a conserved canonical kinase domain with N- and C-terminal flanking regions of unclear functional roles. The N-terminus harbors a rubredoxin-like domain (Rbx), a bacterial protein module characterized by an iron ion coordinated by four cysteine residues. Disruption of the Rbx-metal binding site by simultaneous mutations of all the key cysteine residues significantly impairs PknG activity. This encouraged us to evaluate the effect of a nitro-fatty acid (9- and 10-nitro-octadeca-9-cis-enoic acid; OA-NO2) on PknG activity. Fatty acid nitroalkenes are electrophilic species produced during inflammation and metabolism that react with nucleophilic residues of target proteins (i.e., Cys and His), modulating protein function and subcellular distribution in a reversible manner. Here, we show that OA-NO2 inhibits kinase activity by covalently adducting PknG remote from the catalytic domain. Mass spectrometry-based analysis established that cysteines located at Rbx are the specific targets of the nitroalkene. Cys-nitroalkylation is a Michael addition reaction typically reverted by thiols. However, the reversible OA-NO2-mediated nitroalkylation of the kinase results in an irreversible inhibition of PknG. Cys adduction by OA-NO2 induced iron release from the Rbx domain, revealing a new strategy for the specific inhibition of PknG. These results affirm the relevance of the Rbx domain as a target for PknG inhibition and support that electrophilic lipid reactions of Rbx-Cys may represent a new drug strategy for specific PknG inhibition.

Comprehensive analysis of the HEPN superfamily: identification of novel roles in intra-genomic conflicts, defense, pathogenesis and RNA processing

Background

The major role of enzymatic toxins that target nucleic acids in biological conflicts at all levels has become increasingly apparent thanks in large part to the advances of comparative genomics. Typically, toxins evolve rapidly hampering the identification of these proteins by sequence analysis. Here we analyze an unexpectedly widespread superfamily of toxin domains most of which possess RNase activity.

Results

The HEPN superfamily is comprised of all α-helical domains that were first identified as being associated with DNA polymerase β-type nucleotidyltransferases in prokaryotes and animal Sacsin proteins. Using sensitive sequence and structure comparison methods, we vastly extend the HEPN superfamily by identifying numerous novel families and by detecting diverged HEPN domains in several known protein families. The new HEPN families include the RNase LS and LsoA catalytic domains, KEN domains (e.g. RNaseL and Ire1) and the RNase domains of RloC and PrrC. The majority of HEPN domains contain conserved motifs that constitute a metal-independent endoRNase active site. Some HEPN domains lacking this motif probably function as non-catalytic RNA-binding domains, such as in the case of the mannitol repressor MtlR. Our analysis shows that HEPN domains function as toxins that are shared by numerous systems implicated in intra-genomic, inter-genomic and intra-organismal conflicts across the three domains of cellular life. In prokaryotes HEPN domains are essential components of numerous toxin-antitoxin (TA) and abortive infection (Abi) systems and in addition are tightly associated with many restriction-modification (R-M) and CRISPR-Cas systems, and occasionally with other defense systems such as Pgl and Ter. We present evidence of multiple modes of action of HEPN domains in these systems, which include direct attack on viral RNAs (e.g. LsoA and RNase LS) in conjunction with other RNase domains (e.g. a novel RNase H fold domain, NamA), suicidal or dormancy-inducing attack on self RNAs (RM systems and possibly CRISPR-Cas systems), and suicidal attack coupled with direct interaction with phage components (Abi systems). These findings are compatible with the hypothesis on coupling of pathogen-targeting (immunity) and self-directed (programmed cell death and dormancy induction) responses in the evolution of robust antiviral strategies. We propose that altruistic cell suicide mediated by HEPN domains and other functionally similar RNases was essential for the evolution of kin and group selection and cell cooperation. HEPN domains were repeatedly acquired by eukaryotes and incorporated into several core functions such as endonucleolytic processing of the 5.8S-25S/28S rRNA precursor (Las1), a novel ER membrane-associated RNA degradation system (C6orf70), sensing of unprocessed transcripts at the nuclear periphery (Swt1). Multiple lines of evidence suggest that, similar to prokaryotes, HEPN proteins were recruited to antiviral, antitransposon, apoptotic systems or RNA-level response to unfolded proteins (Sacsin and KEN domains) in several groups of eukaryotes.

Conclusions

Extensive sequence and structure comparisons reveal unexpectedly broad presence of the HEPN domain in an enormous variety of defense and stress response systems across the tree of life. In addition, HEPN domains have been recruited to perform essential functions, in particular in eukaryotic rRNA processing. These findings are expected to stimulate experiments that could shed light on diverse cellular processes across the three domains of life.

Reviewers

This article was reviewed by Martijn Huynen, Igor Zhulin and Nick Grishin

Analysis of two polyhydroxyalkanoate synthases in Bradyrhizobium japonicum USDA 110

Bradyrhizobium japonicum USDA 110 has five polyhydroxyalkanoate (PHA) synthases (PhaC) annotated in its genome: bll4360 (phaC1), bll6073 (phaC2), blr3732 (phaC3), blr2885 (phaC4), and bll4548 (phaC5). All these proteins possess the catalytic triad and conserved amino acid residues of polyester synthases and are distributed into four different PhaC classes. We obtained mutants in each of these paralogs and analyzed phaC gene expression and PHA production in liquid cultures. Despite the genetic redundancy, only phaC1 and phaC2 were expressed at significant rates, while PHA accumulation in stationary-phase cultures was impaired only in the ΔphaC1 mutant. Meanwhile, the ΔphaC2 mutant produced more PHA than the wild type under this condition, and surprisingly, the phaC3 transcript increased in the ΔphaC2 background. A double mutant, the ΔphaC2 ΔphaC3 mutant, consistently accumulated less PHA than the ΔphaC2 mutant. PHA accumulation in nodule bacteroids followed a pattern similar to that seen in liquid cultures, being prevented in the ΔphaC1 mutant and increased in the ΔphaC2 mutant in relation to the level in the wild type. Therefore, we used these mutants, together with a ΔphaC1 ΔphaC2 double mutant, to study the B. japonicum PHA requirements for survival, competition for nodulation, and plant growth promotion. All mutants, as well as the wild type, survived for 60 days in a carbon-free medium, regardless of their initial PHA contents. When competing for nodulation against the wild type in a 1:1 proportion, the ΔphaC1 and ΔphaC1 ΔphaC2 mutants occupied only 13 to 15% of the nodules, while the ΔphaC2 mutant occupied 81%, suggesting that the PHA polymer is required for successful competitiveness. However, the bacteroid content of PHA did not affect the shoot dry weight accumulation.

CLCAs - a family of metalloproteases of intriguing phylogenetic distribution and with cases of substituted catalytic sites

The zinc-dependent metalloproteases with His-Glu-x-x-His (HExxH) active site motif, zincins, are a broad group of proteins involved in many metabolic and regulatory functions, and found in all forms of life. Human genome contains more than 100 genes encoding proteins with known zincin-like domains. A survey of all proteins containing the HExxH motif shows that approximately 52% of HExxH occurrences fall within known protein structural domains (as defined in the Pfam database). Domain families with majority of members possessing a conserved HExxH motif include, not surprisingly, many known and putative metalloproteases. Furthermore, several HExxH-containing protein domains thus identified can be confidently predicted to be putative peptidases of zincin fold. Thus, we predict zincin-like fold for eight uncharacterised Pfam families. Besides the domains with the HExxH motif strictly conserved, and those with sporadic occurrences, intermediate families are identified that contain some members with a conserved HExxH motif, but also many homologues with substitutions at the conserved positions. Such substitutions can be evolutionarily conserved and non-random, yet functional roles of these inactive zincins are not known. The CLCAs are a novel zincin-like protease family with many cases of substituted active sites. We show that this allegedly metazoan family has a number of bacterial and archaeal members. An extremely patchy phylogenetic distribution of CLCAs in prokaryotes and their conserved protein domain composition strongly suggests an evolutionary scenario of horizontal gene transfer (HGT) from multicellular eukaryotes to bacteria, providing an example of eukaryote-derived xenologues in bacterial genomes. Additionally, in a protein family identified here as closely homologous to CLCA, the CLCA_X (CLCA-like) family, a number of proteins is found in phages and plasmids, supporting the HGT scenario.

Bioinformatics analysis identify novel OB fold protein coding genes in C. elegans

BACKGROUND

The C. elegans genome has been extensively annotated by the WormBase consortium that uses state of the art bioinformatics pipelines, functional genomics and manual curation approaches. As a result, the identification of novel genes in silico in this model organism is becoming more challenging requiring new approaches. The Oligonucleotide-oligosaccharide binding (OB) fold is a highly divergent protein family, in which protein sequences, in spite of having the same fold, share very little sequence identity (5-25%). Therefore, evidence from sequence-based annotation may not be sufficient to identify all the members of this family. In C. elegans, the number of OB-fold proteins reported is remarkably low (n=46) compared to other evolutionary-related eukaryotes, such as yeast S. cerevisiae (n=344) or fruit fly D. melanogaster (n=84). Gene loss during evolution or differences in the level of annotation for this protein family, may explain these discrepancies.

METHODOLOGY/PRINCIPAL FINDINGS

This study examines the possibility that novel OB-fold coding genes exist in the worm. We developed a bioinformatics approach that uses the most sensitive sequence-sequence, sequence-profile and profile-profile similarity search methods followed by 3D-structure prediction as a filtering step to eliminate false positive candidate sequences. We have predicted 18 coding genes containing the OB-fold that have remarkably partially been characterized in C. elegans.

CONCLUSIONS/SIGNIFICANCE

This study raises the possibility that the annotation of highly divergent protein fold families can be improved in C. elegans. Similar strategies could be implemented for large scale analysis by the WormBase consortium when novel versions of the genome sequence of C. elegans, or other evolutionary related species are being released. This approach is of general interest to the scientific community since it can be used to annotate any genome.

Protein structure based prediction of catalytic residues

BACKGROUND

Worldwide structural genomics projects continue to release new protein structures at an unprecedented pace, so far nearly 6000, but only about 60% of these proteins have any sort of functional annotation.

RESULTS

We explored a range of features that can be used for the prediction of functional residues given a known three-dimensional structure. These features include various centrality measures of nodes in graphs of interacting residues: closeness, betweenness and page-rank centrality. We also analyzed the distance of functional amino acids to the general center of mass (GCM) of the structure, relative solvent accessibility (RSA), and the use of relative entropy as a measure of sequence conservation. From the selected features, neural networks were trained to identify catalytic residues. We found that using distance to the GCM together with amino acid type provide a good discriminant function, when combined independently with sequence conservation. Using an independent test set of 29 annotated protein structures, the method returned 411 of the initial 9262 residues as the most likely to be involved in function. The output 411 residues contain 70 of the annotated 111 catalytic residues. This represents an approximately 14-fold enrichment of catalytic residues on the entire input set (corresponding to a sensitivity of 63% and a precision of 17%), a performance competitive with that of other state-of-the-art methods.

CONCLUSIONS

We found that several of the graph based measures utilize the same underlying feature of protein structures, which can be simply and more effectively captured with the distance to GCM definition. This also has the added the advantage of simplicity and easy implementation. Meanwhile sequence conservation remains by far the most influential feature in identifying functional residues. We also found that due the rapid changes in size and composition of sequence databases, conservation calculations must be recalibrated for specific reference databases.

Dimerization of Sir3 via its C-terminal winged helix domain is essential for yeast heterochromatin formation

Gene silencing in budding yeast relies on the binding of the Silent Information Regulator (Sir) complex to chromatin, which is mediated by extensive interactions between the Sir proteins and nucleosomes. Sir3, a divergent member of the AAA+ ATPase-like family, contacts both the histone H4 tail and the nucleosome core. Here, we present the structure and function of the conserved C-terminal domain of Sir3, comprising 138 amino acids. This module adopts a variant winged helix-turn-helix (wH) architecture that exists as a stable homodimer in solution. Mutagenesis shows that the self-association mediated by this domain is essential for holo-Sir3 dimerization. Its loss impairs Sir3 loading onto nucleosomes in vitro and eliminates silencing at telomeres and HM loci in vivo. Replacing the Sir3 wH domain with an unrelated bacterial dimerization motif restores both HM and telomeric repression in sir3Δ cells. In contrast, related wH domains of archaeal and human members of the Orc1/Sir3 family are monomeric and have DNA binding activity. We speculate that a dimerization function for the wH evolved with Sir3's ability to facilitate heterochromatin formation.

Characterization of CCDC28B reveals its role in ciliogenesis and provides insight to understand its modifier effect on Bardet-Biedl syndrome

Bardet-Biedl syndrome (BBS) is a genetically heterogeneous disorder that is generally inherited in an autosomal recessive fashion. However, in some families, trans mutant alleles interact with the primary causal locus to modulate the penetrance and/or the expressivity of the phenotype. CCDC28B (MGC1203) was identified as a second site modifier of BBS encoding a protein of unknown function. Here we report the first functional characterization of this protein and show it affects ciliogenesis both in cultured cells and in vivo in zebrafish. Consistent with this biological role, our in silico analysis shows that the presence of CCDC28B homologous sequences is restricted to ciliated metazoa. Depletion of Ccdc28b in zebrafish results in defective ciliogenesis and consequently causes a number of phenotypes that are characteristic of BBS and other ciliopathy mutants including hydrocephalus, left-right axis determination defects and renal function impairment. Thus, this work reports CCDC28B as a novel protein involved in the process of ciliogenesis whilst providing functional insight into the cellular basis of its modifier effect in BBS patients.

The archaeal proteasome is regulated by a network of AAA ATPases

The proteasome is the central machinery for targeted protein degradation in archaea, Actinobacteria, and eukaryotes. In its basic form, it consists of a regulatory ATPase complex and a proteolytic core particle. The interaction between the two is governed by an HbYX motif (where Hb is a hydrophobic residue, Y is tyrosine, and X is any amino acid) at the C terminus of the ATPase subunits, which stimulates gate opening of the proteasomal α-subunits. In archaea, the proteasome-interacting motif is not only found in canonical proteasome-activating nucleotidases of the PAN/ARC/Rpt group, which are absent in major archaeal lineages, but also in proteins of the CDC48/p97/VAT and AMA groups, suggesting a regulatory network of proteasomal ATPases. Indeed, Thermoplasma acidophilum, which lacks PAN, encodes one CDC48 protein that interacts with the 20S proteasome and activates the degradation of model substrates. In contrast, Methanosarcina mazei contains seven AAA proteins, five of which, both PAN proteins, two out of three CDC48 proteins, and the AMA protein, function as proteasomal gatekeepers. The prevalent presence of multiple, distinct proteasomal ATPases in archaea thus results in a network of regulatory ATPases that may widen the substrate spectrum of proteasomal protein degradation.

Accurate multiple sequence alignment of transmembrane proteins with PSI-Coffee

BACKGROUND

Transmembrane proteins (TMPs) constitute about 20~30% of all protein coding genes. The relative lack of experimental structure has so far made it hard to develop specific alignment methods and the current state of the art (PRALINE™) only manages to recapitulate 50% of the positions in the reference alignments available from the BAliBASE2-ref7.

METHODS

We show how homology extension can be adapted and combined with a consistency based approach in order to significantly improve the multiple sequence alignment of alpha-helical TMPs. TM-Coffee is a special mode of PSI-Coffee able to efficiently align TMPs, while using a reduced reference database for homology extension.

RESULTS

Our benchmarking on BAliBASE2-ref7 alpha-helical TMPs shows a significant improvement over the most accurate methods such as MSAProbs, Kalign, PROMALS, MAFFT, ProbCons and PRALINE™. We also estimated the influence of the database used for homology extension and show that highly non-redundant UniRef databases can be used to obtain similar results at a significantly reduced computational cost over full protein databases. TM-Coffee is part of the T-Coffee package, a web server is also available from http://tcoffee.crg.cat/tmcoffee and a freeware open source code can be downloaded from http://www.tcoffee.org/Packages/Stable/Latest.

Crystal structure of a putative isochorismatase hydrolase from Oleispira antarctica

Isochorismatase-like hydrolases (IHL) constitute a large family of enzymes divided into five structural families (by SCOP). IHLs are crucial for siderophore-mediated ferric iron acquisition by cells. Knowledge of the structural characteristics of these molecules will enhance the understanding of the molecular basis of iron transport, and perhaps resolve which of the mechanisms previously proposed in the literature is the correct one. We determined the crystal structure of the apo-form of a putative isochorismatase hydrolase OaIHL (PDB code: 3LQY) from the antarctic γ-proteobacterium Oleispira antarctica, and did comparative sequential and structural analysis of its closest homologs. The characteristic features of all analyzed structures were identified and discussed. We also docked isochorismate to the determined crystal structure by in silico methods, to highlight the interactions of the active center with the substrate. The putative isochorismate hydrolase OaIHL from O. antarctica possesses the typical catalytic triad for IHL proteins. Its active center resembles those IHLs with a D-K-C catalytic triad, rather than those variants with a D-K-X triad. OaIHL shares some structural and sequential features with other members of the IHL superfamily. In silico docking results showed that despite small differences in active site composition, isochorismate binds to in the structure of OaIHL in a similar mode to its binding in phenazine biosynthesis protein PhzD (PDB code 1NF8).

Bioinformatics of the TULIP domain superfamily

Proteins of the BPI (bactericidal/permeability-increasing protein)-like family contain either one or two tandem copies of a fold that usually provides a tubular cavity for the binding of lipids. Bioinformatic analyses show that, in addition to its known members, which include BPI, LBP [LPS (lipopolysaccharide)-binding protein)], CETP (cholesteryl ester-transfer protein), PLTP (phospholipid-transfer protein) and PLUNC (palate, lung and nasal epithelium clone) protein, this family also includes other, more divergent groups containing hypothetical proteins from fungi, nematodes and deep-branching unicellular eukaryotes. More distantly, BPI-like proteins are related to a family of arthropod proteins that includes hormone-binding proteins (Takeout-like; previously described to adopt a BPI-like fold), allergens and several groups of uncharacterized proteins. At even greater evolutionary distance, BPI-like proteins are homologous with the SMP (synaptotagmin-like, mitochondrial and lipid-binding protein) domains, which are found in proteins associated with eukaryotic membrane processes. In particular, SMP domain-containing proteins of yeast form the ERMES [ER (endoplasmic reticulum)-mitochondria encounter structure], required for efficient phospholipid exchange between these organelles. This suggests that SMP domains themselves bind lipids and mediate their exchange between heterologous membranes. The most distant group of homologues we detected consists of uncharacterized animal proteins annotated as TM (transmembrane) 24. We propose to group these families together into one superfamily that we term as the TULIP (tubular lipid-binding) domain superfamily.

Protein domain of unknown function 3233 is a translocation domain of autotransporter secretory mechanism in gamma proteobacteria

Vibrio cholerae, the enteropathogenic gram negative bacteria is one of the main causative agents of waterborne diseases like cholera. About 1/3^rd of the organism's genome is uncharacterised with many protein coding genes lacking structure and functional information. These proteins form significant fraction of the genome and are crucial in understanding the organism's complete functional makeup. In this study we report the general structure and function of a family of hypothetical proteins, Domain of Unknown Function 3233 (DUF3233), which are conserved across gram negative gammaproteobacteria (especially in Vibrio sp. and similar bacteria). Profile and HMM based sequence search methods were used to screen homologues of DUF3233. The I-TASSER fold recognition method was used to build a three dimensional structural model of the domain. The structure resembles the transmembrane beta-barrel with an axial N-terminal helix and twelve antiparallel beta-strands. Using a combination of amphipathy and discrimination analysis we analysed the potential transmembrane beta-barrel forming properties of DUF3233. Sequence, structure and phylogenetic analysis of DUF3233 indicates that this gram negative bacterial hypothetical protein resembles the beta-barrel translocation unit of autotransporter Va secretory mechanism with a gene organisation that differs from the conventional Va system.

Unification of Cas protein families and a simple scenario for the origin and evolution of CRISPR-Cas systems

BACKGROUND

The CRISPR-Cas adaptive immunity systems that are present in most Archaea and many Bacteria function by incorporating fragments of alien genomes into specific genomic loci, transcribing the inserts and using the transcripts as guide RNAs to destroy the genome of the cognate virus or plasmid. This RNA interference-like immune response is mediated by numerous, diverse and rapidly evolving Cas (CRISPR-associated) proteins, several of which form the Cascade complex involved in the processing of CRISPR transcripts and cleavage of the target DNA. Comparative analysis of the Cas protein sequences and structures led to the classification of the CRISPR-Cas systems into three Types (I, II and III).

RESULTS

A detailed comparison of the available sequences and structures of Cas proteins revealed several unnoticed homologous relationships. The Repeat-Associated Mysterious Proteins (RAMPs) containing a distinct form of the RNA Recognition Motif (RRM) domain, which are major components of the CRISPR-Cas systems, were classified into three large groups, Cas5, Cas6 and Cas7. Each of these groups includes many previously uncharacterized proteins now shown to adopt the RAMP structure. Evidence is presented that large subunits contained in most of the CRISPR-Cas systems could be homologous to Cas10 proteins which contain a polymerase-like Palm domain and are predicted to be enzymatically active in Type III CRISPR-Cas systems but inactivated in Type I systems. These findings, the fact that the CRISPR polymerases, RAMPs and Cas2 all contain core RRM domains, and distinct gene arrangements in the three types of CRISPR-Cas systems together provide for a simple scenario for origin and evolution of the CRISPR-Cas machinery. Under this scenario, the CRISPR-Cas system originated in thermophilic Archaea and subsequently spread horizontally among prokaryotes.

CONCLUSIONS

Because of the extreme diversity of CRISPR-Cas systems, in-depth sequence and structure comparison continue to reveal unexpected homologous relationship among Cas proteins. Unification of Cas protein families previously considered unrelated provides for improvement in the classification of CRISPR-Cas systems and a reconstruction of their evolution.

OPEN PEER REVIEW

This article was reviewed by Malcolm White (nominated by Purficacion Lopez-Garcia), Frank Eisenhaber and Igor Zhulin. For the full reviews, see the Reviewers' Comments section.

Yeast homologues of three BLOC-1 subunits highlight KxDL proteins as conserved interactors of BLOC-1

Biogenesis of lysosome-related organelle complex-1 (BLOC-1) is one of the four multi-subunit complexes implicated in sorting cargo to lysosome-related organelles, as loss of function of any of these complexes causes Hermansky-Pudlak syndrome. Eight subunits of BLOC-1 interact with each other and with many other proteins. Identifying new interactors of BLOC-1 will increase understanding of its mechanism of action, and studies in model organisms are useful for finding such interactors. PSI-BLAST searches identify homologues in diverse model organisms, but there are significant gaps for BLOC-1, with none of its eight subunits found in Saccharomyces cerevisiae. Here we use more sensitive searches to identify distant homologues for three BLOC-1 subunits in S. cerevisiae: Blos1, snapin and cappuccino (cno). Published data on protein interactions show that in yeast these are likely to form a complex with three other proteins. One of these is the yeast homologue of the previously uncharacterized KxDL protein, which also interacts with Blos1 and cno in higher eukaryotes, suggesting that KxDL proteins are key interactors with BLOC-1.

Evolution and classification of the CRISPR-Cas systems

The CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) modules are adaptive immunity systems that are present in many archaea and bacteria. These defence systems are encoded by operons that have an extraordinarily diverse architecture and a high rate of evolution for both the cas genes and the unique spacer content. Here, we provide an updated analysis of the evolutionary relationships between CRISPR-Cas systems and Cas proteins. Three major types of CRISPR-Cas system are delineated, with a further division into several subtypes and a few chimeric variants. Given the complexity of the genomic architectures and the extremely dynamic evolution of the CRISPR-Cas systems, a unified classification of these systems should be based on multiple criteria. Accordingly, we propose a 'polythetic' classification that integrates the phylogenies of the most common cas genes, the sequence and organization of the CRISPR repeats and the architecture of the CRISPR-cas loci.

Fold and function of the InlB B-repeat

Host cell invasion by the facultative intracellular pathogen Listeria monocytogenes requires the invasion protein InlB in many cell types. InlB consists of an N-terminal internalin domain that binds the host cell receptor tyrosine kinase Met and C-terminal GW domains that bind to glycosaminoglycans (GAGs). Met binding and activation is required for host cell invasion, while the interaction between GW domains and GAGs enhances this effect. Soluble InlB elicits the same cellular phenotypes as the natural Met ligand hepatocyte growth factor/scatter factor (HGF/SF), e.g. cell scatter. So far, little is known about the central part of InlB, the B-repeat. Here we present a structural and functional characterization of the InlB B-repeat. The crystal structure reveals a variation of the β-grasp fold that is most similar to small ubiquitin-like modifiers (SUMOs). However, structural similarity also suggests a potential evolutionary relation to bacterial mucin-binding proteins. The B-repeat defines the prototype structure of a hitherto uncharacterized domain present in over a thousand bacterial proteins. Generally, this domain probably acts as a spacer or a receptor-binding domain in extracellular multi-domain proteins. In cellular assays the B-repeat acts synergistically with the internalin domain conferring to it the ability to stimulate cell motility. Thus, the B-repeat probably binds a further host cell receptor and thereby enhances signaling downstream of Met.

Cholera- and anthrax-like toxins are among several new ADP-ribosyltransferases

Chelt, a cholera-like toxin from Vibrio cholerae, and Certhrax, an anthrax-like toxin from Bacillus cereus, are among six new bacterial protein toxins we identified and characterized using in silico and cell-based techniques. We also uncovered medically relevant toxins from Mycobacterium avium and Enterococcus faecalis. We found agriculturally relevant toxins in Photorhabdus luminescens and Vibrio splendidus. These toxins belong to the ADP-ribosyltransferase family that has conserved structure despite low sequence identity. Therefore, our search for new toxins combined fold recognition with rules for filtering sequences--including a primary sequence pattern--to reduce reliance on sequence identity and identify toxins using structure. We used computers to build models and analyzed each new toxin to understand features including: structure, secretion, cell entry, activation, NAD+ substrate binding, intracellular target binding and the reaction mechanism. We confirmed activity using a yeast growth test. In this era where an expanding protein structure library complements abundant protein sequence data--and we need high-throughput validation--our approach provides insight into the newest toxin ADP-ribosyltransferases.

Crystal structure of a dimeric archaeal cleavage and polyadenylation specificity factor

Proteins of the metallo-β-lactamase (MβL) fold form a large superfamily of metallo-hydrolase/oxidoreductases. Members of this family are found in all three domains of life and are involved in a variety of biological functions related to hydrolysis, redox processes, DNA repair and uptake, and RNA processing. We classified the archaeal homologs of this superfamily based on sequence similarity and characterized a subfamily of the Cleavage and Polyadenylation Specificity Factor (CPSF) with an uncommon domain composition: in addition to an extended MβL domain, which accommodates the active site for RNA cleavage, this group has two N-terminal KH domains. Here, we present the crystal structure of a member of this group from Methanosarcina mazei. It reveals a dimerization mode of the MβL domain that has not been observed before and suggests that RNA is bound across the dimer interface, recognized by the KH domains of one monomer, and cleaved at the active site of the other.

Homology of SMP domains to the TULIP superfamily of lipid-binding proteins provides a structural basis for lipid exchange between ER and mitochondria

Mitochondria must uptake some phospholipids from the endoplasmic reticulum (ER) for the biogenesis of their membranes. They convert one of these lipids, phosphatidylserine, to phosphatidylethanolamine, which can be re-exported via the ER to all other cellular membranes. The mechanisms underlying these exchanges between ER and mitochondria are poorly understood. Recently, a complex termed ER–mitochondria encounter structure (ERMES) was shown to be necessary for phospholipid exchange in budding yeast. However, it is unclear whether this complex is merely an inter-organelle tether or also the transporter. ERMES consists of four proteins: Mdm10, Mdm34 (Mmm2), Mdm12 and Mmm1, three of which contain the uncharacterized SMP domain common to a number of eukaryotic membrane-associated proteins. Here, we show that the SMP domain belongs to the TULIP superfamily of lipid/hydrophobic ligand-binding domains comprising members of known structure. This relationship suggests that the SMP domains of the ERMES complex mediate lipid exchange between ER and mitochondria.

Contact:andrei.lupas@tuebingen.mpg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Comprehensive analysis of HAMP domains: implications for transmembrane signal transduction

Homodimeric receptors with one or two transmembrane (TM) segments per monomer are universal to life and represent the largest and most diverse group of cellular TM receptors. They frequently share domain types across phyla and, in some cases, have been recombined experimentally into functional chimeras (e.g., the bacterial aspartate chemoreceptor with the human insulin receptor), suggesting that they have a common mechanism. The nature of this mechanism, however, is still being debated. We have proposed a new model for transduction mechanism by axial helix rotation, based on the structure of a widespread domain, HAMP, that frequently occurs in direct continuation of the last TM segment, primarily in histidine kinases and chemoreceptors. Here we show by statistical analysis that HAMP domain sequences have biophysical properties compatible with the two conformations proposed by the model. The analysis also identifies three networks of coevolving residues, which allow the mechanism to subdivide into individual steps. The most extended of these networks is specific for membrane-bound HAMP domains and most likely accepts the signal from the TM helices. In a classification based on sequence clustering, these HAMPs form a central supercluster, surrounded by smaller clusters of divergent HAMPs, which typically combine into arrays of up to 31 consecutive copies and accept conformational input from other HAMP domains. Unexpectedly, the classification shows a division between domains of histidine kinases and those of chemoreceptors; thus, except for a few versatile lineages, HAMP domains are largely specific for one particular output domain. Within proteins using a given output domain, HAMP domains also show extensive coevolution with histidine kinases, but not with chemoreceptors. We attribute the greater capability for recombination among chemoreceptors to their acquisition of a reversible modification system, which acts as a capacitor for the initially deleterious effects of combining domains optimized in different contexts.

Evolution of outer membrane beta-barrels from an ancestral beta beta hairpin

Outer membrane beta-barrels (OMBBs) are the major class of outer membrane proteins from Gram-negative bacteria, mitochondria, and plastids. Their transmembrane domains consist of 8-24 beta-strands forming a closed, barrel-shaped beta-sheet around a central pore. Despite their obvious structural regularity, evidence for an origin by duplication or for a common ancestry has not been found. We use three complementary approaches to show that all OMBBs from Gram-negative bacteria evolved from a single, ancestral beta beta hairpin. First, we link almost all families of known single-chain bacterial OMBBs with each other through transitive profile searches. Second, we identify a clear repeat signature in the sequences of many OMBBs in which the repeating sequence unit coincides with the structural beta beta hairpin repeat. Third, we show that the observed sequence similarity between OMBB hairpins cannot be explained by structural or membrane constraints on their sequences. The third approach addresses a longstanding problem in protein evolution: how to distinguish between a very remotely homologous relationship and the opposing scenario of "sequence convergence." The origin of a diverse group of proteins from a single hairpin module supports the hypothesis that, around the time of transition from the RNA to the protein world, proteins arose by amplification and recombination of short peptide modules that had previously evolved as cofactors of RNAs.