A novel immunity system for bacterial nucleic acid degrading toxins and its recruitment in various eukaryotic and DNA viral systems

The use of nucleases as toxins for defense, offense or addiction of selfish elements is widely encountered across all life forms. Using sensitive sequence profile analysis methods, we characterize a novel superfamily (the SUKH superfamily) that unites a diverse group of proteins including Smi1/Knr4, PGs2, FBXO3, SKIP16, Syd, herpesviral US22, IRS1 and TRS1, and their bacterial homologs. Using contextual analysis we present evidence that the bacterial members of this superfamily are potential immunity proteins for a variety of toxin systems that also include the recently characterized contact-dependent inhibition (CDI) systems of proteobacteria. By analyzing the toxin proteins encoded in the neighborhood of the SUKH superfamily we predict that they possess domains belonging to diverse nuclease and nucleic acid deaminase families. These include at least eight distinct types of DNases belonging to HNH/EndoVII- and restriction endonuclease-fold, and RNases of the EndoU-like and colicin E3-like cytotoxic RNases-folds. The N-terminal domains of these toxins indicate that they are extruded by several distinct secretory mechanisms such as the two-partner system (shared with the CDI systems) in proteobacteria, ESAT-6/WXG-like ATP-dependent secretory systems in Gram-positive bacteria and the conventional Sec-dependent system in several bacterial lineages. The hedgehog-intein domain might also release a subset of toxic nuclease domains through auto-proteolytic action. Unlike classical colicin-like nuclease toxins, the overwhelming majority of toxin systems with the SUKH superfamily is chromosomally encoded and appears to have diversified through a recombination process combining different C-terminal nuclease domains to N-terminal secretion-related domains. Across the bacterial superkingdom these systems might participate in discriminating `self' or kin from `non-self' or non-kin strains. Using structural analysis we demonstrate that the SUKH domain possesses a versatile scaffold that can be used to bind a wide range of protein partners. In eukaryotes it appears to have been recruited as an adaptor to regulate modification of proteins by ubiquitination or polyglutamylation. Similarly, another widespread immunity protein from these toxin systems, namely the suppressor of fused (SuFu) superfamily has been recruited for comparable roles in eukaryotes. In animal DNA viruses, such as herpesviruses, poxviruses, iridoviruses and adenoviruses, the ability of the SUKH domain to bind diverse targets has been deployed to counter diverse anti-viral responses by interacting with specific host proteins.     

Polymorphic toxin systems: comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

ABSTRACT: BACKGROUND: Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflict. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis. RESULTS: Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses. CONCLUSIONS: Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.     

The RAGNYA fold: a novel fold with multiple topological variants found in functionally diverse nucleic acid, nucleotide and peptide-binding proteins

Using sensitive structure similarity searches, we identify a shared alpha+beta fold, RAGNYA, principally involved in nucleic acid, nucleotide or peptide interactions in a diverse group of proteins. These include the Ribosomal proteins L3 and L1, ATP-grasp modules, the GYF domain, DNA-recombination proteins of the NinB family from caudate bacteriophages, the C-terminal DNA-interacting domain of the Y-family DNA polymerases, the uncharacterized enzyme AMMECR1, the siRNA silencing repressor of tombusviruses, tRNA Wybutosine biosynthesis enzyme Tyw3p, DNA/RNA ligases and related nucleotidyltransferases and the Enhancer of rudimentary proteins. This fold exhibits three distinct circularly permuted versions and is composed of an internal repeat of a unit with two-strands and a helix. We show that despite considerable structural diversity in the fold, its representatives show a common mode of nucleic acid or nucleotide interaction via the exposed face of the sheet. Using this information and sensitive profile-based sequence searches: (1) we predict the active site, and mode of substrate interaction of the Wybutosine biosynthesis enzyme, Tyw3p, and a potential catalytic role for AMMECR1. (2) We provide insights regarding the mode of nucleic acid interaction of the NinB proteins, and the evolution of the active site of classical ATP-grasp enzymes and DNA/RNA ligases. (3) We also present evidence for a bacterial origin of the GYF domain and propose how this version of the fold might have been utilized in peptide interactions in the context of nucleoprotein complexes.     

The 1.14 A crystal structure of yeast cytosine deaminase: evolution of nucleotide salvage enzymes and implications for genetic chemotherapy

