Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 521.
Published online 2017 September 12. doi:  10.1038/s41467-017-00599-0
PMCID: PMC5595980

Discovery of a proteolytic flagellin family in diverse bacterial phyla that assembles enzymatically active flagella


Bacterial flagella are cell locomotion and occasional adhesion organelles composed primarily of the polymeric protein flagellin, but to date have not been associated with any enzymatic function. Here, we report the bioinformatics-driven discovery of a class of enzymatic flagellins that assemble to form proteolytically active flagella. Originating by a metallopeptidase insertion into the central flagellin hypervariable region, this flagellin family has expanded to at least 74 bacterial species. In the pathogen, Clostridium haemolyticum, metallopeptidase-containing flagellin (which we termed flagellinolysin) is the second most abundant protein in the flagella and is localized to the extracellular flagellar surface. Purified flagellar filaments and recombinant flagellin exhibit proteolytic activity, cleaving nearly 1000 different peptides. With ~ 20,000 flagellin copies per  ~ 10-μm flagella this assembles the largest proteolytic complex known. Flagellum-mediated extracellular proteolysis expands our understanding of the functional plasticity of bacterial flagella, revealing this family as enzymatic biopolymers that mediate interactions with diverse peptide substrates.


One of the most remarkable molecular structures in nature is the bacterial flagellum: a self-assembling, 10-μm long molecular nanomachine responsible for cell motility1, 2. Although primarily associated with locomotion, the roles of flagella have expanded to include a diverse range of biological phenomena including adhesion2, mechanosensing3, biofilm interactions35, and virulence6, 7. This functional diversity is mirrored at the molecular level, as the genes encoding flagellar proteins and the structure of flagellar machinery can vary considerably between species8.

The primary structural component of bacterial flagella is the protein flagellin. Up to 20,000 flagellin monomers assemble to form a helical, hollow filament ~ 20 nm in diameter and  ~ 10 μm in length9, which rotates via a proton/sodium motive force to drive cell motility10. Thus, any function encoded by a single flagellin monomer can be massively amplified upon filament polymerization. Flagellin monomers are composed of three structural domains: slowly evolving N-terminal and C-terminal coiled-coil domains that interact in cis to form the core of the filament, and a central hypervariable region of extreme sequence variation, which forms the filament surface8, 11 (Fig. 1a) and thus, is largely responsible for interfacing with the environment6. These interactions play an important role in adherence and colonization of host cells in a diverse range of mammalian1214 and non-mammalian bacterial pathogens15, 16. As a result, flagellins also constitute an important class of pathogen-associated molecular patterns (PAMPs), and are recognized by toll-like receptor 5 in the innate immune system of mammals, flies and plants17.

Fig. 1
Structural model of the flagellar filament and identification of uncharacterized surface-exposed flagellin domains. a Structural model of the flagellar filament and constituent flagellin proteins (left), highlighting the interior flagellin N-terminal ...

Given the extreme sequence diversity of the hypervariable region and its potential as a scaffold for surface-localized functions, we hypothesized that there are flagellins with novel functionality that can be predicted from available bacterial genomes. Here, by computational identification of unexpected domain fusions18, we report the discovery and experimental validation of a family of enzymatically active flagellins present in the genomes of 74 bacterial species including the pathogenic clostridia Clostridium haemolyticum, and strains of Clostridium novyi and Clostridium botulinum. We find that these flagellins harbor a catalytically active zinc-metallopeptidase domain that is localized to extracellular flagellar filaments, so resulting in flagella-embedded protease activity in structures of up to 10 µm. Flagellin-mediated proteolysis expands our understanding of the functional plasticity of flagellar filaments as enzymatic biopolymers, with potential for numerous roles in saprophytic bacteria and in pathogens including biofilm interactions and colonization, tissue colonization and virulence.


Computational prediction of a proteolytic flagellin family

To explore novel functionality in bacterial flagellins, we gathered all 26,587 predicted flagellin sequences from the NCBI non-redundant database and examined them for novel domain architectures. Sixty-one different domains were identified within the flagellin hypervariable region, revealing tremendous domain diversity associated with surface-exposed structures on the flagellar filament (Fig. 1a). Bacterial flagellins are not known to possess enzymatic activity. Therefore, it was unexpected to identify several putative enzymatic domains within the flagellin hypervariable regions (Table 1, Supplementary Table 1). The most frequently occurring flagellin enzymatic domain identified was the Peptidase M9 family domain, a thermolysin-like gluzincin metallopeptidase domain found in bacterial collagenases (Fig. 1b) (Table 1). Sequences annotated as Peptidase M2619 family domains, which include metallopeptidases that cleave mammalian IgA, were also identified within some flagellins (Table 1). After detection of additional homologs via PSI-BLAST, a family of 86 metallopeptidase-containing flagellin sequences were identified from 74 bacterial species and 35 genera (Supplementary Table 2). These species are phylogenetically diverse, indicative of lateral gene transfer, and most commonly occur in the Firmicutes, Betaproteobacteria and Gammaproteobacteria (Fig. 2). In addition, metallopeptidase-containing flagellins are present in several animal pathogens of both medical and agricultural importance including C. haemolyticum, C. novyi A str. 4570, and C. botulinum C/D str. DC5 (Supplementary Table 2). Despite the presence of flagellin and protease domains in all identified members of this flagellin family, most have been labeled as “hypothetical proteins” in the database (Supplementary Table 2), which may have contributed to their previous lack of identification and characterization.

Table 1
Top 10 most frequent protein domains detected within flagellin hypervariable regions according to the NCBI conserved domain database
Fig. 2
Phylogenomic distribution of proteolytic flagellins on bacterial tree of life. Genera with predicted proteolytic flagellin genes are indicated by red lines. The presence/absence pattern is highly scattered indicative of extensive lateral gene transfer. ...

To verify homology of this flagellin hypervariable region family to zinc metallopeptidases, we performed in-depth sequence analysis (Fig. 3a, b) and structural homology modeling (Fig. 3c–e). Reciprocal PSI-BLAST searches20 and three independent structure prediction methods2123 all identified the peptidase domain of clostridial collagenases as the closest homolog. Notwithstanding a low overall sequence identity of ~13%, we constructed a high-confidence structural model based on the available crystal structures of clostridial collagenases (Fig. 3c), which further guided the identification of key proteolytic motifs. Remarkably, all of the identified flagellins conserve the critical Zn++ ion-binding histidine residues in the HExxH catalytic motif—the glutamate being the general base, as well as additional key residues in close proximity to the active site, including the more distal third zinc ion ligating glutamate and Ala265 forming the metalloprotease hydrophobic “basement” (Fig. 3b–d). Thus, members of the domain family are putative gluzincin metallopeptidases19. Their functional importance in these flagellins was further indicated by conservation of the enzymatic motifs rather than their loss through pseudogenization. Genomic context analysis also identified flagellar genes and other peptidases as gene neighbors of metallopeptidase-containing flagellins (Supplementary Fig. 1), suggesting that they form both part of the flagellar machinery and are involved in extracellular proteolytic pathways.

