Structural Study of Metal Binding and Coordination in Ancient Metallo-β-Lactamase PNGM-1 Variants.

The increasing incidence of community- and hospital-acquired infections with multidrug-resistant (MDR) bacteria poses a critical threat to public health and the healthcare system. Although β-lactam antibiotics are effective against most bacterial infections, some bacteria are resistant to β-lactam antibiotics by producing β-lactamases. Among β-lactamases, metallo-β-lactamases (MBLs) are especially worrisome as only a few inhibitors have been developed against them. In MBLs, the metal ions play an important role as they coordinate a catalytic water molecule that hydrolyzes β-lactam rings. We determined the crystal structures of different variants of PNGM-1, an ancient MBL with additional tRNase Z activity. The variants were generated by site-directed mutagenesis targeting metal-coordinating residues. In PNGM-1, both zinc ions are coordinated by six coordination partners in an octahedral geometry, and the zinc-centered octahedrons share a common face. Structures of the PNGM-1 variants confirm that the substitution of a metal-coordinating residue causes the loss of metal binding and β-lactamase activity. Compared with PNGM-1, subclass B3 MBLs lack one metal-coordinating residue, leading to a shift in the metal-coordination geometry from an octahedral to tetrahedral geometry. Our results imply that a subtle change in the metal-binding site of MBLs can markedly change their metal-coordination geometry and catalytic activity.

1. Introduction

The increasing prevalence of multidrug-resistant (MDR) bacteria in the community and hospitals represents a worldwide concern [1,2] and poses a serious threat to the healthcare system [3]. The World Health Organization (WHO) warns that unless the necessary actions are taken against MDR bacteria, the estimated death toll from MDR bacteria could exceed that from n class="Disease">cancer by 2050 [4].

β-lactam antibiotics, including penicillins, n class="Chemical">cephalosporins, carbapenems, and monobactams, have been the most effective antibiotics for almost 80 years since the first use of penicillin [5]. β-lactam antibiotics bind to bacterial penicillin-binding proteins (PBPs) via their four-membered β-lactam ring and interrupt the protein’s transpeptidase activity, which is required for synthesis of the peptidoglycan layer. This β-lactam-mediated inhibition of cell wall synthesis causes bacterial cell lysis. However, the emergence of β-lactamase genes and acquisition of these resistance genes by bacteria have led to a rapidly growing number of MDR bacteria that are resistant to β-lactam antibiotics [6]. β-lactamases can be divided into two groups according to their main catalytic motifs: serine β-lactamases and metallo-β-lactamases (MBLs). MBLs have a broad spectrum of hydrolytic activity and are resilient against most serine β-lactamase inhibitors; only monobactams exhibit some stability against MBLs [7].

Unlike serine β-lactamases in which a n class="Chemical">serine residue nucleophilically attacks the carbonyl carbon atom in the β-lactam ring, MBLs use a hydroxide ion to cleave the C–N bond, wherein the hydroxide ion is derived from a catalytic water molecule and coordinated by metal ion(s), usually Zn2+ [8,9] (Scheme 1). MBLs are classified into three subclasses: B1, B2, and B3. In general, B1 and B3 MBLs have two metal-binding sites (MBSs), whereas B2 MBLs have only one MBS [10,11]. As the substrate-attacking hydroxide ion is directly coordinated by metal ion(s) in the active site, the coordination geometry of metal binding plays an essential role in the structural and catalytic understanding of MBLs, which is critical for the development of a specific inhibitor against MBLs. The two MBSs in B1 and B3 MBLs, which are referred to as Zn1 and Zn2, have a different affinity for zinc ions [10]. B2 MBLs possess only Zn2.

Scheme 1

Mechanism of β-lactam hydrolysis catalyzed by metallo-β-lactamases (MBLs). (A) Hydrolysis mechanism catalyzed by MBLs containing two zinc ions. (B) Hydrolysis mechanism catalyzed by MBLs containing one zinc ion. The hydrolytic mechanism is divided into three main steps: First, MBL and β-lactam form an enzyme-substrate complex. Then, the hydroxide ion coordinated by zinc ion(s) attacks the β-lactam ring via a nucleophilic attack. Finally, a water molecule donates a proton to complete the reaction cycle.