Cytosine deaminase (CD) catalyzes the deamination of cytosine and is only present in prokaryotes and fungi, where it is a member of the pyrimidine salvage pathway. The enzyme is of interest both for antimicrobial drug design and gene therapy applications against tumors. The structure of Saccharomyces cerevisiae CD has been determined in the presence and absence of a mechanism-based inhibitor, at 1.14 and 1.43 A resolution, respectively. The enzyme forms an alpha/beta fold similar to bacterial cytidine deaminase, but with no similarity to the alpha/beta barrel fold used by bacterial cytosine deaminase or mammalian adenosine deaminase. The structures observed for bacterial, fungal, and mammalian nucleic acid deaminases represent an example of the parallel evolution of two unique protein folds to carry out the same reaction on a diverse array of substrates.     

Genetic determination of essential residues of the Vibrio cholerae actin cross-linking domain reveals functional similarity with glutamine synthetases

Actin cross-linking domains (ACDs) are distinct domains found in several bacterial toxins, including the Vibrio cholerae MARTX toxin. The ACD of V. cholerae (ACD_Vc) catalyses the formation of an irreversible iso-peptide bond between lysine 50 and glutamic acid 270 on two actin molecules in an ATP- and Mg/Mn²⁺-dependent manner. In vivo , cross-linking depletes the cellular pool of G-actin leading to actin cytoskeleton depolymerization. While the actin cross-linking reaction performed by these effector domains has been significantly characterized, the ACD_Vc catalytic site has remained elusive due to lack of significant homology to known proteins. Using multiple genetic approaches, we have identified regions and amino acids of ACD_Vc required for full actin cross-linking activity. Then, using these functional data and structural homology predictions, it was determined that several residues demonstrated to be important for ACD_Vc activity are conserved with active-site residues of the glutamine synthetase family of enzymes. Thus, the ACDs are a family of bacterial toxin effectors that may be evolutionarily related to ligases involved in amino acid biosynthesis.     

Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous cofactor participating in numerous redox reactions. It is also a substrate for regulatory modifications of proteins and nucleic acids via the addition of ADP-ribose moieties or removal of acyl groups by transfer to ADP-ribose. In this study, we use in-depth sequence, structure and genomic context analysis to uncover new enzymes and substrate-binding proteins in NAD-utilizing metabolic and macromolecular modification systems. We predict that Escherichia coli YbiA and related families of domains from diverse bacteria, eukaryotes, large DNA viruses and single strand RNA viruses are previously unrecognized components of NAD-utilizing pathways that probably operate on ADP-ribose derivatives. Using contextual analysis we show that some of these proteins potentially act in RNA repair, where NAD is used to remove 2'-3' cyclic phosphodiester linkages. Likewise, we predict that another family of YbiA-related enzymes is likely to comprise a novel NAD-dependent ADP-ribosylation system for proteins, in conjunction with a previously unrecognized ADP-ribosyltransferase. A similar ADP-ribosyltransferase is also coupled with MACRO or ADP-ribosylglycohydrolase domain proteins in other related systems, suggesting that all these novel systems are likely to comprise pairs of ADP-ribosylation and ribosylglycohydrolase enzymes analogous to the DraG-DraT system, and a novel group of bacterial polymorphic toxins. We present evidence that some of these coupled ADP-ribosyltransferases/ribosylglycohydrolases are likely to regulate certain restriction modification enzymes in bacteria. The ADP-ribosyltransferases found in these, the bacterial polymorphic toxin and host-directed toxin systems of bacteria such as Waddlia also throw light on the evolution of this fold and the origin of eukaryotic polyADP-ribosyltransferases and NEURL4-like ARTs, which might be involved in centrosomal assembly. We also infer a novel biosynthetic pathway that might be involved in the synthesis of a nicotinate-derived compound in conjunction with an asparagine synthetase and AMPylating peptide ligase. We use the data derived from this analysis to understand the origin and early evolutionary trajectories of key NAD-utilizing enzymes and present targets for future biochemical investigations.     

Nucleotide sequence of the staphylococcal enterotoxin C3 gene sequence comparison of all three type C staphylococcal enterotoxins