Fig. 3
Sequence analysis and structural modeling of the proteolytic flagellin family. a Multiple alignment of 86 proteolytic flagellin sequences from 74 species, with key functional motifs depicted as sequence logos b including a putative Ca++ ion-binding site ...

Flagellin from C. haemolyticum is an active protease

We selected the flagellin metallopeptidase protein FliA(H) from the animal pathogen C. haemolyticum for biochemical characterization. Notably, C. haemolyticum is associated with collagenolytic activity and host-tissue degradation24. FliA(H) forms a subfamily of flagellins containing a metallopeptidase domain within the hypervariable region, and the absence of this domain in related sequences implies an ancestral insertion event (Supplementary Fig. 2). To focus on the putative gluzincin metallopeptidase domain in the absence of flagellin polymerization, we designed a recombinant cDNA construct containing the FliA(H)-hypervariable region only lacking the N- terminal and C-terminal coiled-coil domain. We expressed and purified the putative metallopeptidase domain to high purity (Fig. 4a), without other protease contaminants as verified by liquid chromatography−tandem mass spectrometry (LC−MS/MS) (Supplementary Fig. 3, Supplementary Table 3).

Fig. 4
Peptidase activity assays of purified recombinant FliA(H)-hypervariable region and flagellar filaments from C. haemolyticum. a 12% SDS-PAGE of purified FliA(H)-hypervariable region (band “i”), electrophoresing close to that predicted for ...

We utilized a high throughput proteomic approach known as proteomic identification of protease cleavage sites (PICS)25, 26 to confirm and characterize the metallopeptidase activity of the FliA(H) hypervariable region and peptide substrate specificity. PICS profiles both the prime (P′) and non-prime (P) amino acid preferences of cleavage sites using biological-derived database searchable peptide libraries. PICS confirmed FliA(H)-hypervariable region as an active peptidase, cleaving 391 peptides in a trypsin-generated peptide library, with each peptide ending in either an arginine or a lysine residue (Fig. 4b). To exclude potential library bias, we repeated the PICS analysis using a GluC-generated peptide library, with each peptide ending at glutamate or aspartate residues. Again, the FliA(H)-hypervariable region cleaved 498 peptides with a highly similar substrate preference (Supplementary Fig. 4a). Employing both peptide libraries is advantageous as each includes within the peptide sequences either the acidic or basic residues, respectively, for full characterization of the substrate specificity. We therefore propose the name “flagellinolysin” and the symbol “FlaMP” for FliA(H) and related metallopeptidase flagellins.

Amino acid occurrences at each substrate subsite were calculated and displayed as heat maps after alignment of the cleaved sequences at the P1↓P1′ scissile bond, where P and P′ residues lie to the N- terminal and C-terminal sides of the scissile bond ↓, respectively (Fig. 4b and Supplementary Fig. 4b). From these, the active site specificity profile from P6 to P6′ was derived (Fig. 4c, Supplementary Fig. 4a for complementary and confirmatory GluC-generated peptide library profiles). Despite the homology with clostridial collagenase, the substrate specificity preference of C. haemolyticum flagellinolysin (Fig. 4b, c) was more similar to their mammalian counterparts, the matrix metalloproteinases (MMPs) (Supplementary Discussion), which cleave a diverse range of proteins in extracellular matrices. Noteworthy, this is in agreement with our structural model (Supplementary Fig. 5). That is: (i) complete absence in the flagellinolysin family of the characteristic “double glycine” motif, which is a key determinant of P1′ substrate specificity in clostridial collagenases27, 28, and (ii) the presence of a deep pocket S1′ subsite typically found in MMPs.

To further confirm protease activity and substrate specificity, we performed quenched fluorescence (QF) peptide cleavage assays using three peptides: ALG↓L, which closely matches the cleavage specificity determined by PICS; the classical MMP substrate QF24 (PLG↓L); and PLG↓V to test the dependency for leucine in P1′ (Fig. 4d). As demonstrated by the fluorescence emission progress curves, C. haemolyticum flagellinolysin readily cleaved both ALG↓L and QF24 (PLG↓L), but did not cleave PLG↓V. This result coincides perfectly with the active site specificity profile determined by PICS and revealed the essential character of the S1′ subsite for bond scission with the absence of PLG↓V cleavage confirming the strong preference for Leu in P1′ similar to MMPs.

Confirming that flagellinolysin was a metalloprotease, we found that peptidase activity was completely abolished in the presence of the metal chelator EDTA (Fig. 4d, red filled circles), indicating divalent cation-dependent proteolytic activity. In contrast, assays with inhibitors of all other protease classes had no effect on activity. Indeed, in our structural model, Asp206 (which substitutes for Glu430 in clostridial collagenase H), forms a Ca++ ion-binding site in close proximity to the active site together with the backbone oxygens of Asn241, Leu244, and Asp247 (Fig. 3d). This structural feature was recently identified in clostridial collagenases to stabilize the active site cleft and was indispensable for full proteolytic activity29. Notably, we mutated the catalytic residue glutamate 234 to alanine, which resulted in complete loss of cleavage of the ALG↓L quenched fluorescent peptide by the recombinant Glu234Ala flagellinolysin (Supplementary Fig. 6).

To provide a structural context for these active site specificity results, we performed high-resolution peptide docking30 using the flagellinolysin homology model and a nonapeptide mimicking the consensus specificity motif identified in PICS (AVTYY↓LVIA) (Fig. 3e, Supplementary Fig. 5). Indeed, the peptide proved an excellent fit and further highlighted the deep hydrophobic S2 and S1′ pockets in the active site.

Flagellinolysins localize to flagellar filament surfaces

To confirm that proteolytic flagellins form components of assembled flagellar filaments, we grew C. haemolyticum anaerobically, and purified intact flagellar filaments by ultracentrifugation after shearing. Transmission electron microscopy (TEM) imaging showed flagellar filaments to be highly pure and intact (Fig. 5a). SDS-PAGE and liquid chromatography−tandem mass spectrometry (LC−MS/MS) verified that flagellar filaments were composed of two different flagellin proteins: the non-proteolytic “structural” flagellin (NCBI accession WP_039229459), which was the dominant band iii, with 82% protein sequence coverage by LC−MS/MS; and the second most abundant protein band ii, identified as the full-length proteolytic flagellin protein (annotated as FliA(H) in the database), and having a predicted molecular weight of 60.9 kDa (Fig. 4e) with 79% protein coverage by LC-MS/MS, including the HExxH active site catalytic motif (Supplementary Fig. 7).

