Crystal structure of BinB: A receptor binding component of the binary toxin from Lysinibacillus sphaericus

The binary toxin (Bin), produced by Lysinibacillus sphaericus, is composed of BinA (42 kDa) and BinB (51 kDa) proteins, which are both required for full toxicity against Culex and Anopheles mosquito larvae. Specificity of Bin toxin is determined by the binding of BinB component to a receptor present on the midgut epithelial membranes, while BinA is proposed to be a toxic component. Here, we determined the first crystal structure of the active form of BinB at a resolution of 1.75 Å. BinB possesses two distinct structural domains in its N‐ and C‐termini. The globular N‐terminal domain has a β‐trefoil scaffold which is a highly conserved architecture of some sugar binding proteins or lectins, suggesting a role of this domain in receptor‐binding. The BinB β‐rich C‐terminal domain shares similar three‐dimensional folding with aerolysin type β‐pore forming toxins, despite a low sequence identity. The BinB structure, therefore, is a new member of the aerolysin‐like toxin family, with probably similarities in the cytolytic mechanism that takes place via pore formation. Proteins 2014; 82:2703–2712. © 2014 Wiley Periodicals, Inc.     

Thermodynamic properties of the effector domains of MARTX toxins suggest their unfolding for translocation across the host membrane

MARTX (multifunctional autoprocessing repeats‐in‐toxin) family toxins are produced by Vibrio cholerae, Vibrio vulnificus, Aeromonas hydrophila and other Gram‐negative bacteria. Effector domains of MARTX toxins cross the cytoplasmic membrane of a host cell through a putative pore formed by the toxin's glycine‐rich repeats. The structure of the pore is unknown and the translocation mechanism of the effector domains is poorly understood. We examined the thermodynamic stability of the effector domains of V. cholerae and A. hydrophila MARTX toxins to elucidate the mechanism of their translocation. We found that all but one domain in each toxin are thermodynamically unstable and several acquire a molten globule state near human physiological temperatures. Fusion of the most stable cysteine protease domain to the adjacent effector domain reduces its thermodynamic stability ～ 1.4‐fold (from 21.8 to 16.1 kJ mol^?1). Precipitation of several individual domains due to thermal denaturation is reduced upon their fusion into multi‐domain constructs. We speculate that low thermostability of the MARTX effector domains correlates with that of many other membrane‐penetrating toxins and implies their unfolding for cell entry. This study extends the list of thermolabile bacterial toxins, suggesting that this quality is essential and could be susceptible for selective targeting of pathogenic toxins.     

Nuclear Import of Chromatin Remodeler Isw1 Is Mediated by Atypical Bipartite cNLS and Classical Import Pathway

The protein Isw1 of Saccharomyces cerevisiae is an imitation‐switch chromatin‐remodeling factor. We studied the mechanisms of its nuclear import and found that the nuclear localization signal (NLS) mediating the transport of Isw1 into the nucleus is located at the end of the C‐terminus of the protein (aa1079–1105). We show that it is an atypical bipartite signal with an unconventional linker of 19 aa (KRIR X₁₉ KKAK) and the only nuclear targeting signal within the Isw1 molecule. The efficiency of Isw1 nuclear import was found to be modulated by changes to the amino acid composition in the vicinity of the KRIR motif, but not by the linker length. Live‐cell imaging of various karyopherin mutants and in vitro binding assays of Isw1NLS to importin‐α revealed that the nuclear translocation of Isw1 is mediated by the classical import pathway. Analogous motifs to Isw1NLS are highly conserved in Isw1 homologues of other yeast species, and putative bipartite cNLS were identified in silico at the end of the C‐termini of imitation switch (ISWI) proteins from higher eukaryotes. We suggest that the C‐termini of the ISWI family proteins play an important role in their nuclear import.     

Crystal structure of the variable domain of the Streptococcus gordonii surface protein SspB