Summary The structural gene entC3, which encodes staphylococcal enterotoxin C3 was cloned from the genome of Staphylococcus aureus FRI-913 and sequenced. The primary amino acid sequence of the toxin was deduced from the nucleotide sequence data. entC3 contains 801 by and encodes a precursor protein of 266 amino acids. Glutamic acid was found to be the N-terminus of mature enterotoxin C3. Thus, the first 27 residues of the toxin precursor comprise the signal peptide, and the mature toxin contains 239 amino acids with a molecular weight of 27 563 daltons. Enterotoxin C3 differs from enterotoxin C2 by four amino acids and from enterotoxin C1 by nine residues. The 167 C-terminal residues of the three toxins are identical, except for one conservative amino acid substitution in enterotoxin C3. The degree of immunological relatedness among the three Type C enterotoxins is proportional to their molecular relatedness. This study also provides evidence that the N-termini of Type C enterotoxins determine subtype-specific antigenic epitopes, while more conserved C-terminal regions determine biological properties and cross-reactive antigenic epitopes shared with other pyrogenic toxins.     

Deduced amino acid sequence of Escherichia coli adenosine deaminase reveals evolutionarily conserved amino acid residues: implications for catalytic function

The goal of the research reported here is to identify evolutionarily conserved amino acid residues associated with enzymatic deamination of adenosine. To do this, we isolated molecular clones of the Escherichia coli adenosine deaminase gene by functional complementation of adenosine deaminase deficient bacteria and deduced the amino acid sequence of the enzyme from the nucleotide sequence of the gene. Nucleotide sequence analysis revealed the presence of a 996-nucleotide open reading frame encoding a protein of 332 amino acids having a molecular weight of 36,345. The deduced amino acid sequence of the E. coli enzyme has approximately 33% identity with those of the mammalian adenosine deaminases. With conservative amino acid substitutions the overall sequence homology approaches 50%, suggesting that the structures and functions of the mammalian and bacterial enzymes are similar. Additional amino acid sequence analysis revealed specific residues that are conserved among all three adenosine deaminases and four AMP deaminases for which sequence information is currently available. In view of previously published enzymological data and the conserved amino acid residues identified in this study, we propose a model to account for the enzyme-catalyzed hydrolytic deamination of adenosine. Potential catalytic roles are assigned to the conserved His 214, Cys 262, Asp 295, and Asp 296 residues of mammalian adenosine deaminases and the corresponding conserved amino acid residues in bacterial adenosine deaminase and the eukaryotic AMP deaminases.     

Diversity of bacterial type II toxin-antitoxin systems: a comprehensive search and functional analysis of novel families

Type II toxin-antitoxin (TA) systems are generally composed of two genes organized in an operon, encoding a labile antitoxin and a stable toxin. They were first discovered on plasmids where they contribute to plasmid stability by a phenomenon denoted as 'addiction', and subsequently in bacterial chromosomes. To discover novel families of antitoxins and toxins, we developed a bioinformatics approach based on the 'guilt by association' principle. Extensive experimental validation in Escherichia coli of predicted antitoxins and toxins increased significantly the number of validated systems and defined novel toxin and antitoxin families. Our data suggest that toxin families as well as antitoxin families originate from distinct ancestors that were assembled multiple times during evolution. Toxin and antitoxin families found on plasmids tend to be promiscuous and widespread, indicating that TA systems move through horizontal gene transfer. We propose that due to their addictive properties, TA systems are likely to be maintained in chromosomes even though they do not necessarily confer an advantage to their bacterial hosts. Therefore, addiction might play a major role in the evolutionary success of TA systems both on mobile genetic elements and in bacterial chromosomes.     

Exchange of a single amino acid switches the substrate properties of RhoA and RhoD toward glucosylating and transglutaminating toxins

Rho GTPases are the preferred targets of various bacterial cytotoxins, including Clostridium difficile toxins A and B, Clostridium sordellii lethal toxin, the cytotoxic necrotizing factors (CNF1) from Escherichia coli, and the dermonecrotizing toxin (DNT) from Bordetella species. The toxins inactivate or activate specific sets of Rho GTPases by mono-O-glucosylation and deamidation/transglutamination, respectively. Here we studied the structural basis of the recognition of RhoA, which is modified by toxin B, CNF1, and DNT, in comparison with RhoD, which is solely a substrate for lethal toxin. We found that a single amino acid residue in RhoA and RhoD defines the substrate specificity for toxin B and lethal toxin. Change of serine 73 to phenylalanine in RhoA turned RhoA into a substrate for lethal toxin. Accordingly, change of the equivalently positioned phenylalanine 85 in RhoD with serine allowed glucosylation by toxin B. Comparable results were achieved with the Rho-activating and transglutaminating enzymes CNF1 and DNT. Here, amino acid glutamate 64 of RhoA and the equivalent aspartate 76 of RhoD define substrate specificity for CNF1 and DNT, respectively. These data indicate that single amino acid residues located in the switch II region of Rho proteins determine enzyme specificity for diverse bacterial toxins.     