Fig. 5
Immunogold electron microscopy of intact flagellar filaments from C. haemolyticum. a Flagellar filaments (control lacking antibody). b Flagellar filaments labeled with gold 10-nm nanoparticle-conjugated affinity purified polyclonal antibodies raised against ...

To visualize the localization of proteolytic flagellins on assembled flagellar filaments, immunogold labeling electron microscopy was performed using affinity purified polyclonal antibodies raised against the proteolytic domain of flagellinolysin (Fig. 5a–d). The specificity of the antibody was verified by immunostaining using both the recombinant protease domain and full-length proteolytic flagellin from purified sheared filaments (Supplementary Fig. 8). Immunogold TEM imaging revealed highly specific localization of the flagellin protease domain to the surface of individual filaments (Fig. 5b, c), where it was periodically and uniformly distributed, as well as on the surface of flagellar bundles (Fig. 5d).

C. haemolyticum flagella demonstrate proteolytic activity

Having verified proteolytic activity of the isolated flagellinolysin peptidase domain and that flagellinolysin was displayed on the exposed surface of flagella, we next confirmed proteolytic activity in the context of whole flagella. To do so, we purified flagellar filaments sheared from cells and analyzed these by both PICS (Fig. 4f, g, Supplementary Fig. 4c) and quenched fluorescent peptide cleavage assays (Fig. 4h). PICS again revealed robust peptidase activity with the purified flagellar filaments cleaving 269 different peptides, and with a substrate specificity matching that observed for the recombinant monomeric protease domain (Fig. 4g). Quenched fluorescence peptide cleavage assays also demonstrated activity of filaments on the peptide ALG↓L, which was abolished in all controls (Fig. 4h). Peptidase activity was significantly reduced in the supernatant following ultracentrifugation, demonstrating that the observed peptidase activity was filament-associated and not due to free protease monomers in solution. The minor residual (12.8%) activity in the supernatant probably resulted from fragmented filaments that did not sediment in these conditions. Finally, we observed minor flagellar filament-associated autocatalytic activity and modest cleavage of bovine native fibronectin, which was also inhibited by EDTA (Supplementary Fig. 9). The loss of the ~440 kDa fibronectin dimer indicates that a C-terminal cleavage occurs, potentially at 2403EY↓LGA2407, which matches the PICS determined cleavage specificity (Fig. 4c). This site lies N-terminal to the C-terminal interchain crosslink sites between Cys2459 and Cys2463 on bovine fibronectin. As plasma fibronectin is a major secretory protein of liver hepatocytes and cellular fibronectin is an important cell adhesive glycoprotein in connective tissues, fibronectin degradation and depolymerization of the multimeric 440-kDa crosslinked form may facilitate C. haemolyticum colonization, tissue spreading and necrosis. Taken together, these results demonstrated that C. haemolyticum proteolytic flagellins are abundant constituents of assembled flagella, confer peptidase activity to intact filaments, and degrade biological targets of potential physiological relevance to C. haemolyticum virulence.


We report the discovery of bacterial flagella with enzymatic activity. Through bioinformatic analysis of the flagellin hypervariable region, we uncovered hidden functional diversity of flagellin-associated domains. These included peptidase M9 and M2619 enzymatic domains embedded within surface-exposed flagellin hypervariable regions of bacterial flagella. The acquisition of peptidase function may augment its sensory function by generating small peptides and amino acids. Although a remarkable repertoire of domains exists within the rapidly growing collection of bacterial genomes, surprisingly few have been previously characterized. The most common of these, we showed is a catalytically active gluzincin-like metallopeptidase domain present in 74 bacterial species to date, which are scattered across the bacterial tree of life. Despite their abundance in different species, which exceeds that of other characterized flagellin domains such as D2/D3, this proteolytic flagellin family has gone mostly unnoticed to date31.

Indicative of lateral gene transfer, the proteolytic flagellins are highly dispersed on the bacterial tree of life, even at the inter-phylum level. Interestingly, we observed evidence for two separate mechanisms of lateral transfer: one involving possible transfer of a full-length proteolytic flagellin gene (e.g., flagellin gene insertion into Comamonas aquatica, Supplementary Fig. 1); and a second involving a partial insertion of only the protease domain into the flagellin hypervariable region (Supplementary Fig. 2). Consistent with the latter mechanism, C. haemolyticum and several closely related genes possess the metallopeptidase domain, but it is completely absent in the larger flagellin family of clostridial flagellins from which C. haemolyticum flagellin is derived (Supplementary Fig. 2). This further implies an ancestral protease insertion event in the ancestor of the C. haemolyticum flagellin subfamily. Thus, both full and partial gene insertion mechanisms have been utilized.

As flagellinolysin domains form a discrete family that is distant from M9, as well as other gluzincin families, e.g., M2 and M419, the exact identity of the ancestral protease insertion is undetermined. If the proteolytic flagellin hypervariable region evolved by acquisition of a collagenase-related domain, a subsequent shift in substrate specificity must have occurred, in part through substitution of the key glycine–glycine motif, perhaps driven by selection for saprophytic degradation of environmental or biofilm peptides rather than collagen. Indeed, we were unable to demonstrate native or denatured type I collagen cleavage by flagellinolysin that lacks collagen binding exosite domains essential for collagen triple helicase activity28, 32.

To test our bioinformatic prediction of flagellin-associated proteolytic function, we focussed on C. haemolyticum, an important animal pathogen and cause of bacillary hemoglobinuria33. In this organism, proteolytic flagellins assembled as abundant components of flagellar filaments yielding extracellular peptidase activity. Considering a relative abundance of up to 10% in assembled filaments, proteolytically active flagella would consist of 1000–2000 active sites per filament. Thereby, proteolytic flagellins form the largest protease assemblies known to date with up to  ~ 20,000 flagellin copies co-polymerized in a single macromolecular structure  ~ 10 µm long. Especially when considering the extra cell surface area provided by flagellar filaments, the immobilization of enzymes on flagella results in highly localized and amplified enzymatic gradients in controlled environments around a single cell. When directed by cell receptors, flagellated bacteria turn and undergo directed motion toward ascending gradients nutrients and amino acids, which are potentially generated by the flagella itself in cis, to focally improve the nutrient environment for increased pathogenic potential.

