Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

New mini- zincin structures provide a minimal scaffold for members of this metallopeptidase superfamily

BMC Bioinformatics volume 15, Article number: 1 (2014) Cite this article

Abstract

Background

The Acel_2062 protein from Acidothermus cellulolyticus is a protein of unknown function. Initial sequence analysis predicted that it was a metallopeptidase from the presence of a motif conserved amongst the Asp-zincins, which are peptidases that contain a single, catalytic zinc ion ligated by the histidines and aspartic acid within the motif (HEXXHXXGXXD). The Acel_2062 protein was chosen by the Joint Center for Structural Genomics for crystal structure determination to explore novel protein sequence space and structure-based function annotation.

Results

The crystal structure confirmed that the Acel_2062 protein consisted of a single, zincin-like metallopeptidase-like domain. The Met-turn, a structural feature thought to be important for a Met-zincin because it stabilizes the active site, is absent, and its stabilizing role may have been conferred to the C-terminal Tyr113. In our crystallographic model there are two molecules in the asymmetric unit and from size-exclusion chromatography, the protein dimerizes in solution. A water molecule is present in the putative zinc-binding site in one monomer, which is replaced by one of two observed conformations of His95 in the other.

Conclusions

The Acel_2062 protein is structurally related to the zincins. It contains the minimum structural features of a member of this protein superfamily, and can be described as a “mini- zincin”. There is a striking parallel with the structure of a mini-Glu-zincin, which represents the minimum structure of a Glu-zincin (a metallopeptidase in which the third zinc ligand is a glutamic acid). Rather than being an ancestral state, phylogenetic analysis suggests that the mini-zincins are derived from larger proteins.

Background

A metallopeptidase is a proteolytic enzyme that has one or two metal ions as an integral part of its catalytic machinery located within the active site. There are many families of metallopeptidases that bind a single zinc ion required for catalysis. The zinc ion is tetrahedrally co-ordinated by three residues from the peptidase and a water molecule that becomes activated to be the nucleophile in the catalytic reaction. In many of these families, but not all, the residues that ligate the zinc ion (referred to here as “ligands”) are two histidines within an HEXXH motif, and a third coordinating residue that is C-terminal to this motif, which can be a glutamic acid, a histidine or an aspartic acid. Metallopeptidases with an HEXXH motif are known as zincins. In a metallopeptidase such as thermolysin, the third zinc ligand is usually a glutamic acid, and these peptidases are known as Glu-zincins. In a metallopeptidase such as matrix metallopeptidase 1 (MMP1), the zinc ligands are the three histidines within an HEXXHXXGXXH motif. In MMP1 there is also an important region known as the Met-turn, in which there is a conserved methionine that structurally supports the active site. For this reason, metallopeptidases such as MMP1 are known as Met-zincins [1]. In some zincins the third zinc ligand may be an aspartic acid (within the motif HEXXHXXGXXD), and there is no Met-turn; these peptidases are known as Asp-zincins [2]. All of the zincins share structural similarities and in the MEROPS classification and database, the different families that can be recognized by sequence similarities are all included in clan MA. This clan is subdivided into three subclans. Subclan MA(E) containing the Glu-zincins, subclan MA(M) containing the Met-zincins and subclan MA(D) containing the Asp-zincins [3].

A zincin structure contains at least two subdomains: an N-terminal subdomain which includes the HEXXH motif, and a C-terminal subdomain that includes the third zinc ligand. The active site is therefore between the two subdomains. Within the zincins, the size of the C-terminal subdomain varies enormously, from being large in the matrix metallopeptidases (peptidase family M10), to being just a single helix, which is the case for snapalysin from Streptomyces lividans (peptidase family M7), which has one of the smallest sequences and structures in the clan [4].

Very recently, a minimal structure for a Glu-zincin has been determined [5], representing the smallest known member of this subclan so far discovered. No catalytic activity could be detected. The family has been provisionally assigned the name M95, but will not appear in the MEROPS database until peptidase activity has been experimentally confirmed. Lenart et al. [6] identified a number of protein families in which an HEXXH motif was conserved. One of these families was Pfam family PF06262 (DUF1025)[Pfam:PF06262], which includes atleast 400 sequences from bacteria. Members of this family have the Asp-zincin-like motif HEXXHXXGXXD. In this paper, we report the structure of a member of this family: the Acel_2062 protein from Acidothermus cellulolyticus, a cellulolytic thermophile found in hot springs [7]. The domain architecture represents a minimal structure for a zincin, with very little sequence beyond the third zinc ligand and an absence of the Met-turn.

Methods

Protein expression and purification

