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

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.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Mutational Effects on Carbapenem Hydrolysis of YEM-1, a New Subclass B2 Metallo-β-Lactamase from Yersinia mollaretii

Analysis of the genome sequence of Yersinia mollaretii ATCC 43969 identified the bla _YEM gene, encoding YEM-1, a putative subclass B2 metallo-β-lactamase. The objectives of our work were to produce and purify YEM-1 and to complete its kinetic characterization. YEM-1 displayed the narrowest substrate range among known subclass B2 metallo-β-lactamases, since it can hydrolyze imipenem, but not other carbapenems, such as biapenem, meropenem, doripenem, and ertapenem, with high catalytic efficiency. A possible explanation of this activity profile is the presence of tyrosine at residue 67 (loop L1), threonine at residue 156 (loop L2), and serine at residue 236 (loop L3). We showed that replacement of Y67 broadened the activity profile of the enzyme for all carbapenems but still resulted in poor activity toward the other β-lactam classes.

Recognition and interaction of CphA from Aeromonas hydrophila with imipenem and biapenem by spectroscopic analysis in combination with molecular docking

The molecular recognition and interaction of CphA from Aeromonas hydrophila with imipenem (Imip) and biapenem (Biap) were studied by means of the combined use of fluorescence spectra and molecular docking. The results showed that both the fluorescence quenching of CphA by Imip and Biap were caused through the combined dynamic and static quenching, and the latter was dominating in the process; the microenvironment and conformational of CphA were altered upon the addition of Imip and Biap from synchronous and three-dimensional fluorescence. The binding of CphA with Imip or Biap caused a conformational change in the loop of CphA, and through the conformational change, the loop opened the binding pocket of CphA to allow for an induced fit of the newly introduced ligand. In the binding of CphA with Imip, the whole molecule entered into the active pocket of CphA. The binding was driven by enthalpy change, and the binding force between them was mainly hydrogen bonding and Van der Waals force; whereas in the binding of CphA with Biap, only the beta-lactam ring of Biap entered into the binding pocket of CphA while the side chain was located outside the active pocket. The binding was driven by the enthalpy change and entropy change together, and the binding force between them was mainly electrostatic interaction. This study provided an insight into the recognition and binding of CphA with antibiotics, which may be helpful for designing new substrate for beta-lactamase and developing new antibiotics resistant to superbugs.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

Shaping Substrate Selectivity in a Broad-Spectrum Metallo-β-Lactamase

Metallo-β-lactamases (MBLs) are the major group of carbapenemases produced by bacterial pathogens. The design of MBL inhibitors has been limited by, among other issues, incomplete knowledge about how these enzymes modulate substrate recognition. While most MBLs are broad-spectrum enzymes, B2 MBLs are exclusive carbapenemases. This narrower substrate profile has been attributed to a sequence insertion present in B2 enzymes that limits accessibility to the active site. In this work, we evaluate the role of sequence insertions naturally occurring in the B2 enzyme Sfh-I and in the broad-spectrum B1 enzyme SPM-1. We engineered a chimeric protein in which the sequence insertion of SPM-1 was replaced by the one present in Sfh-I. The chimeric variant is a selective cephalosporinase, revealing that the substrate profile of MBLs can be further tuned depending on the protein context. These results also show that the stable scaffold of MBLs allows a modular engineering much richer than the one observed in nature.

Theoretical studies of the hydrolysis of antibiotics catalyzed by a metallo-β-lactamase

In this paper, hybrid QM/MM molecular dynamics (MD) simulations have been performed to explore the mechanisms of hydrolysis of two antibiotics, Imipenen (IMI), an antibiotic belonging to the subgroup of carbapenems, and the Cefotaxime (CEF), a third-generation cephalosporin antibiotic, in the active site of a mono-nuclear β-lactamase, CphA from Aeromonas hydrophila. Significant different transition state structures are obtained for the hydrolysis of both antibiotics: while the TS of the CEF is an ionic species with negative charge on nitrogen, the IMI TS presents a tetrahedral-like character with negative charge on oxygen atom of the carbonyl group of the lactam ring. Thus, dramatic conformational changes can take place in the cavity of CphA to accommodate different substrates, which would be the origin of its substrate promiscuity. Since CphA shows only activity against carbapenem antibiotic, this study sheds some light into the origin of the selectivity of the different MbL and, as a consequence, into the discovery of specific and potent MβL inhibitors against a broad spectrum of bacterial pathogens. We have finally probed that a re-parametrization of semiempirical methods should be done to properly describe the behavior the metal cation in active site, Zn(2+), when used in QM/MM calculations.

QM/MM simulations of amyloid-β 42 degradation by IDE in the presence and absence of ATP

The ability of the insulin-degrading enzyme (IDE) to degrade amyloid-β 42 (Aβ42), a process regulated by ATP, has been studied as an alternative path in the development of drugs against Alzheimer's disease. In this study, we calculated the potential of mean force for the degradation of Aβ42 by IDE in the presence and absence of ATP by umbrella sampling with hybrid quantum mechanics and molecular mechanics (QM/MM) calculations, using the SCC-DFTB QM Hamiltonian and Amber ff99SB force field. Results indicate that the reaction occurs in two steps: The first step is characterized by the formation of the intermediate. The second step is characterized by breaking the peptide bond of the substrate, the latter being the rate-determining step. In our simulations, the activation energy barrier in the absence of ATP is 15 ± 2 kcal mol(-1), which is 7 kcal mol(-1) lower than in the presence of ATP, indicating that the presence of the nucleotide decreases the reaction rate by about 10(5) times.

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.

A variety of roles for versatile zinc in metallo-β-lactamases