The Antigen I/II (AgI/II) family of proteins are cell wall anchored adhesins expressed on the surface of oral streptococci. The AgI/II proteins interact with molecules on other bacteria, on the surface of host cells, and with salivary proteins. Streptococcus gordonii is a commensal bacterium, and one of the primary colonizers that initiate the formation of the oral biofilm. S. gordonii expresses two AgI/II proteins, SspA and SspB that are closely related. One of the domains of SspB, called the variable (V‐) domain, is significantly different from corresponding domains in SspA and all other AgI/II proteins. As a first step to elucidate the differences among these proteins, we have determined the crystal structure of the V‐domain from S. gordonii SspB at 2.3 Å resolution. The domain comprises a β‐supersandwich with a putative binding cleft stabilized by a metal ion. The overall structure of the SspB V‐domain is similar to the previously reported V‐domain of the Streptococcus mutans protein SpaP, despite their low sequence similarity. In spite of the conserved architecture of the binding cleft, the cavity is significantly smaller in SspB, which may provide clues about the difference in ligand specificity. We also verified that the metal in the binding cleft is a calcium ion, in concurrence with previous biological data. It was previously suggested that AgI/II V‐domains are carbohydrate binding. However, we tested that hypothesis by screening the SspB V‐domain for binding to over 400 glycoconjucates and found that the domain does not interact with any of the carbohydrates.     

Model peptide studies of sequence regions in the elastomeric biomineralization protein, Lustrin A. I. The C-domain consensus-PG-, -NVNCT-motif

The lustrin superfamily represents a unique group of biomineralization proteins localized between layered aragonite mineral plates (i.e., nacre layer) in mollusk shell. Recent atomic force microscopy (AFM) pulling studies have demonstrated that the lustrin‐containing organic nacre layer in the abalone, Haliotis rufescens, exhibits a typical sawtooth force‐extension curve with hysteretic recovery. This force extension behavior is reminiscent of reversible unfolding and refolding in elastomeric proteins such as titin and tenascin. Since secondary structure plays an important role in force‐induced protein unfolding and refolding, the question is, What secondary structure(s) exist within the major domains of Lustrin A? Using a model peptide (FPGKNVNCTSGE) representing the 12‐residue consensus sequence found near the N‐termini of the first eight cysteine‐rich domains (C‐domains) within the Lustrin A protein, we employed CD, NMR spectroscopy, and simulated annealing/minimization to determine the secondary structure preferences for this sequence. At pH 7.4, we find that the 12‐mer sequence adopts a loop conformation, consisting of a “bend” or “turn” involving residues G3–K4 and N7–C8–T9, with extended conformations arising at F1–G3; K4–V6; T9–S10–G11 in the sequence. Minor pH‐dependent conformational effects were noted for this peptide; however, there is no evidence for a salt‐bridge interaction between the K4 and E12 side chains. The presence of a loop conformation within the highly conserved —PG—, —NVNCT— sequence of C1–C8 domains may have important structural and mechanistic implications for the Lustrin A protein with regard to elastic behavior. © 2002 Wiley Periodicals, Inc. Biopolymers 63: 358–369, 2002     

The membrane localization domains of two distinct bacterial toxins form a 4‐helix‐bundle in solution

Membrane localization domain (MLD) was first proposed for a 4‐helix‐bundle motif in the crystal structure of the C1 domain of Pasteurella multocida toxin (PMT). This structure motif is also found in the crystal structures of several clostridial glycosylating toxins (TcdA, TcdB, TcsL, and TcnA). The Ras/Rap1‐specific endopeptidase (RRSP) module of the multifunctional autoprocessing repeats‐in‐toxins (MARTX) toxin produced by Vibrio vulnificus has sequence homology to the C1‐C2 domains of PMT, including a putative MLD. We have determined the solution structure for the MLDs in PMT and in RRSP using solution state NMR. We conclude that the MLDs in these two toxins assume a 4‐helix‐bundle structure in solution.     

