Molecular dynamics simulations of the mononuclear zinc-beta-lactamase from Bacillus cereus complexed with benzylpenicillin and a quantum chemical study of the reaction mechanism

Metallo-β-lactamases in the Age of Multidrug Resistance: From Structure and Mechanism to Evolution, Dissemination, and Inhibitor Design

Antimicrobial resistance is one of the major problems in current practical medicine. The spread of genes coding for resistance determinants among bacteria challenges the use of approved antibiotics, narrowing the options for treatment. Resistance to carbapenems, last resort antibiotics, is a major concern. Metallo-β-lactamases (MBLs) hydrolyze carbapenems, penicillins, and cephalosporins, becoming central to this problem. These enzymes diverge with respect to serine-β-lactamases by exhibiting a different fold, active site, and catalytic features. Elucidating their catalytic mechanism has been a big challenge in the field that has limited the development of useful inhibitors. This review covers exhaustively the details of the active-site chemistries, the diversity of MBL alleles, the catalytic mechanism against different substrates, and how this information has helped developing inhibitors. We also discuss here different aspects critical to understand the success of MBLs in conferring resistance: the molecular determinants of their dissemination, their cell physiology, from the biogenesis to the processing involved in the transit to the periplasm, and the uptake of the Zn(II) ions upon metal starvation conditions, such as those encountered during an infection. In this regard, the chemical, biochemical and microbiological aspects provide an integrative view of the current knowledge of MBLs.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Peptide Hydrolysis in Thermolysin: Ab Initio QM/MM Investigation of the Glu143-Assisted Water Addition Mechanism

Thermolysin (TLN) is one of the best-studied zinc metalloproteases. Yet the mechanism of action is still under debate. In order to investigate the energetic feasibility of the currently most favored mechanism, we have docked a tripeptide to the active site of TLN and computed the free energy profile at the quantum mechanics/molecular mechanics level of theory. The mechanism consists of three distinct steps:  (i) a Zn-bound water molecule is deprotonated by Glu143 and attacks the carbonyl bond of the substrate; (ii) Glu143 transfers the proton to the amide nitrogen atom; (iii) the nitrogen atom is protonated and the peptide bond is irreversibly broken. The free energy barriers for steps i and iii have almost equal heights, 14.8 and 14.7 kcal/mol, respectively, and are in good agreement with the effective experimental activation barrier obtained for similar substrates, 12.1-13.6 kcal/mol. Transition state stabilization for nucleophilic attack is achieved by formation of a weak coordination bond between the substrate carbonyl oxygen atom and the Zn ion and of three strong hydrogen bonds between the substrate and protonated His231 and two solvent molecules. The transition state for the nucleophilic attack (step i) is more tightly bonded than the enzyme-substrate complex, implying that TLN complies with Pauling's hypothesis regarding transition-state stabilization. Glu143, at first unfavorably oriented for protonation of the amide nitrogen atom, displayed large structural fluctuations that facilitated reorganization of the local hydrogen-bond network and transport of the proton to the leaving group on the nanosecond time scale. The present simulations give further evidence that Glu143 is a highly effective proton shuttle which should be assigned a key role in any reaction mechanism proposed for TLN.

QM/MM molecular dynamics studies of metal binding proteins

Mixed quantum-classical (quantum mechanical/molecular mechanical (QM/MM)) simulations have strongly contributed to providing insights into the understanding of several structural and mechanistic aspects of biological molecules. They played a particularly important role in metal binding proteins, where the electronic effects of transition metals have to be explicitly taken into account for the correct representation of the underlying biochemical process. In this review, after a brief description of the basic concepts of the QM/MM method, we provide an overview of its capabilities using selected examples taken from our work. Specifically, we will focus on heme peroxidases, metallo-β-lactamases, α-synuclein and ligase ribozymes to show how this approach is capable of describing the catalytic and/or structural role played by transition (Fe, Zn or Cu) and main group (Mg) metals. Applications will reveal how metal ions influence the formation and reduction of high redox intermediates in catalytic cycles and enhance drug metabolism, amyloidogenic aggregate formation and nucleic acid synthesis. In turn, it will become manifest that the protein frame directs and modulates the properties and reactivity of the metal ions.

Asp120Asn mutation impairs the catalytic activity of NDM-1 metallo-β-lactamase: experimental and computational study