Differential effects of cholera toxin on guanine nucleotide regulation of beta-adrenergic agonist high affinity binding and adenylate cyclase activation in frog erythrocyte membranes

The guanine nucleotide regulatory protein(s) regulates both adenylate cyclase activity and the affinity of adenylate cyclase-coupled receptors for hormones or agonist drugs. Cholera toxin catalyzes the covalent modification of the nucleotide regulatory protein of adenylate cyclase systems. Incubation of frog erythrocyte membranes with cholera toxin and NAD+ did not substantially alter the dose dependency for guanine nucleotide activation of adenylate cyclase activity. In contrast, toxin treated membranes demonstrated a 10 fold increase in the concentrations of guanine nucleotide required for a half maximal effect in regulating beta-adrenergic receptor affinity for the agonist (+/-) [3H]hydroxybenzylisoproterenol. The data emphasize the bifunctional nature of the guanine nucleotide regulatory protein and suggest that distinct structural domains of the guanine nucleotide regulatory protein may mediate the distinct regulatory effects on adenylate cyclase and receptor affinity for agonists.     

Functional diversification of the RING finger and other binuclear treble clef domains in prokaryotes and the early evolution of the ubiquitin system

Recent studies point to a diverse assemblage of prokaryotic cognates of the eukaryotic ubiquitin (Ub) system. These systems span an entire spectrum, ranging from those catalyzing cofactor and amino acid biosynthesis, with only adenylating E1-like enzymes and ubiquitin-like proteins (Ubls), to those that are closer to eukaryotic systems by virtue of possessing E2 enzymes. Until recently E3 enzymes were unknown in such prokaryotic systems. Using contextual information from comparative genomics, we uncover a diverse group of RING finger E3s in prokaryotes that are likely to function with E1s, E2s, JAB domain peptidases and Ubls. These E1s, E2s and RING fingers suggest that features hitherto believed to be unique to eukaryotic versions of these proteins emerged progressively in such prokaryotic systems. These include the specific configuration of residues associated with oxyanion-hole formation in E2s and the C-terminal UFD in the E1 enzyme, which presents the E2 to its active site. Our study suggests for the first time that YukD-like Ubls might be conjugated by some of these systems in a manner similar to eukaryotic Ubls. We also show that prokaryotic RING fingers possess considerable functional diversity and that not all of them are involved in Ub-related functions. In eukaryotes, other than RING fingers, a number of distinct binuclear (chelating two Zn atoms) and mononuclear (chelating one zinc atom) treble clef domains are involved in Ub-related functions. Through detailed structural analysis we delineated the higher order relationships and interaction modes of binuclear treble clef domains. This indicated that the FYVE domain acquired the binuclear state independently of the other binuclear forms and that different treble clef domains have convergently acquired Ub-related functions independently of the RING finger. Among these, we uncover evidence for notable prokaryotic radiations of the ZF-UBP, B-box, AN1 and LIM clades of treble clef domains and present contextual evidence to support their role in functions unrelated to the Ub-system in prokaryotes. In particular, we show that bacterial ZF-UBP domains are part of a novel cyclic nucleotide-dependent redox signaling system, whereas prokaryotic B-box, AN1 and LIM domains have related functions as partners of diverse membrane-associated peptidases in processing proteins. This information, in conjunction with structural analysis, suggests that these treble clef domains might have been independently recruited to the eukaryotic Ub-system due to an ancient conserved mode of interaction with peptides.     

Molecular weaponry: diverse effectors delivered by the Type VI secretion system

The Type VI secretion system is a widespread bacterial nanomachine, used to deliver toxins directly into eukaryotic or prokaryotic target cells. These secreted toxins, or effectors, act on diverse cellular targets, and their action provides the attacking bacterial cell with a significant fitness advantage, either against rival bacteria or eukaryotic host organisms. In this review, we discuss the delivery of diverse effectors by the Type VI secretion system, the modes of action of the so‐called ‘anti‐bacterial’ and ‘anti‐eukaryotic’ effectors, the mechanism of self‐resistance against anti‐bacterial effectors and the evolutionary implications of horizontal transfer of Type VI secretion system‐associated toxins. Whilst it is likely that many more effectors remain to be identified, it is already clear that toxins delivered by this secretion system represent efficient weapons against both bacteria and eukaryotes.     