Induced autoprocessing of the cytopathic Makes caterpillars floppy‐like effector domain of the Vibrio vulnificus MARTX toxin

The multifunctional‐autoprocessing repeats‐in‐toxin (MARTX_Vv) toxin that harbours a varied repertoire of effector domains is the primary virulence factor of Vibrio vulnificus. Although ubiquitously present among Biotype I toxin variants, the ‘Makes caterpillars floppy‐like’ effector domain (MCF_Vv) is previously unstudied. Using transient expression and protein delivery, MCF_Vv and MCF_Ah from the Aeromonas hydrophila MARTX_Ah toxin are shown for the first time to induce cell rounding. Alanine mutagenesis across the C‐terminal subdomain of MCF_Vv identified an Arg‐Cys‐Asp (RCD) tripeptide motif shown to comprise a cysteine protease catalytic site essential for autoprocessing of MCF_Vv. The autoprocessing could be recapitulated in vitro by the addition of host cell lysate to recombinant MCF_Vv, indicating induced autoprocessing by cellular factors. The RCD motif is also essential for cytopathicity, suggesting autoprocessing is essential first to activate the toxin and then to process a cellular target protein resulting in cell rounding. Sequence homology places MCF_Vv within the C58 cysteine protease family that includes the type III secretion effectors YopT from Yersinia spp. and AvrPphB from Pseudomonas syringae. However, the catalytic site RCD motif is unique compared with other C58 peptidases and is here proposed to represent a new subgroup of autopeptidase found within a number of putative large bacterial toxins.     

Autoprocessing of the Vibrio cholerae RTX toxin by the cysteine protease domain

Vibrio cholerae RTX is a large multifunctional bacterial toxin that causes actin crosslinking. Due to its size, it was predicted to undergo proteolytic cleavage during translocation into host cells to deliver activity domains to the cytosol. In this study, we identified a domain within the RTX toxin that is conserved in large clostridial glucosylating toxins TcdB, TcdA, TcnA, and TcsL; putative toxins from V. vulnificus, Yersinia sp., Photorhabdus sp., and Xenorhabdus sp.; and a filamentous/hemagglutinin-like protein FhaL from Bordetella sp. In vivo transfection studies and in vitro characterization of purified recombinant protein revealed that this domain from the V. cholerae RTX toxin is an autoprocessing cysteine protease whose activity is stimulated by the intracellular environment. A cysteine point mutation within the RTX holotoxin attenuated actin crosslinking activity suggesting that processing of the toxin is an important step in toxin translocation. Overall, we have uncovered a new mechanism by which large bacterial toxins and proteins deliver catalytic activities to the eukaryotic cell cytosol by autoprocessing after translocation.     

Identification of elicitor‐responsive proteins in rice leaves by a proteomic approach

Experiments were conducted to identify the differentially expressed proteins in rice (Oryza sativa L.) plants after treatment with the glycoprotein elicitor CSB I, purified from ZC₁₃, a race of the rice blast fungus Magnaporthe grisea. The interactions of two near isogenic lines of rice, C101A51 and CO39, with ZC₁₃ resulted in completely incompatible and compatible types, respectively. Proteins were extracted from rice leaves at 12 and 24 h after treatment with CSB I. Temporal changes in total proteins were examined using 2‐DE. Among more than 900 protein spots reproducibly detected on each gel, 11 were up‐regulated, three were down‐regulated and seven were newly induced during, at a minimum, one time point. Twenty‐one differentially expressed proteins were identified by linear ion trap quadrupole (LTQ)‐MS/MS. The identified proteins were classified into six categories based on their putative function reported: (i) defense proteins (PR‐10a, PR‐5 and putative salt‐induced protein), (ii) signal transduction (nucleoside diphosphate kinase and putative profilin), (iii) ROS (Mn‐SOD, Cu/Zn‐SOD, GST and CAT), (iv) programmed cell death (translationally controlled tumor protein), (v) molecule biosynthesis (putative ribosomal protein S5, putative ribosomal protein L12, putative translational elongation factor Tu and putative chaperonin 21 precursor) and (vi) metabolism (putative fructose‐bisphosphate aldolase class‐I, putative malate dehydrogenase, cytoplasmic malate dehydrogenase, putative acid phosphatase, putative transketolase1 and gamma hydroxybutyrate dehydrogenase‐like protein). All of these proteins (except Cu/Zn‐SOD, putative acid phosphatase and translationally controlled tumor protein) were induced faster and to a higher degree in C101A51 than in CO39. These data suggest that the incompatible rice line may possess a more sensitive recognition system that can identify and react to specific chemical, biological or physical triggers in a more efficient manner, thus eliciting an early and fast defense response.     