Previous studies have exploited the functional potential of the flagellin hypervariable region in vitro by replacing it with adhesive peptides34, 35, fluorescence domains36, and glycosyl hydrolases37. However, prior to our discovery, only structural and adhesive functions for flagellins have been observed in nature. Hence, the newly discovered metallopeptidase flagellinolysin family (FlaMP) we report here highlights an enzymatic role for bacterial flagellar filaments. It also raises the possibility of enzymatic or other unknown functionality in existing flagellin hypervariable region domains of unknown function, e.g., the D3 domain found in Salmonella flagellins. Future work to explore proteolytic flagellins in additional species will examine the numerous potential biological pathways impacted by flagellum-mediated extracellular proteolysis, including biofilm or host-tissue adherence/colonization, virulence, bacterial defense, or extracellular cell–cell communication.

Lastly, our findings provide a fascinating and multilayered story of molecular evolution, involving not only protein domain recombination, but also lateral gene transfer. First, a collagenase-related gene appears to have inserted into a flagellin hypervariable region, presumably within a collagenase-containing lineage such as Clostridium. This is consistent with numerous studies that have documented intragenic recombination in flagellins, which serves as a mechanism for antigenic diversification3843. It is reasonable to assume that this domain insertion likely happened after the evolution of microbial collagenase and MMP-like proteins. The proteolytic flagellin gene then spread to other microorganisms through lateral transfer (Fig. 2, Supplementary Fig. 1), possibly within host-associated and aquatic microbiomes. The remarkable degree of lateral transfer combined with the observed enzymatic activity and clear evidence of selection at the sequence level, suggest that these transfers confer a functional benefit to their recipients.


Computational methods