Bovine interferon alpha genes. Structure and expression

The bovine genome contains a gene family of interferon-alpha s (bIFN-alpha) that consists of at least five distinct members. Four of the bIFN-alpha genes isolated show a high degree of homology (97% in the nucleotide sequence and 93% in amino acid sequence). The overall homology in amino acid sequence of bIFN-alpha to human, murine, and rat IFNs-alpha is approximately 60%. Yet there are amino acid clusters (positions 28-41 and 118-146) which are highly conserved throughout the mammalian evolution and in which the overall homology can be as high as 86%. Within the C terminus conserved cluster there is a sequence containing 9 amino acids completely conserved in 16 mammalian IFNs-alpha and of these, 7 are also shared with a similar domain in some bacterial toxins, implying a common functional role for these domains. One of the genes, IFN-alpha C, was expressed in Escherichia coli. The purified bacterial IFN (specific activity, 2 X 10(8) units/mg) exhibited antiviral activity on bovine cells but no detectable activity was demonstrated on human and simian cells.     

Relative character-state space, amount of potential phylogenetic information, and heterogeneity of nucleotide and amino acid characters

We examined a broad selection of protein-coding loci from a diverse array of clades and genomes to quantify three factors that determine whether nucleotide or amino acid characters should be preferred for phylogenetic inference. First, we quantified the difference in observed character-state space between nucleotides and amino acids. Second, we quantified the loss of potential phylogenetic signal from silent substitutions when amino acids are used. Third, we used the disparity index to quantify the relative compositional heterogeneity of nucleotides and amino acids and then determined how commonly convergent (rather than unique) shifts in nucleotide and amino acid composition occur in a phylogenetic context. The greater potential phylogenetic signal for nucleotide characters was found to be enormous (on average 440% that of amino acids), whereas the greater observed character-state space for amino acids was less impressive (on average 150.4% that of nucleotides). While matrices of amino acid sequences had less compositional heterogeneity than their corresponding nucleotide sequences, heterogeneity in amino acid composition may be more homoplasious than heterogeneity in nucleotide composition. Given the ability of increased taxon sampling to better utilize the greater potential phylogenetic signal of nucleotide characters and decrease the potential for artifacts caused by heterogeneous nucleotide composition among taxa, we suggest that increased taxon sampling be performed whenever possible instead of restricting analyses to amino acid characters.     

Comparative effects of cholera and Bordetella pertussis toxins on cyclic AMP and GTP levels and on lipolysis in rat adipocytes incubated in vitro

The respective effects of cholera and Bordetella pertussis toxins were studied in time and concentration dependent experiments, following glycerol and fatty acid release, GTP and cAMP levels. Cholera toxin, after a lag time of 30 min, stimulated linearly GTP and cAMP accumulation and lipolysis (maximal effect: 2-fold increase at 5 micrograms/ml). Pertussis toxin presented a biphasic effect both in time and concentration dependent studies. Up to a maximum reached after 2 h with 1.4 units LPF/ml the stimulation affected GTP (3 fold) and cAMP (7 fold) levels, glycerol and fatty acid release (15 fold). Beyond this, an inhibition occurred, yielding a decrease towards basal values of GTP and cAMP content whereas the glycerol and fatty acid release was stopped. These results, which are the first reporting the fluctuation of the GTP content of intact cells challenged with bacterial toxins, show a close relationship between GTP and cyclic AMP levels and lipolytic activity.     

Modeling nucleic acids

Nucleic acids are an important class of biological macromolecules that carry out a variety of cellular roles. For many functions, naturally occurring DNA and RNA molecules need to fold into precise three-dimensional structures. Due to their self-assembling characteristics, nucleic acids have also been widely studied in the field of nanotechnology, and a diverse range of intricate three-dimensional nanostructures have been designed and synthesized. Different physical terms such as base-pairing and stacking interactions, tertiary contacts, electrostatic interactions and entropy all affect nucleic acid folding and structure. Here we review general computational approaches developed to model nucleic acid systems. We focus on four key areas of nucleic acid modeling: molecular representation, potential energy function, degrees of freedom and sampling algorithm. Appropriate choices in each of these key areas in nucleic acid modeling can effectively combine to aid interpretation of experimental data and facilitate prediction of nucleic acid structure.     

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.