PNGM-1 is derived from a metagenomic library from deep-sea sediments that predate the antibiotic era [12]. MBLs belong to the MBL superfamily, which also includes tRNase Zs [13]. PNGM-1 has both MBL and tRNase Z activities and is considered to be the evolutionary origin of B3 MBLs [14]. Like B1 and B3 MBLs, PNGM-1 has two MBSs.

MBLs hydrolyze almost all clinically used β-lactam antibiotics [15]; however, there is no effective inhibitor against MBLs yet. Compared to the catalytic serine residue in n class="Chemical">serine β-lactamases, the metal ions in the active site of MBLs do not allow the covalent binding of an irreversible or suicidal inhibitor. In addition, inhibitors with a high affinity for zinc ion(s) are toxic and not suitable for clinical applications due to putative nonspecific inhibition of ubiquitous zinc-binding proteins. We aimed to understand how metal-coordinating residues in the active site of MBLs affect the enzyme’s affinity for metal ions on a structural basis.

In this study, we determined the crystal structures of five Pn class="Chemical">NGM-1 variants and confirmed that each variant exhibits a different metal binding affinity. The present structural study improves the understanding of metal binding in MBLs.

2. Results

2.1. Two Metal-Binding Sites in PNGM-1

The wild-type of PNGM-1 has two MBSs in its substrate-binding pocket (Figure 1A), and each is occupied by a zinc ion (Figure S1). Each zinc ion is coordinated by four residues and two n class="Chemical">water molecules in an octahedral geometry, whereby Asp210 is involved in the binding of both zinc ions. The zinc ion in the first MBS (Zn1) is coordinated by the residues His91, His93, His188, and Asp210; the zinc ion in the second MBS (Zn2), by Asp95, His96, His279, and Asp210. The zinc ion at Zn1 is more tightly bound with shorter coordination distances than that at Zn2 (Figure S2 and Table S1). The coordination distances between the zinc ion at Zn1 and the coordinating residues range from 2.0 to 2.9 Å, and between that at Zn2 and the residues from 2.1 to 3.0 Å. The B-factor of the zinc ion at Zn1 is lower than that of Zn2 (Table S2).

Figure 1

The metal-binding site (MBS) of wild-type PNGM-1. (A) The MBSs of Zn1 (right) and Zn2 (left) are displayed. A zinc ion (grey sphere) is bound to each MBS and coordinates to six coordination partners (shown in stick mode) with coordination distances of 2.0–3.1 Å. Water molecules are displayed as red spheres. The catalytic and non-catalytic water molecules are labeled as Wcat and Wnoncat, respectively. Coordinations are indicated by dashed lines. (B) Scheme of the metal-binding geometry. The coordination partners and metal ion in each MBS exhibit an octahedron geometry. The binding site consists of two combined zinc-centered octahedrons.

2.2. Metal Coordination Geometry of the MBSs

In each MBS, a zinc ion is located at the center surrounded by six coordinating moieties, namely, four amino acid residues and two n class="Chemical">water molecules (or one hydroxide ion during the catalytic reaction) in an octahedral geometry (Figure 1B). The two octahedrons (Zn1 and Zn2) are attached to each other by sharing a plane formed by the two coordinating water molecules and Asp210. All zinc-coordinating atom distances are relatively consistent (2.0–2.3 Å) except for the long coordination distance (3.0 Å) between the zinc ions and the non-catalytic water molecule.

2.3. Metal Occupancy in Two MBSs

We determined the crystal structures of five PNGM-1 variants: H91A, n class="CellLine">H93A, H96A, H257A, and H279A (Figure 2 and Figure S3). His91 and His93 coordinate the zinc ion at Zn1, whereas His96 and His279 coordinate that at Zn2. His257 is not directly involved in metal coordination but plays a role in the tRNase Z activity of PNGM-1. In the structure of the H91A variant, the electron density for the zinc ion at Zn1 completely disappears compared to that in the wild-type structure, indicating the loss of the zinc ion in this variant. In the H93A variant, the electron density for the zinc ion at Zn1 is weak to allow only a partial occupancy of a zinc ion or a full occupancy for a water molecule. In the crystal structures of the H96A and H279A variants, the electron density of the zinc ion at Zn2 is markedly lower compared to that in the wild-type structure, indicating that either a water molecule with full occupancy or a zinc ion with half occupancy is present in this site. As His257 is distant from the MBSs and not involved in metal binding, the Zn1 and Zn2 in the H257A variant are occupied by metal ions.