New Delhi metallo-β-lactamase-1 (NDM-1) has attracted extensive attention in recent years for its high activity for hydrolyzing almost all β-lactam antibiotics. Like other metallo-β-lactamases (MβLs), NDM-1 features an invariant Asp120 that ligates the zinc ion (ZN2) in the active site. Previous studies showed that substitutions of Asp120 with residues such as Ala, Ser, Asn and Glu dramatically impaired the MβL (BcII, IMP-1, L1) activity, but no consensus about the exact role of Asp120 has reached. Here we constructed D120N mutant of NDM-1 by site-directed mutagenesis. The replacement of Asp120 with Asn, which has much weaker metal ligating capabilities than Asp, severely impaired the lactamase activity without abolishing the ZN2 site. Molecular dynamics simulations suggested that the ZN1-ZN2 distance increased because of mutation, leading to a rearrangement of the active site, including the bridging OH(-). Thereby, the Mulliken charges of ZN1 and ZN2 redistributed, especially for ZN2, which might be the major cause of the impaired activity. Reducing the point charges of Asp120 carboxyl oxygens weakened the ionic interactions between Asp120 and ZN2, and the positions of the zinc ions were also changed as a result. It is proposed that Asp120 acts as a strong ZN2 ligand, positioning ZN2 for catalytically important interactions with the substrate, stabilizing the negatively charged amide nitrogen of the hydrolyzed intermediate, and more importantly, orienting the ZN-bound OH(-) for nucleophilic attacks and protonation. These functions are of general importance for catalyzing β-lactam antibiotics by NDM-1 as well as other MβLs.

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.

On the active site of mononuclear B1 metallo β-lactamases: a computational study

Metallo-β-lactamases (MβLs) are Zn(II)-based bacterial enzymes that hydrolyze β-lactam antibiotics, hampering their beneficial effects. In the most relevant subclass (B1), X-ray crystallography studies on the enzyme from Bacillus Cereus point to either two zinc ions in two metal sites (the so-called '3H' and 'DCH' sites) or a single Zn(II) ion in the 3H site, where the ion is coordinated by Asp120, Cys221 and His263 residues. However, spectroscopic studies on the B1 enzyme from B. Cereus in the mono-zinc form suggested the presence of the Zn(II) ion also in the DCH site, where it is bound to an aspartate, a cysteine, a histidine and a water molecule. A structural model of this enzyme in its DCH mononuclear form, so far lacking, is therefore required for inhibitor design and mechanistic studies. By using force field based and mixed quantum-classical (QM/MM) molecular dynamics (MD) simulations of the protein in aqueous solution we constructed such structural model. The geometry and the H-bond network at the catalytic site of this model, in the free form and in complex with two common β-lactam drugs, is compared with experimental and theoretical findings of CphA and the recently solved crystal structure of new B2 MβL from Serratia fonticola (Sfh-I). These are MβLs from the B2 subclass, which features an experimentally well established mono-zinc form, in which the Zn(II) is located in the DCH site. From our simulations the εεδ and δεδ protomers emerge as possible DCH mono-zinc reactive species, giving a novel contribution to the discussion on the MβL reactivity and to the drug design process.

DIFFERENT PROTONATION STATES OF THE BACILLUS CEREUS BINUCLEAR ZINC METALLO-β-LACTAMASE ACTIVE SITE STUDIED BY COMBINED QUANTUM MECHANICAL AND MOLECULAR MECHANICAL SIMULATIONS

Three different protonation states of the active site of the Bacillus cereus zinc-β-lactamase in its binuclear form are studied using combined quantum mechanics/molecular mechanics molecular dynamics simulations. The reliability of the quantum mechanical model, the self-consistent-charge density-functional-based tight binding method, in describing the zinc centers are tested through comparisons with ab initio quantum mechanical results. We found that this model gave relatively accurate results for structures and performed much better than the MNDO type semi-empirical method for the particular systems. The enzyme simulations suggested that when the overall charge of the active site is +1, i.e., both Asp90 and Wat1 (a water molecule coordinated with the first zinc ion) deprotonated, the second zinc ion is coordinated with Asp90 and Wat1, and a second water molecule cannot coordinate with the second zinc ion. When the overall charge is +2, i.e., either Asp90 or Wat1 protonated, Asp90 and Wat1 form a stable hydrogen bond. Depending on the proton being on Asp90 or on Wat1, the active site structure produced by the simulations is either similar to molecule A or to molecule B, both contained in the same crystal structure that has two enzyme molecules in a single asymmetric unit. The simulations of the +2 charge states also reproduced the experimentally observed "loose" coordination around the second zinc for the Bacillus Cereus enzyme. Based on the simulations and a gas phase potential energy surface scanning using ab initio model, we argue that the penta-coordination around the second zinc ion is not a stable arrangement. Mechanistic implications of these results are discussed.

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.

