Quantum mechanical and molecular dynamics simulations of ureases and Zn beta-lactamases

Molecular Engineering of Bacillus paralicheniformis Acid Urease To Degrade Urea and Ethyl Carbamate in Model Chinese Rice Wine

Bacillus paralicheniformis urease (BpUrease) has been shown to be a promising biocatalyst for degrading the carcinogenic chemical ethyl carbamate (EC or urethane) in rice wine. However, low EC affinity and catalytic efficiency limit the practical application of BpUrease. In this study, we improved the EC degradation capability of BpUrease by site-saturation mutagenesis (SSM). The best variant L253P/L287N showed a 49% increase in EC affinity, 1027% increase in catalytic efficiency ( k_cat/ K_m), and 583% increase in half-life ( t_1/2) at 70 °C. Homology modeling analysis suggest that mutation of Leu253 to Pro increased the BpUrease EC specificity by affecting the interaction between Arg339 with the catalytic residue His323, while Leu287Asn mutation benefits EC specificity and affinity by changing the interaction networks among the residues in the catalytic pocket. Our results show that the L253P/L287N variant efficiently degraded urea and EC in a model rice wine, making it a good candidate for practical application in the food industry.

Zn Coordination Chemistry: Development of Benchmark Suites for Geometries, Dipole Moments, and Bond Dissociation Energies and Their Use To Test and Validate Density Functionals and Molecular Orbital Theory

We present nonrelativistic and relativistic benchmark databases (obtained by coupled cluster calculations) of 10 Zn-ligand bond distances, 8 dipole moments, and 12 bond dissociation energies in Zn coordination compounds with O, S, NH3, H2O, OH, SCH3, and H ligands. These are used to test the predictions of 39 density functionals, Hartree-Fock theory, and seven more approximate molecular orbital theories. In the nonrelativisitic case, the M05-2X, B97-2, and mPW1PW functionals emerge as the most accurate ones for this test data, with unitless balanced mean unsigned errors (BMUEs) of 0.33, 0.38, and 0.43, respectively. The best local functionals (i.e., functionals with no Hartree-Fock exchange) are M06-L and τ-HCTH with BMUEs of 0.54 and 0.60, respectively. The popular B3LYP functional has a BMUE of 0.51, only slightly better than the value of 0.54 for the best local functional, which is less expensive. Hartree-Fock theory itself has a BMUE of 1.22. The M05-2X functional has a mean unsigned error of 0.008 Å for bond lengths, 0.19 D for dipole moments, and 4.30 kcal/mol for bond energies. The X3LYP functional has a smaller mean unsigned error (0.007 Å) for bond lengths but has mean unsigned errors of 0.43 D for dipole moments and 5.6 kcal/mol for bond energies. The M06-2X functional has a smaller mean unsigned error (3.3 kcal/mol) for bond energies but has mean unsigned errors of 0.017 Å for bond lengths and 0.37 D for dipole moments. The best of the semiempirical molecular orbital theories are PM3 and PM6, with BMUEs of 1.96 and 2.02, respectively. The ten most accurate functionals from the nonrelativistic benchmark analysis are then tested in relativistic calculations against new benchmarks obtained with coupled-cluster calculations and a relativistic effective core potential, resulting in M05-2X (BMUE = 0.895), PW6B95 (BMUE = 0.90), and B97-2 (BMUE = 0.93) as the top three functionals. We find significant relativistic effects (∼0.01 Å in bond lengths, ∼0.2 D in dipole moments, and ∼4 kcal/mol in Zn-ligand bond energies) that cannot be neglected for accurate modeling, but the same density functionals that do well in all-electron nonrelativistic calculations do well with relativistic effective core potentials. Although most tests are carried out with augmented polarized triple-ζ basis sets, we also carried out some tests with an augmented polarized double-ζ basis set, and we found, on average, that with the smaller basis set DFT has no loss in accuracy for dipole moments and only ∼10% less accurate bond lengths.

Monoligand Zn(II) Complexes: Ab Initio Benchmark Calculations and Comparison with Density Functional Theory Methodologies