Figure 2

The metal-binding sites of PNGM-1 wild-type and variants. The 2FoFc electron density maps contoured at 1.0 e/Å3 are shown (blue mesh) for (A) PNGM-1 wild-type, (B) H91A variant, (C) H93A variant, (D) H96A variant, (E) H257A variant, and (F) H279A variant.

Among all the variants, H91A and H96A exhibit the lowest metal occupancy. The n class="Disease">MBSs in PNGM-1 are solvent exposed (Figure S4). The residues His93, Aps95, His279, and Asp210 occupy equatorial positions along the enzyme’s surface. The residues His91 and His96 are buried in the substrate-binding pocket of PNGM-1.

2.4. Catalytic Water Molecule between the Two Zinc Ions

The catalytic mechanism of MBLs requires a substrate-attacking hydroxide ion bound by zinc ion(s) (Scheme 1). Tn class="Chemical">his hydroxide ion could be directly recruited from the bulk solvent or generated by the deprotonation of a water molecule located at the same position. In the PNGM-1 structures, two water molecules are involved in metal coordination, where one (catalytic water) attacks the C–N bond of the β-lactam ring (Figure 1). In the structures of all PNGM-1 variants, except in that of H96A, the electron density corresponding to the catalytic water remains despite the complete or partial loss of a zinc ion (Figure 2).

Both water molecules in the n class="Disease">MBSs are coordinated by four or five coordinating partners, including each other (Figure 1A). The catalytic water molecule lies on a plane formed by the four zinc-coordinating residues His93, His188, His279, and Asp95 and interacts with the non-catalytic water molecule in an apical position. The binding of the non-catalytic water with the coordination partner of the two zinc ions, Asp210, and the catalytic water is almost linear and unidirectional. The binding of the non-catalytic water molecule is weaker than that of the catalytic water. The B-factor of the catalytic water molecule is lower than that of the non-catalytic water molecule (Table S2).

2.5. Structural Comparisons between PNGM-1 and Other B3 MBLs

We structurally compared the MBSs of Pn class="Chemical">NGM-1 with those of tRNase Z from Bacillus subtilis (Bs-tRNase Z) [16] and B3 MBL from Elizabethkingia meningoseptica (GOB-18) [17]. In both Bs-tRNase Z and GOB-18, two metal ions are essential for enzyme activity. Although the positions of the two metal ions are well conserved among the compared enzymes, the metal coordination of Bs-tRNase Z and GOB-18 differs from that of PNGM-1. As shown for PNGM-1, Bs-tRNase Z contains two MBSs that exhibit octahedral geometries and share a common face. In contrast, GOB-18 lacks a metal-coordinating residue that corresponds to Asp210 in PNGM-1; therefore, its MBSs exhibit tetrahedral geometries share only one intersection point (Figure 3).

Figure 3

Comparison of the MBSs between PNGM-1 and Bs-tRNase Z and B3 MBL GOB-18. (A) The MBS of PNGM-1 (green) superimposed on that of Bs-tRNase Z (magenta). The metal ions (grey spheres) and water molecules (red spheres) belong to the structure of wild-type PNGM-1. (B) The MBS of PNGM-1 (green) superimposed on that of B3 MBL GOB-18 (cyan). The metal ions (grey spheres) and red water molecules belong to the structure of wild-type PNGM-1, whereas the cyan water molecule belongs to that of GOB-18. Metal coordination schemes are shown at the bottom.