Catalytic role of the metal ion in the metallo-beta-lactamase GOB

Metallo-beta-lactamases (MbetaLs) stand as one of the main mechanisms of bacterial resistance toward carbapenems. The rational design of an inhibitor for MbetaLs has been limited by an incomplete knowledge of their catalytic mechanism and by the structural diversity of their active sites. Here we show that the MbetaL GOB from Elizabethkingia meningoseptica is active as a monometallic enzyme by using different divalent transition metal ions as surrogates of the native Zn(II) ion. Of the metal derivatives in which Zn(II) is replaced, Co(II) and Cd(II) give rise to the most active enzymes and are shown to occupy the same binding site as the native ion. However, Zn(II) is the only metal ion capable of stabilizing an anionic intermediate that accumulates during nitrocefin hydrolysis, in which the C-N bond has already been cleaved. This finding demonstrates that the catalytic role of the metal ion in GOB is to stabilize the formation of this intermediate prior to nitrogen protonation. This role may be general to all MbetaLs, whereas nucleophile activation by a Zn(II) ion is not a conserved mechanistic feature.

Common mechanistic features among metallo-beta-lactamases: a computational study of Aeromonas hydrophila CphA enzyme

Metallo-beta-lactamases (MbetaLs) constitute an increasingly serious clinical threat by giving rise to beta-lactam antibiotic resistance. They accommodate in their catalytic pocket one or two zinc ions, which are responsible for the hydrolysis of beta-lactams. Recent x-ray studies on a member of the mono-zinc B2 MbetaLs, CphA from Aeromonas hydrophila, have paved the way to mechanistic studies of this important subclass, which is selective for carbapenems. Here we have used hybrid quantum mechanical/molecular mechanical methods to investigate the enzymatic hydrolysis by CphA of the antibiotic biapenem. Our calculations describe the entire reaction and point to a new mechanistic description, which is in agreement with the available experimental evidence. Within our proposal, the zinc ion properly orients the antibiotic while directly activating a second catalytic water molecule for the completion of the hydrolytic cycle. This mechanism provides an explanation for a variety of mutagenesis experiments and points to common functional facets across B2 and B1 MbetaLs.

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.

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.

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

Herein we briefly review theoretical contributions that have increased our understanding of the structure and function of metallo-beta-lactamases and ureases. Both are bimetallic metalloenzymes, with the former containing two zinc ions and the latter containing two nickel ions. We describe the use of several different methodologies, including quantum chemical calculations, molecular dynamic simulations, as well as mixed QM/MM approaches and how they have impacted our understanding of the structure and function of metallo-beta-lactamases and ureases.

Inhibitor Binding by Metallo-β-lactamase IMP-1 from Pseudomonas aeruginosa: Quantum Mechanical/Molecular Mechanical Simulations

The dynamics of the IMP-1 enzyme complexed with three prototypical inhibitors are investigated using a quantum mechanical/molecular mechanical (QM/MM) method based on the self-consistent-charge density-functional tight-binding model. The binding patterns of the inhibitors observed in X-ray diffraction experiments are well reproduced in 600 ps molecular dynamics simulations at room temperature. These inhibitors anchor themselves in the enzyme active site by direct coordination with the two zinc ions, displacing the hydroxide nucleophile that bridges the two zinc ions. In addition, they also interact with several active-site residues and those in two mobile loops. The excellent agreement with experimental structural data validates the QM/MM treatment used in our simulations.

Antibiotic resistance: mono- and dinuclear zinc complexes as metallo-beta-lactamase mimics