A systematic theoretical study on several models of Zn(II) complexes has been carried out employing both ab initio correlated wave function and density functional methods. The performance of five different functionals namely PW91, PBE, B3LYP, MPWLYP1M, and TPSS in the prediction of metal-ligand bond distances, binding energies, and proton affinities has been assessed comparing the results to those obtained with the MP2 and CCSD(T) wave function methodologies. Several basis sets ranging from double-ζ up to quintuple-ζ quality have been used, including the recently developed all-electron correlation consistent basis sets for zinc. It is shown that all the tested functionals overestimate both the metal-ligand bond distances and the binding energies, being that the B3LYP and TPSS functionals are the ones that perform the best. An analysis of the metal-ligand interaction energy shows that induction and charge-transfer effects play a prominent role in the bonding of these systems, even for those complexes with the less polarizable ligands. This finding highlights the importance of a correct description of the polarization of the monomers' charge densities by any theoretical method which aims to be applied to the study of Zn(II) complexes.

Overcoming differences: The catalytic mechanism of metallo-β-lactamases

Metallo-β-lactamases are the latest resistance mechanism of pathogenic and opportunistic bacteria against carbapenems, considered as last resort drugs. The worldwide spread of genes coding for these enzymes, together with the lack of a clinically useful inhibitor, have raised a sign of alarm. Inhibitor design has been mostly impeded by the structural diversity of these enzymes. Here we provide a critical review of mechanistic studies of the three known subclasses of metallo-β-lactamases, analyzed at the light of structural and mutagenesis investigations. We propose that these enzymes present a modular structure in their active sites that can be dissected into two halves: one providing the attacking nucleophile, and the second one stabilizing a negatively charged reaction intermediate. These are common mechanistic elements in all metallo-β-lactamases. Nucleophile activation does not necessarily requires a Zn(II) ion, but a Zn(II) center is essential for stabilization of the anionic intermediate. Design of a common inhibitor could be therefore approached based in these convergent mechanistic features despite the structural differences.

Catalytic mechanisms of metallohydrolases containing two metal ions

At least one-third of enzymes contain metal ions as cofactors necessary for a diverse range of catalytic activities. In the case of polymetallic enzymes (i.e., two or more metal ions involved in catalysis), the presence of two (or more) closely spaced metal ions gives an additional advantage in terms of (i) charge delocalisation, (ii) smaller activation barriers, (iii) the ability to bind larger substrates, (iv) enhanced electrostatic activation of substrates, and (v) decreased transition-state energies. Among this group of proteins, enzymes that catalyze the hydrolysis of ester and amide bonds form a very prominent family, the metallohydrolases. These enzymes are involved in a multitude of biological functions, and an increasing number of them gain attention for translational research in medicine and biotechnology. Their functional versatility and catalytic proficiency are largely due to the presence of metal ions in their active sites. In this chapter, we thus discuss and compare the reaction mechanisms of several closely related enzymes with a view to highlighting the functional diversity bestowed upon them by their metal ion cofactors.

Mysteries of metals in metalloenzymes

Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn(2+) in the active site but becomes much more active with Fe(2+). For β-lactamases, the replacement of the native Zn(2+) with Ni(2+) was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe(2+)- and Mn(2+)-dependent ribonucleotide reductases and W(4+)- and Mo(4+)-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe(2+)) by creating protein structures that utilize another metal (e.g., Mn(2+)) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.

Selective recognition of parallel and anti-parallel thrombin-binding aptamer G-quadruplexes by different fluorescent dyes

G-quadruplexes (G4) have been found increasing potential in applications, such as molecular therapeutics, diagnostics and sensing. Both Thioflavin T (ThT) and N-Methyl mesoporphyrin IX (NMM) become fluorescent in the presence of most G4, but thrombin-binding aptamer (TBA) has been reported as the only exception of the known G4-forming oligonucleotides when ThT is used as a high-throughput assay to identify G4 formation. Here, we investigate the interactions between ThT/NMM and TBA through fluorescence spectroscopy, circular dichroism and molecular docking simulation experiments in the absence or presence of cations. The results display that a large ThT fluorescence enhancement can be observed only when ThT bind to the parallel TBA quadruplex, which is induced to form by ThT in the absence of cations. On the other hand, great promotion in NMM fluorescence can be obtained only in the presence of anti-parallel TBA quadruplex, which is induced to fold by K+ or thrombin. The highly selective recognition of TBA quadruplex with different topologies by the two probes may be useful to investigate the interactions between conformation-specific G4 and the associated proteins, and could also be applied in label-free fluorescent sensing of other biomolecules.

Molecular Dynamics Study of Helicobacter pylori Urease