Characterisation of a novel gene family of putative cyclic nucleotide- and calmodulin-regulated ion channels in Arabidopsis thaliana

In plants, cyclic GMP is involved in signal transduction in response to light and gibberellic acid. For cyclic AMP, a potential role during the plant cell cycle was recently reported. However, cellular targets for cyclic nucleotides in plants are largely unknown. Here we report on the identification and characterisation of a new gene family in Arabidopsis, which share features with cyclic nucleotide-gated channels from animals and inward-rectifying K⁺ channels from plants. The identified gene family comprises six members (Arabidopsis thaliana cyclic nucleotide-gated channels, AtCNGC1–6) with significant homology among the deduced proteins. Hydrophobicity analysis predicted six membrane-spanning domains flanked by hydrophilic amino and carboxy termini. A putative cyclic nucleotide binding domain (CNBD) which contains several residues that are invariant in other CNBDs was located in the carboxy terminus. This domain overlaps with a predicted calmodulin (CaM) binding site, suggesting interaction between cyclic nucleotide and CaM regulation. We demonstrated interaction of the carboxy termini of AtCNGC1 and AtCNGC2 with CaM in yeast, indicating that the CaM binding sites are functional. Furthermore, it was shown that both AtCNGC1 and AtCNGC2 can partly complement the K⁺-uptake-deficient yeast mutant CY162. Therefore, we propose that the identified genes constitute a family of plant cyclic nucleotide- and CaM-regulated ion channels.     

Surface plasmon resonance analysis of the binding of high‐risk mucosal HPV E6 oncoproteins to the PDZ1 domain of the tight junction protein MAGI‐1

The E6 oncoproteins from high‐risk mucosal human papillomavirus (HPV) induce cervical cancer via two major activities, the binding and the degradation of the p53 protein and PDZ domain‐containing proteins. Human MAGI‐1 is a multi‐PDZ domain protein implicated into protein complex assembly at cell–cell contacts. High‐risk mucosal HPV E6 proteins interact with the PDZ1 domain of MAGI‐1 via a C‐terminal consensus binding motif. Here, we developed a medium throughput protocol to accurately measure by surface plasmon resonance affinity constants of protein domains binding to peptidic sequences produced as recombinant fusions to the glutathione‐S‐transferase (GST). This approach was applied to measure the binding of MAGI‐1 PDZ1 to the C‐termini of viral or cellular proteins. Both high‐risk mucosal HPV E6 C‐terminal peptides and cellular partners of MAGI‐1 PDZ1 bind to MAGI‐1 PDZ1 with comparable dissociation constants in the micromolar range. MAGI‐1 PDZ1 shows a preference for C‐termini with a valine at position 0 and a negative charge at position ?3, confirming previous studies performed with HPV18 E6. A detailed combined analysis via site‐directed mutagenesis of the HPV16 C‐terminal peptide and PDZ1 indicated that interactions mediated by charged residues upstream the PDZ‐binding motif strongly contribute to binding selectivity of this interaction. In addition, our work highlighted the K₄₉₉ residue of MAGI‐1 as a novel determinant of binding specificity. Finally, we showed that MAGI‐1 PDZ1 also binds to the C‐termini of LPP and Tax proteins, which were already known to bind to PDZ proteins but not to MAGI‐1. Copyright © 2010 John Wiley & Sons, Ltd.     