According to structural superimposition, the zinc positions are almost identical in the structures of PNGM-1 and Bs-tRn class="Chemical">Nase Z; the distances between the zinc ions in the Zn1 and Zn2 are only 0.4 and 0.2 Å, respectively. The position of the phosphate group of the bound tRNA substrate in the Bs-tRNase Z structure superimposes well with that of the non-catalytic water molecule in the PNGM-1 structure (phosphate–water distance, 0.9 Å). This position is suitable for attack by the catalytic water molecule in PNGM-1. The catalytic water molecule is not present in the structure of the tRNA–Bs-tRNAse Z complex.

The differences in metal coordination between B3 MBL GOB-18 and Pn class="Chemical">NGM-1 cause a shift in the position of the catalytic water molecule within the active site. According to structural superimposition, the position of the catalytic water molecule in the GOB-18 structure is shifted by 2.4 Å toward residue Asp210 (only present in PNGM1) compared to that in the PNGM-1 structure. The GOB-18 structure does not contain a non-catalytic water molecule.

3. Discussion

MBLs hydrolyze a broad spectrum of β-lactam antibiotics and are promising targets for the development of antibacterial agents, especially those against MDR bacteria. The metal ions in MBLs are essential for catalytic activity in which the catalytic n class="Chemical">water molecule, which is coordinated to the metal ions, plays a key role by attacking the β-lactam ring. We scrutinized the metal binding and coordination in wild-type PNGM-1 and in different variants with substituted residues within the MBS. Each MBS of PNGM-1 is formed by four metal-coordinating residues, and the substitution of only one of these residues dramatically reduces the metal binding affinity of the enzyme (Figure S3 and Table S2). The PNGM-1 variants lose their catalytic activity for both RNAs and β-lactam rings [14].

The metal coordination geometries in the two n class="Disease">MBSs of PNGM-1 are almost identical. The metal coordination of Zn1 is slightly stronger than that of Zn2, as evidenced by shorter coordination distances. From a structural viewpoint, the substitution of His91 and His96 by alanine markedly weakens the metal binding affinity for Zn1 and Zn2, respectively, compared to that of the wild-type and the remaining variants. No electron density was observed for the zinc ion at Zn1 in the variant H91A. In H96A, even though electron density exists for the zinc ion at Zn2, the density is weak and noisy, implying weak binding of the metal ion. Both H91 and H96 are located deep in the catalytic pocket of PNGM-1; thus, the substitution of a bulky His residue to a sterically smaller Ala residue could lower the overall stability of the MBS, in addition to the loss of one metal coordination (Figure S2).

Geometrically, the MBS of Pn class="Chemical">NGM-1 consists of two zinc-centered octahedrons sharing a common face. The common face is formed by two water molecules and Asp210. In B3 MBL GOB-18, the residue corresponding to Asp210 in PNGM-1 is missing; therefore, the MBSs of GOB-18 exhibit tetrahedral metal coordination geometries, where the two metal-centered tetrahedrons share one common point, which is the catalytic water molecule (Figure S5). The change in the metal coordination geometry shifts the position of the catalytic water molecule, altering the direction of the nucleophilic attack by the catalytic hydroxide molecule in tRNase Z and MBL.

Structural inspection of PNGM-1 structures and comparison with those of Bs-tRn class="Chemical">Nase Z and B3 MBL GOB-18 showed that a variation in the metal coordination by substituting only one metal-coordinating residue could already affect the enzyme activity and substrate specificity. Currently, MBLs represent a serious threat to our community and healthcare system, and specific inhibitors against MBLs are being systematically pursued. For development of MBL inhibitors, interruption of the canonical metal coordination in MBLs could be an effective strategy. The present structural study of the metal binding and coordination in MBLs is useful for development of antibacterial agents against MDR bacteria containing MBLs.

4. Materials and Methods

4.1. Reagents

The expression vector, pET-28a(+), was purchased from n class="Chemical">Novagen (San Diego, CA, USA). The expression host, Escherichia coli BL21(DE3), and all restriction enzymes were purchased from New England Biolabs (Ipswich, MA, USA). Luria-Bertani (LB) medium was purchased from BD Biosciences (San Jose, CA, USA). Pre-stained protein markers for SDS-PAGE and a gel filtration calibration kit were purchased from MBI Fermentas (Hanover, MD, USA) and GE Healthcare (Piscataway, NJ, USA), respectively.

