- Research article
- Open Access
Pepsin homologues in bacteria
BMC Genomicsvolume 10, Article number: 437 (2009)
Peptidase family A1, to which pepsin belongs, had been assumed to be restricted to eukaryotes. The tertiary structure of pepsin shows two lobes with similar folds and it has been suggested that the gene has arisen from an ancient duplication and fusion event. The only sequence similarity between the lobes is restricted to the motif around the active site aspartate and a hydrophobic-hydrophobic-Gly motif. Together, these contribute to an essential structural feature known as a psi-loop. There is one such psi-loop in each lobe, and so each lobe presents an active Asp. The human immunodeficiency virus peptidase, retropepsin, from peptidase family A2 also has a similar fold but consists of one lobe only and has to dimerize to be active. All known members of family A1 show the bilobed structure, but it is unclear if the ancestor of family A1 was similar to an A2 peptidase, or if the ancestral retropepsin was derived from a half-pepsin gene. The presence of a pepsin homologue in a prokaryote might give insights into the evolution of the pepsin family.
Homologues of the aspartic peptidase pepsin have been found in the completed genomic sequences from seven species of bacteria. The bacterial homologues, unlike those from eukaryotes, do not possess signal peptides, and would therefore be intracellular acting at neutral pH. The bacterial homologues have Thr218 replaced by Asp, a change which in renin has been shown to confer activity at neutral pH. No pepsin homologues could be detected in any archaean genome.
The peptidase family A1 is found in some species of bacteria as well as eukaryotes. The bacterial homologues fall into two groups, one from oceanic bacteria and one from plant symbionts. The bacterial homologues are all predicted to be intracellular proteins, unlike the eukaryotic enzymes. The bacterial homologues are bilobed like pepsin, implying that if no horizontal gene transfer has occurred the duplication and fusion event might be very ancient indeed, preceding the divergence of bacteria and eukaryotes. It is unclear whether all the bacterial homologues are derived from horizontal gene transfer, but those from the plant symbionts probably are. The homologues from oceanic bacteria are most closely related to memapsins (or BACE-1 and BACE-2), but are so divergent that they are close to the root of the phylogenetic tree and to the division of the A1 family into two subfamilies.
Peptidases are widespread enzymes that catalyse the hydrolysis of peptide bonds. Not only do peptidases hydrolyse proteins to amino acids and peptides for nutrition and recycling, but they also perform some of the most important post-translational processing events leading to the activation (or inactivation) of many other proteins, including other enzymes and peptide hormones. Over 2% of the protein coding genes in a genome encode peptidases, and there are over 500 peptidase genes in the human genome. Peptidases exist in at least six catalytic types, depending on the nature of the nucleophile in the catalytic reaction (either the hydroxyl of a serine or threonine, the thiol of a cysteine, or activated water bound either to a metal ion (metallopeptidases) or aspartate or glutamate residues). Peptidases are grouped into over 250 different evolutionary families .
The pepsin family (MEROPS family A1 ) contains aspartic-type endopeptidases. Activated water is bound by two aspartate residues each of which is located in a motif Xaa-Xaa-Asp-Xbb-Gly-Xbb where Xaa is a hydrophobic amino acid and Xbb is either Ser or Thr. A third important residue is a tyrosine which interacts with the substrate and is located on a beta-hairpin loop known as the "flap" which forms part of the roof to the active site. A second motif, part of a structure known as the psi-loop  consists of two hydrophobic residues followed by a glycine. Two of these motifs are conserved in all sequences known or predicted to be active within family A1. The psi-loops are important in stabilizing the structure, especially the active centre. Many members of family A1 have three conserved disulphide bridges. Homologues are known from animals, plants, fungi and protozoa, but no homologues have been reported from prokaryotes. Family A1 has recently been divided into two subfamilies, with the majority of characterized peptidases in subfamily A1A. Peptidases within this subfamily include: pepsin A (EC 184.108.40.206), gastricsin (EC 220.127.116.11) and chymosin (EC 18.104.22.168), which are digestive enzymes in the stomach; cathepsin D (EC 22.214.171.124), which is a lysosomal enzyme that degrades phagocytozed peptides and proteins; renin (EC 126.96.36.199), which is a peptidase that processes angiotensinogen; and memapsins 1 (or BACE-2, EC 188.8.131.52) and 2 (BACE-1, EC 184.108.40.206), widely distributed peptidases of unknown function but which may be involved in Alzheimer's disease (memapsin-2 is currently under investigation as a drug target ). Subfamily A1B contains predominantly peptidase homologues from plants including nucellin  and nepenthesin (EC 220.127.116.11; ). Besides peptidases, family A1 also contains a number of proteins which cannot be peptidases because the active site aspartates are not conserved. These include pregnancy-associated glycoproteins  and a xylanase inhibitor .
The tertiary structures of several members of this family have been determined, and each shows a bilobed protein. Each lobe contains one of the active site aspartates in a psi-loop and the two lobes are structurally similar. This similarity can only be observed in the tertiary structure because, apart from motifs contributing to each psi-loop, there is no observable similarity in the protein sequence. This led to the hypothesis that an ancestral gene duplication had taken place, followed by gene fusion [9, 9]. This raises the question of how did the original gene product function, being only one half of a modern-day pepsin? The discovery of a peptidase, retropepsin (EC 18.104.22.168), within the genome of the human immunodeficiency virus (HIV) that also represents a single lobe of an extremely divergent pepsin (so divergent that it is the representative of a different family, A2, in MEROPS, but within the same clan or superfamily AA) and which has to dimerize to be active, provided evidence to answer this question . Proteins with a single pepsin lobe also exist in eukaryotes. The structure of the Saccharomyces cerevisiae protein Ddi1 (DNA-damage inducible protein 1), which is not a peptidase, has been resolved and shown to contain a domain with a retropepsin-like structure, and the protein also has to dimerize to be active . This domain was probably derived from a retroelement and represents an even more ancient divergence than that of the A1 and A2 families, as has been shown in a recent phylogeny of clan AA peptidases . It has been assumed that the HIV peptidase represents the ancestral state and that the ancestral half-pepsin also had to dimerize to be active. However, it is also possible that the HIV peptidase was derived from a normal, bilobed pepsin when a virus captured half of its host's gene. This "chicken and egg" situation might be resolvable if pepsin homologues could be found in prokaryotes, because this would represent an extremely distant sequence divergence. If a prokaryote pepsin homologue contained only one lobe, then it would support the hypothesis of the ancient gene duplication and fusion event. However, if the prokaryote pepsin homologue were bilobed, then either a) this hypothesis is incorrect; b) the gene duplication precedes the most recent common ancestor of prokaryotes and eukaryotes; or c) the bacterial homologue genes were derived by horizontal gene transfer.
We report the discovery of the first prokaryote pepsin homologues in the completed genome sequences of several proteobacteria.
Results and Discussion
Over 960 completely sequenced prokaryote genomes were analysed (914 bacterial and 55 archaean), and 953 of these contained no detectable homologue of pepsin. This included all of the archaean genomes. However, pepsin homologues were detected in seven bacterial genomes, and these are shown in Table 1. Characteristics of the predicted proteins of these homologues are shown in Table 2. All the bacteria containing pepsin homologues are members of the class Gammaproteobacteria, and all except Marinomonas are members of the order Alteromonadales; Marinomonas is a member of the order Oceanospirillales. Despite the genomes having been completely sequenced, no pepsin homologue was detected in the proteomes of the following Shewanella species: S. baltica, S. frigidimarina, S. halifaxensis, S. oneidensis, S. pealeana, S. putrefaciens, S. sp. ANA-3, S. sp. MR-4, S. sp. MR-7, S. sp. W3-18-1 or S. woodyi, nor in Sinorhizobium meliloti.
BlastP searches of the NCBI non-redundant database using the bacterial pepsin homologues as queries returned only one other bacterial sequence. This was the A5A_A0203 gene product from Vibrio cholerae strain MZO-2 (UniProt accession A6A7Y6). This sequence consists of 382 residues with the peptidase unit predicted to be residues 14-257, and the active site residues to be Asp41, Phe87 and Asp249 (numbered according to the translated coding sequence). The genome of this particular strain is incomplete and no pepsin homologue is present in any of the other thirteen strains of Vibrio cholerae for which the genome sequences are complete. Because this sequence might be a contaminant, it has not been included in the analysis. A search of the UniProt sequence database using HMMER also failed to return any other sequences from bacteria. Two Psi-Blast searches of the NCBI non-redundant protein sequence database were undertaken using the amino acid sequences of the pepsin homologues from Shewanella amazonesis and Marinomonas as the seed sequences for 19 and 14 iterations respectively (returning 2888 known pepsin homologues); these also failed to find any new bacterial homologues.
Four homologues were found among the environmental samples nucleotide sequence database (see Table 3), and all were from marine environments. The predicted protein sequences of AACY023056044 and AACY024067159 were virtually identical. None of these fragments were identical to any known pepsin homologue, but each was most closely related to a homologue from Shewanella and are likely to be derived from at least three other bacterial species. Because these are fragments only, they are not considered further.
Several of the bacteria with pepsin homologues are marine and psychrophilic, including Colwellia psychrerythraea (formerly Vibrio psychroerythus)  and Shewanella detrificans. Shewanella loihica is also marine but associated with deep-sea, hydrothermal vents. Shewanella amazonensis was isolated from marine sediments from the shallow waters of the Amazon river delta, and Shewanella sediminis from sediments in Halifax Harbour, Canada . Marinomonas sp. MWYL1 is a salt marsh species initially isolated from the root surface of the grass Spartina anglica. Exceptionally, Sinorhizobium medicae is found associated with root nodules of Medicago species . There are completed genome for several other Shewanella species that are also marine (S. baltica and S. frigidimarina, for example), but these do not contain pepsin homologues, so there is no apparent correlation between environment and presence of a pepsin homologue.
An alignment containing only bacterial homologues, human pepsin A and memapsins 1 and 2, generated by extracting the sequences from the MUSCLE alignment and removing all the gap-only columns, is shown in Fig. 1. As can be seen from this figure, both active site aspartates (Asp32 and Asp215) are conserved, indicating that the bacterial homologues are bilobed. This means that each homologue has the same structure for the peptidase unit as pepsin, and each would be active in the monomeric form.
None of the proteins were predicted to possess a signal peptide. This means that unlike the majority of members of peptidase family A1, these bacterial proteins are not secreted (the aspartic peptidase BcAP1 from the plant pathogenic fungus Botrytis cinerea also does not possess a signal peptide nor disulfide bonds ). Many members of family A1 are active at acidic pH, being secreted in the stomach, or to the lysosome or plant and fungal vacuoles, but a cytoplasmic peptidase would presumably be active at neutral pH. Amongst the mammalian pepsin homologues, renin is secreted into the blood and is also active at neutral pH. This change in pH optimum has been explained by the replacement of Thr218 by Ala, thereby preventing a hydrogen bond forming which affects the acidity of the active site residue Asp215. That this replacement is also found in the HIV retropepsin, which is also active at neutral pH, supports this hypothesis . Intriguingly, six of the seven bacterial homologues also have Ala218 (see Fig. 1). In pepsin A, Thr218 is hydrogen-bonded to Asp303, but in renin both residues are replaced by Ala and site-directed mutagenesis of Ala303 for Asp lowered the pH optimum . In the sequences from Shewanella and Colwellia, Asp303 is replaced by Leu, which is an isosteric replacement. In the Sinorhizobium sequence Asp303 is replaced by Arg, and in the Marinomonas sequence it is replaced by Ser. There would be no problem accommodating the smaller Ser in the available space, but the larger Arg could be more problematic.
Fig. 1 shows that of the cysteines forming the three disulphide bridges found in many of the eukaryotic homologues, only those forming the first are conserved in the Sinorhizobium and Marinomonas homologues. Cys282 from the third disulphide bridge is retained in the sequences from Colwellia, Shewanella loihica and S. sediminis, which may mean that the proteins are thiol-dependent. Disulphide bridges would not be expected in intracellular proteins. The Sinorhizobium and Marinomonas homologues also differ in having Phe75 instead of Tyr (although it must be acknowledged that the alignment here is by no means certain); although this is a common replacement, none of the homologues with Phe75 has ever been biochemically characterized or shown to be catalytically active. Replacement of Tyr75 for Phe in rhizopuspepsin by site-directed mutagenesis led to some, weakened activity . The hydrophobic-hydrophobic-Gly motifs in the psi-loops are conserved in all the bacterial homologues except Sinorhizobium and Marinomonas . Because these bacterial homologues have inserts here, the alignment is uncertain. Although it is possible that the inserts might compensate for the loss of the hydrophobic-hydrophobic-Gly motif, which would imply a different structure in this region, it is much more likely that the two proteins lacking these motifs are not active as peptidases.
From the solved tertiary structure the residues forming the S1 and S1' substrate-binding pockets for human pepsin A are known. Pepsin A prefers large hydrophobic residues (Phe and Leu) in substrates for both P1 and P1', and the substrate-binding pockets are correspondingly lined with hydrophobic residues . The bacterial homologues, with the exception of that from Marinomonas, have most of the substrate-binding residues conserved except Thr77, which is replaced by another hydrophobic residue, and Thr218, which is replaced by Ala (see Fig. 1).
The phylogenetic tree (Fig. 2) shows that all bar two of the bacterial sequences are close to the origin of the division between subfamilies on the tree, implying that horizontal transfer of genes from a recent eukaryote species to a bacterium is unlikely (the key is found in Additional File 1). This is confirmed by the results from Alien Hunter (see Table 2) which finds horizontal transfer of genes unlikely in all species except Marinomonas. Considering that Marinomonas is a commensal organism living on the root surface of a grass, the origin of a pepsin homologue via horizontal transfer of a gene from the host is not unexpected. The Marinomonas homologue is most closely related to that from Sinorhizobium. These sequences do not cluster with the other bacterial homologues on the tree. Sinorhizobium is also a commensal organism, living in the root nodules of legumes, but the region containing the gene was not predicted to be the result of horizontal transfer. There are therefore two groups of bacterial pepsins, one containing sequences from Shewanella and Colwellia and one containing sequences from Sinorhizobium and Marinomonas.
To further investigate whether the pepsin homologues in Shewanella species were derived from horizontal gene transfer, we attempted to estimate the rate of mutation. We hypothesized that if the bacterial pepsins are of ancient origin then their rate of mutation must be very low otherwise we would not be able to see the sequence similarity to eukaryotic proteins. We calculated the percentage identities between the pepsin sequences from Shewanella species and compared these with the percentage identities from other peptidase families. The families of signal peptidase 2 (peptidase family A8) and the ClpP subunit of endopeptidase Clp (family S16) were chosen because these are well conserved in bacteria, and all the Shewanella species have only one homologue per family. The percentage identities with respect to the homologues from Shewanella denitrificans are shown in Table 4, which shows that the ClpP subunit is the most stable, followed by signal peptidase 2 and then pepsin. The pepsin homologues are changing at twice the rate of the other two families. The closest homologues to S. denitrificans signal peptidase 2 and ClpP are the respective proteins from S. loihica. However, of the Shewanella pepsin homologues, that from S. loihica is most distantly related to S. denitrificans, even more distantly related than the sequence from a species in a different genus, Colwellia. To put the percentage identities into a geological time-frame, percentage identities were calculated for cathepsin D, a pepsin homologue found only in animals where divergence times can be estimated from the fossil record. The percentage identity between cathepsin D sequences from human and pig is 88% (human and pig diverged around 65 million years ago ), and human and Xenopus tropicalis is 73% (the species diverged around 350 million years ago ); while that of human cathepsin D and nemepsin-2 from the nematode Ancylostoma caninum is 52% (the species diverged around 660 million years ago ). The mutation rate for pepsin homologues amongst bacteria must therefore be considerably higher than eukaryote members of the family, and this rapid mutation rate might explain their placing near the root of the phylogenetic tree.
The evidence as to whether the ancestral pepsin gene in bacteria originated from a horizontal gene transfer from a eukaryote or not can be summed up as follows. The positioning of most of the bacterial pepsins close to the divergence of the two subfamilies on the phylogenetic tree, the cytoplasmic location and lack of disulfide bridges and the scores from Alien Hunter imply no lateral transference; the absence of homologues in other closely related bacterial species and Archaea and the apparent fast mutation rate are points in favour of horizontal gene transfer.
If the bacterial genes are the result of lateral gene transference from an ancient eukaryote gene, then there is no known current eukaryote gene that resembles it, because nearly all modern pepsin homologues are secreted proteins and the ancient gene product would presumably be cytoplasmic. The eukaryotic peptidases that are most closely related to the bacterial homologues are the memapsins, as can be seen in Fig. 2. The memapsins and bacterial pepsins cluster with high confidence. It is possible that memapsins, which are widely distributed in mammalian tissues, might be closer to the ancestral peptidase from which all other vertebrate pepsin homologues are derived.
Aspartic-type peptidases unrelated to pepsin are known in bacteria, including the gpr peptidase from Bacillus megaterium (peptidase family A25)  and omptin from Escherichia coli (peptidase family A26) . Neither peptidase is inhibited by pepstatin. Sporulation factor SpoIIGA (peptidase family U4)  has also been claimed to be an aspartic peptidase because of the presence of a single Asp-Ser-Gly motif (the protein is assumed to be active as a dimer), but this alone is not sufficient because a same motif occurs around the active site Asp of the serine-type peptidase subtilisin. Pepstatin-sensitive aspartic peptidases have previously been found in Escherichia coli and Haemophilus influenzae, but the sequences were not homologous to that of pepsin but homologues of gluconate permease . A recent publication  has reported an acidic peptidase from a Synergistes species isolated from the anaerobic digester used for treatment of tannery solid waste. This peptidase is inhibited 75% by 0.01 mM pepstatin, but not by inhibitors of other catalytic types. No sequence is available, but this may be the first characterized pepsin homologue from a bacterium. Unlike the homologues from the marine bacteria, this is presumably a secreted protein. No complete genome sequence for any Synergistes species is publicly available, so we were unable to test for the presence of a pepsin homologue. The genomes of Dethiosulfovibrio peptidovorans and Thermanaerovibrio acidaminovorans, both of which are members of the same bacterial phylum Synergistetes, have been partially sequenced (genome projects 20741 and 29531, respectively). These were searched with the amino acid sequences of pepsin homologues from Shewanella amazonensis and Sinorhizobium medicae, but no pepsin homologues were found, even at an E value of 10.
The pepsin family of peptidases (family A1) is not just confined to eukaryotes, but is also present in some bacteria. These bacterial homologues are predicted to be intracellular, whereas nearly all other members of family A1 enter the secretory pathway. The bacterial homologues are structurally similar to pepsin, consisting of two lobes each of which bears one active site Asp. Two of the genes for bacterial pepsins might be derived from horizontal gene transfer from a eukaryote, but whether all the bacterial genes for the bacterial pepsins are so derived is debatable, and evidence is presented for and against the case. If no horizontal gene transfer has occurred, then the hypothetical gene duplication and fusion event that gave rise to the ancestral pepsin gene would therefore have occurred before the divergence of bacteria and eukaryotes. If so, this has implications for the origin of the viral and retrotransposon peptidases in family A2, which were either derived from the single-lobed ancestral sequence, or secondarily from the horizontal transfer of half a pepsin gene from a host. For the retropepsin gene to be derived by vertical transfer the ancestral gene would have to have been present in a prokaryote. However, retroviruses infect only eukaryotes and retrotransposons are unknown in prokaryotes, making a more recent horizontal transfer the more likely origin of the retropepsin gene.
Data sources and analyses
All the complete, predicted proteomes (in FastA format) from bacterial genome sequencing projects were downloaded from the National Center for Biotechnology Information (NCBI) FTP site http://www.ncbi.nlm.nih.gov/Ftp/. Each proteome was submitted to the MEROPS batch Blast Internet utility . This service performs a BlastP search  against a representative set of peptidase sequences for every sequence in the proteome library. For each sequence that is detected as homologous to a peptidase sequence (where the expect value is less than 0.001) the extent of the peptidase unit is determined and the potential active site residues are predicted. The peptidase unit is defined as the subdomains determined from the tertiary structure of the type example peptidase (in this case human pepsin A) that bear the known active site residues.
The MEROPS batch Blast software determines to which peptidase family a sequence belongs. To further characterize the sequences, each sequence returned in the MEROPS batch Blast output was submitted to a further BlastP search against the complete library of peptidase unit sequences. During this second round of searching, homologous fragments are distinguished from possible false positive fragments by using a more rigorous expect value threshold (< e-10). A sequence was defined as a fragment if in the BlastP alignment between the query and hit sequences there was no overlap in the regions of sequence where active site residues were predicted to be.
Any pepsin homologue detected was used in a further BlastP search of the NCBI non-redundant protein sequence database to try to find more distantly related homologues. The protein sequence database at NCBI from environmental samples was also searched for homologues using BlastP. Additionally, TBlastN searches  were done against the environmental samples nucleotide sequence database and against all nucleotide sequences in the NCBI non-redundant nucleotide sequence database. The UniProt database was searched using the HMMER software http://hmmer.janelia.org/ and a Hidden Markov Model built from any pepsin homologues found. Psi-Blast  searches of the NCBI non-redundant database were also undertaken, accumulating hits with an E-value of 0.005 or less.
Sequences were aligned using MUSCLE  to holotypes of the pepsin family. In the MEROPS database a holotype sequence is chosen to be a representative of each different peptidase characterized biochemically. Holotypes have also been selected for every different human, mouse and Arabidopsis thaliana gene product predicted to be a peptidase homologue. Any sequences where the active site residues were not conserved were removed, leaving 135 sequences.
A phylogenetic tree derived from the MUSCLE alignment was generated using the PROML (protein maximum likelihood) program from the Phylip package . A consensus tree was also generated using the BOOTSEQ, PROML and CONSENSE programs from the Phylip package for a subset of the holotype library, retaining all holotypes from subfamily A1A plus the sequences of nepenthesin and nucellin from subfamily A1B. The number of bootstraps was 100, otherwise default values were used.
Signal peptide prediction
Prediction of horizontally transfer DNA
In order to find regions where DNA might have been horizontally transferred between species, rather than vertically transferred from a common ancestor, the program Alien Hunter was used . This program uses compositional biases (including CG content and third codon usage) to detect potential horizontally transferred genetic material. The default value window of 5000 bases was used. Any region scoring higher than the calculated threshold for the genome is potentially derived by horizontal transfer.
Rawlings ND, Morton FR, Kok CY, Kong J, Barrett AJ: MEROPS: the peptidase database. Nucleic Acids Res. 2008, 36: D320-D325. 10.1093/nar/gkm954.
Rawlings ND, Barrett AJ: Evolutionary families of peptidases. Biochem J. 1993, 290: 205-218.
Cooper JB, Khan G, Taylor G, Tickle IJ, Blundell TL: X-ray analyses of aspartic proteinases. II. Three-dimensional structure of the hexagonal crystal form of porcine pepsin at 2.3 A resolution. J Mol Biol. 1990, 214: 199-222. 10.1016/0022-2836(90)90156-G.
Ghosh AK, Gemma S, Tang J: beta-Secretase as a therapeutic target for Alzheimer's disease. Neurotherapeutics. 2008, 5: 399-408. 10.1016/j.nurt.2008.05.007.
Bi X, Khush GS, Bennett J: The rice nucellin gene ortholog OsAsp1 encodes an active aspartic protease without a plant-specific insert and is strongly expressed in early embryo. Plant Cell Physiol. 2005, 46: 87-98. 10.1093/pcp/pci002.
Athauda SB, Matsumoto K, Rajapakshe S, Kuribayashi M, Kojima M, Kubomura-Yoshida N, Iwamatsu A, Shibata C, Inoue H, Takahashi K: Enzymic and structural characterization of nepenthesin, a unique member of a novel subfamily of aspartic proteinases. Biochem J. 2004, 381: 295-306. 10.1042/BJ20031575.
Hughes AL, Green JA, Piontkivska H, Roberts RM: Aspartic proteinase phylogeny and the origin of pregnancy-associated glycoproteins. Mol Biol Evol. 2003, 20: 1940-1945. 10.1093/molbev/msg217.
Sansen S, De Ranter CJ, Gebruers K, Brijs K, Courtin CM, Delcour JA, Rabijns A: Structural basis for inhibition of Aspergillus niger xylanase by triticum aestivum xylanase inhibitor-I. J Biol Chem. 2004, 279: 36022-36028. 10.1074/jbc.M404212200.
Tang J, James MNG, Hsu IN, Jenkins JA, Blundell TL: Structural evidence for gene duplication in the evolution of the acid proteases. Nature. 1978, 271: 618-621. 10.1038/271618a0.
Lapatto R, Blundell T, Hemmings A, Overington J, Wilderspin A, Wood S, Merson JR, Whittle PJ, Danley DE, Geoghegan KF, et al: X-ray analysis of HIV-1 proteinase at 2.7 Å resolution confirms structural homology among retroviral enzymes. Nature. 1989, 342: 299-302. 10.1038/342299a0.
Sirkis R, Gerst JE, Fass D: Ddi1, a eukaryotic protein with the retroviral protease fold. J Mol Biol. 2006, 364: 376-387. 10.1016/j.jmb.2006.08.086.
Llorens C, Futami R, Renaud G, Moya A: Bioinformatic flowchart and database to investigate the origins and diversity of Clan AA peptidases. Biol Direct. 2009, 4: 3-10.1186/1745-6150-4-3.
Methe BA, Nelson KE, Deming JW, Momen B, Melamud E, Zhang X, Moult J, Madupu R, Nelson WC, Dodson RJ, et al: The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses. Proc Natl Acad Sci USA. 2005, 102: 10913-10918. 10.1073/pnas.0504766102.
Brettar I, Christen R, Hofle MG: Shewanella denitrificans sp. nov., a vigorously denitrifying bacterium isolated from the oxic-anoxic interface of the Gotland Deep in the central Baltic Sea. Int J Syst Evol Microbiol. 2002, 52: 2211-2217. 10.1099/ijs.0.02255-0.
Hau HH, Gralnick JA: Ecology and biotechnology of the genus Shewanella. Annu Rev Microbiol. 2007, 61: 237-258. 10.1146/annurev.micro.61.080706.093257.
Van Landschoot A, de Ley J: Intra- and intergeneric similarities of the rRNA cistrons of Alteromonas, Marinomonas (gen. nov.) and some other gram-negative bacteria. J Gen Microbiol. 1983, 129: 3057-3074.
Rome S, Fernandez MP, Brunel B, Normand P, Cleyet-Marel JC: Sinorhizobium medicae sp. nov., isolated from annual Medicago spp. Int J Syst Bacteriol. 1996, 46: 972-980.
Ten Have A, Dekkers E, Kay J, Phylip LH, van Kan JA: An aspartic proteinase gene family in the filamentous fungus Botrytis cinerea contains members with novel features. Microbiology. 2004, 150: 2475-2489. 10.1099/mic.0.27058-0.
Ido E, Han HP, Kezdy FJ, Tang J: Kinetic studies of human immunodificiency virus type 1 protease and its active-site hydrogen bond mutant A28S. J Biol Chem. 1991, 266: 24359-24366.
Yamauchi T, Nagahama M, Hori H, Murakami K: Functional characterization of Asp-317 mutant of human renin expressed in COS cells. FEBS Lett. 1988, 230: 205-208. 10.1016/0014-5793(88)80672-0.
Park YN, Aikawa J, Nishiyama M, Horinouchi S, Beppu T: Involvement of a residue at position 75 in the catalytic mechanism of a fungal aspartic proteinase, Rhizomucor pusillus pepsin. Replacement of tyrosine 75 on the flap by asparagine enhances catalytic efficiency. Protein Eng. 1996, 9: 869-875. 10.1093/protein/9.10.869.
Abad-Zapatero C, Rydel TJ, Neidhart DJ, Luly J, Erickson JW: Inhibitor binding induces structural changes in porcine pepsin. Adv Exp Med Biol. 1991, 306: 9-21.
Lee MS: Molecular clock calibrations and metazoan divergence dates. J Mol Evol. 1999, 49: 385-391. 10.1007/PL00006562.
Hugall AF, Foster R, Lee MS: Calibration choice, rate smoothing, and the pattern of tetrapod diversification according to the long nuclear gene RAG-1. Syst Biol. 2007, 56: 543-563. 10.1080/10635150701477825.
Carroll TM, Setlow P: Site-directed mutagenesis and structural studies suggest that the germination protease, GPR, in spores of Bacillus species is an atypical aspartic acid protease. J Bacteriol. 2005, 187: 7119-7125. 10.1128/JB.187.20.7119-7125.2005.
Vandeputte-Rutten L, Kramer RA, Kroon J, Dekker N, Egmond MR, Gros P: Crystal structure of the outer membrane protease OmpT from Escherichia coli suggests a novel catalytic site. EMBO J. 2001, 20: 5033-5039. 10.1093/emboj/20.18.5033.
Imamura D, Zhou R, Feig M, Kroos L: Evidence that the Bacillus subtilis SpoIIGA protein is a novel type of signal-transducing aspartic protease. J Biol Chem. 2008, 283: 15287-15299. 10.1074/jbc.M708962200.
Hill J, Phylip LH: Bacterial aspartic proteinases. FEBS Lett. 1997, 409: 357-360. 10.1016/S0014-5793(97)00547-4.
Kumar AG, Nagesh N, Prabhakar TG, Sekaran G: Purification of extracellular acid protease and analysis of fermentation metabolites by Synergistes sp. utilizing proteinaceous solid waste from tanneries. Bioresour Technol. 2008, 99: 2364-2372. 10.1016/j.biortech.2007.05.001.
Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2008, 36: D13-D21. 10.1093/nar/gkm1000.
Rawlings ND, Morton FR: The MEROPS batch BLAST: a tool to detect peptidases and their non-peptidase homologues in a genome. Biochimie. 2008, 90: 243-259. 10.1016/j.biochi.2007.09.014.
Altschul SF, Madden TL, Schöffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.
Felsenstein J: PHYLIP - phylogeny inference package. Cladistics. 1989, 5: 164-166.
Emanuelsson O, Brunak S, von Heijne G, Nielsen H: Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007, 2: 953-971. 10.1038/nprot.2007.131.
Vernikos GS, Parkhill J: Interpolated variable order motifs for identification of horizontally acquired DNA: revisiting the Salmonella pathogenicity islands. Bioinformatics. 2006, 22: 2196-2203. 10.1093/bioinformatics/btl369.
Goodstadt L, Ponting CP: CHROMA: consensus-based colouring of multiple alignments for publication. Bioinformatics. 2001, 17: 845-846. 10.1093/bioinformatics/17.9.845.
We would like to thank Dr K. Woodwark, Dr P. Gardner and Dr A. J. Barrett for helpful suggestions and discussions. We would also like to thank the diligent reviewers of this manuscript who made many constructive suggestions. This work was supported by the Wellcome Trust [grant number WT077044/Z/05/Z].
NDR conceived the study, performed the sequence analyses and drafted the manuscript. AB participated in the design and helped to draft the manuscript. Both authors read and approved the final manuscript.