Identification of the amino acid residues and domains in the cysteine‐rich protein of Chinese wheat mosaic virus that are important for RNA silencing suppression and subcellular localization

Cysteine‐rich proteins (CRPs) encoded by some plant viruses in diverse genera function as RNA silencing suppressors. Within the N‐terminal portion of CRPs encoded by furoviruses, there are six conserved cysteine residues and a Cys–Gly–X–X–His motif (Cys, cysteine; Gly, glycine; His, histidine; X, any amino acid residue) with unknown function. The central domains contain coiled‐coil heptad amino acid repeats that usually mediate protein dimerization. Here, we present evidence that the conserved cysteine residues and Cys–Gly–X–X–His motif in the CRP of Chinese wheat mosaic virus (CWMV) are critical for protein stability and silencing suppression activity. Mutation of a leucine residue in the third coiled‐coil heptad impaired CWMV CRP activity for suppression of local silencing, but not for the promotion of cell‐to‐cell movement of Potato virus X (PVX). In planta and in vitro analysis of wild‐type and mutant proteins indicated that the ability of the CRP to self‐interact was correlated with its suppression activity. Deletion of up to 40 amino acids at the C‐terminus did not abolish suppression activity, but disrupted the association of CRP with endoplasmic reticulum (ER), and reduced its activity in the enhancement of PVX symptom severity. Interestingly, a short region in the C‐terminal domain, predicted to form an amphipathic α‐helical structure, was responsible for the association of CWMV CRP with ER. Overall, our results demonstrate that the N‐terminal and central regions are the functional domains for suppression activity, whereas the C‐terminal region primarily functions to target CWMV CRP to the ER.     

Research note: Analysis of expressed sequence tags from the chrysophycean alga Ochromonas danica (Heterokontophyta)

A cDNA library was constructed from the chrysophycean alga, Ochromonas danica E. G. Pringsheim. 5′‐end sequencing of about 600 cDNA clones yielded 476 authentic expressed sequence tags (EST) of which 275 showed significant matches (E‐value <10^?4) to sequences in a public database. The annotation of these ESTs was carried out to assess subcellular localization of the putative proteins using several internet‐accessible prediction programs for subcellular localization. These analyses revealed that putative plastid proteins in Ochromonas possess N‐terminal bipartite presequences with a conserved phenylalanine at the N‐terminus of the predicted transit peptide‐like domains, similar to other ‘red‐lineage’ secondary symbiotic organisms. The examination of sequences of 3′‐UTR revealed that, similarly to chlorophyte algae, UGUAA may represent a putative polyadenylation signal in O. danica.     

A binding motif for Hsp90 in the A chains of ADP‐ribosylating toxins that move from the endoplasmic reticulum to the cytosol