Biomimetic systems containing one or two zinc(II) ions supported by phenolate ligands were developed as functional mimics of metallo-beta-lactamase. These complexes were shown to catalytically hydrolyze beta-lactam substrates, such as oxacillin and penicillin G. The dinuclear zinc complex 1, which has a coordinated water molecule, exhibits high beta-lactamase activity, whereas the dinuclear zinc complex 2, which has no water molecules, but labile chloride ligands, shows a much lower activity. The high beta-lactamase activity of complex 1 can be ascribed to the presence of a zinc-bound water molecule that is activated by being hydrogen bonded to acetate substituents. The kinetics of the hydrolysis of oxacillin by complex 1 and the effect of pH on the reaction rates are reported in detail. In addition, the kinetic parameters obtained for the synthetic analogues are compared with those of the natural metallo-beta-lactamase from Bacillus cereus (BcII). To understand the role of the second metal ion in hydrolysis, the syntheses and catalytic activities of two mononuclear complexes (3 and 4) that include coordinated water molecules are described. Interestingly, the mononuclear zinc complexes 3 and 4 also exhibit high activity, supporting the assumption that the second zinc ion is not crucial for the beta-lactamase activity.

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.

Computational studies of the farnesyltransferase ternary complex part I: substrate binding

Farnesyltransferase (FTase) catalyzes the transfer of farnesyl from farnesyl diphosphate (FPP) to cysteine residues at or near the C-terminus of protein acceptors with a CaaX motif (a, aliphatic; X, Met). Farnesylation is a critical modification to many switch proteins involved in the extracellular signal transduction pathway, which facilitates their fixation on the cell membrane where the extracellular signal is processed. The malfunction caused by mutations in these proteins often causes uncontrolled cell reproduction and leads to tumor formation. FTase inhibitors have been extensively studied as potential anticancer agents in recent years, several of which have advanced to different phases of clinical trials. However, the mechanism of this biologically important enzyme has not been firmly established. Understanding how FTase recruits the FPP substrate is the first and foremost step toward further mechanistic investigations and the design of effective FTase inhibitors. Molecular dynamic simulations were carried out on the ternary structure of FTase complexed with the FPP substrate and an acetyl-capped tetrapeptide (acetyl-CVIM), which revealed that the FPP substrate maintains an inactive conformation and the binding of the diphosphate group can be largely attributed to residues R291beta, K164alpha, K294beta, and H248beta. The FPP substrate assumes an extended conformation in the binding site with restricted rotation of the backbone dihedral angles; however, it does not have a well-defined conformation when unbound in solution. This is evident from multinanosecond MD simulations of the FPP substrate in a vacuum and solution. Our conclusion is further supported by theoretical J coupling calculations. Our results on the FPP binding are in good agreement with previous experimental kinetic studies on FTase mutants. The hypothetical conformational activation of the FPP substrate is currently under investigation.

Site-selective binding of Zn(II) to metallo-beta-lactamase L1 from Stenotrophomonas maltophilia

Extended X-ray absorption fine structure studies of the metallo-beta-lactamase L1 from Stenotrophomonas maltophilia containing 1 and 2 equiv of Zn(II) and containing 2 equiv of Zn(II) plus hydrolyzed nitrocefin are presented. The data indicate that the first, catalytically dominant metal ion is bound by L1 at the consensus Zn1 site. The data further suggest that binding of the first metal helps preorganize the ligands for binding of the second metal ion. The di-Zn enzyme displays a well-defined metal-metal interaction at 3.42 A. Reaction with the beta-lactam antibiotic nitrocefin results in a product-bound species, in which the ring-opened lactam rotates in the active site to present the S1 sulfur atom of nitrocefin to one of the metal ions for coordination. The product bridges the two metal ions, with a concomitant lengthening of the Zn-Zn interaction to 3.62 A.

Complexes of thiomandelate and captopril mercaptocarboxylate inhibitors to metallo-β-lactamase by polarizable molecular mechanics. Validation on model binding sites by quantum chemistry