Helicobacter pylori have been implicated in an array of gastrointestinal disorders including, but not limited to, gastric and duodenal ulcers and adenocarcinoma. This bacterium utilizes an enzyme, urease, to produce copious amounts of ammonia through urea hydrolysis in order to survive the harsh acidic conditions of the stomach. Molecular dynamics (MD) studies on the H. pylori urease enzyme have been employed in order to study structural features of this enzyme that may shed light on the hydrolysis mechanism. A total of 400 ns of MD simulation time were collected and analyzed in this study. A wide-open flap state previously observed in MD simulations on Klebsiella aerogenes [Roberts et al. J. Am. Chem. Soc.2012, 134, 9934] urease has been identified in the H. pylori enzyme that has yet to be experimentally observed. Critical distances between residues on the flap, contact points in the closed state, and the separation between the active site Ni²⁺ ions and the critical histidine α322 residue were used to characterize flap motion. An additional flap in the active site was elaborated upon that we postulate may serve as an exit conduit for hydrolysis products. Finally we discuss the internal hollow cavity and present analysis of the distribution of sodium ions over the course of the simulation.

Biapenem inactivation by B2 metallo β-lactamases: energy landscape of the hydrolysis reaction

Background

A general mechanism has been proposed for metallo β-lactamases (MβLs), in which deprotonation of a water molecule near the Zn ion(s) results in the formation of a hydroxide ion that attacks the carbonyl oxygen of the β-lactam ring. However, because of the absence of X-ray structures that show the exact position of the antibiotic in the reactant state (RS) it has been difficult to obtain a definitive validation of this mechanism.

Methodology/Principal Findings

We have employed a strategy to identify the RS, which does not rely on substrate docking and/or molecular dynamics. Starting from the X-ray structure of the enzyme:product complex (the product state, PS), a QM/MM scan was used to drive the reaction uphill from product back to reactant. Since in this process also the enzyme changes from PS to RS, we actually generate the enzyme:substrate complex from product and avoid the uncertainties associated with models of the reactant state. We used this strategy to study the reaction of biapenem hydrolysis by B2 MβL CphA. QM/MM simulations were carried out under 14 different ionization states of the active site, in order to generate potential energy surfaces (PESs) corresponding to a variety of possible reaction paths.

Conclusions/Significance

The calculations support a model for biapenem hydrolysis by CphA, in which the nucleophile that attacks the β-lactam ring is not the water molecule located in proximity of the active site Zn, but a second water molecule, hydrogen bonded to the first one, which is used up in the reaction, and thus is not visible in the X-ray structure of the enzyme:product complex.

Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF)

Currently the Protein Data Bank (PDB) contains over 18,000 structures that contain a metal ion including Na, Mg, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Pd, Ag, Cd, Ir, Pt, Au, and Hg. In general, carrying out classical molecular dynamics (MD) simulations of metalloproteins is a convoluted and time consuming process. Herein, we describe MCPB (Metal Center Parameter Builder), which allows one, to conveniently and rapidly incorporate metal ions using the bonded plus electrostatics model (Hoops et al., J. Am. Chem. Soc. 1991, 113, 8262-8270) into the AMBER Force Field (FF). MCPB was used to develop a Zinc FF, ZAFF, which is compatible with the existing AMBER FFs. The PDB was mined for all Zn containing structures with most being tetrahedrally bound. The most abundant primary shell ligand combinations were extracted and FFs were created. These include Zn bound to CCCC, CCCH, CCHH, CHHH, HHHH, HHHO, HHOO, HOOO, HHHD, and HHDD (O = water and the remaining are 1 letter amino acid codes). Bond and angle force constants and RESP charges were obtained from B3LYP/6-31G* calculations of model structures from the various primary shell combinations. MCPB and ZAFF can be used to create FFs for MD simulations of metalloproteins to study enzyme catalysis, drug design and metalloprotein crystal refinement.

Identification of novel urease inhibitors by high-throughput virtual and in vitro screening

Ureases are important in both agriculture and human health. Bacterial ureases are directly involved in many farm-field problems and pathological conditions. Here, we report a structure-based virtual screening of an in-house compound bank of about 6000 molecular entities by computational docking and binding free energy calculations followed by in vitro screening. Applied protocol leads to the identification of novel urease inhibitors, which can serve as starting points for structural optimization.

Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field

The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.

Energies, Geometries, and Charge Distributions of Zn Molecules, Clusters, and Biocenters from Coupled Cluster, Density Functional, and Neglect of Diatomic Differential Overlap Models