Cholera toxin (Ctx) is an AB‐type protein toxin that acts as an adenosine diphosphate (ADP)‐ribosyltransferase to disrupt intracellular signalling in the target cell. It moves by vesicle carriers from the cell surface to the endoplasmic reticulum (ER) of an intoxicated cell. The catalytic CtxA1 subunit then dissociates from the rest of the toxin, unfolds, and activates the ER‐associated degradation system for export to the cytosol. Translocation occurs through an unusual ratchet mechanism in which the cytosolic chaperone Hsp90 couples CtxA1 refolding with CtxA1 extraction from the ER. Here, we report that Hsp90 recognises two peptide sequences from CtxA1: an N‐terminal RPPDEI sequence (residues 11–16) and an LDIAPA sequence in the C‐terminal region (residues 153–158) of the 192 amino acid protein. Peptides containing either sequence effectively blocked Hsp90 binding to full‐length CtxA1. Both sequences were necessary for the ER‐to‐cytosol export of CtxA1. Mutagenesis studies further demonstrated that the RPP residues in the RPPDEI motif are required for CtxA1 translocation to the cytosol. The LDIAPA sequence is unique to CtxA1, but we identified an RPPDEI‐like motif at the N‐ or C‐termini of the A chains from four other ER‐translocating toxins that act as ADP‐ribosyltransferases: pertussis toxin, Escherichia coli heat‐labile toxin, Pseudomonas aeruginosa exotoxin A, and Salmonella enterica serovar Typhimurium ADP‐ribosylating toxin. Hsp90 plays a functional role in the intoxication process for most, if not all, of these toxins. Our work has established a defined RPPDEI binding motif for Hsp90 that is required for the ER‐to‐cytosol export of CtxA1 and possibly other toxin A chains as well.     

Disulfide bonds and disorder in granulin‐3: An unusual handshake between structural stability and plasticity

Granulins (GRNs) are a family of small (～6 kDa) proteins generated by the proteolytic processing of their precursor, progranulin (PGRN), in many cell types. Both PGRN and GRNs are implicated in a plethora of biological functions, often in opposing roles to each other. Lately, GRNs have generated significant attention due to their implicated roles in neurodegenerative disorders. Despite their physiological and pathological significance, the structure‐function relationships of GRNs are poorly defined. GRNs contain 12 conserved cysteines forming six intramolecular disulfide bonds, making them rather exceptional, even among a few proteins with high disulfide bond density. Solution NMR investigations in the past have revealed a unique structure containing putative interdigitated disulfide bonds for several GRNs, but GRN‐3 was unsolvable due to its heterogeneity and disorder. In our previous report, we showed that abrogation of disulfide bonds in GRN‐3 renders the protein completely disordered (Ghag et al., Prot Eng Des Sel 2016). In this study, we report the cellular expression and biophysical analysis of fully oxidized, native GRN‐3. Our results indicate that both E. coli and human embryonic kidney (HEK) cells do not exclusively make GRN‐3 with homogenous disulfide bonds, likely due to the high cysteine density within the protein. Biophysical analysis suggests that GRN‐3 structure is dominated by irregular loops held together only by disulfide bonds, which induced remarkable thermal stability to the protein despite the lack of regular secondary structure. This unusual handshake between disulfide bonds and disorder within GRN‐3 could suggest a unique adaptation of intrinsically disordered proteins towards structural stability.     

Role of AIP56 disulphide bond and its reduction by cytosolic redox systems for efficient intoxication

Apoptosis‐inducing protein of 56 kDa (AIP56) is a major virulence factor of Photobacterium damselae subsp. piscicida, a gram‐negative pathogen that infects warm water fish species worldwide and causes serious economic losses in aquacultures. AIP56 is a single‐chain AB toxin composed by two domains connected by an unstructured linker peptide flanked by two cysteine residues that form a disulphide bond. The A domain comprises a zinc‐metalloprotease moiety that cleaves the NF‐kB p65, and the B domain is involved in binding and internalisation of the toxin into susceptible cells. Previous experiments suggested that disruption of AIP56 disulphide bond partially compromised toxicity, but conclusive evidences supporting the importance of that bond in intoxication were lacking. Here, we show that although the disulphide bond of AIP56 is dispensable for receptor recognition, endocytosis, and membrane interaction, it needs to be intact for efficient translocation of the toxin into the cytosol. We also show that the host cell thioredoxin reductase‐thioredoxin system is involved in AIP56 intoxication by reducing the disulphide bond of the toxin at the cytosol. The present study contributes to a better understanding of the molecular mechanisms operating during AIP56 intoxication and reveals common features shared with other AB toxins.     

Identification and characterization of two chitin-binding proteins from the peritrophic membrane of the silkworm, Bombyx mori L