Using the polarizable molecular mechanics method SIBFA, we have performed a search for the most stable binding modes of D- and L-thiomandelate to a 104-residue model of the metallo-beta-lactamase from B. fragilis, an enzyme involved in the acquired resistance of bacteria to antibiotics. Energy balances taking into account solvation effects computed with a continuum reaction field procedure indicated the D-isomer to be more stably bound than the L-one, conform to the experimental result. The most stably bound complex has the S(-) ligand bridging monodentately the two Zn(II) cations and one carboxylate O(-) H-bonded to the Asn193 side chain. We have validated the SIBFA energy results by performing additional SIBFA as well as quantum chemical (QC) calculations on small (88 atoms) model complexes extracted from the 104-residue complexes, which include the residues involved in inhibitor binding. Computations were done in parallel using uncorrelated (HF) as well as correlated (DFT, LMP2, MP2) computations, and the comparisons extended to corresponding captopril complexes (Antony et al., J Comput Chem 2002, 23, 1281). The magnitudes of the SIBFA intermolecular interaction energies were found to correctly reproduce their QC counterparts and their trends for a total of twenty complexes.

Reactions of Pyrazolylborate−Zinc−Hydroxide Complexes Related to β-Lactamase Activity

Simple beta-lactams and their hydrolysis products, the beta-amino acids, react with TpZn-OH under deprotonation. The latter become semibidentate carboxylate ligands with a NH...O hydrogen bond, and the former become N-bound beta-lactamide ligands. Likewise the antibiotic derivatives 6-aminopenicillanic acid and 7-aminocephalosporanic acid are incorporated as carboxylate ligands. beta-Lactams bearing nitrophenyl or acyl substituents at the nitrogen atoms are opened hydrolytically by TpZn-OH, and the resulting N-substituted beta-amino acids are attached to zinc by their carboxylate functions. Only with trifluoroacetyl as the N-substituent does the hydrolytic cleavage occur at the external amide bond, yielding the free beta-lactam and TpZn-trifluoroacetate. The kinetic investigation of the opening reactions has shown them to be of second order like all other TpZn-OH-induced hydrolytic cleavages, thereby supporting the four-center mechanism for the monozinc beta-lactamases.

Structural determinants of substrate binding to Bacillus cereus metallo-beta-lactamase

Binding and hydrolysis of the beta-lactams cefotaxime, cephapirin, imipenem, and benzylpenicillin by the metallo-beta-lactamase from Bacillus cereus were studied by presteady state kinetic measurements. In all cases, the substrate was unmodified in the most populated reaction intermediate, and no chemically modified substrate species accumulated to a detectable amount. The cephalosporins tested showed similar formation rate constants for this intermediate, and they differed mostly in their decay rates. Formation of a non-productive enzyme.substrate complex was detected for imipenem. The substrate binding differences can be accounted for by considering the structural features of each substrate. The apoenzyme could not bind any of the substrates, but binding was restored when the apoenzyme was reconstituted with Zn(II), revealing that the metal ions are the main determinants of substrate binding. This evidence is in line with the lack of an optimized substrate recognition patch in B1 and B3 metallo-beta-lactamases that provides a broad substrate spectrum.

Substrate binding to mononuclear metallo-β-lactamase from Bacillus cereus

Structure and dynamics of substrate binding (cefotaxime) to the catalytic pocket of the mononuclear zinc-beta-lactamase from Bacillus cereus are investigated by molecular dynamics simulations. The calculations, which are based on the hydrogen-bond pattern recently proposed by Dal Peraro et al. (J Biol Inorg Chem 2002; 7:704-712), are carried out for both the free and the complexed enzyme. In the resting state, active site pattern and temperature B-factors are in agreement with crystallographic data. In the complexed form, cefotaxime is accommodated into a stable orientation in the catalytic pocket within the nanosecond timescale, interacting with the enzyme zinc-bound hydroxide and the surrounding loops. The beta-lactam ring remains stable and very close to the hydroxide nucleophile agent, giving a stable representation of the productive enzyme-substrate complex.

Molecular modeling of penicilloate anions: an RHF-SCF analysis

