Skip to main content

Bioinformatic analysis of the neprilysin (M13) family of peptidases reveals complex evolutionary and functional relationships



The neprilysin (M13) family of endopeptidases are zinc-metalloenzymes, the majority of which are type II integral membrane proteins. The best characterised of this family is neprilysin, which has important roles in inactivating signalling peptides involved in modulating neuronal activity, blood pressure and the immune system. Other family members include the endothelin converting enzymes (ECE-1 and ECE-2), which are responsible for the final step in the synthesis of potent vasoconstrictor endothelins. The ECEs, as well as neprilysin, are considered valuable therapeutic targets for treating cardiovascular disease. Other members of the M13 family have not been functionally characterised, but are also likely to have biological roles regulating peptide signalling. The recent sequencing of animal genomes has greatly increased the number of M13 family members in protein databases, information which can be used to reveal evolutionary relationships and to gain insight into conserved biological roles.


The phylogenetic analysis successfully resolved vertebrate M13 peptidases into seven classes, one of which appears to be specific to mammals, and insect genes into five functional classes and a series of expansions, which may include inactive peptidases. Nematode genes primarily resolved into groups containing no other taxa, bar the two nematode genes associated with Drosophila DmeNEP1 and DmeNEP4. This analysis reconstructed only one relationship between chordate and invertebrate clusters, that of the ECE sub-group and the DmeNEP3 related genes. Analysis of amino acid utilisation in the active site of M13 peptidases reveals a basis for their biochemical properties. A relatively invariant S1' subsite gives the majority of M13 peptidases their strong preference for hydrophobic residues in P1' position. The greater variation in the S2' subsite may be instrumental in determining the specificity of M13 peptidases for their substrates and thus allows M13 peptidases to fulfil a broad range of physiological roles.


The M13 family of peptidases have diversified extensively in all species examined, indicating wide ranging roles in numerous physiological processes. It is predicted that differences in the S2' subsite are fundamental to determining the substrate specificities that facilitate this functional diversity.


The neprilysin (M13) family of zinc-metallopeptidases is a large group of medically and developmentally important enzymes of which mammalian neprilysin (EC, neutral endopeptidase, NEP) was the first member to be biochemically characterised [1]. They metabolise bioactive peptides and are involved in a number of biological processes in mammals including modulation of neurotransmitter levels, reproduction, control of blood pressure and cancer progression. The majority of M13 peptidases described so far are endopeptidases with a strong preference for cleaving the amino-terminal bond of hydrophobic residues [2]. Many are selectively inhibited by phosphoramidon and, typically, the substrates are small to medium sized peptides, including tachykinins, opioid peptides, big-endothelins and bombesin. Recently, it has been shown that two members of the neprilysin family cleave the Alzheimer's amyloid β-peptide (Aβ) in the mammalian brain and that ex-vivo expression of neprilysin reduces amyloid plaque burden in a mouse model [3]. Neprilysin and the neprilysin-like peptidases are typically type II integral membrane proteins with their active sites facing the extracellular environment [2]. Soluble neprilysin-like enzymes however do occur in mammalian (human MMEL2 and rodent SEP/NL1 and NEPII) and insect (Drosophila melanogaster) tissues [4, 5]. Both mammalian and insect soluble enzymes are strongly expressed in the testes suggesting a physiological role in reproduction. Indeed, female mice mated with males lacking SEP/NL1 have smaller litters confirming the important role for this enzyme in mammalian reproduction.

The endothelin converting enzymes (ECE-1 and ECE-2) and the ECE-like group of enzymes all have distinctive biology. ECE-1 exists as four isoforms and has a physiological role in the metabolism of endothelins by generating mature endothelin from the inactive precursor big-endothelin [6]. ECE-1 knockout mice show a fatal developmental phenotype with severe disruption to several developmental processes, including craniofacial development [7]. ECE-2 is predominantly neurally expressed and is thought to be involved in the processing of peptides prior to secretion [8].

ECEL-1 and its rodent homologue, damage-induced neuronal endopeptidase (DINE), are two of the least well-characterised members of the M13 family [9]. DINE was identified due to its up-regulation after neuronal damage and has been shown to be neuro-protective. DINE-knockout mice develop normally, but die immediately after birth due to an inability to inflate their lungs [10]. No physiological substrate for ECEL-1/DINE has been identified, which has hindered understanding the neuro-protective mechanism and its role in respiratory control.

PHEX is a gene that has been shown to be deficient in patients suffering X-linked hypophosphataemic rickets (XLH) [11]. PHEX has no known natural substrate, although it has been reported that it may cleave FGF-23, a member of the fibroblast growth factor family that inhibits renal tubular phosphate transport [12]. PHEX is unusual for an M13 peptidase in that it appears to have a preference for acidic residues at its S1' site [13]. Kell is an important blood group antigen but is also found in the sertoli cells in the testes [14, 15]. Kell can convert big endothelin-3 to biologically active endothelin-3 from its precursor protein [14]. Kell is an atypical M13 peptidase having no transmembrane domain, but is instead normally covalently anchored to the membrane protein XK [16]. Deficiency in Kell does not cause disease, but XK null patients (who also lack Kell) suffer from McLeod's syndrome, which leads to acanthocytic anaemia [16].