β-Lactamases inactivate the important β-lactam antibiotics by catalysing the hydrolysis of the β-lactam ring, thus. One class of these enzymes, the metallo-β-lactamases, bind two zinc ions at the active site and these play important roles in the catalytic mechanism.

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.

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

BACKGROUND

The first line of defense by bacteria against β-lactam antibiotics is the expression of β-lactamases, which cleave the amide bond of the β-lactam ring. In the reaction of biapenem inactivation by B2 metallo β-lactamases (MβLs), after the β-lactam ring is opened, the carboxyl group generated by the hydrolytic process and the hydroxyethyl group (common to all carbapenems) rotate around the C5-C6 bond, assuming a new position that allows a proton transfer from the hydroxyethyl group to C2, and a nucleophilic attack on C3 by the oxygen atom of the same side-chain. This process leads to the formation of a bicyclic compound, as originally observed in the X-ray structure of the metallo β-lactamase CphA in complex with product.

METHODOLOGY/PRINCIPAL FINDINGS

QM/MM and metadynamics simulations of the post-hydrolysis steps in solution and in the enzyme reveal that while the rotation of the hydroxyethyl group can occur in solution or in the enzyme active site, formation of the bicyclic compound occurs primarily in solution, after which the final product binds back to the enzyme. The calculations also suggest that the rotation and cyclization steps can occur at a rate comparable to that observed experimentally for the enzymatic inactivation of biapenem only if the hydrolysis reaction leaves the N4 nitrogen of the β-lactam ring unprotonated.

CONCLUSIONS/SIGNIFICANCE

The calculations support the existence of a common mechanism (in which ionized N4 is the leaving group) for carbapenems hydrolysis in all MβLs, and suggest a possible revision of mechanisms for B2 MβLs in which the cleavage of the β-lactam ring is associated with or immediately followed by protonation of N4. The study also indicates that the bicyclic derivative of biapenem has significant affinity for B2 MβLs, and that it may be possible to obtain clinically effective inhibitors of these enzymes by modification of this lead compound.

Carbapenems: past, present, and future

In this review, we summarize the current "state of the art" of carbapenem antibiotics and their role in our antimicrobial armamentarium. Among the β-lactams currently available, carbapenems are unique because they are relatively resistant to hydrolysis by most β-lactamases, in some cases act as "slow substrates" or inhibitors of β-lactamases, and still target penicillin binding proteins. This "value-added feature" of inhibiting β-lactamases serves as a major rationale for expansion of this class of β-lactams. We describe the initial discovery and development of the carbapenem family of β-lactams. Of the early carbapenems evaluated, thienamycin demonstrated the greatest antimicrobial activity and became the parent compound for all subsequent carbapenems. To date, more than 80 compounds with mostly improved antimicrobial properties, compared to those of thienamycin, are described in the literature. We also highlight important features of the carbapenems that are presently in clinical use: imipenem-cilastatin, meropenem, ertapenem, doripenem, panipenem-betamipron, and biapenem. In closing, we emphasize some major challenges and urge the medicinal chemist to continue development of these versatile and potent compounds, as they have served us well for more than 3 decades.

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.

Engineered mononuclear variants in Bacillus cereus metallo-beta-lactamase BcII are inactive

Metallo-beta-lactamases (MbetaLs) are zinc enzymes able to hydrolyze almost all beta-lactam antibiotics, rendering them inactive, at the same time endowing bacteria high levels of resistance. The design of inhibitors active against all classes of MbetaLs has been hampered by their structural diversity and by the heterogeneity in metal content in enzymes from different sources. BcII is the metallo-beta-lactamase from Bacillus cereus, which is found in both the mononuclear and dinuclear forms. Despite extensive studies, there is still controversy about the nature of the active BcII species. Here we have designed two mutant enzymes in which each one of the metal binding sites was selectively removed. Both mutants were almost inactive, despite preserving most of the structural features of each metal site. These results reveal that neither site isolated in the MbetaL scaffold is sufficient to render a fully active enzyme. This suggests that only the dinuclear species is active or that the mononuclear variants can be active only if aided by other residues that would be metal ligands in the dinuclear species.

Mutational analysis of the zinc- and substrate-binding sites in the CphA metallo-beta-lactamase from Aeromonas hydrophila

The subclass B2 CphA (Carbapenemase hydrolysing Aeromonas) beta-lactamase from Aeromonas hydrophila is a Zn(2+)-containing enzyme that specifically hydrolyses carbapenems. In an effort to evaluate residues potentially involved in metal binding and/or catalysis (His(118), Asp(120), His(196) and His(263)) and in substrate specificity (Val(67), Thr(157), Lys(224) and Lys(226)), site-directed mutants of CphA were generated and characterized. Our results confirm that the first zinc ion is in interaction with Asp(120) and His(263), and thus is located in the 'cysteine' zinc-binding site. His(118) and His(196) residues seem to be interacting with the second zinc ion, as their replacement by alanine residues has a negative effect on the affinity for this second metal ion. Val(67) plays a significant role in the binding of biapenem and benzylpenicillin. The properties of a mutant with a five residue (LFKHV) insertion just after Val(67) also reveals the importance of this region for substrate binding. This latter mutant has a higher affinity for the second zinc ion than wild-type CphA. The T157A mutant exhibits a significantly modified activity spectrum. Analysis of the K224Q and N116H/N220G/K224Q mutants suggests a significant role for Lys(224) in the binding of substrate. Lys(226) is not essential for the binding and hydrolysis of substrates. Thus the present paper helps to elucidate the position of the second zinc ion, which was controversial, and to identify residues important for substrate binding.