The CDART tool44 was used to identify all proteins in the NCBI non-redundant protein database containing a predicted flagellin N-hypervariable region-C domain architecture (search performed on 31 January 2017). Identified sequences possessing a central putative gluzincin domain within the flagellin hypervariable region were retrieved and aligned using MUSCLE version v3.8.3145 with default parameters, and visualized with Jalview version 2.80b146. A structural model of the FliA(H)-hypervariable region was generated using the modeling servers SWISS-MODEL47 and iTASSER22 with the crystal structures of clostridial collagenase G (PDB entry 2y3u), H (PDB entry 4ar1), and T (PDB entry 4ar8) as templates28, 29. Phyre 2.0’s intensive modeling procedure was also used and yielded highly similar results21. Energy minimization and secondary structure assignment were performed within the UCSF Chimera package48. Peptide docking was initiated with models based on available collagenase complex structures28, 29 and the PICS peptide cleavage specificity results using the Rosetta FlexPepDock web server30, allowing full peptide flexibility and side-chain flexibility to the protease. Molecular graphic figures were made using the molecular visualization system PyMOL ( Phylogenetic analysis was performed by mapping proteolytic flagellin presence/absence on to a recently constructed tree of life49.

Protein expression and purification

FliA(H)-hypervariable region (Lys149–Ile448, NCBI RefSeq protein sequence accession number WP_039229452.1) cDNA was constructed by gene synthesis and subcloned into the KpnI and XhoI sites in pET-32a, and expressed in Escherichia coli BL21(DE3) cells (Bio Basic Inc.). Protein expression was induced using 0.1 mM IPTG and continued for 3 h at 30 °C. After induction, cells were collected by centrifugation (4000 r.p.m., 10 min). The pellet was re-suspended in 50 mM Tris, 300 mM NaCl, pH 8.0. The cell suspension was sonicated for 30 min in ice bath (sonication power 450 W, a sonication cycle: working time 2 s and pause time 6 s). Additional sonication cycles were repeated until the solution became translucent. The cell debris was removed by centrifugation at 12,000 r.p.m. for 20 min. After ultrasonication, supernatant was purified using HisTrap column (Ni2+ sepharose column, GE healthcare). The sample was applied to the column and washed with at least five column volumes of the washing buffer (50 mM Tris, 300 mM NaCl, pH 8.0). The protein was subsequently eluted with 20 mM Tris-HCl (pH 8.0), 300 mM imidazole, 300 mM NaCl. The protein was then applied to a Superdex 75 column (GE Healthcare) equilibrated in 1X PBS, pH 7.0 while monitoring OD at 280 nm, and the eluted fractions were analyzed via SDS-PAGE. Typically, 2.5 mg of highly pure protein were obtained per 500 ml of E. coli culture. Uncropped gel images used in the preparation of figures are shown in Supplementary Fig. 10.

Proteomic identification of protease cleavage sites assay

PICS assays using whole proteome peptide libraries were performed as described in detail previously25, 26, 5052. To produce E. coli proteome-derived peptide libraries, cell pellets collected from E. coli K12 JM109 cultures (New England Biolabs) were lysed in the presence of protease inhibitors and cell debris was removed by centrifugation. A concentration of 4 M guanidine hydrochloride was used to denature soluble proteins, and cysteine side chains were alkylated using iodoacetamide. After chloroform/methanol precipitation, the pelleted protein interphase was re-suspended, and trypsin (TPCK treated, Sigma-Aldrich) digestion was applied at a protease to proteome ratio of 1:100 (w/w). Reductive dimethylation was used to block primary amines and samples were desalted using size exclusion chromatography and further purified using reversed-phase chromatography. Peptide eluates were concentrated, re-suspended in water, and stored as 200 µg peptide libraries at –80 °C. PICS cleavage assays were performed by incubating 200 µg of peptide library with recombinant FliA(H)-HVR protein. Cleaved peptides were selectively biotinylated, affinity purified and desalted using reversed-phase solid phase extraction. Eluates were vacuum dried to near dryness and stored at –80 °C until LC−MS/MS analysis on a high-resolution quadrupole Time-Of-Flight mass spectrometer (Impact II, Bruker Daltonics), as detailed below.

Spectrum to peptide matching and identification

Peptides were identified at a  1% false discovery rate from the UniProt E. coli K12 database (November 2013) using two search engines, Mascot v2.4.1 (Matrix Science, London, UK) and X! Tandem53, in conjunction with PeptideProphet54 as implemented in the Trans-Proteomic Pipeline v.4.6. Search parameters included a mass tolerance of 15 p.p.m. for the parental ions and 0.1 Da for fragment ions, and allowed up to two missed cleavages. The following fixed peptide modifications were used: carbamidomethylation of cysteine residues (+ 57.02 Da) and dimethylation of lysine ε-amines (+ 28.03 Da). N-terminal dimethylation of uncleaved library carryover peptides (+ 28.03 Da), methionine oxidation (+ 15.99 Da), and thioacylation of protease-generated neo N-termini (+ 88.00 Da) were set as variable modifications. Peptide lists from both search engines were combined within the Trans-Proteomic Pipeline for further analysis. A web-based bioinformatics tool, WebPICS55 was used to reconstruct the non-prime side of each identified unique cleavage site28, 29. Obtained cleavage sites were aligned along the scissile peptide bond and visualized as heat maps in GnuPlot ( and iceLogos56. The MS raw data associated with the present paper are available upon request.

Quenched fluorescence protease activity assays

Stock solutions (1.0 mM) of synthetic quenched fluorescence (QF) peptide substrates (GenScript Inc.) were dissolved in DMSO and working stocks (100 μM) were prepared using the molar extinction coefficient of the conjugated quencher, (2,4)-dinitrophenyl, of 6.985 cm–1 mM–1 at 400 nm. Purified recombinant FliA(H)-hypervariable region protease was incubated at a final concentration of 0.5 µM with 10 µM QF-peptide substrate in 100 µl of 150 mM NaCl, 10 mM CaCl2, 50 mM HEPES, pH 7.5, in the presence of HALT protease inhibitor cocktail (Life Technologies) plus/minus 20 mM EDTA at 37 °C. Three different quenched fluorescent peptides (PLG↓L), Mca-Pro-Leu-Gly-Leu-Dpa-Ala-Arg; ALG↓L, Mca-Ala-Leu-Gly-Leu-Dpa-Ala-Arg; and PLG↓V, Mca-Pro-Leu-Gly-Val-Dpa-Ala-Arg)51 were tested using a multi-wavelength fluorescence scanner (POLARstar OPTIMA, BMG Lab technologies). The excitation and emission wavelengths were 320 and 405 nm, respectively. Fluorescence was measured in arbitrary units over 30 min at 45 s intervals. Values are means of triplicate measurements.

Purification and analysis of flagellar filaments

Cultures of C. haemolyticum were grown anaerobically (in an atmosphere of 80% N2, 10% CO2 and 10% H2) overnight at 30 °C and collected by centrifugation. The cell pellet was re-suspended in 0.1 M Tris-HCl, pH 7.0 in 1/20th of the culture original volume. Flagella were removed from cells by mechanical shearing using a Teflon tissue homogenizer. The flagellar filaments were then purified after removing whole cells by two centrifugation cycles at 5000×g, 15 min each, followed by ultracentrifugation at 130,000×g for 1 h (70 Ti rotor, Beckman Coulter Canada Inc.). The pellets containing the flagella were then re-suspended in ultrapure water followed by a second round of ultracentrifugation (130,000×g, 1 h). The purified flagellin was quantified using BCA protein assay and re-suspended in 50% glycerol and stored at –20 °C57. Fidelity of purification was confirmed by SDS-PAGE (12%). Protein bands of 62 kDa (proteolytic flagellin) and 46 kDa (structural flagellin) were excised and analyzed by in-gel trypsin digestion and sequenced by LC−MS/MS (Advanced Analysis Centre, University of Guelph, Guelph, Ontario, Canada).

Immunogold transmission electron microscopy

The concentration of affinity purified primary rabbit anti-FliA(H) antibody against isolated flagellar proteins used for immunogold TEM was optimized and guided by enhanced chemiluminescent (ECL)-western blotting of the purified recombinant FliA(H) protein, expressed with a thioredoxin tag (Trx), and isolated flagellar proteins. The secondary antibody used was goat HRP-conjugated anti-rabbit IgG.

Isolated flagella (50 ng) from C. haemolyticum were directly pipetted on formvar-coated nickel grids (Ted Pella Inc. Lot# 221014; 200 mesh, Ni; Prod No. 01800 N) and air-dried for 15 min. The grids were blocked in 100 μl of blocking buffer (1X PBS + 0.1% Triton X100 + 1% BSA) for 2 h before incubation with 100 μl of anti-FliA(H) primary antibody diluted 1:50 from stock concentration of 2 mg ml−1 in wash buffer (1X PBS + 1% BSA). The primary antibody incubation was performed overnight in a wet chamber at 4 °C after which the grid was washed three times with wash buffer. For immunogold labeling, grids were incubated in goat anti-rabbit IgG-Gold (10 nm) secondary antibody, diluted 1:100 in the wash buffer, for 1 h. The grids were negatively stained with 1% (w/v) ammonium molybdate. TEM was performed using a Philips CM10 microscope at the Electron and Confocal Microscopy Facility, Department of Biology, University of Waterloo.

Statistical methods

The statistical approaches used in PICS are standardized and involve calculation of position-specific amino frequencies in an experimental versus reference set of peptides, and visualization as sequence logos25, 56 (e.g., Fig. 4). Residue frequencies were calculated based on iterative sampling of peptides from the experimental and reference sets, and P-values were estimated using standard t-tests (distributions are inferred to be normal based on central limit theorem). Use of a large cohort (hundreds to thousands) of individual cleaved sequences provides robust statistical confidence for determination of subsite specificity25.

Statistical confidence for proteins and peptides identified by LC−MS/MS was quantified with PeptideProphet54 using false discovery rate calculations, which estimate the proportion of total predictions that are false due to chance.

Homologs of proteolytic flagellins were identified using NCBI’s CDART tool and PSI-BLAST20. C. haemolyticum flagellinolysin (UniProt entry Q8RR94) was used as an initial query in two rounds of PSI-BLAST, and an E-value threshold of 0.001 to narrow high-confidence matches.

Mass spectrometry

LC–MS/MS analysis was performed on a nano-LC system (EASY-nLC 1000, Thermo Scientific, USA) coupled to a high-resolution quadrupole Time-Of-Flight mass spectrometer (Impact II, Bruker Daltonics) using the CaptiveSpray ionization source (Bruker Daltonics), a 2-cm long, 75-μm inner diameter fused silica trap column, and a 20 cm long, 50-μm inner diameter fused silica fritted analytical column. The trap column was packed with 5-μm diameter Aqua C18 beads (Phenomenex, USA), whereas the analytical column was packed with 1.9-μm diameter Reprosil-Pur C18-AQ beads (Dr. Maisch, Ammerbuch, Germany). Peptides were eluted using a 0 – 80% gradient of organic phase over 90 min. Buffer A was 0.1% formic acid and buffer B was 100% (v/v) acetonitrile with 0.1% formic acid. MS/MS data were acquired automatically using the otofControl software (Bruker Daltonics) for information-dependent acquisition. Error of mass measurement was usually within 5 p.p.m. and was not allowed to exceed 10 p.p.m. Peak lists and mzXML-files of the acquired high resolution Time-Of-Flight mass spectrometer data were created using the Compass DataAnalysis software 4.2 (Bruker Daltonics).

LC−MS/MS of recombinant FliA(H)-hypervariable region

Tandem mass spectrometry was performed essentially as described in detail by Gundry et al.58. In brief, free sulfhydryl groups of 10 μg of FliA(H)-hypervariable region protease in 1X phosphate buffer saline (PBS; 10 mM Na2HPO4, 1.8 mM KH2PO4, 137 mM NaCl, 2.7 mM KCl, pH 7.4) were alkylated with iodoacetamide (20 mM iodoacetamide, 1 h, 20 °C) after initial cysteine reduction with 10 mM dithiothreitol (DTT; 30 min, 20 °C). Adding further DTT (10 mM, 15 min, 20 °C) ended labeling. Trypsin Gold (Promega) was added in a 1:50 w/w ratio and the sample was digested overnight at 37 °C. The reaction was stopped by adding formic acid to a final concentration of 2.5% and desalted by C18-STAGE tip59 prior mass spectrometry analysis.

LC–MS/MS analysis was performed as described below. Peptides were identified at a 1% false discovery rate from a combined UniProt database, containing both human and E. coli proteins and supplemented by our target protein, thioredoxin-tagged FliA(H)-hypervariable region (flagellinolysin). Two search engines, Mascot v2.4.1 (Matrix Science, London, UK) and X! Tandem53 were used in conjunction with PeptideProphet54 as implemented in the Trans-Proteomic Pipeline v.4.6. Search parameters included a mass tolerance of 15 p.p.m. for the parental ions and 0.1 Da for fragment ions, and allowed for one missed cleavage. Carbamidomethylation of cysteine residues (+ 57.02 Da) was used as fixed modification, whereas methionine oxidation (+ 15.99 Da) was set as a variable. Peptide lists from both search engines were combined within the Trans-Proteomic Pipeline for further analysis. The MS raw data associated with the present paper are available upon request.

In-gel digestion and mass spectrometry

In-gel digests were performed as described by Huesgen et al.60. Recombinant flagellinolysin was resolved by 12% SDS-PAGE, stained with Coomassie G-250, and nine bands were excised, destained with 60% acetonitrile, 20 mM ammonium bicarbonate, and then washed in 100% acetonitrile before lyophilization. Gel bands were rehydrated with 10 μl Trypsin Gold (Promega) in 20 mM ammonium bicarbonate (12 ng/μl) by passive diffusion for 1 h at 4 °C. Excess solution was removed and 20 μl 20 mM ammonium bicarbonate was added to the gel plugs and digested overnight at 37 °C. Supernatants were removed and three rounds of active extraction were performed using 20 μl of 1% formic acid, and two times 20 μl of 30% acetonitrile, 1% formic acid. All supernatants were pooled and concentrated to 10 μl using a SpeedVac, and 1–2 μl were analyzed by LC–MS/MS using the same set-up as above, but 30 min gradients. Search parameters included 15 p.p.m. tolerance for MS and 0.10 Da for MS/MS, variable methionine oxidation and propionamide cysteine, and a maximum of one missed cleavage.

Data availability

The authors declare that all data supporting the findings of this study are available from the corresponding authors upon reasonable request.

Electronic supplementary material


This work was supported by the Natural Sciences and Engineering Research Council of Canada, through a Discovery Grant to A.C.D., a postdoctoral fellowship from the Michael Smith Foundation of Health Research (MSFHR) to U.E., a Canada Research Chair in Metalloprotease Proteomics and Systems Biology (C.M.O.), and project grants from the Canadian Institutes of Health Research (CIHR), as well as with infrastructure grants from the MSFHR and the Canada Foundation for Innovation to C.M.O. We would like to thank Glen Armstrong, Klaus Winzer, and Brendan McConkey for helpful discussions. We also thank Simon Chuong and Mishi Groh for technical assistance on immunogold electron microscopy. We gratefully acknowledge Nestor Solis and Theo Klein for invaluable assistance with the mass spectrometer.

Author contributions

Author contributions

U.E., H.B., G.M., J.C., I.W. and J.A. performed experiments. U.E., M.J.M., T.H., T.C.C., J.A., C.M.O. and A.C.D. performed data analysis and interpretation. U.E., H.B., M.J.M., C.M.O. and A.C.D. wrote and prepared the manuscript.


Competing interests

The authors declare no competing financial interests.


Ulrich Eckhard, Hina Bandukwala and Michael J. Mansfield contributed equally to this work.

Christopher M. Overall and Andrew C. Doxey jointly supervised this work.

Electronic supplementary material

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-00599-0.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Christopher M. Overall, ac.cbu@llarevo.sirhc.

Andrew C. Doxey, ac.oolretawu@yexodca.


1. Macnab RM. How bacteria assemble flagella. Annu. Rev. Microbiol. 2003;57:77–100. doi: 10.1146/annurev.micro.57.030502.090832. [PubMed] [Cross Ref]
2. Moens S, Vanderleyden J. Functions of bacterial flagella. Crit. Rev. Microbiol. 1996;22:67–100. doi: 10.3109/10408419609106456. [PubMed] [Cross Ref]
3. Belas R. Biofilms, flagella, and mechanosensing of surfaces by bacteria. Trends Microbiol. 2014;22:517–527. doi: 10.1016/j.tim.2014.05.002. [PubMed] [Cross Ref]
4. Houry A, Briandet R, Aymerich S, Gohar M. Involvement of motility and flagella in Bacillus cereus biofilm formation. Microbiology. 2010;156:1009–1018. doi: 10.1099/mic.0.034827-0. [PubMed] [Cross Ref]
5. Guttenplan SB, Kearns DB. Regulation of flagellar motility during biofilm formation. FEMS Microbiol. Rev. 2013;37:849–871. doi: 10.1111/1574-6976.12018. [PMC free article] [PubMed] [Cross Ref]
6. Ramos HC, Rumbo M, Sirard JC. Bacterial flagellins: Mediators of pathogenicity and host immune responses in mucosa. Trends Microbiol. 2004;12:509–517. doi: 10.1016/j.tim.2004.09.002. [PubMed] [Cross Ref]
7. Duan Q, Zhou M, Zhu L, Zhu G. Flagella and bacterial pathogenicity. J. Basic Microbiol. 2013;53:1–8. doi: 10.1002/jobm.201100335. [PubMed] [Cross Ref]
8. Beatson SA, Minamino T, Pallen MJ. Variation in bacterial flagellins: from sequence to structure. Trends Microbiol. 2006;14:151–155. doi: 10.1016/j.tim.2006.02.008. [PubMed] [Cross Ref]
9. Macnab RM. The bacterial flagellum: reversible rotary propellor and type III export apparatus. J. Bacteriol. 1999;181:7149–7153. [PMC free article] [PubMed]
10. Chevance FFV, Hughes KT. Coordinating assembly of a bacterial macromolecular machine. Nat. Rev. Microbiol. 2008;6:455–465. doi: 10.1038/nrmicro1887. [PubMed] [Cross Ref]
11. Samatey FA, et al. Structure of the bacterial flagellar protofilament and implications for a switch for supercoiling. Nature. 2001;410:331–337. doi: 10.1038/35066504. [PubMed] [Cross Ref]
12. Allen-Vercoe E, Woodward MJ. The role of flagella, but not fimbriae, in the adherence of Salmonella enterica serotype Enteritidis to chick gut explant. J. Med. Microbiol. 1999;48:771–780. doi: 10.1099/00222615-48-8-771. [PubMed] [Cross Ref]
13. Attridge SR, Rowley D. The role of the flagellum in the adherence of Vibrio cholerae. J. Infect. Dis. 1983;147:864–872. doi: 10.1093/infdis/147.5.864. [PubMed] [Cross Ref]
14. Morooka T, Umeda A, Amako K. Motility as an intestinal colonization factor for Campylobacter jejuni. J. Gen. Microbiol. 1985;131:1973–1980. [PubMed]
15. De Weger LA, et al. Flagella of a plant-growth-stimulating Pseudomonas fluorescens strain are required for colonization of potato roots. J. Bacteriol. 1987;169:2769–2773. doi: 10.1128/jb.169.6.2769-2773.1987. [PMC free article] [PubMed] [Cross Ref]
16. Meron D, et al. Role of flagella in virulence of the coral pathogen Vibrio coralliilyticus. Appl. Environ. Microbiol. 2009;75:5704–5707. doi: 10.1128/AEM.00198-09. [PMC free article] [PubMed] [Cross Ref]
17. Hayashi F, et al. The innate immune response to bacterial flagellin is mediated by Toll-like receptor 5. Nature. 2001;410:1099–1103. doi: 10.1038/35074106. [PubMed] [Cross Ref]
18. Lobb B, Doxey AC. Novel function discovery through sequence and structural data mining. Curr. Opin. Struct. Biol. 2016;38:53–61. doi: 10.1016/ [PubMed] [Cross Ref]
19. Rawlings ND, Barrett AJ, Bateman A. MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2012;40:D343–D350. doi: 10.1093/nar/gkr987. [PMC free article] [PubMed] [Cross Ref]
20. Altschul SF, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–3402. doi: 10.1093/nar/25.17.3389. [PMC free article] [PubMed] [Cross Ref]
21. Kelley LA, Sternberg MJE. Protein structure prediction on the Web: a case study using the Phyre server. Nat. Protoc. 2009;4:363–371. doi: 10.1038/nprot.2009.2. [PubMed] [Cross Ref]
22. Yang J, et al. The I-TASSER Suite: protein structure and function prediction. Nat. Methods. 2014;12:7–8. doi: 10.1038/nmeth.3213. [PMC free article] [PubMed] [Cross Ref]
23. Jaroszewski L, Rychlewski L, Li Z, Li W, Godzik A. FFAS03: a server for profile-profile sequence alignments. Nucleic Acids Res. 2005;33:W284–W288. doi: 10.1093/nar/gki418. [PMC free article] [PubMed] [Cross Ref]
24. MacLennan JD. The histotoxic clostridial infections of man. Bacteriol. Rev. 1962;26:177–276. [PMC free article] [PubMed]
25. Schilling O, Overall CM. Proteome-derived, database-searchable peptide libraries for identifying protease cleavage sites. Nat. Biotechnol. 2008;26:685–694. doi: 10.1038/nbt1408. [PubMed] [Cross Ref]
26. Schilling O, Huesgen PF, Barré O, Auf dem Keller U, Overall CM. Characterization of the prime and non-prime active site specificities of proteases by proteome-derived peptide libraries and tandem mass spectrometry. Nat. Protoc. 2011;6:111–120. doi: 10.1038/nprot.2010.178. [PubMed] [Cross Ref]
27. Eckhard U, et al. Biochemical characterization of the catalytic domains of three different clostridial collagenases. Biol. Chem. 2009;390:11–18. doi: 10.1515/BC.2009.004. [PubMed] [Cross Ref]
28. Eckhard U, Schönauer E, Nüss D, Brandstetter H. Structure of collagenase G reveals a chew-and-digest mechanism of bacterial collagenolysis. Nat. Struct. Mol. Biol. 2011;18:1109–1114. doi: 10.1038/nsmb.2127. [PMC free article] [PubMed] [Cross Ref]
29. Eckhard U, Schönauer E, Brandstetter H. Structural basis for activity regulation and substrate preference of clostridial collagenases G, H, and T. J. Biol. Chem. 2013;288:20184–20194. doi: 10.1074/jbc.M112.448548. [PMC free article] [PubMed] [Cross Ref]
30. London N, Raveh B, Cohen E, Fathi G, Schueler-Furman O. Rosetta FlexPepDock web server - High resolution modeling of peptide-protein interactions. Nucleic Acids Res. 2011;39:W249–W253. doi: 10.1093/nar/gkr431. [PMC free article] [PubMed] [Cross Ref]
31. Doxey AC, Lynch MD, Müller KM, Meiering EM, McConkey BJ. Insights into the evolutionary origins of clostridial neurotoxins from analysis of the Clostridium botulinum strain A neurotoxin gene cluster. BMC Evol. Biol. 2008;8:316. doi: 10.1186/1471-2148-8-316. [PMC free article] [PubMed] [Cross Ref]
32. Pelman GR, Morrison CJ, Overall CM. Pivotal molecular determinants of peptidic and collagen triple helicase activities reside in the S3’ subsite of matrix metalloproteinase 8 (MMP-8): the role of hydrogen bonding potential of ASN188 and TYR189 and the connecting cis bond. J. Biol. Chem. 2005;280:2370–2377. doi: 10.1074/jbc.M409603200. [PubMed] [Cross Ref]
33. Uzal, F. A., Songer, J. G., Prescott, J. F. & Popoff, M. R. Clostridial Diseases of Animals. doi:10.1002/9781118728291 (John Wiley & Sons, Inc, 2016).
34. Westerlund-Wikström B, et al. Functional expression of adhesive peptides as fusions to Escherichia coli flagellin. Protein Eng. 1997;10:1319–1326. doi: 10.1093/protein/10.11.1319. [PubMed] [Cross Ref]
35. Tanskanen J, Korhonen TK, Westerlund-Wikstrom B. Construction of a multihybrid display system: flagellar filaments carrying two foreign adhesive peptides. Appl. Environ. Microbiol. 2000;66:4152–4156. doi: 10.1128/AEM.66.9.4152-4156.2000. [PMC free article] [PubMed] [Cross Ref]
36. Klein A, Tóth B, Jankovics H, Muskotál A, Vonderviszt F. A polymerizable GFP variant. Protein Eng. Des. Sel. 2012;25:153–157. doi: 10.1093/protein/gzs003. [PubMed] [Cross Ref]
37. Szabó V, Muskotál A, Tóth B, Mihovilovic MD, Vonderviszt F. Construction of a xylanase A variant capable of polymerization. PLoS ONE. 2011;6:e25388. doi: 10.1371/journal.pone.0025388. [PMC free article] [PubMed] [Cross Ref]
38. Frankel G, Newton SM, Schoolnik GK, Stocker BA. Intragenic recombination in a flagellin gene: characterization of the H1-j gene of Salmonella typhi. EMBO J. 1989;8:3149–3152. [PubMed]
39. Smith NH, Beltran P, Selander RK. recombination of Salmonella phase 1 flagellin genes generates new serovars. J. Bacteriol. 1990;172:2209–2216. doi: 10.1128/jb.172.5.2209-2216.1990. [PMC free article] [PubMed] [Cross Ref]
40. Harrington CS, Thomson-Carter FM, Carter PE. Evidence for recombination in the flagellin locus of Campylobacter jejuni: Implications for the flagellin gene typing scheme. J. Clin. Microbiol. 1997;35:2386–2392. [PMC free article] [PubMed]
41. Reid SD, Selander RK, Whittam TS. Sequence diversity of flagellin (fliC) alleles in pathogenic Escherichia coli. J. Bacteriol. 1999;181:153–160. [PMC free article] [PubMed]
42. Fitzgerald C, Jones K, Fields P. Evidence for inter-species recombination of the flagellin (flaA) locus between Campylobacter coli and Campylobacter jejuni. Int. J. Med. Microbiol. 2001;291:66–67.
43. Mortimer CKB, Gharbia SE, Logan JMJ, Peters TM, Arnold C. Flagellin gene sequence evolution in Salmonella. Infect. Genet. Evol. 2007;7:411–415. doi: 10.1016/j.meegid.2006.12.001. [PubMed] [Cross Ref]
44. Geer LY, Domrachev M, Lipman DJ, Bryant SH. CDART: protein homology by domain architecture. Genome Res. 2002;12:1619–1623. doi: 10.1101/gr.278202. [PubMed] [Cross Ref]
45. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–1797. doi: 10.1093/nar/gkh340. [PMC free article] [PubMed] [Cross Ref]
46. Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ. Jalview Version 2-A multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–1191. doi: 10.1093/bioinformatics/btp033. [PMC free article] [PubMed] [Cross Ref]
47. Bordoli L, Schwede T. Automated protein structure modeling with swiss-model workspace and the protein model portal. Methods Mol. Biol. 2012;857:107–136. doi: 10.1007/978-1-61779-588-6_5. [PubMed] [Cross Ref]
48. Yang Z, et al. UCSF Chimera, MODELLER, and IMP: an integrated modeling system. J. Struct. Biol. 2012;179:269–278. doi: 10.1016/j.jsb.2011.09.006. [PMC free article] [PubMed] [Cross Ref]
49. Hug LA, et al. A new view of the tree of life. Nat. Microbiol. 2016;1:16048. doi: 10.1038/nmicrobiol.2016.48. [PubMed] [Cross Ref]
50. Eckhard U, Huesgen PF, Brandstetter H, Overall CM. Proteomic protease specificity profiling of clostridial collagenases reveals their intrinsic nature as dedicated degraders of collagen. J. Proteom. 2014;100:102–114. doi: 10.1016/j.jprot.2013.10.004. [PMC free article] [PubMed] [Cross Ref]
51. Eckhard U, et al. Active site specificity profiling of the matrix metalloproteinase family: proteomic identification of 4300 cleavage sites by nine MMPs explored with structural and synthetic peptide cleavage analyses. Matrix Biol. 2016;49:37–60. doi: 10.1016/j.matbio.2015.09.003. [PubMed] [Cross Ref]
52. Barré O, et al. Cleavage specificity analysis of six type II transmembrane serine proteases (TTSPs) using PICS with proteome-derived peptide libraries. PLoS ONE. 2014;9:e105984. doi: 10.1371/journal.pone.0105984. [PMC free article] [PubMed] [Cross Ref]
53. Craig R, Beavis RC. TANDEM: matching proteins with tandem mass spectra. Bioinformatics. 2004;20:1466–1467. doi: 10.1093/bioinformatics/bth092. [PubMed] [Cross Ref]
54. Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal. Chem. 2002;74:5383–5392. doi: 10.1021/ac025747h. [PubMed] [Cross Ref]
55. Schilling O, auf dem Keller U, Overall CM. Factor Xa subsite mapping by proteome-derived peptide libraries improved using WebPICS, a resource for proteomic identification of cleavage sites. Biol. Chem. 2011;392:1031–1037. doi: 10.1515/BC.2011.158. [PubMed] [Cross Ref]
56. Colaert N, Helsens K, Martens L, Vandekerckhove J, Gevaert K. Improved visualization of protein consensus sequences by iceLogo. Nat. Methods. 2009;6:786–787. doi: 10.1038/nmeth1109-786. [PubMed] [Cross Ref]
57. Paul CJ, et al. Flagellin diversity in Clostridium botulinum groups I and II: a new strategy for strain identification. Appl. Environ. Microbiol. 2007;73:2963–2975. doi: 10.1128/AEM.02623-06. [PMC free article] [PubMed] [Cross Ref]
58. Gundry RL, et al. Preparation of proteins and peptides for mass spectrometry analysis in a bottom-up proteomics workflow. Curr. Protoc. Mol. Biol. 2009;Supp. 10:Unit10.25. [PMC free article] [PubMed]
59. Rappsilber J, Ishihama Y, Mann M. Stop and go extraction tips for matrix-assisted laser desorption/ionization, nanoelectrospray, and LC/MS sample pretreatment in proteomics. Anal. Chem. 2003;75:663–670. doi: 10.1021/ac026117i. [PubMed] [Cross Ref]
60. Huesgen PF, et al. LysargiNase mirrors trypsin for protein C-terminal and methylation-site identification. Nat. Methods. 2015;12:55–58. doi: 10.1038/nmeth.3177. [PubMed] [Cross Ref]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group