Invertebrate M13 peptidases have been found in organisms ranging from Hydra vulgaris through to highly derived insects such as Drosophila melanogaster. The sea snail Aplysia californica and the mussels, Mytilus edulis and Mytilus galloprovincialis, have neprilysin-like peptidases [2]. These activities are involved in control of feeding in the snail and have been implicated in modulating the response of mussel immune effector cells [17, 18]. Free-living and parasitic nematodes also possess neprilysin-like peptidases that cleave peptide bonds N-terminal to hydrophobic residues and are inhibited by phosphoramidon [1922]. A role for neprilysin in nematode locomotion and reproduction has been established in a study of a deletion mutant of Caenorhabditis elegans NEP1 [23]. Neprilysin-like activity is enriched in the brain neuropil and in isolated synaptic membranes of relatively basal insects, the locusts,Schistocerca gregaria and Locusta migratoria and the cockroach Leucophea maderae, indicating an evolutionarily conserved role for M13 peptidases in the functioning of nervous systems that use neuropeptides extensively as neurotransmitters/modulators [2, 24]. Recently an ECE-like gene was identified in L. migratoria and was shown to be highly expressed in the central nervous system and the midgut [25]. Insect M13 peptidases are associated with metamorphosis [26, 27] and immunity to bacterial, fungal and protozoan infections [28, 29]. The most thoroughly characterised insect M13 peptidase is D. melanogaster NEP2, which is expressed in the stellate cells of Malpighian (renal) tubules and in the testes of adult male flies [4, 30], where expression is strongest in the elongating cyst cells.

Genome sequencing projects and individual gene studies have provided a large and expanding set of protein sequence data for comparative genomic and phylogenetic studies. Phylogenetics has traditionally concentrated on deducing the evolutionary relationships of various taxa being examined, ranging from whole reconstruction of phyla level relationships to finer resolution studies of individual species groups. However, phylogenetic techniques can also be used to attempt to unravel the functional and evolutionary relationships of a set of paralogous genes. The work presented here provides an extensive examination of the functional relationships of M13 peptidases and provides new insights into the evolution of this medically important family of peptidases.

Results and Discussion

Identification and alignment of M13 protein sequences

A combination of methods (including, BLASTP, PSI-BLAST and HMMs) were employed to identify M13 sequences from the genomes of seventeen organisms, eleven of which were from fully sequenced genomes and six of which were from individual entries in GenBank (Additional file 1). Simple BLASTP analysis against NCBI, ENSEMBL and various species specific databases identified the majority of sequences used in this study (Additional file 1); human neprilysin [NP_009220] was used as the query sequence. The remainder were identified from the raw genomic data of D. pseudoobscura and Ap. mellifera using the SHARKhunt software [31], This was particularly valuable for the analysis of the Ap. mellifera genome, as we identified two unannotated genes (Ame1 and Ame6).

The large dataset of protein sequences described above was used to produce a multiple sequence alignment (Figure 1, Additional file 2). The program, MUSCLE was used to align these sequences [32, 33]. This program was chosen as it aligns areas of strong local homology, which is advantageous because of the strong conservation of catalytic regions, but high variability in other regions of M13 proteins, in particular, the N-terminal region. This alignment was edited to remove large uninformative insertions and was then analysed manually to determine the conservation of key catalytic residues (Figure 1) as defined by mutagenesis studies and by the crystal structure of human neprilysin [34, 35].

Figure 1
figure 1

A section of a multiple sequence alignment of M13 peptidases. A multiple sequence alignment of 111 protein sequences was generated using MUSCLE [32,33] and was as the basis for the analysis presented. A highly conserved section of the alignment, representing residues 541 to 652 of human neprilysin, contains important catalytic residues. These residues include the HExxH zinc binding motif and the catalytically important GENIAD and VNAFY motifs which are coloured blue. For full alignment see additional file 2.

Phylogenetic analysis

Three methods of phylogenetic reconstruction were compared in this study and were resolved into a single consensus tree (Figure 2; for bootstrap values see Additional file 3). The three methods were able to resolve clades at the extremes of the branches of the tree with a strong consistency across all three methods. Deeper resolution was less complete with many branches originating from the same node (Figure 2).

Figure 2
figure 2

Phylogenetic analysis of M13 peptidases. Majority consensus tree of all three methods of phylogenetic reconstruction. The tree was generated using CTree [65]. See additional file 3 for bootstrap values. Key to species; Dme Drosophila melanogaster, Aga Anopheles gambiae, Dps Drosophila pseudoobscura, Lmi Locusta migratoria, Bmo Bombyx mori, Ame Apis mellifera, Hsa Homo sapiens, Rno Rattus norvegicus, Mmu Mus musculus, Fru Fugu rubripes, Pfl Perca flavescens, Xla Xenopus laevus, Ocu Oryctolagus cuniculus, Bta Bos taurus, Cpo Cavia porcelllus, Cin Ciona intestinalis, Cel Caenorhabditis elegans and Cbr Caenorhabditis briggsae.

Insect neprilysins form five main clades

Insect neprilysins from An. gambiae, Ap. mellifera and D. pseudoobscura cluster in five groups of strongly conserved genes that delineate with the five D. melanogaster genes that show the greatest similarity to human neprilysin (Figure 2 III, VI, VIII, IX). These five clades are present in all phylogenetic reconstructions and each D. melanogaster gene can be seen to have a D. pseudoobscura, An. gambiae, and Ap. mellifera homologue, indicating that these genes evolved before the divergence of the Hymenoptera and the Diptera, approximately 35 and 260 million years ago, respectively [36, 37].

The majority of the five clades described above contain only insect genes. Interestingly, several insect NEP genes, including DmeNEP1 and DmeNEP4, form a cluster with the C. elegans gene CelTO5A8.4 and its C. briggsae homologue CbrPO4554 (Figure 2 VI). DmeNEP1, in particular, also shares an extremely high conservation of active site residues with CelTO5A8.4, suggesting they are functionally very similar. This is the only example of Drosophila genes clustering with nematode genes. Neighbour-joining places CelTO5A8.4 and CbrPO4554 outside the DmeNEP1/DmeNEP4 cluster. However, maximum likelihood places them within the cluster with a strong relationship to DmeNEP1. Therefore, it is unclear what is the true position of these nematode genes within this cluster and whether DmeNEP1 and DmeNEP4 were formed from a duplication before or after divergence of the Nematoda and Insecta