The peritrophic membrane (PM) is a semi‐permeable lining of the insect midgut, broadly analogous to the mucous lining of vertebrate gut. The PM proteins are important achievements for the function of the PM. In this study, two chitin‐binding proteins (BmPM‐P43 and BmPM‐P41) from the PM of the silkworm, Bombyx mori, were identified and cloned. These proteins showed the molecular mass of 43 and 41 kDa, respectively. The deduced amino acid sequences codes for a protein of 381 amino acid residues and 364 amino acid residues, containing 12 and 14 cysteine residues followed by similar domain, both of them have 5 cysteine residues in similar position in the C‐terminal. The confirmation of these proteins was performed by western blot analysis of recombinant BmPM‐P43 and BmPM‐P41. The chitin‐binding activity analysis showed that the BmPM‐P43 and BmPM‐P41 could bind to chitin strongly. It is concluded that BmPM‐P43 and BmPM‐P41 contains a polysaccharide deacetylase domain instead of peritrophin domain, indicated that these two proteins may belong to a new chitin‐binding protein family. © 2010 Wiley Periodicals, Inc.     

Structural and biochemical characterization of the nucleoside hydrolase from C. elegans reveals the role of two active site cysteine residues in catalysis

Nucleoside hydrolases (NHs) catalyze the hydrolysis of the N‐glycoside bond in ribonucleosides and are found in all three domains of life. Although in parasitic protozoa a role in purine salvage has been well established, their precise function in bacteria and higher eukaryotes is still largely unknown. NHs have been classified into three homology groups based on the conservation of active site residues. While many structures are available of representatives of group I and II, structural information for group III NHs is lacking. Here, we report the first crystal structure of a purine‐specific nucleoside hydrolase belonging to homology group III from the nematode Caenorhabditis elegans (CeNH) to 1.65Å resolution. In contrast to dimeric purine‐specific NHs from group II, CeNH is a homotetramer. A cysteine residue that characterizes group III NHs (Cys253) structurally aligns with the catalytic histidine and tryptophan residues of group I and group II enzymes, respectively. Moreover, a second cysteine (Cys42) points into the active site of CeNH. Substrate docking shows that both cysteine residues are appropriately positioned to interact with the purine ring. Site‐directed mutagenesis and kinetic analysis proposes a catalytic role for both cysteines residues, with Cys253 playing the most prominent role in leaving group activation.     

Analysis of the peroxiredoxin family: using active-site structure and sequence information for global classification and residue analysis

Peroxiredoxins (Prxs) are a widespread and highly expressed family of cysteine‐based peroxidases that react very rapidly with H₂O₂, organic peroxides, and peroxynitrite. Correct subfamily classification has been problematic because Prx subfamilies are frequently not correlated with phylogenetic distribution and diverge in their preferred reductant, oligomerization state, and tendency toward overoxidation. We have developed a method that uses the Deacon Active Site Profiler (DASP) tool to extract functional‐site profiles from structurally characterized proteins to computationally define subfamilies and to identify new Prx subfamily members from GenBank(nr). For the 58 literature‐defined Prx test proteins, 57 were correctly assigned, and none were assigned to the incorrect subfamily. The >3500 putative Prx sequences identified were then used to analyze residue conservation in the active site of each Prx subfamily. Our results indicate that the existence and location of the resolving cysteine vary in some subfamilies (e.g., Prx5) to a greater degree than previously appreciated and that interactions at the A interface (common to Prx5, Tpx, and higher order AhpC/Prx1 structures) are important for stabilization of the correct active‐site geometry. Interestingly, this method also allows us to further divide the AhpC/Prx1 into four groups that are correlated with functional characteristics. The DASP method provides more accurate subfamily classification than PSI‐BLAST for members of the Prx family and can now readily be applied to other large protein families. Proteins 2011. © 2010 Wiley‐Liss, Inc.