- Research article
- Open Access
Pepsin homologues in bacteria
© Rawlings and Bateman; licensee BioMed Central Ltd. 2009
- Received: 28 November 2008
- Accepted: 16 September 2009
- Published: 16 September 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.
- Horizontal Gene Transfer
- Horizontal Transfer
- Active Site Residue
- Aspartic Peptidase
- Disulphide Bridge
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 126.96.36.199), gastricsin (EC 188.8.131.52) and chymosin (EC 184.108.40.206), which are digestive enzymes in the stomach; cathepsin D (EC 220.127.116.11), which is a lysosomal enzyme that degrades phagocytozed peptides and proteins; renin (EC 18.104.22.168), which is a peptidase that processes angiotensinogen; and memapsins 1 (or BACE-2, EC 22.214.171.124) and 2 (BACE-1, EC 126.96.36.199), 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 188.8.131.52; ). 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 184.108.40.206), 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.
Homologues of pepsin from completely sequenced bacterial genomes
Marinomonas sp. MWYL1
Predicted characteristics of bacterial homologues of pepsin
Predicted active site residues
Alien_Hunter score (threshold)
D36, Y68, D300
Marinomonas sp. MWYL1
D49, F95, D274
D37, Y69, D282
D37, Y69, D287
D41, Y73, D399
D41, Y88, D364
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.
Bacterial pepsin homologues in environmental samples
closest pepsin homologue
Eastern Pacific ocean surface water
Eastern Pacific ocean surface water
Eastern Pacific ocean surface water
microbial mat from gray whale carcass in the Santa Cruz Basin, Pacific Ocean (depth 1674 metres)
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.
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).
Comparison of replacements among bacteria with pepsin homologues
A8 (signal peptidase 2)
S16 (endpeptidase Clp)
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.
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].
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Rawlings ND, Barrett AJ: Evolutionary families of peptidases. Biochem J. 1993, 290: 205-218.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Lee MS: Molecular clock calibrations and metazoan divergence dates. J Mol Evol. 1999, 49: 385-391. 10.1007/PL00006562.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Hill J, Phylip LH: Bacterial aspartic proteinases. FEBS Lett. 1997, 409: 357-360. 10.1016/S0014-5793(97)00547-4.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP - phylogeny inference package. Cladistics. 1989, 5: 164-166.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Goodstadt L, Ponting CP: CHROMA: consensus-based colouring of multiple alignments for publication. Bioinformatics. 2001, 17: 845-846. 10.1093/bioinformatics/17.9.845.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.