4.2. Gene Cloning and Site-Directed Mutagenesis

Gene cloning and site-directed mutagenesis were performed as described previously [14]. Briefly, the PblaPNGM-1 gene (GenBank ID: MF445022) was obtained from the functional metagenomic library from the deep-sea sediments of Edison Seamount in Papua n class="Chemical">New Guinea [18] and amplified by polymerase chain reaction (PCR) with the following primer pairs: forward primer, 5′-ATACCATGGGCCACCATCATCATCATCATGACGACGACGACAAGGCAGGTGGAAAAGTAACCTC-3′; backward primer, 5′-GAGAAGCTTTTAGCTTCCCACTCCCAAATC-3′ (restriction sites are written in bold, the hexahistidine (His6)-tag site is underlined, and the enterokinase (EK) recognition site is highlighted in italics). The amplified DNA and pET-28a (+) vector (Novagen, San Diego, CA, USA) were double-digested with NcoI and HindIII, and the digested DNA was then cloned into the digested pET-28a (+) vector.

Site-directed mutagenesis to generate the PNGM-1 variants H91A, n class="CellLine">H93A, H96A, H257A, and H279A was performed using a QuikChange II Site-Directed Mutagenesis Kit (Stratagene, Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s instructions. After verifying the DNA sequences, the resulting plasmids expressing the corresponding PNGM-1 variants were individually transformed into E. coli BL21(DE3) cells.

4.3. Expression and Purification of PNGM-1 Variants

The transformed E. coli BL21(DE3) cells containing the n class="Species">pET-28a (+) vectors with the PNGM-1 variant genes were grown in LB medium containing 50 mg/mL kanamycin at 303 K until an OD600nm of 0.6 was reached. Then, 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) was added to the culture to induce protein expression. After 16 h of cultivation at 289 K, the cells were harvested by centrifugation at 4700× g for 10 min at 277 K and resuspended in ice-cold 20 mM sodium phosphate pH 7.0. The cells were disrupted by sonication and centrifuged at 20,000× g for 60 min at 277 K. The clarified supernatant was loaded onto a His-Bind column (Novagen, San Diego, CA, USA) equilibrated with binding buffer (20 mM sodium phosphate pH 7.9, 10 mM imidazole, 500 mM NaCl). The His-tagged PNGM-1 variants were eluted with the same buffer containing 500 mM imidazole. For further purification, the His6-tag was removed from the PNGM-1 variants using enterokinase according to the manufacturer’s instructions (Novagen). The reaction mixture was desalted and concentrated using a Fast Desalting column (Amersham Biosciences, Buckinghamshire, UK) and then loaded onto a Mono S column (Amersham Biosciences, Buckinghamshire, UK) pre-equilibrated with 10 mM sodium phosphate pH 7.0. PNGM-1 variants were eluted with a linear gradient of NaCl (0–0.5 M) in phosphate buffer. The proteins were further purified by size-exclusion chromatography on a Superdex 200 (16/60) column (GE Healthcare) equilibrated with 10 mM MES, pH 6.8, at a flow rate of 1 mL/min.

4.4. Crystallization

The purified wild-type PNGM-1 and its variants H91A, H93A, H96A, n class="Mutation">H257A, and H279A were crystallized in 96-well Intelli-Plates (Art Robbins) and Hampton research buffer sets at 287 K using a Hydra II e-drop automated pipetting system (Matrix). Crystallization drops consisted of 0.5 µL protein solution and 0.5 µL reservoir solution, which were equilibrated against 50 μL reservoir solution at 287 K. After obtaining initial crystals, we optimized the crystallization conditions and refined it in terms of reliability. Each crystallization condition was as described in Table S3. For crystal optimization in 24 wells plates, hanging drops were prepared by mixing 1 μL protein solution with the same volume of mother liquor. Fully-grown single crystals were transferred to a mixture of mother liquor and cryoprotectant solution (20% glycerol, 20% ethylene glycol, or 15% MPD). Crystals were flash-cooled in liquid nitrogen at 100 K for data collection.

4.5. Data Collection, Structure Determination, and Refinement