DmeNEP3 contains the ECE-like VNAYY motif (Figure 1) and was seen to cluster with the extended clade containing the ECE sub-group of enzymes (Figure 2 III). This is the only example of insect and vertebrate genes clustering in this analysis and is seen with all three types of analyses. There is however some variation as to the position of the C. intestinalis gene Cin2. C. intestinalis is a urochordate and shares a very ancient common ancestor with vertebrates, having diverged from that lineage before the evolution of the Craniata and therefore provides the best insight into an intermediate between vertebrates and invertebrates [38]. The neighbour-joining method places Cin2 outside the cluster of insect and vertebrate genes. Character based methods however, return trees which better reflect the evolutionary relationships of these species with Cin2 inside the clade between the insect and vertebrate clusters (Figure 2 III). It seems reasonable to assume that Cin2 evolved from the same ancestral ECE from which the vertebrate ECEs also evolved. This provides strong evidence that these genes arose early in the metazoa before the evolution of the majority of organisms used in this analysis.

DmeNEP2 is a soluble secreted enzyme found in the renal tubules and testes of D. melanogaster. Interestingly, the mammalian group of soluble peptidases (MmuSEP and HsaMMELII) appears to have evolved recently (Figure 2 IV), after the split from the urochordates and is therefore un-related to DmeNEP2. DmeNEP2 expression is modulated by dietary phosphate, leading to the suggestion that this enzyme is functionally related to HsaPHEX [39]. However, the fly enzyme lacks the highly conserved PHEX motifs [40] and presumably does not share the preference of HsaPHEX for acidic residues in the P1' position [4, 12, 30]. Therefore, it is not surprising that in the current analysis DmeNEP2 shows no association with the PHEX cluster.

Vertebrate M13 genes form distinct functionally related clades

The majority of the vertebrate M13 peptidases delineate into functionally related clusters representing the PHEX, Kell, ECE-1, ECE-2, ECEL1, neprilysin and MMELLII groups of peptidases (Figure 2 I-V). All of these clusters contain human and rodent sequences, with the majority also containing Fugu rubripes NEPs. Most of these groups have one or more C. intestinalis sequence forming a root at the base of the cluster, indicating that these genes evolved before the divergence of these lineages. The Kell group of genes has neither a C. intestinalis nor a F. rubripes sequence associated with it, suggesting that Kell is a recently evolved member of the M13 family, restricted to mammals (Figure 2 II). A single C. intestinalis sequence roots the clade containing NEP/MMELLII clusters (Figure 2 IV). As mentioned previously, this implies that the neprilysin and soluble members of the family arose during the evolution of vertebrates. Both clusters contain fish sequences indicating that, although these sequences evolved after the evolution of the Craniata, they must have arisen early in vertebrate evolution. There are two F. rubripes sequences which cluster with neprilysin, indicating either a gene duplication event in fish, or a gene loss in mammals. Similarly, it appears that ECE-1 and ECE-2 arose from a common ancestral gene after the divergence of the vertebrates and urochordates. Interestingly, no F. rubripes sequence is found in the ECE-1 sub-cluster (Figure 2 III), which may indicate that ECE-2 is the prototypical member of this group as ECE-2 is seen to cluster with a F. rubripes sequence. Though, it is possible that F. rubripes may have lost an ECE-1 gene, the drastic consequences of losing ECE-1 in mammals makes this unlikely. Also the genome of another fish, Tetraodon nigroviridis, contains two ECE-like genes and the frog, Xenopus laevis, has an ECE-1 protein, and therefore it is conceivable that the ECE-1 sequence of F. rubripes has yet to be determined. Maximum likelihood analysis suggests the possibility that Kell is derived from the ECE subgroup, which is consistent with the report that Kell efficiently cleaves big-endothelin-3, although the physiological relevance of this in vitro study is unclear [16]. The placing of Kell, however, may simply be due to its unusual characteristics masking true evolutionary relationships.

Two C. intestinalis (Cin3 and Cin4) genes are seen to cluster with the PHEX group of peptidases indicating a potential recent, gene duplication in Ciona (Figure 2 V). The predicted proteins are very similar to each other with only minor changes to the side chains that form the substrate binding pockets, suggesting possible substrate specificity differences.

Invertebrate gene expansions

The majority of the nematode genes analysed in this work form a series of extended clades (Figure 2 VII, X-XII). The majority of C. elegans genes cluster with a C. briggsae homologue indicating that these gene expansions arose before these species split, approximately 100 million years ago (mya) [41]. However some species-specific expansions have also occurred, consistent with these families continuing to evolve in these lineages. The most thoroughly described C. elegans M13 peptidase ZK20.6 (NEP-1) clusters with several worm peptidases, but is quite distinct from sequences from outside the Nematoda.

Outside of the five clades described above, the majority of Drosophila sequences form a series of Drosophila-specific expansions (Figure 2 XIII). Many of these expansions contain at least one D. pseudoobscura sequence again suggesting that the expansion occurred before the divergence of these species approximately 35 mya [37]. The expansion of the NEP family is larger in D. melanogaster, which may have resulted from different selection pressures on these species. These expansions appear to have occurred fairly early in the diversification of the old world Drosophila as D. melanogaster and D. pseudoobscura represent two of the major species complexes of the Drosophila genus. The fact that the clusters containing these expansions contain no An. gambiae or Ap. mellifera sequences is confirmatory evidence for these expansions being specific to Drosophila and not due to gene loss in other organisms. All of the genes belonging to these groups are expressed (as full length ESTs), but some lack key catalytic residues, indicating they probably encode catalytically inactive proteins with possibly novel non-peptidase roles. Non-enzymatic functions for metallopeptidase-like proteins have been shown to be important for development in other organisms [4244], but there is currently no evidence for physiological roles for these non-peptidase members of the M13 family. Nevertheless, the fact that these proteins have been conserved as non-peptidases over 35 million years of evolution suggests that they do indeed have important functions.

Conservation of catalytically and structurally important residues