An ab initio restricted Hartree-Fock self-consistent field (RHF-SCF) analysis of penicilloate anions was performed at the TZV level with GAMESS. Geometry optimization was initialized by the semi-empirical AM1 method followed by optimization at the 6-31++G** level. The total energy obtained was -1116.0997 a.u. for the penicilloate amine, -1115.3164 a.u. for the imine, -1115.2969 a.u. for the enamine and -1115.2017 a.u. for the amine that was deprotonated at the thiazolidine nitrogen. Formation of the free thiolate in the imine and enamine anions by deprotonation of the penicilloate amine is associated with: (1) an increase in total energy (2) an increase in the energy of the highest occupied molecular orbital (HOMO) to that of anti-bonding (3) a decrease in chemical hardness (4) an increase in the chemical potential (5) a more negative Mulliken net charge on the sulfur atom and (6) an increase in the Mulliken atomic population on the former thiazolidine sulfur atom in the HOMO. The RHF-SCF analysis presented here suggests a potential role for the thiolate sulfur of penicilloate anions, especially of the imine, as a chemically reactive soft nucleophile.

The 1.5-A structure of Chryseobacterium meningosepticum zinc beta-lactamase in complex with the inhibitor, D-captopril

The crystal structure of the class-B beta-lactamase, BlaB, from the pathogenic bacterium, Chryseobacterium meningosepticum, in complex with the inhibitor, d-captopril, has been solved at 1.5-A resolution. The enzyme has the typical alphabeta/betaalpha metallo-beta-lactamase fold and the characteristic two metal binding sites of members of the subclass B1, in which two Zn2+ ions were identified. d-Captopril, a diastereoisomer of the commercial drug, captopril, acts as an inhibitor by displacing the catalytic hydroxyl ion required for antibiotic hydrolysis and intercalating its sulfhydryl group between the two Zn2+ ions. Interestingly, d-captopril is located on one side of the active site cleft. The x-ray structure of the complex of the closely related enzyme, IMP-1, with a mercaptocarboxylate inhibitor, which also contains a sulfhydryl group bound to the two Zn2+ ions, shows the ligand to be located on the opposite side of the active site cleft. A molecule generated by fusion of these two inhibitors would cover the entire cleft, suggesting an interesting approach to the design of highly specific inhibitors.

A substrate variant as a high-affinity, reversible inhibitor: insight from the X-ray structure of cilastatin bound to membrane dipeptidase

An analysis of the X-ray structure of cilastatin bound to membrane dipeptidase, together with docking studies, is presented here to reveal how a simple amide may act as a high-affinity, reversible, amidase inhibitor. Cilastatin binds as a normal substrate and is orientated in a perfect near-attack conformer for formation of a tetrahedral intermediate with the zinc-bound water/hydroxide. This intermediate is fated, however, only to revert to its starting components as scission of the amide bond is prevented by the precise fit of cilastatin within the active site. The cilastatin alkyl end groups that are tightly buttressed against amino acid residues on opposite sides of the active site, are aligned along the C-N reaction coordinate axis thereby preventing collapse of the intermediate via rupture of the C-N bond. Such a feature could have more general applicability in the explicit design of substrate variants as selective, tight-binding, and reversible inhibitors.

Molecular dynamics simulations of the dinuclear zinc-beta-lactamase from Bacteroides fragilis complexed with imipenem

Herein, we present results from MD simulations of the Michaelis complex formed between the dizinc beta-lactamase from B. fragilis and imipenem. We considered two catalytically important configurations, which differ in the presence or absence of a hydroxide bridge connecting the two zinc ions in the active site. The structural and dynamical effects induced by substrate binding, the specific roles of the conserved residues and the zinc-bound water molecules, the near attack conformers of the Michaelis complex, and so forth, are discussed in detail. The relative stability of the two configurations was estimated from QM linear scaling calculations on the enzyme-substrate complex combined with Poisson-Boltzmann electrostatic calculations and normal mode calculations. Importantly, we find that the two configurations have similar energies, indicating that these two structures could readily be interchanged, thereby facilitating catalysis. The configuration with the hydroxide bound to the two zinc ions is predicted to be the resting form of the enzyme, while the configuration without the bridge is the reactive form that was found to place the hydroxide in position to attack the carbonyl of the beta-lactam ring. Thus, we propose that the enzyme initiates catalysis by converting from the hydroxide bridge form into the configuration that lacks the hydroxide bridge. This interconversion increases the nucleophilicity of the hydroxide ion and exposes it to the beta-lactam carbonyl, which ultimately facilitates nucleophilic attack. The implications of the observed modes of binding, the possible influence of mutating the Lys184 and Asn193 residues on substrate binding, and the reaction mechanism are also discussed in detail.