X-ray data were collected at the beamlines BL-5C and BL-11C at the Pohang Light Source in South Korea. Data were integrated and scaled with the DENZO and SCALEPACK algorithms [19]. Initial phases of the structures were determined by molecular replacement (MR) with the phaser program of the CCP4 software package [20]. The structure of Pn class="Chemical">NGM-1 (PDB ID: 6j4n) was used as the search model for MR. Manual model building into the electron density map was performed using the COOT program [21]. Multiple cycles of structure refinement were performed with Refmac5 of the CCP4 software package [22]. The coordinates and geometries of all structures were validated by WHATIF [23] and SFCheck [24]. Graphic presentations were created with PyMOL [25]. The diffraction data and structure refinement statistics are provided in Table 1.

Table 1

Data collection and refinement statistics.

	H91A	H93A	H96A	H257A	H279A
Data collection
Wavelength (Å)	0.9796	0.9796	0.9796	0.9796	0.9796
Resolution range	48.0–1.9(2.01–1.94)	40.3–1.6(1.67–1.61)	45.8–2.5(2.54–2.45)	48.6–2.0(2.07–2.0)	48.0–2.2(2.23–2.16)
Space group	P 1 21 1	P 1 21 1	C 1 2 1	P 1 21 1	C 1 2 1
Unit cell a, b, c (Å)α, β, γ (°)	121.4, 83.1, 163.690, 110.7, 90	121.8, 82.8, 164.090, 111.1, 90	121.6, 82.5, 163.790, 110.8, 90	79.8, 143.7, 79.890, 111.8, 90	121.9, 83.0, 164.090, 110.9, 90
Total reflections	1082322	2483252	255443	524840	397308
Unique reflections	214810 (16994)	388130 (37693)	50474 (4821)	109298 (10188)	79596 (7627)
Multiplicity	5.1 (4.0)	6.4 (5.5)	5.1 (4.7)	4.8 (3.5)	5.0 (3.7)
Completeness (%)	95.0 (75.9)	98.9 (96.5)	90.2 (87.0)	97.3 (91.5)	96.5 (92.6)
Mean I/sigma (I)	14.5 (2.3)	28.0 (4.4)	16.9 (3.8)	11.3 (2.3)	9.6(2.2)
Rmerge (%)	12.2 (39.7)	8.6 (31.0)	12.4 (41.4)	13.9 (38.2)	16.5 (43.7)
Refinement
Reflections used in refinement	214406 (16993)	387739 (37693)	50455 (4819)	109291 (10188)	79574 (7627)
Reflections used for Rfree	10759 (868)	19516 (1870)	2481 (270)	5527 (483)	3906 (344)
Rwork	0.22 (0.35)	0.20 (0.21)	0.19 (0.25)	0.24 (0.28)	0.21 (0.27)
Rfree	0.26 (0.36)	0.23 (0.24)	0.26 (0.35)	0.29 (0.33)	0.26 (0.32)
Number of non-hydrogen atoms	23700	24548	11853	11315	11670
macromolecules	22738	23404	11540	10640	11349c
ligands	8	8	4	8	4
solvent	954	1136	309	667	317
Protein residues	2887	2976	1460	1347	1444
RMS (bonds)	0.014	0.015	0.016	0.022	0.015
RMS (angles)	1.84	1.92	1.97	2.09	2.00
Ramachandran favored (%)	95.5	96.7	93.1	92.7	93.8
Ramachandran allowed (%)	4.1	2.7	6.0	5.6	5.2
Ramachandran outliers (%)	0.4	0.6	0.9	1.7	1.0
Rotamer outliers (%)	4.2	3.0	5.3	5.0	4.0
Clash score	7.3	4.3	6.7	9.0	5.9
Average B-factor	28.1	14.3	31.3	26.9	20.4
macromolecules	28.1	14.2	31.48	26.7	20.5
ligands	21.4	9.4	26.21	19.8	18.0
solvent	27.9	17.6	24.74	30.6	16.7

Values in parentheses are for the shell with the highest resolution. where is the mean intensity of the ith observation of symmetry-related reflections hkl. where is the calculated protein structure factor from the atomic model (Rfree was calculated as Rwork with a randomly selected 5% of the reflections).