Site-directed mutagenesis studies together with the elucidation of the crystal structure of human neprilysin complexed with specific inhibitors have led to a detailed understanding of the structure of the active site and the catalytic mechanism [34]. In human neprilysin, this involves co-ordination of the zinc ion by His583, His587 and Glu646, and the involvement of Glu584 in polarising the water molecule that attacks the peptide bond of the substrate. These residues and the substrate co-ordinating asparagine (Asn542) were generally conserved across the protein sequences examined, apart from the aforementioned invertebrate non-peptidases (Figure 1, Additional file 2).

There is also good conservation of the cysteine residues involved in the formation of the disulphide bonds of neprilysin which is indicative of a general conservation of the tertiary structure of M13 peptidases. However, a cysteine bridge between Cys233 and Cys241 of human neprilysin is only present in mammalian neprilysins and no other M13 peptidases (Additional file 2). The absence of this disulphide bond in HsaECE-1 and RnoNEPII shows that it is not essential for M13 peptidase activity per se. Several genes code for proteins that also lack an equivalent to Cys142 of human neprilysin that forms a bridge with Cys410. Interestingly, this bond is also not required for peptidase activity of DmeNEP2, which lacks an equivalent to Cys142 [4]. Conservation of Cys410 is much more common, which may suggest an alternative disulphide bridge or function for this residue in other proteins. Further structural data on M13 peptidases will hopefully resolve these issues. Thirty sequences had significant deletions or substitutions in important catalytic or structural positions and were all invertebrate sequences except for Rattus norvegicus Kell, which has a lysine instead of the water-activating Glu584 and a serine rather than the zinc co-ordinating Glu646 of human neprilysin (Additional file 2).

Of particular importance in determining cleavage site specificity are the S1' and S2' sub-sites that interact with the side-chains of the two residues immediately C-terminal to the scissile peptide bond. Multiple alignment of the neprilysin sequences indicated that there was a greater degree of variation in the S2' than in the S1' subsite (Figure 3). The majority of S1' substitutions are relatively conservative, maintaining a hydrophobic environment (Figure 3). The S1' pocket of the Ap. mellifera protein Ame7 also contains a glutamate for isoleucine substitution, as well as the substitution of tyrosine and serine for the hydrophobic phenylalanine and valine, respectively resulting in a pocket that would presumably be more accommodating for charged or polar residues.

Figure 3
figure 3

Utilisation of residues by ligand binding subsites of M13 peptidases. The utilisation of the twenty amino acids in both the S1' and S2' subsites was examined. The percentage contribution of each amino acid to either binding site was calculated and residues were placed in the order of descending frequency for the S1' subsite.

The S2' subsite shows a greater spectrum of properties, from being predominantly charged, as in human neprilysin, to being predominantly polar as in D. melanogaster NEP2, to being almost entirely hydrophobic as in Ap. mellifera Ame4. Interestingly, the S2' subsite of PHEX is highly polar, which is consistent with PHEX's strong P2' preference for polar residues [12]. The conservation of a hydrophobic S1' subsite provides a likely explanation why M13 peptidases from distantly related organisms retain the ability to cleave similar substrates, whereas the range of side chains found at the S2' subsite might determine the different peptide bond preferences seen between enzymes [2, 4, 25, 30]. Variation in the properties of the S2' subsite may be particularly informative about evolutionary pressures on diversifying the physiological substrates of these peptidases [4, 45].

The 541VNAFY545 motif of human neprilysin is important for the orientation of the peptide bond between the P1 and P1' positions of the substrate [34]. This motif is conserved, although there is some variation which is frequently associated with specific functional classes. The ECE subgroup of enzymes have VNAYY rather than the VNAFY motif [25, 46], which is important for determining the specificity of ECE-1 for big-endothelin [6]. Interestingly, the ECE-like group of enzymes (ECEL-1, DINE, etc), for which there are no known physiological substrates, are unique in that all have a LNAYY motif, which may influence substrate specificity. The VNAYY motif is present in a number of other sequences including the Kell proteins as well as various insect M13 peptidases. Both DmeNEP3 and DmeNEP4 contain the VNAYY motif, whereas DmeNEP1 and DmeNEP2 contain the neprilysin-like VNAFY. All the sequences that cluster with DmeNEP5, apart from Ame6, contain a unique version of the VNAFY motif that has a histidine in place of the tyrosine. The tyrosine in the VNAFY motif of human neprilysin is part of a hydrogen bonding network, which is conserved in other M13 peptidases [46, 47]. It is not clear what effect this substitution or some of the other substitutions described here would have; however it is likely that structural changes occur, with possible knock-on effects on the properties of the enzymes.

The GENIAD motif contains the zinc-co-ordinating glutamate and is generally conserved amongst the M13 peptidases (Figure 1). In some of the more divergent sequences significant variation is seen in this motif and will probably indicate a loss of activity, as residues in this sequence have been shown to be critical for endopeptidase activity [35]. However, some sequences show small changes which may not disrupt peptidase activity. For example, the majority of F. rubripes genes contain a glycine to alanine substitution producing an AENIAD consensus. Interestingly the Kell proteins have a distinct motif (LENAAD) that is conserved across all Kell proteins examined. All substitutions in this motif are conservative, with similar overall side chain properties.

The analysis presented here highlights the functional diversity and complicated evolutionary relationships of M13 peptidases. The phylogenetic analysis successfully resolved vertebrate M13 peptidases into seven classes, one of which appears to be specific to mammals, and insect genes into five functional classes and a series of expansions, which may include inactive peptidases. Nematode genes primarily resolved into groups containing no other taxa, bar the two nematode genes associated with DmeNEP1 and DmeNEP4. This analysis reconstructed only one relationship between chordate and invertebrate clusters, that of the ECE sub-group and the DmeNEP3 related genes. This may be because sequences have diverged so far that reconstructing an evolutionary relationship is impossible. Evolution of these classes could be quite cryptic due to the high rate of evolution and gene loss in the invertebrates [4851]. If this is the case then it is possible that these proteins are also no longer functionally similar. However, as the S1' and S2' subsites are the main factors in defining M13 substrate specificity, active site analysis may still give insights into functional similarities.