The American Type Culture Collection (ATCC) provided the genomic DNA used to clone Acel_2062 (ATCC Number: ATCC 43068). Protein production and crystallization of the Acel_2062 protein was carried out by standard JCSG protocols [8]. Clones were generated using the Polymerase Incomplete Primer Extension (PIPE) cloning method [9]. The gene encoding Acel_2062 (GenBank: YP_873820[GenBank:YP_873820]; UniProtKB: A0LWM4[UniProtKB:A0LWM4]) was synthesized with codons optimized for Escherichia coli expression (Codon Devices, Cambridge, MA) and cloned into plasmid pSpeedET, which encodes an expression and purification tag followed by a tobacco etch virus (TEV) protease cleavage site (MGSDKIHHHHHHENLYFQ/G) at the amino terminus of the full-length protein. Escherichia coli GeneHogs (Invitrogen) competent cells were transformed and dispensed on selective LB-agar plates. The cloning junctions were confirmed by DNA sequencing. Expression was performed in a selenomethionine-containing medium at 37°C. Selenomethionine was incorporated via inhibition of methionine biosynthesis [10], which does not require a methionine auxotrophic strain. At the end of fermentation, lysozyme was added to the culture to a final concentration of 250 μg/ml, and the cells were harvested and frozen. After one freeze/thaw cycle the cells were homogenized in lysis buffer [50 mM HEPES, 50 mM NaCl, 10 mM imidazole, 1 mM Tris(2-carboxyethyl)phosphine-HCl (TCEP), pH 8.0] and passed through a Microfluidizer (Microfluidics). The lysate was clarified by centrifugation at 32,500 x g for 30 minutes and loaded onto a nickel-chelating resin (GE Healthcare) pre-equilibrated with lysis buffer, the resin was washed with wash buffer [50 mM HEPES, 300 mM NaCl, 40 mM imidazole, 10% (v/v) glycerol, 1 mM TCEP, pH 8.0], and the protein was eluted with elution buffer [20 mM HEPES, 300 mM imidazole, 10% (v/v) glycerol, 1 mM TCEP, pH 8.0]. The eluate was buffer exchanged with TEV buffer [20 mM HEPES, 200 mM NaCl, 40 mM imidazole, 1 mM TCEP, pH 8.0] using a PD-10 column (GE Healthcare), and incubated with 1 mg of TEV protease per 15 mg of eluted protein for 2 hours at 20°–25°C followed by overnight at 4°C. The protease-treated eluate was passed over nickel-chelating resin (GE Healthcare) pre-equilibrated with HEPES crystallization buffer [20 mM HEPES, 200 mM NaCl, 40 mM imidazole, 1 mM TCEP, pH 8.0] and the resin was washed with the same buffer. The flow-through and wash fractions were combined and concentrated to 15.6 mg/ml by centrifugal ultrafiltration (Millipore) for crystallization trials.

Protein crystallization