A new zinc-protein coordination site in intracellular metal trafficking: solution structure of the Apo and Zn(II) forms of ZntA(46-118)

Zinc, a metal ion that functions in a wide variety of catalytic and structural sites in metalloproteins, is shown here to adopt a novel coordination environment in the Escherichia coli transport protein ZntA. The ZntA protein is a P-type ATPase that pumps zinc out of the cytoplasm and into the periplasm. It is physiologically selective for Zn(II) and functions with metalloregulatory proteins in the cell to keep the zinc quota within strict limits. Yet, the N-terminal cytoplasmic domain contains a region that is highly homologous to the yeast Cu(I) metallochaperone Atx1. To investigate how the structure of this region may influence its function, this fragment, containing residues 46-118, has been cloned out of the gene and overexpressed. We report here the solution structure of this fragment as determined by NMR. Both the apo and Zn(II)-ZntA(46-118) structures have been determined. It contains a previously unknown protein coordination site for zinc that includes two cysteine residues, Cys59 and Cys62, and a carboxylate residue, Asp58. The solvent accessibility of this site is also remarkably high, a feature that increasingly appears to be a characteristic of domains of heavy metal ion transport proteins. The participation of Asp58 in this ZntA metal ion binding site may play an important role in modulating the relative affinities and metal exchange rates for Zn(II)/Pb(II)/Cd(II) as compared with other P-type ATPases, which are selective for Cu(I) or Ag(I).

Binding of D- and L-captopril inhibitors to metallo-beta-lactamase studied by polarizable molecular mechanics and quantum mechanics

The bacterial Zn2+ metallo-beta-lactamase from B. fragilis is a zinc-enzyme with two potential metal ion binding sites. It cleaves the lactam ring of antibiotics, thus contributing to the acquired resistance of bacteria against antibiotics. The present study bears on the binuclear form of the enzyme. We compare several possible binding modes of captopril, a mercaptocarboxamide inhibitor of several zinc-metalloenzymes. Two diastereoisomers of captopril were considered, with either a D- or an L-proline residue. We have used the polarizable molecular mechanics procedure SIBFA (Sum of Interactions Between Fragments ab initio computed). Two beta-lactamase models were considered, encompassing 104 and 188 residues, respectively. The energy balances included the inter and intramolecular interaction energies as well as the contribution from solvation computed using a continuum reaction field procedure. The thiolate ion of the inhibitor is binding to both metal ions, expelling the bridging solvent molecule from the uncomplexed enzyme. Different competing binding modes of captopril were considered, either where the inhibitor binds in a monodentate mode to the zinc cations only with its thiolate ion, or in bidentate modes involving additional zinc binding by its carboxylate or ketone carbonyl groups. The additional coordination by the inhibitor's carboxylate or carbonyl group always occurs at the zinc ion, which is bound by a histidine, a cysteine, and an aspartate side chain. For both diastereomers, the energy balances favor monodentate binding of captopril via S-. The preference over bidentate binding is small. The interaction energies were recomputed in model sites restricted to captopril, the Zn2+ cations, and their coordinating end side chains from beta-lactamase (98 atoms). The interaction energies and their ranking among competing arrangements were consistent with those computed by ab initio HF and DFT procedures.

IMP-1 metallo-beta-lactamase: effect of chelators and assessment of metal requirement by electrospray mass spectrometry

Metallo-beta-lactamases have attracted considerable attention due to their role in microbial resistance to beta-lactam antibiotics. IMP-1, the binuclear Zn-dependent beta-lactamase produced by Pseudomonas aeruginosa and other microorganisms, is of particular interest in view of its increasing prevalence. An examination of the susceptibility of IMP-1 to inactivation by six different divalent metal ion chelators has revealed that all except Zincon cause inhibition by forming a complex with the holoenzyme. Exposure of the enzyme to dipicolinic acid (DPA), the most potent inhibitor, results in the production of the mononuclear Zn form of the protein as determined by electrospray ionization mass spectrometry (ESI-MS) under nondenaturing conditions. This mononuclear Zn species was found to be catalytically competent. Studies with the chromophoric chelator 4-(2-pyridylazo)resorcinol (PAR) show that the two zinc centers in IMP-1 differ in their accessibility, a feature that could be overcome in the presence of guanidine hydrochloride (GdnHCl, 1.5 M).