Our analysis shows the M13 family of peptidases to have diversified extensively in all species examined, indicating wide ranging roles in numerous physiological processes. It is predicted that differences in the S2' subsite are fundamental to determining the substrate specificities that facilitate this functional diversity. The work presented here provides the most thorough and sophisticated analysis to date of the phylogenetic relationships of these peptidases and will provide a strong framework for the study of these genes across model systems and in humans.


In silico identification of M13 genes

M13 gene sequences were identified by probing sequence repositories using the BLASTP program [52]. Human neprilysin [NP_009220] was used as the query sequence. Protein sequences from: Drosophila melanogaster, Caenorhabditis briggsae, C. elegans, Anopheles gambiae, Ciona intestinalis, Fugu rubripes, Mus musculus, and Rattus norvegicus were obtained from NCBI, ENSEMBL and species specific databases (Additional file 1) [5358]. Human and all other mammalian sequences were identified by the BLASTP program from the sequence data at NCBI. Rat Kell was produced by concatenating the two partial sequences found in the R. norvegicus genome. M13 genes from D. pseudoobscura and Ap. mellifera were identified from raw genomic data using the SHARKhunt program [31]. SHARKhunt uses a search protocol that combines PSI-BLAST with profile Hidden Markov Model (HMM) techniques. The gene model employed was generated from fifteen M13 protein sequences from across the taxa examined here (Additional file 1). Some Apis mellifera sequences were also obtained from annotated entries in NCBI (Additional file 1).

Phylogenetic reconstruction of M13 genes

Protein sequences determined using the methods described above, were aligned using the program MUSCLE with default settings [32, 33]. To facilitate analysis, this alignment was edited using BioEdit software [59] to remove gaps and uninformative insertions. Three methods were used to reconstruct phylogenetic relationships of 111 NEP-like proteins. These were: neighbour-joining, maximum parsimony and maximum likelihood. Neighbour-joining [60] analysis was performed on the alignment described above using PAUP 4.0 [61] and was set to bootstrap the tree 1000 times (Additional file 4). Maximum parsimony analysis was performed using PAUP 4.0 [61]. Trees were generated by the random addition (ten replicates) of sequences from the alignment described above. After completion of trees, branches were swapped using the "tree-bisection and reconnection" method and the most parsimonious trees were saved. From this a consensus tree was generated and bootstrapped 1000 times (Additional file 5). Maximum likelihood analysis was carried out using the PROML function of PHYLIP [62]. Trees were generated using the PAM250 substitution score matrix and a gamma distribution of 1.78 as determined using ProtTest [63]. This model was used to generate 100 bootstrap replicates from which a consensus tree was generated. This tree was re-rooted using NOTUNG [64] to produce a tree that more closely resembled distance and parsimony trees (Additional file 6). These trees were used to generate a majority consensus tree using the consense function of PHYLIP [62]. For clarity an unrooted circular tree omitting bootstrap values (Figure 2) was produced using CTree and a rectangular cladogram incorporating bootstrap values was also produced (Additional file 3) [65].

Site-specific analysis of M13 sequences