The Acel_2062 protein was crystallized using the nanodroplet vapor diffusion method [11] with standard JCSG crystallization protocols [12]. Sitting drops composed of 100 nl protein solution mixed with 100 nl crystallization solution in a sitting drop format were equilibrated against a 50 μl reservoir at 277 K for 72 days prior to harvest. The crystallization reagent consisted of 24% polyethylene glycol 8000, 0.167 M calcium acetate, 0.1 M MES pH 6.17. Glycerol was added to a final concentration of 20% (v/v) as a cryoprotectant. Initial screening for diffraction was carried out using the Stanford Automated Mounting system (SAM) [13] at the Stanford Synchrotron Radiation Lightsource (SSRL, Menlo Park, CA). Data were collected at 3 wavelengths corresponding to the inflection(l1), high remote(l2) and peak energy(l3) of a selenium MAD (multi-wavelength anomalous dispersion) experiment at 100 K using a MARCCD 325 detector (Rayonix) at Stanford Synchrotron Radiation Lightsource (SSRL) beamline 9_2. Data processing was carried out using XDS [14] and the statistics are presented in Table 1. The structure was determined by the MAD method using programs SHELX [15] and autoSHARP [16], and refinement was carried out using REFMAC5 [17]. The structure was validated using the JCSG Quality Control server (http://smb.slac.stanford.edu/jcsg/QC).

Table 1 Summary of crystal parameters, data collection, and refinement statistics for Acel_2062 from acidothermus cellulolyticus 11b [PDB:3e11]
Full size table

Determination of oligomeric state

The oligomeric state of the Acel_2062 protein in solution was determined using a 0.8 × 30 cm2 Shodex Protein KW-803 size exclusion column (Thomson Instruments) [9] pre-calibrated with gel filtration standards (Bio-Rad). The mobile phase consisted of 20 mM Tris pH 8.0, 150 mM NaCl, and 0.02% (w/v) sodium azide. The apparent molecular weight was calculated using the Bio-Rad Gel Filtration Standard set (#151-1901) and a linear regression of log10MW.

Bioinformatics

To find homologues of the Acel_2062 protein, a Blastp search was conducted against the non-redundant protein sequence database at NCBI [21], using standard parameters. Structure diagrams were prepared using PyMol. Domain diagrams were taken from Pfam release 27 [22]. Secondary structure topology diagrams were generated by HERA [23] and downloaded from PDBSum website (http://www.ebi.ac.uk/pdbsum/). The alignment was prepared using MAFFT [24] and ESPript 2.2 (http://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi). PISA analysis [25] of the dimer interface was performed using the PDBe server at the European Bioinformatics Institute (http://www.ebi.ac.uk/msd-srv/prot_int/). The electrostatic surface was displayed using PyMol (http://www.pymol.org/) and a Delphi [26] embedded script kindly provided by Qingping Xu. Coot [27] was used to superimpose structures from the Protein Data Bank (PDB). Molecular graphics and analyses were performed with the UCSF Chimera package [28]. The theoretical pI was calculated using the Expasy website (http://web.expasy.org/compute_pi).

Figure 1
figure 1

Crystal structure. The crystal structure of the Acel_2062 protein from Acidothermus cellulolyticus, shown in magenta, reveals a single zincin-like domain. The putative zinc ligands (His95, His99 and Asp105), the putative catalytic Glu96, and Tyr113 are shown as sticks.

Full size image

Results and discussion

Structure description

The crystal structure of the Acel_2062 protein was determined to 1.8 Å resolution using the MAD phasing method. Atomic coordinates and experimental structure factors to 1.8 Å resolution [PDB:3e11] have been deposited in the Protein Data Bank (http://www.wwpdb.org, [29]). Data-collection, model and refinement statistics are summarized in Table 1. The final model includes two protein molecules (residues 1–113), two acetate molecules, two (presumed to be structural) calcium ions and 193 water molecules in the asymmetric unit. The calcium ions are near the centre of the dimer interface, and may be important for dimerization. The calcium ions are coordinated by Asp18, Glu38 and via waters by Glu15. No zinc was found in the structure, either because little zinc was present during purification and crystallization, the enzyme is in latent state or the protein is not an enzyme. The Matthews coefficient (V M ; [30]) is 2.07 Å3/Da and the estimated solvent content is 40.54%. The Ramachandran plot produced by MolProbity [31] shows that 98.0% of the residues are in favoured regions, with no outliers. The side-chain atoms of Glu22, Asp37, Glu43 and Glu106 on chain A and Glu15, Glu43 and Glu110 on chain B had poor electron density and were omitted from the model. For monomer A, the structure is composed of four helices, two 310 helices and three beta strands (see Figure 1). In monomer B, the N-terminus forms a fourth strand and the dimer is formed from this strand inserting into the beta sheet of monomer A.

The Acel_2062 protein was crystallized as a dimer, with the putative active sites at opposite ends of the dimer. PISA analysis of the structure indicates that the solvent-excluded surface area for the proposed dimer is ~639 Å2. From size exclusion chromatography, the molecular weight of the Acel_2062 protein in solution is estimated to be 26,824, which is a ratio of 2.09 over the expected molecular weight of the monomer (12,833) indicating that the protein exists as a dimer (Additional file 1: Figure S1).

The crystal structure of Acel_2062 protein has a minimal zincin-like fold because it retains a three-stranded mixed beta-sheet (rather than the five-stranded beta-sheet of many zincins), the loops are much shorter and the overall sequence length of all members of this family (~110 aa) is significantly shorter than the average for a matrix metalloprotease (MMP)-like domain (~160 aa). The distant homology prediction program FFAS [32] recognizes this similarity, with a marginal statistical significance (Z-score of −8.9 as compared to −9.5 as the significant threshold), suggesting that DUF1025 family is distantly related to metalloproteases.

There is no signal peptide, and the Acel_2062 protein is presumably intracellular.

Putative active site

From the presence of the HEXXHXXGXXD motif, the potential zinc ligands in the Acel_2062 protein are predicted to be His95, His99 and Asp105, and Glu96 is predicted to be a catalytic residue. In the crystal structure (PDB:3E11), the Glu96 is hydrogen-bonded to five water molecules in both monomers. Conservation of the active site residues is shown in Figure 2. Of the two active sites in both monomers A and B, the one in B is empty and the other in A is occupied by a single water molecule which is coordinated by His95 (at 2.7 Å), His99 (at 3.3 Å) and Asp105 (at 2.9 Å), exactly where the zinc ion would be expected to be. Delphi calculations show that the entire pocket is very acidic. There are five negatively charged residues, plus the carboxyl group from the C-terminus in that area. The surface electrostatic potential is shown in Figure 3. The molecule is negatively charged overall (with 15 Glu and 11 Asp compared to one Lys and eight Arg), with a theoretical pI calculated to be 4.29. This may be an adaptation to the acidic hot spring environment in which Acidothermus cellulolyticus lives.

Figure 2
figure 2

Sequence alignment of the Acel_2062 protein and a selection of its homologues. The UniProt accession and the range of the peptidase domain are shown on the left. The zincin motif is boxed in red. Conserved residues shown in white text highlighted in red. Key to sequences (ordered locus name, species): A0LWM4[UniProt:A0LWM4] (Acel_2062, Acidothermus cellulolyticus), F2NA85[UniProt:F2NA85] (Corgl_0144, Coriobacterium glomerans), C0ACH6[UniProt:C0ACH6] (ObacDRAFT_6101, Diplosphaera colitermitum), E6JBL2[UniProt:E6JBL2] (ES5_13138, Dietzia cinnamea), A5P6T1[UniProt:A5P6T1] (ED21_26213, Erythrobacter sp. SD-21), A1B0F0[UniProt:A1B0F0] (Pden_0883, Paracoccus denitrificans), E8L598[UniProt:E8L598] (Met49242DRAFT_2641, Methylocystis sp. ATCC 49242), Q1PUY8[UniProt:Q1PUY8] (kustc0300, Candidatus Kuenenia stuttgartiensis), F8CQR3[UniProt:F8CQR3] (LILAB_30480, Myxococcus fulvus), D0LZA2[UniProt:D0LZA2] (Hoch_3865, Haliangium ochraceu), H8MHL2[UniProt:H8MHL2] (COCOR_02006, Corallococcus coralloides), Q08VI2[UniProt: Q08VI2] (STAUR_2801, Stigmatella aurantiaca), E2T6Y2[UniProt:E2T6Y2] (TMBG_02375, Mycobacterium tuberculosis), Q5SLR6[UniProt:Q5SLR6] (TTHA0227, Thermus thermophilus), F2NR95[UniProt:F2NR95] (Marky_2224, Marinithermus hydrothermalis), F9MVY9[UniProt:F9MVY9] (HMPREF9130_1347, Peptoniphilus sp. oral taxon 375 str. F0436), D3PNK3[UniProt:D3PNK3] (Mrub_0627, Meiothermus ruber), B0TE39[UniProt:B0TE39] (Helmi_16090, Heliobacterium modesticaldum), E4U5S7[UniProt:E4U5S7] (Ocepr_0229, Oceanithermus profundus).

Full size image
Figure 3
figure 3

Electrostatic surface of the Acel_2062 protein. The electrostatic surface potential representation of the active pocket was made with a Delphi embedded JCSG script, (kindly provided by Qingping Xu) and displayed using Pymol. Negative electrostatic potential predominates over the entire surface of Acel_2062. The putative active site residues (Glu96 and Tyr113) and zinc ligands are (His95, His99 and Asp105) are shown in sticks. The color scale is in units of kT/e ranging from −10 to +10.

Full size image

In the Acel_2062 protein, and other members of the family, the HEXXH motif is very close to the C-terminus. Like Glu- and Asp-zincins, there is no Met-turn. The first His (His95) of the HEXXH motif exists in two conformations in monomer B. In the more stable conformation, His95 is electrostatically bound to the carboxyl group of the C-terminal residue (Tyr113). This is the situation also found in monomer A, so it is possible that the water (and presumably also the zinc) only bind when His95 and Tyr113 interact. The structure of the Acel_2062 protein was superimposed upon that of the reprolysin BaP1 peptidase from the snake Bothrops asper [PDB:2w14] and this clearly shows that Tyr113 occupies a similar position to the methionine of the Met-turn in the BaP1 peptidase (see Figure 4). So although there is no Met-turn, Tyr113 may compensate for it.

Figure 4
figure 4

Superposition of the active sites of the Acel_2062 protein and Bap1 peptidase from Bothrops asper . The comparison between the active centres containing the HEEXH motif; Coot ssm superposition of the two PDB entries (3E11 and 2 W14), has been displayed using Chimera; 3E11-A monomer is shown in turquoise, B in pink, 2 W14 is shown in magenta.

Full size image

Sequence similarities

Over 500 homologues of the Acel_2062 protein were found from the Blastp search. In 80 of these proteins, the HEXXHXXGXXD motif is not conserved. The third zinc ligand has been replaced in many of these homologues, often with glutamic acid, which is the third zinc ligand in Glu-zincins.

All of the homologues are from bacteria belonging to seven different phyla. Most homologues come from species in the phylum Firmicutes (294), which are Gram-positive bacteria. There are 185 homologues from species in the phylum Proteobacteria, fourteen from Chloroflexi, five from Planctomycetes, three from Verrucomicrobia, and one each from species in the phyla Caldiserica and Nitrospirae. Most species have only one homologue, but Stigmatella aurantiaca has two, though only one has the third zinc ligand conserved.

The different Pfam domain architectures for members of this family are shown in Figure 5. The vast majority of homologues have the simple domain architecture of the Acel_2062 protein. Seventeen homologues include tetratricopeptide repeats (TPR), which mediate protein-protein interactions and the assembly of multiprotein complexes [33]. A TPR repeat motif consists of several tandem repeats of a 34-residue sequence. Eight proteins with TPR repeats are predicted to have signal peptides and are presumably secreted. A homologue from Mycobacterium tuberculosis possesses an N-terminal transmembrane domain and a domain found in proteins known to be important for septum formation during spore formation [34].

Figure 5
figure 5

Domain architectures for the Acel_2062 protein and its homologues. The different domain architectures for proteins containing an M94 metallopeptidase domain is shown. The UniProt identifier and source organism are given for an example of each domain architecture. Key to domains: DUF1025 is the metallopeptidase domain;TPR, tetratricopeptide repeats; TM, transmembrane region; Septum form, a domain found in proteins predicted to play a role in septum formation during cell division; S, signal peptide.

Full size image

Structural similarities

Amongst known metallopeptidases, DALI analysis shows that the Acel_2062 protein structure is most similar to that of a Met-zincin from the archaean Methanocorpusculum labreanum (peptidase family M54; archaelysins or archaemetzincins [PDB:3lmc]). The Acel_2062 protein structure is also similar to the carboxy-terminal domain (residues 511 to 624) of Escherischia coli HtpG/Hsp90 protein, which is a chaperone protein. This C-terminal domain is important for dimerization, but the mechanism of dimerization via the C-terminal helices is completely different to that of the Acel_2062 dimer. Two of the helices and the beta sheet can be superimposed, and the beta strands run in the same direction. The relationship between Hsp90 and a zincin has not previously been recognized in either the SCOP [35] or CATH [36] databases, and suggests a common evolutionary origin. These structural relationships cover the entire Acel_2062 protein sequence.

The Met-turn is important in all Met-zincins because the methionine is crucial for structurally stabilizing the active site. Even in snapalysins from Streptomyces, which also have short sequences, there is a conserved methionine approximately ten residues C-terminal to the third zinc ligand. In the Acel_2062 protein although there are no residues that correspond to the C-terminal 36 residues of the archaemetzincin, which includes the essential Met168, the C-terminal Tyr113 occupies a similar position to the methionine. Although it is tempting to suggest that this tyrosine performs a similar role, it should be noted that the tyrosine is present in only 75 homologues in the family and replaced by tryptophan in a further 220 homologues and by phenylalanine in a homologue from Rhodococcus sp. AW25M09. In a further 221 homologues, the full-length sequence falls short of Tyr113. In astacin, a tyrosine (Tyr149) has been shown to be an important residue, contributing to transition state binding, which can be replaced with much reduced efficiency by phenylalanine [37]. If proven to be a peptidase, the Acel_2062 protein and its homologues would form family M94 in MEROPS.

The structure of another putative metallopeptidase, the TTHA0227 protein from Thermus thermophilus ([PDB:2ejq]; unpublished), is similar. The TTHA0227 protein is also a dimer but has been crystallized with a single magnesium ion, and the structure also lacks any zinc ions. This is also an acidic protein with a theoretical pI calculated to be 4.53. The TTHA0227 protein also contains the Asp-zincin metal-binding motif, and the potential zinc coordinating residues are His95, His99 and Asp109. Chain A is shorter at the C-terminus, lacking residues 109–130, which means that the third zinc ligand is missing. The putative catalytic residue is Glu96. There is a four-residue insert preceding the essential glycine (Gly106, Additional file 2: Figure S2). In chain A, the final helix is continuous, whereas in chain B, there is a turn and there are two, opposing helices. This second helix does not superimpose with the final helix in the 3E11 structure, because it is pointing in the opposite direction so that the faces of the helices that oppose each other in 2EJQ are different from those in 3E11 (Figure 6).

Figure 6
figure 6

Superposition of structures for the Acel_2062 and TTHA0227 proteins. The structure for the Acel_2062 protein (chain A from PDB entry 3E11) is shown in magenta, and the structure for the TTHA0227 protein (chain B from 2EJQ) is shown in cyan. The putative zinc ligands (His95, His99, Asp105 and Asp108) and active site residues (Glu96 and Tyr113) are shown as sticks. Note that the terminal helices do not superimpose, and that Asp108 from the 2EJQ structure is not interacting with the other putative zinc ligands.

Full size image

There are several orthologues of the TTHA0227 protein which have the insert within the Asp-zincin motif. In the homologues from Oceanithermus profundus and Meiothermus ruber the glycine (Gly106), which in Met-zincins is important for the turn that permits the zinc ligands to face one another, is replaced by tryptophan. It had been thought that only a glycine was acceptable in this position [38], although several bacterial homologues of pappalysin, family M43, have asparagine at this position, including a homologue from Methanosarcina acetivorans for which the crystal structure has been solved [39], and a homologue from Cytophaga hutchinsonii (Chut1718 gene product) has threonine at this position. Comparisons of the structural topologies of various Met-zincins and the mini-Glu-zincin are shown in Figure 7.

Figure 7
figure 7

Secondary structure topologies for the Acel_2062 protein and other metallopeptidases. Helices are shown as red barrels, strands are shown as pink arrows and connecting regions of random coil and turns are shown as thin blue arrows. “N” and “C” indicate the N- and C-terminus. Numbers indicate the residue range of each structural element. The positions of the zinc ligands and catalytic Glu are indicated by the letters “H”, “D” and “E”. Some images have been rotated and/or flipped to that these active site residues are in equivalent positions. Key to structures: 3E11[PDB:3e11], mini- zincin, Acidothermus cellulolyticus; 4JIU[PDB:4jiu], mini-Glu-zincin, Pyrococcus abyssi; 3LMC[PDB:3lmc], archaelysin, Methanocorpusculum labreanum; 1KUH[PDB:1kuh], snapalysin, Streptomyces caespitosus; 2EJQ[PDB:2ejq], mini- zincin, Thermus thermophilus; 1AST[PDB:1ast], astacin, Astacus astacus.

Full size image

The structure of the Acel_2062 protein represents the minimum sequence known for a zincin domain. The question remains: is this a situation that has developed within this family, or is it a relic of the ancestral zincin gene? One way to answer this question is to look at the species distribution of peptidases from the various families in the clan. Table 2 shows the phyletic distribution of all families within clan MA. The number of phyla within each of the three superkingdoms (Archaea, Bacteria, Eukaryota) containing at least one homologue within each family is shown. Amongst the Glu-zincins, families M1, M3, M32 and M48 are widely distributed in phyla from all three superkingdoms (M41 is also widely distributed in bacteria and eukaryotes, but is absent from archaea). The mini-Glu-zincins, from family M95, have a much narrower distribution and are absent from eukaryotes. These observations imply that the last common universal ancestor most likely possessed a homologue from each of these families, and that the larger Glu-zincin structure is the ancestral state. The distribution of Asp- and Met-zincins is more restricted in all families, and the only family that is well represented in all domains of life is family M54, the archaelysins. There are homologues from archaea (93 species), bacteria (47 species) and eukaryotes (69 species), though these are found in less than half of the bacterial and eukaryote phyla. It is likely that an archaelysin most closely represents the ancestral Met-zincin structure. Archaelysin possesses the Met-turn [40, 41], so the implication is that the Met-turn has been lost from an ancestor of family M94 and functionally replaced by a C-terminal aromatic residue (tyrosine or phenylalanine). The much narrower distribution of members of M94 supports the hypothesis that the family is a more recent development.

Table 2 Phyletic distribution of experimentally confirmed and hypothetical peptidase families in MEROPS clan MA
Full size table

Conclusions

The Acel_2062 protein from Acidothermus cellulolyticus is a protein of unknown function, but was predicted to be a metallopeptidase from the presence of a motif (HEXXHXXGXXD) conserved amongst the Asp-zincins, which contain a single, catalytic zinc ion ligated by the histidines and aspartic acid within the motif. The tertiary structure of the Acel_2062 protein was determined by the Joint Center for Structural Genomics, and confirmed the presence of a single, zincin-like metallopeptidase-like domain. In our crystallographic model there are two molecules in the asymmetric unit and from size-exclusion chromatography, the protein dimerizes in solution. A water molecule is present in the putative zinc-binding site in one monomer, which is replaced by one of two observed conformations of His95 in the other. The C-terminal Tyr113 may be important for stabilizing the putative active site. Additional experimentation would be required to prove that the Acel_2062 protein is a metallopeptidase.

Although the Acel_2062 protein is structurally related to the zincins, it contains the minimum structural features of a member of this protein superfamily, and can be described as a “mini- zincin”. There is a striking parallel with the structure of a mini-Glu-zincin, which represents the minimum structure of a Glu-zincin (a metallopeptidase in which the third zinc ligand is a glutamic acid). Rather than being an ancestral state, phylogenetic analysis suggests that the mini-zincins are derived from larger proteins.

References

  1. Stöcker W, Grams F, Baumann U, Reinemer P, Gomis-Rüth FX, McKay DB, Bode W: The metzincins–topological and sequential relations between the astacins, adamalysins, serralysins, and matrixins (collagenases) define a superfamily of zinc-peptidases. Protein Sci. 1995, 4: 823-840.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Fushimi N, Ee CE, Nakajima T, Ichishima E: Aspzincin, a family of metalloendopeptidases with a new zinc-binding motif. Identification of new zinc-binding sites (His(128), His(132), and Asp(164)) and three catalytically crucial residues (Glu(129), Asp(143), and Tyr(106)) of deuterolysin from Aspergillus oryzae by site-directed mutagenesis. J Biol Chem. 1999, 274: 24195-24201. 10.1074/jbc.274.34.24195.

    Article  PubMed  CAS  Google Scholar 

  3. Rawlings ND, Barrett AJ: Evolutionary families of metallopeptidases. Methods Enzymol. 1995, 248: 183-228.

    Article  PubMed  CAS  Google Scholar 

  4. Kurisu G, Kinoshita T, Sugimoto A, Nagara A, Kai Y, Kasai N, Harada S: Structure of the zinc endoprotease from Streptomyces caespitosus. J Biochem. 1997, 121 (2): 304-308. 10.1093/oxfordjournals.jbchem.a021587.

    Article  PubMed  CAS  Google Scholar 

  5. Lopéz-Pelegrín M, Cerdà-Costa N, Martínez-Jiménez F, Cintas-Pedrola A, Canals A, Peinado JR, Marti-Renom MA, Lopéz-Otín C, Arolas JL, Gomis-Rüth FX: A novel family of soluble minimal scaffolds provides structural insight into the catalytic domains of integral-membrane metallopeptidases. J Biol Chem. 2013,  -in press

    Google Scholar 

  6. Lenart A, Dudkiewicz M, Grynberg M, Pawlowski K: CLCAs - a family of metalloproteases of intriguing phylogenetic distribution and with cases of substituted catalytic sites. PLoS One. 2013, 8: e62272-10.1371/journal.pone.0062272.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  7. Barabote RD, Xie G, Leu DH, Normand P, Necsulea A, Daubin V, Médigue C, Adney WS, Xu XC, Lapidus A, Parales RE, Detter C, Pujic P, Bruce D, Lavire C, Challacombe JF, Brettin TS, Berry AM: Complete genome of the cellulolytic thermophile Acidothermus cellulolyticus 11B provides insights into its ecophysiological and evolutionary adaptations. Genome Res. 2009, 19 (6): 1033-1043. 10.1101/gr.084848.108.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  8. Elsliger MA, Deacon AM, Godzik A, Lesley SA, Wooley J, Wüthrich K, Wilson IA: The JCSG high-throughput structural biology pipeline. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2010, 66: 1137-1142. 10.1107/S1744309110038212.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  9. Klock HE, Koesema EJ, Knuth MW, Lesley SA: Combining the polymerase incomplete primer extension method for cloning and mutagenesis with microscreening to accelerate structural genomics efforts. Proteins. 2008, 71: 982-994. 10.1002/prot.21786.

    Article  PubMed  CAS  Google Scholar 

  10. Van Duyne GD, Standaert RF, Karplus PA, Schreiber SL, Clardy J: Atomic structures of the human immunophilin FKBP-12 complexes with FK506 and rapamycin. J Mol Biol. 1993, 229: 105-124. 10.1006/jmbi.1993.1012.

    Article  PubMed  CAS  Google Scholar 

  11. Santarsiero BD, Yegian DT, Lee CC, Spraggon G, Gu J, Scheibe D, Uber DC, Cornell EW, Nordmeyer RA, Kolbe WF, Jin J, Jones AL, Jaklevic JM, Schultz PG, Stevens RC: An approach to rapid protein crystallization using nanodroplets. J Appl Crystallogr. 2002, 35: 278-281. 10.1107/S0021889802001474.

    Article  CAS  Google Scholar 

  12. Lesley SA, Kuhn P, Godzik A, Deacon AM, Mathews I, Kreusch A, Spraggon G, Klock HE, McMullan D, Shin T, Vincent J, Robb A, Brinen LS, Miller MD, McPhillips TM, Miller MA, Scheibe D, Canaves JM, Guda C, Jaroszewski L, Selby TL, Elsliger MA, Wooley J, Taylor SS, Hodgson KO, Wilson IA, Schultz PG, Stevens RC: Structural genomics of the Thermotoga maritima proteome implemented in a high-throughput structure determination pipeline. Proc Natl Acad Sci U S A. 2002, 99: 11664-11669. 10.1073/pnas.142413399.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  13. Cohen AE, Ellis PJ, Miller MD, Deacon AM, Phizackerley RP: An automated system to mount cryo-cooled protein crystals on a synchrotron beamline, using compact sample cassettes and a small-scale robot. J Appl Crystallogr. 2002, 35: 720-726. 10.1107/S0021889802016709.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  14. Kabsch W: XDS. Acta Crystallogr Sect D Biol Crystallogr. 2010, 66: 125-132. 10.1107/S0907444909047337.

    Article  CAS  Google Scholar 

  15. Sheldrick GM: A short history of SHELX. Acta Crystallogr Sect A Found Crystallogr. 2008, 64: 112-122. 10.1107/S0108767307043930.

    Article  CAS  Google Scholar 

  16. Vonrhein C, Bricogne G: Exploiting structure similarity in refinement: automated NCS and target-structure restraints in BUSTER. Acta Crystallogr Sect D Biol Crystallogr. 2012, 68: 368-380. 10.1107/S0907444911056058.

    Article  Google Scholar 

  17. Winn MD, Isupov MN, Murshudov GN: Use of TLS parameters to model anisotropic displacements in macromolecular refinement. Acta Crystallogr Sect D Biol Crystallogr. 2001, 57: 122-133. 10.1107/S0907444900014736.

    Article  CAS  Google Scholar 

  18. Diederichs K, Karplus PA: Improved R-factors for diffraction data analysis in macromolecular crystallography. Nature Struct Biol. 1997, 4: 269-275. 10.1038/nsb0497-269.

    Article  PubMed  CAS  Google Scholar 

  19. Weiss MS: Global indicators of X-ray data quality. J Appl Cryst. 2001, 34: 130-135. 10.1107/S0021889800018227.

    Article  CAS  Google Scholar 

  20. Cruickshank DW: Remarks about protein structure precision. Acta Crystal Sect D, Biol Cryst. 1999, 55: 583-601. 10.1107/S0907444998012645.

    Article  CAS  Google Scholar 

  21. NCBI Resource Coordinators: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2013, 41: D8-D20.

    Article  PubMed Central  Google Scholar 

  22. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer EL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40: D290-D301. 10.1093/nar/gkr1065.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  23. Hutchinson EG, Thornton JM: HERA–a program to draw schematic diagrams of protein secondary structures. Proteins. 1990, 8: 203-212. 10.1002/prot.340080303.

    Article  PubMed  CAS  Google Scholar 

  24. Katoh K, Standley DM: MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013, 30: 772-780. 10.1093/molbev/mst010.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  25. Krissinel E, Henrick K: Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007, 372: 774-797. 10.1016/j.jmb.2007.05.022.

    Article  PubMed  CAS  Google Scholar 

  26. Li L, Li C, Sarkar S, Zhang J, Witham S, Zhang Z, Wang L, Smith N, Petukh M, Alexov E: DelPhi: a comprehensive suite for DelPhi software and associated resources. BMC Biophys. 2012, 4: 9-

    Article  Google Scholar 

  27. Emsley P, Lohkamp B, Scott WG, Cowtan K: Features and development of Coot. Acta Crystallogr D Biol Crystallogr. 2010, 66: 486-501. 10.1107/S0907444910007493.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  28. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera–a visualization system for exploratory research and analysis. J Comput Chem. 2004, 25: 1605-1612. 10.1002/jcc.20084.

    Article  PubMed  CAS  Google Scholar 

  29. Rose PW, Bi C, Bluhm WF, Christie CH, Dimitropoulos D, Dutta S, Green RK, Goodsell DS, Prlic A, Quesada M, Quinn GB, Ramos AG, Westbrook JD, Young J, Zardecki C, Berman HM, Bourne PE: The RCSB protein data bank: new resources for research and education. Nucleic Acids Res. 2013, 41: D475-D482. 10.1093/nar/gks1200.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  30. Matthews BW: Solvent content of protein crystals. J Mol Biol. 1968, 33: 491-497. 10.1016/0022-2836(68)90205-2.

    Article  PubMed  CAS  Google Scholar 

  31. Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, Murray LW, Richardson JS, Richardson DC: MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallographica Sect D Biol Crystallogr. 2010, 66: 12-21. 10.1107/S0907444909042073.

    Article  CAS  Google Scholar 

  32. Jaroszewski L, Li Z, Cai XH, Weber C, Godzik A: FFAS server: novel features and applications. Nucleic Acids Res. 2011, 39: W38-W44. 10.1093/nar/gkr441.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  33. D’Andrea LD, Regan L: TPR proteins: the versatile helix. Trends Biochem Sci. 2003, 28: 655-662. 10.1016/j.tibs.2003.10.007.

    Article  PubMed  Google Scholar 

  34. Slayden RA, Knudson DL, Belisle JT: Identification of cell cycle regulators in Mycobacterium tuberculosis by inhibition of septum formation and global transcriptional analysis. Microbiology. 2006, 152: 1789-1797. 10.1099/mic.0.28762-0.

    Article  PubMed  CAS  Google Scholar 

  35. Andreeva A, Howorth D, Chandonia JM, Brenner SE, Hubbard TJ, Chothia C, Murzin AG: Data growth and its impact on the SCOP database: new developments. Nucleic Acids Res. 2008, 36: D419-D425.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  36. Sillitoe I, Cuff AL, Dessailly BH, Dawson NL, Furnham N, Lee D, Lees JG, Lewis TE, Studer RA, Rentzsch R, Yeats C, Thornton JM, Orengo CA: New functional families (FunFams) in CATH to improve the mapping of conserved functional sites to 3D structures. Nucleic Acids Res. 2013, 41: D499-D507. 10.1093/nar/gks1266.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Yiallouros I, Grosse Berkhoff E, Stöcker W: The roles of Glu93 and Tyr149 in astacin-like zinc peptidases. FEBS Lett. 2000, 484: 224-228. 10.1016/S0014-5793(00)02163-3.

    Article  PubMed  CAS  Google Scholar 

  38. Gomis-Rüth FX: Structural aspects of the metzincin clan of metalloendopeptidases. Mol Biotechnol. 2003, 24: 157-202. 10.1385/MB:24:2:157.

    Article  PubMed  Google Scholar 

  39. Garcia-Castellanos R, Tallant C, Marrero A, Sola M, Baumann U, Gomis-Rüth FX: Substrate specificity of a metalloprotease of the pappalysin family revealed by an inhibitor and a product complex. Arch Biochem Biophys. 2007, 457: 57-72. 10.1016/j.abb.2006.10.004.

    Article  PubMed  CAS  Google Scholar 

  40. Waltersperger SM, Widmer C, Baumann U: Crystal structure of archaemetzincin AmzA from Methanopyrus kandleri at 1.5 A resolution. Proteins. 2010, 78: 2720-

    PubMed  CAS  Google Scholar 

  41. Graef C, Schacherl M, Waltersperger S, Baumann U: Crystal structures of archaemetzincin reveal a moldable substrate-binding site. PLoS One. 2012, 7: 43863-10.1371/journal.pone.0043863.

    Article  Google Scholar 

Download references

Acknowledgements

We are grateful to the Sanford Burnham Medical Research Institute for hosting the DUF Annotation Jamboree in June 2013 that allowed the authors to collaborate on this work. We would like to thank all the participants of this workshop for their intellectual contributions to this work: L. Aravind, Alex Bateman, Debanu Das, Robert D. Finn, Adam Godzik, William Hwang, Lukasz Jaroszewski, Alexey Murzin, Padmaja Natarajan, Daniel Rigden, Mayya Sedova, Anna Sheydina, John Wooley. We thank the members of the JCSG high-throughput structural biology pipeline for their contribution to this work. This work was supported in part by National Institutes of Health grant U54 GM094586 from the NIGMS Protein Structure Initiative to the Joint Center for Structural Genomics; intramural funds of the National Library of Medicine, USA, to LA; NIH grant R01GM101457 to AG; Howard Hughes Medical Institute to RDF; and Wellcome Trust grant WT077044/Z/05/Z for funding for open access charges. Portions of this research were carried out at the Stanford Synchrotron Radiation Lightsource, a Directorate of SLAC National Accelerator Laboratory and an Office of Science User Facility operated for the U.S. Department of Energy Office of Science by Stanford University. The SSRL Structural Molecular Biology Program is supported by the DOE Office of Biological and Environmental Research, and by the National Institutes of Health, National Institute of General Medical Sciences (including P41GM103393). The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of NIGMS, NCRR or NIH.

Author information

Authors and Affiliations

  1. Joint Center for Structural Genomics, La Jolla, CA, 92037, USA

    Christine B Trame

  2. Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, 94025, USA

    Christine B Trame & Herbert L Axelrod

  3. Sandford-Burnham Institute, La Jolla, CA, 92037, USA

    Yuanyuan Chang

  4. Wellcome Trust Sanger Institute, Hinxton, Cambridgeshire, CB10 1SA, UK

    Ruth Y Eberhardt, Penelope Coggill, Marco Punta & Neil D Rawlings

  5. European Molecular Biology Laboratory, European Bioinformatics Institute, Hinxton, Cambridgeshire, CB10 1SD, UK

    Ruth Y Eberhardt, Penelope Coggill, Marco Punta & Neil D Rawlings

Authors
  1. Christine B Trame
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yuanyuan Chang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Herbert L Axelrod
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ruth Y Eberhardt
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Penelope Coggill
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Marco Punta
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Neil D Rawlings
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Neil D Rawlings.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CBT performed X-ray structure determination and prepared some of the figures; NDR, YC, HLA, RYE, PC and MP analysed the sequence-structure-function relationships and prepared the manuscript, tables and the other figures. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2013_6243_MOESM1_ESM.png

Additional file 1: Figure S1.: Molecular weight determination by size-exclusion chromatography. A). Elution profile for Bio-Rad Gel Filtration Standard set (#151-1901) comprising vitamin B12 (1,350 Da), horse myoglobin (17 kDa), chicken ovalbumin (44 kDa), bovine gamma-globulin (158 kDa) and bovine thyroglobulin (670 kDa). B) Elution profile for the Acel_2062 protein. (PNG 57 KB)

12859_2013_6243_MOESM2_ESM.png

Additional file 2: Figure S2.: Crystal structure of the TTHA0227 protein from Thermus thermophilus. The PDB entry 2EJQ[PDB:2ejq] is a dimer, and because different elements are missing from both monomers, the structure of both monomers is shown. Chain A is shown in beige, chain B in cyan. Residues described below and in the text are labelled and shown as sticks. The monomers differ in terms of the residues that cannot be resolved. Chain A lacks the C-terminal residues Asp109-Gly130, and chain B lacks residues Pro54-Leu64 as well as the C-terminal Gly-Glu-Gly residues. Also in 2EJQ, Asp109 is too far away from the other potential zinc ligands to be a ligand itself. There are other potential zinc ligands in residues 110–130, including Glu113, Asp114, and Asp119. Only Asp119 is close enough to the imidazolium rings of the histidines to act as the third zinc ligand. Unfortunately, this Asp119 is poorly conserved, whereas Asp109 is well conserved. Because Asp109 and the final helix are close to the dimer interface, the structure here may be distorted because of the dimerization. (PNG 33 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), 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

Trame, C.B., Chang, Y., Axelrod, H.L. et al. New mini- zincin structures provide a minimal scaffold for members of this metallopeptidase superfamily. BMC Bioinformatics 15, 1 (2014). https://doi.org/10.1186/1471-2105-15-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-15-1

Keywords

BMC Bioinformatics

ISSN: 1471-2105

Contact us