We present benchmark databases of Zn-ligand bond distances, bond angles, dipole moments, and bond dissociation energies for Zn-containing small molecules and Zn coordination compounds with H, CH3, C2H5, NH3, O, OH, H2O, F, Cl, S, and SCH3 ligands. The test set also includes clusters with Zn-Zn bonds. In addition, we calculated dipole moments and binding energies for Zn centers in coordination environments taken from zinc metalloenzyme X-ray structures, representing both structural and catalytic zinc centers. The benchmark values are based on relativistic-core coupled cluster calculations. These benchmark calculations are used to test the predictions of four density functionals, namely B3LYP and the more recently developed M05-2X, M06, and M06-2X levels of theory, and six semiempirical methods, including neglect of diatomic differential overlap (NDDO) calculations incorporating the new PM3 parameter set for Zn called ZnB, developed by Brothers and co-workers, and the recent PM6 parametrization of Stewart. We found that the best DFT method to reproduce dipole moments and dissociation energies of our Zn compound database is M05-2X, which is consistent with a previous study employing a much smaller and less diverse database and a much larger set of density functionals. Here we show that M05-2X geometries and single-point coupled cluster calculations with M05-2X geometries can also be used as benchmarks for larger compounds, where coupled cluster optimization is impractical, and in particular we use this strategy to extend the geometry, binding energy, and dipole moment databases to additional molecules, and we extend the tests involving crystal-site coordination compounds to two additional proteins. We find that the most predictive NDDO methods for our training set are PM3 and MNDO/d. Notably, we also find large errors in B3LYP for the coordination compounds based on experimental X-ray geometries.

A comparative study of semiempirical, ab initio, and DFT methods in evaluating metal-ligand bond strength, proton affinity, and interactions between first and second shell ligands in Zn-biomimetic complexes

Although theoretical methods are now available which give very accurate results, often comparable to the experimental ones, modeling chemical or biological interesting systems often requires less demanding and less accurate theoretical methods, mainly due to computer limitations. Therefore, it is crucial to know the precision of such less reliable methods for relevant models and data. This has been done in this work for small zinc-active site models including O- (H(2)O and OH(-)) and N-donor (NH(3) and imidazole) ligands. Calculations using a number of quantum mechanical methods were carried out to determine their precision for geometries, coordination number relative stability, metal-ligand bond strengths, proton affinities, and interaction energies between first and second shell ligands. We have found that obtaining chemical accuracy can be as straightforward as HF geometry optimization with a double-zeta plus polarization basis followed by a B3LYP energy calculation with a triple-zeta quality basis set including diffuse and polarization functions. The use of levels as low as PM3 geometry optimization followed by a B3LYP single-point energy calculation with a double-zeta quality basis including polarization functions already yields useful trends in bond length, proton affinities or bond dissociation energies, provided that appropriate caution is taken with the optimized structures. The reliability of these levels of calculation has been successfully demonstrated for real biomimetic cases.

Insights into the role and structure of plant ureases

The broad distribution of ureases in leguminous seeds, as well as the accumulation pattern of the protein during seed maturation, are suggestive of an important physiological role for this enzyme. Since the isolation and characterization of jack bean urease by Sumner in 1926, many investigations have been dedicated to the structural and biological features of this enzyme; nevertheless, many questions still remain. It has been reported that ureases from plants (jack bean and soybean seeds) display biological properties unrelated to their ureolytic activity, notably a high insecticidal activity against Coleoptera (beetles) and Hemiptera (bugs), suggesting that ureases might be involved in plant defense. Besides the insecticidal activity, canatoxin, a jack bean urease isoform, causes convulsions and death in mice and rats, induces indirect hemagglutination (hemilectin activity) and promotes exocytosis in several cell types. Not only plant ureases but also some microbial ureases (found in Bacillus pasteurii and Helicobacter pylori) are able to induce activation of platelets in a process mediated by lipoxygenase-derived metabolites. This review summarizes the biological and structural properties of plant ureases, compares them with those displayed by bacterial ureases, and discusses the significance of these findings.

Protonation state and substrate binding to B2 metallo-beta-lactamase CphA from Aeromonas hydrofila

The zinc enzymes metallo beta-lactamases counteract the beneficial action of beta-lactam antibiotics against bacterial infections, by hydrolyzing their beta-lactam rings. To understand structure/function relationships on a representative member of this class, the B2 M beta L CphA, we have investigated the H-bond pattern at the Zn enzymatic active site and substrate binding mode by molecular simulation methods. Extensive QM calculations at the DFT-BLYP level on eleven models of the protein active site, along with MD simulations of the protein in aqueous solution, allow us to propose two plausible protonation states for the unbound enzyme, which are probably in equilibrium. Docking procedures along with MD simulations and QM calculations suggest that in the complex between the enzyme and its substrate (biapenem), the latter is stable in only one of the two protonation states, in addition it exhibits two different binding modes, of which only one agrees with previous proposals. In both cases, the substrate is polarized as in aqueous solution. We conclude that addressing mechanistic issues on this class of enzymes requires a careful procedure to assign protonation states and substrate docking modes.

Antibiotic Deactivation by a Dizinc β-Lactamase: Mechanistic Insights from QM/MM and DFT Studies

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods and density functional theory (DFT) were used to investigate the initial ring-opening step in the hydrolysis of moxalactam catalyzed by the dizinc L1 beta-lactamase from Stenotrophomonas maltophilia. Anchored at the enzyme active site via direct metal binding as suggested by a recent X-ray structure of an enzyme-product complex (Spencer, J.; et al. J. Am. Chem. Soc. 2005, 127, 14439), the substrate is well aligned with the nucleophilic hydroxide that bridges the two zinc ions. Both QM/MM and DFT results indicate that the addition of the hydroxide nucleophile to the carbonyl carbon in the substrate lactam ring leads to a metastable intermediate via a dominant nucleophilic addition barrier. The potential of mean force obtained by SCC-DFTB/MM simulations and corrected by DFT/MM calculations yields a reaction free energy barrier of 23.5 kcal/mol, in reasonable agreement with the experimental value of 18.5 kcal/mol derived from kcat of 0.15 s(-1). It is further shown that zinc-bound Asp120 plays an important role in aligning the nucleophile, but accepts the hydroxide proton only after the nucleophilic addition. The two zinc ions are found to participate intimately in the catalysis, consistent with the proposed mechanism. In particular, the Zn(1) ion is likely to serve as an "oxyanion hole" in stabilizing the carbonyl oxygen, while the Zn(2) ion acts as an electrophilic catalyst to stabilize the anionic nitrogen leaving group.

Antibiotic Binding to Dizinc β-Lactamase L1 from Stenotrophomonas maltophilia: SCC-DFTB/CHARMM and DFT Studies

A dizinc beta-lactamase (L1 from Stenotrophomonas maltophilia) complexed with an antibiotic compound (moxalactam) has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach. The QM region is described by the self-consistent charge-density functional tight binding (SCC-DFTB) model while the MM by CHARMM. The Michaelis complex, which is constructed from a recent X-ray structure of the L1 enzyme with the hydrolyzed moxalactam, is simulated by molecular dynamics. The simulation yields valuable insights into substrate-enzyme interaction, whose implications in the enzyme catalysis are discussed. Finally, the QM/MM results are compared with a high-level density functional theory study of a truncated active-site model and the agreement provides strong support for the SCC-DFTB treatment of the QM region.

Role of zinc content on the catalytic efficiency of B1 metallo beta-lactamases

Metallo beta-lactamases (MbetaL) are enzymes naturally evolved by bacterial strains under the evolutionary pressure of beta-lactam antibiotic clinical use. They have a broad substrate spectrum and are resistant to all the clinically useful inhibitors, representing a potential risk of infection if massively disseminated. The MbetaL scaffold is designed to accommodate one or two zinc ions able to activate a nucleophilic hydroxide for the hydrolysis of the beta-lactam ring. The role of zinc content on the binding and reactive mechanism of action has been the subject of debate and still remains an open issue despite the large amount of data acquired. We report herein a study of the reaction pathway for binuclear CcrA from Bacteroides fragilis using density functional theory based quantum mechanics-molecular mechanics dynamical modeling. CcrA is the prototypical binuclear enzyme belonging to the B1 MbetaL family, which includes several harmful chromosomally encoded and transferable enzymes. The involvement of a second zinc ion in the catalytic mechanism lowers the energetic barrier for beta-lactam hydrolysis, preserving the essential binding features found in mononuclear B1 enzymes (BcII from Bacillus cereus) while providing a more efficient single-step mechanism. Overall, this study suggests that uptake of a second equivalent zinc ion is evolutionary favored.