Specific sites in M13 protein sequences were compared using the multiple sequence alignment described above. Sites of interest comprised residues previously identified as being important for catalysis and substrate binding by site-directed mutagenesis studies and from the crystal structure of human neprilysin [34, 35]. In the absence of structural data, only residues directly aligned with catalytically important regions of human neprilysin were considered and not adjacent residues that may form portions of binding sites in other enzymes.


  1. Turner AJ: Neprilysin. Handbook of Proteolytic Enzymes. Edited by: Barret AJ, Rawlings ND, Woessner JF. 2004, London , Elsevier academic press, 1: 419-425. 2nd

    Chapter  Google Scholar 

  2. Turner AJ, Isaac RE, Coates D: The neprilysin (NEP) family of zinc metalloendopeptidases: genomics and function. Bioessays. 2001, 23 (3): 261-269. 10.1002/1521-1878(200103)23:3<261::AID-BIES1036>3.0.CO;2-K.

    Article  CAS  PubMed  Google Scholar 

  3. Hemming ML, Patterson M, Reske-Nielsen C, Lin L, Isacson O, Selkoe DJ: Reducing amyloid plaque burden via ex vivo gene delivery of an Abeta-degrading protease: a novel therapeutic approach to Alzheimer disease. PLoS Med. 2007, 4 (8): e262-10.1371/journal.pmed.0040262.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Thomas JE, Rylett CM, Carhan A, Bland ND, Bingham RJ, Shirras AD, Turner AJ, Isaac RE: Drosophila melanogaster NEP2 is a new soluble member of the neprilysin family of endopeptidases with implications for reproduction and renal function. Biochem J. 2005, 386 (Pt 2): 357-366.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Carpentier M, Desgroseillers L, Boileau G: Neprilysin-2. Handbook of Proteolytic Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Esevier academic press, 1: 426-429. 2nd

    Chapter  Google Scholar 

  6. Ahn K, Johnson GD: Endothelin-converting enzyme 1. Handbook of Proteolytic Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Elsevier Academic Press, 1: 429-434. 2nd

    Chapter  Google Scholar 

  7. Yanagisawa H, Hammer RE, Richardson JA, Emoto N, Williams SC, Takeda S, Clouthier DE, Yanagisawa M: Disruption of ECE-1 and ECE-2 reveals a role for endothelin-converting enzyme-2 in murine cardiac development. J Clin Invest. 2000, 105 (10): 1373-1382. 10.1172/JCI7447.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Ahn K: Endothelin-converting enzyme 2. Handbook of Proteolytic Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Elsevier Academic Press, 1: 434-436. 2nd

    Chapter  Google Scholar 

  9. Valdenaire O: Endothelin-converting enzyme-like peptidase 1. Handbook of Proteolytic Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Elsevier Academic Press, 1: 436-437. 2nd

    Chapter  Google Scholar 

  10. Schweizer A, Valdenaire O, Koster A, Lang Y, Schmitt G, Lenz B, Bluethmann H, Rohrer J: Neonatal lethality in mice deficient in XCE, a novel member of the endothelin-converting enzyme and neutral endopeptidase family. J Biol Chem. 1999, 274 (29): 20450-20456. 10.1074/jbc.274.29.20450.

    Article  CAS  PubMed  Google Scholar 

  11. HypConsortium: A gene (PEX) with homologies to endopeptidases is mutated in patients with X-linked hypophosphatemic rickets. The HYP Consortium. Nat Genet. 1995, 11 (2): 130-136. 10.1038/ng1095-130.

    Article  Google Scholar 

  12. Boileau G: PHEX endopeptidase. Handbook of Proteolytic Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Elsevier Academic press, 1: 438-440. 2nd

    Chapter  Google Scholar 

  13. Campos M, Couture C, Hirata IY, Juliano MA, Loisel TP, Crine P, Juliano L, Boileau G, Carmona AK: Human recombinant endopeptidase PHEX has a strict S1' specificity for acidic residues and cleaves peptides derived from fibroblast growth factor-23 and matrix extracellular phosphoglycoprotein. Biochem J. 2003, 373 (Pt 1): 271-279. 10.1042/BJ20030287.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Lee S, Lin M, Mele A, Cao Y, Farmar J, Russo D, Redman C: Proteolytic processing of big endothelin-3 by the kell blood group protein. Blood. 1999, 94 (4): 1440-1450.

    CAS  PubMed  Google Scholar 

  15. Camara-Clayette V, Rahuel C, Lopez C, Hattab C, Verkarre V, Bertrand O, Cartron JP: Transcriptional regulation of the KEL gene and Kell protein expression in erythroid and non-erythroid cells. Biochem J. 2001, 356 (Pt 1): 171-180. 10.1042/0264-6021:3560171.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Redman C, Lee S: Kell blood-group protein. Handbook of Proteolytc Enzymes. Edited by: Barret A, Rawlings ND, Woessner JF. 2004, London , Elsevier Academic Press, 1: 438-440. 2nd

    Google Scholar 

  17. Shipp MA, Stefano GB, D'Adamio L, Switzer SN, Howard FD, Sinisterra J, Scharrer B, Reinherz EL: Downregulation of enkephalin-mediated inflammatory responses by CD10/neutral endopeptidase 24.11. Nature. 1990, 347 (6291): 394-396. 10.1038/347394a0.

    Article  CAS  PubMed  Google Scholar 

  18. Caselgrandi E, Kletsas D, Ottaviani E: Neutral endopeptidase-24.11 (NEP) deactivates PDGF- and TGF-beta-induced cell shape changes in invertebrate immunocytes. Cell Biol Int. 2000, 24 (2): 85-90. 10.1006/cbir.1999.0453.

    Article  CAS  PubMed  Google Scholar 

  19. Masler EP: In vitro metabolism of an insect neuropeptide by homogenates of the nematode Caenorhabditis elegans. J Helminthol. 2003, 77 (1): 43-48. 10.1079/JOH2002152.

    Article  CAS  PubMed  Google Scholar 

  20. Redmond DL, Knox DP, Newlands G, Smith WD: Molecular cloning and characterisation of a developmentally regulated putative metallopeptidase present in a host protective extract of Haemonchus contortus. Mol Biochem Parasitol. 1997, 85 (1): 77-87. 10.1016/S0166-6851(96)02812-5.

    Article  CAS  PubMed  Google Scholar 

  21. Jones BF, Hotez PJ: Molecular cloning and characterization of Ac-mep-1, a developmentally regulated gut luminal metalloendopeptidase from adult Ancylostoma caninum hookworms. Mol Biochem Parasitol. 2002, 119 (1): 107-116. 10.1016/S0166-6851(01)00409-1.

    Article  CAS  PubMed  Google Scholar 

  22. Sajid M, Isaac RE: Identification and properties of a neuropeptide-degrading endopeptidase (neprilysin) of Ascaris suum muscle. Parasitology. 1995, 111 ( Pt 5): 599-608.

    Article  CAS  Google Scholar 

  23. Spanier B, Sturzenbaum SR, Holden-Dye LM, Baumeister R: Caenorhabditis elegans Neprilysin NEP-1: an Effector of Locomotion and Pharyngeal Pumping. J Mol Biol. 2005

    Google Scholar 

  24. Isaac RE, Nassel DR: Identification and localization of a neprilysin-like activity that degrades tachykinin-related peptides in the brain of the cockroach, Leucophaea maderae, and locust, Locusta migratoria. J Comp Neurol. 2003, 457 (1): 57-66. 10.1002/cne.10561.

    Article  CAS  PubMed  Google Scholar 

  25. Macours N, Poels J, Hens K, Luciani N, De Loof A, Huybrechts R: An endothelin-converting enzyme homologue in the locust, Locusta migratoria: functional activity, molecular cloning and tissue distribution. Insect Mol Biol. 2003, 12 (3): 233-240. 10.1046/j.1365-2583.2003.00406.x.

    Article  CAS  PubMed  Google Scholar 

  26. Zhao X, Mita K, Shimada T, Okano K, Quan GX, Kanke E, Kawasaki H: Isolation and expression of an ecdysteroid-inducible neutral endopeptidase 24.11-like gene in wing discs of Bombyx mori. Insect Biochem Mol Biol. 2001, 31 (12): 1213-1219. 10.1016/S0965-1748(01)00069-8.

    Article  CAS  PubMed  Google Scholar 

  27. Wilson CL, Shirras AD, Isaac RE: Extracellular peptidases of imaginal discs of Drosophila melanogaster. Peptides. 2002, 23 (11): 2007-2014. 10.1016/S0196-9781(02)00188-2.

    Article  CAS  PubMed  Google Scholar 

  28. Zhu Y, Johnson TJ, Myers AA, Kanost MR: Identification by subtractive suppression hybridization of bacteria-induced genes expressed in Manduca sexta fat body. Insect Biochem Mol Biol. 2003, 33 (5): 541-559. 10.1016/S0965-1748(03)00028-6.

    Article  CAS  PubMed  Google Scholar 

  29. Aguilar R, Jedlicka AE, Mintz M, Mahairaki V, Scott AL, Dimopoulos G: Global gene expression analysis of Anopheles gambiae responses to microbial challenge. Insect Biochem Mol Biol. 2005, 35 (7): 709-719. 10.1016/j.ibmb.2005.02.019.

    Article  CAS  PubMed  Google Scholar 

  30. Bland ND, Thomas JE, Audsley N, Shirras AD, Turner AJ, Isaac RE: Expression of NEP2, a soluble neprilysin-like endopeptidase, during embryogenesis in Drosophila melanogaster. Peptides. 2007, 28 (1): 127-135. 10.1016/j.peptides.2006.08.032.

    Article  CAS  PubMed  Google Scholar 

  31. Pinney JW, Shirley MW, McConkey GA, Westhead DR: metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella. Nucl Acids Res. 2005, 33 (4): 1399-1409. 10.1093/nar/gki285.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5 (1): 113-10.1186/1471-2105-5-113.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Oefner C, D'Arcy A, Hennig M, Winkler FK, Dale GE: Structure of human neutral endopeptidase (neprilysin) complexed with phosphoramidon. Journal of Molecular Biology. 2000, 296 (2): 341-349. 10.1006/jmbi.1999.3492.

    Article  CAS  PubMed  Google Scholar 

  35. Crine P, Dion N, Boileau G: Endopeptidase-24.11. Cell-Surface Peptidases in health and disease. Edited by: Kenny AJ, Boustead CM. 1997, Oxford , BIOS Scientific Publishers LTD, 79-98.

    Google Scholar 

  36. Gaunt MW, Miles MA: An insect molecular clock dates the origin of the insects and accords with palaeontological and biogeographic landmarks. Mol Biol Evol. 2002, 19 (5): 748-761.

    Article  CAS  PubMed  Google Scholar 

  37. FlyBaseConsortium: The FlyBase database of the Drosophila genome projects and community literature. Nucl Acids Res. 2003, 31 (1): 172-175. 10.1093/nar/gkg094.

    Article  Google Scholar 

  38. Satoh N, Satou Y, Davidson B, Levine M: Ciona intestinalis: an emerging model for whole-genome analyses. Trends Genet. 2003, 19 (7): 376-381. 10.1016/S0168-9525(03)00144-6.

    Article  CAS  PubMed  Google Scholar 

  39. Ito M, Akai E, Izuka M, Segawa H, Kuwahata M, Miyamoto K: Cloning and characterization of three PHEX homologues in Drosophila. J Bone Miner Metab. 2004, 22 (1): 3-11. 10.1007/s00774-003-0440-8.

    Article  CAS  PubMed  Google Scholar 

  40. Bianchetti L, Oudet C, Poch O: M13 endopeptidases: New conserved motifs correlated with structure, and simultaneous phylogenetic occurrence of PHEX and the bony fish. Proteins. 2002, 47 (4): 481-488. 10.1002/prot.10075.

    Article  CAS  PubMed  Google Scholar 

  41. Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chinwalla A, Clarke L, Clee C, Coghlan A, Coulson A, D'Eustachio P, Fitch DHA, Fulton LA, Fulton RE, Griffiths-Jones S, Harris TW, Hillier LDW, Kamath R, Kuwabara PE, Mardis ER, Marra MA, Miner TL, Minx P, Mullikin JC, Plumb RW, Rogers J, Schein JE, Sohrmann M, Spieth J, Stajich JE, Wei C, Willey D, Wilson RK, Durbin R, Waterston RH: The Genome Sequence of Caenorhabditis briggsae: A Platform for Comparative Genomics. PLoS Biology. 2003, 1 (2): e45-10.1371/journal.pbio.0000045.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Brooks DR, Appleford PJ, Murray L, Isaac RE: An essential role in molting and morphogenesis of Caenorhabditis elegans for ACN-1, a novel member of the angiotensin-converting enzyme family that lacks a metallopeptidase active site. J Biol Chem. 2003, 278 (52): 52340-52346. 10.1074/jbc.M308858200.

    Article  CAS  PubMed  Google Scholar 

  43. Guo R, Rowe PS, Liu S, Simpson LG, Xiao ZS, Darryl Quarles LD: Inhibition of MEPE cleavage by Phex. Biochem Biophys Res Commun. 2002, 297 (1): 38-45. 10.1016/S0006-291X(02)02125-3.

    Article  CAS  PubMed  Google Scholar 

  44. Rowe PS, Garrett IR, Schwarz PM, Carnes DL, Lafer EM, Mundy GR, Gutierrez GE: Surface plasmon resonance (SPR) confirms that MEPE binds to PHEX via the MEPE-ASARM motif: a model for impaired mineralization in X-linked rickets (HYP). Bone. 2005, 36 (1): 33-46. 10.1016/j.bone.2004.09.015.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Voisin S, Rognan D, Gros C, Ouimet T: A Three-dimensional Model of the Neprilysin 2 Active Site Based on the X-ray Structure of Neprilysin: IDENTIFICATION OF RESIDUES INVOLVED IN SUBSTRATE HYDROLYSIS AND INHIBITOR BINDING OF NEPRILYSIN 2. J Biol Chem. 2004, 279 (44): 46172-46181. 10.1074/jbc.M407333200.

    Article  CAS  PubMed  Google Scholar 

  46. Sansom CE, Hoang MV, Turner AJ: Molecular modelling and site-directed mutagenesis of the active site of endothelin-converting enzyme. Protein Eng. 1998, 11 (12): 1235-1241. 10.1093/protein/11.12.1235.

    Article  CAS  PubMed  Google Scholar 

  47. Manzetti S: Computer modeling and nanosecond simulation of the enzyme-substrate complex of the common lymphoblastic leukemia antigen (neprilysin) indicates shared residues at the primary specificity pocket (S1') with matrix metalloproteases. Journal of Molecular Modeling. 2003, 9 (5): 348-354. 10.1007/s00894-003-0158-5.

    Article  CAS  PubMed  Google Scholar 

  48. Kortschak RD, Samuel G, Saint R, Miller DJ: EST Analysis of the Cnidarian Acropora millepora Reveals Extensive Gene Loss and Rapid Sequence Divergence in the Model Invertebrates. Current Biology. 2003, 13 (24): 2190-2195. 10.1016/j.cub.2003.11.030.

    Article  CAS  PubMed  Google Scholar 

  49. Coghlan A, Wolfe KH: Fourfold faster rate of genome rearrangement in nematodes than in Drosophila. Genome Res. 2002, 12 (6): 857-867. 10.1101/gr.172702.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Ranz JM, Casals F, Ruiz A: How malleable is the eukaryotic genome? Extreme rate of chromosomal rearrangement in the genus Drosophila. Genome Res. 2001, 11 (2): 230-239. 10.1101/gr.162901.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Hughes AL, R F: Differential loss of ancestral gene families as a source of genomic divergence in animals. PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES B-BIOLOGICAL SCIENCES. 2004, 271 (S107-S109):

  52. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. NCBI. []

  54. ENSEMBL. []

  55. Ciona genome project . []

  56. Harris TW, Stein LD: WormBase: methods for data mining and comparative genomics. Methods Mol Biol. 2006, 351: 31-50.

    PubMed  Google Scholar 

  57. Crosby MA, Goodman JL, Strelets VB, Zhang P, Gelbart WM: FlyBase: genomes by the dozen. Nucleic Acids Res. 2007, 35 (Database issue): D486-91. 10.1093/nar/gkl827.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Dehal P, Satou Y, Campbell RK, Chapman J, Degnan B, De Tomaso A, Davidson B, Di Gregorio A, Gelpke M, Goodstein DM, Harafuji N, Hastings KE, Ho I, Hotta K, Huang W, Kawashima T, Lemaire P, Martinez D, Meinertzhagen IA, Necula S, Nonaka M, Putnam N, Rash S, Saiga H, Satake M, Terry A, Yamada L, Wang HG, Awazu S, Azumi K, Boore J, Branno M, Chin-Bow S, DeSantis R, Doyle S, Francino P, Keys DN, Haga S, Hayashi H, Hino K, Imai KS, Inaba K, Kano S, Kobayashi K, Kobayashi M, Lee BI, Makabe KW, Manohar C, Matassi G, Medina M, Mochizuki Y, Mount S, Morishita T, Miura S, Nakayama A, Nishizaka S, Nomoto H, Ohta F, Oishi K, Rigoutsos I, Sano M, Sasaki A, Sasakura Y, Shoguchi E, Shin-i T, Spagnuolo A, Stainier D, Suzuki MM, Tassy O, Takatori N, Tokuoka M, Yagi K, Yoshizaki F, Wada S, Zhang C, Hyatt PD, Larimer F, Detter C, Doggett N, Glavina T, Hawkins T, Richardson P, Lucas S, Kohara Y, Levine M, Satoh N, Rokhsar DS: The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. Science. 2002, 298 (5601): 2157-2167. 10.1126/science.1080049.

    Article  CAS  PubMed  Google Scholar 

  59. BioEdit. []

  60. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.

    CAS  PubMed  Google Scholar 

  61. Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (and Other Methods) 4.0 B10. Sinauer associates, Inc. Sunderland, MA, USA.

  62. PHYLIP. []

  63. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.

    Article  CAS  PubMed  Google Scholar 

  64. Chen K, Durand D, Farach-Colton M: NOTUNG: a program for dating gene duplications and optimizing gene family trees. J Comput Biol. 2000, 7 (3-4): 429-447. 10.1089/106652700750050871.

    Article  CAS  PubMed  Google Scholar 

  65. Archer J, Robertson DL: CTree: comparison of clusters between phylogenetic trees made easy. Bioinformatics. 2007

    Google Scholar 

Download references


This work was supported by a Biotechnology and Biological Sciences Research Council studentship to NB and a Medical Research Council studentship to JP. Thanks to Dr. R Simpson for advice on phylogenetic reconstruction.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Nicholas D Bland.

Additional information

Authors' contributions

NB, JP and JET carried out the analyses. NB, AJT and REI planned the study and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Proteins used in this study. A comprehensive list of all proteins used in this study including accession numbers, source of the sequence and abbreviations used in the manuscript. Underlined sequences were used to generate the SHARKhunt gene model. (PDF 27 KB)


Additional file 2: Multiple sequence alignment of M13 protein sequences. A multiple sequence alignment of M13 proteins generated using MUSCLE and edited to remove uninformative regions. (FAS 83 KB)


Additional file 3: Consensus tree of phylogenetic analyses of M13 peptidases. A majority consensus cladogram of three methods of phylogenetic reconstruction of M13 proteins, including percentage bootstrap values in the order: neighbour-joining/maximum parsimony/maximum likelihood. (PDF 13 KB)


Additional file 4: Neighbour-joining analysis of M13 peptidases. Cladogram of neighbour-joining reconstruction of M13 proteins, including percentage bootstrap values. (PDF 10 KB)


Additional file 5: Maximum parsimony analysis of M13 peptidases. Cladogram of maximum parsimony reconstruction of M13 proteins, including percentage bootstrap values. (PDF 10 KB)


Additional file 6: Maximum likelihood analysis of M13 peptidases. Cladogram of maximum likelihood reconstruction of M13 proteins, including percentage bootstrap values. (PDF 11 KB)

Authors’ original submitted files for images

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bland, N.D., Pinney, J.W., Thomas, J.E. et al. Bioinformatic analysis of the neprilysin (M13) family of peptidases reveals complex evolutionary and functional relationships. BMC Evol Biol 8, 16 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: