Class A β-lactamases and inhibitors: In silico analysis of the binding mode and the relationship with resistance

In silico analysis of luteolin derivatives as antibacterial agents targeting DNA gyrase and CTX-M-15 extended-spectrum β-lactamase of Escherichia coli

Luteolin exhibited antibacterial activity against Escherichia coli and its chemical structure similar to that of ciprofloxacin (CPF) which works by inhibiting DNA gyrase. Filtrate from passion fruit extract containing luteolin and its derivatives could inhibit extended-spectrum β-lactamase (ESBL)-producing E. coli. Antibacterial compounds that can also inhibit ESBL will be valuable compounds to overcome the problem of resistant bacteria. This study aimed to ensure the potency of luteolin and luteolin derivatives targeting DNA gyrase and ESBL by in silico approach. Docking simulation of ligands L1-L14 was performed using AutoDock Vina, and pharmacokinetics and toxicity (absorption, distribution, metabolism, excretion, and toxicity) profiles were predicted by pKCSM online. The docking result revealed higher binding affinity on DNA gyrase (PDB.1KZN) of 12 luteolin derivatives (energy <-7.6 kcal/mol) compared to CPF and higher affinity (energy <-6.27 kcal/mol) of all compounds than clavulanic acid against ESBL CTX-M-15 (PDB.4HBU). The compounds could be absorbed through the human intestine moderately, which showed low permeability to blood-brain barrier, nontoxic and nonhepatotoxic. The most active luteolin glycoside (L6) is capable to inhibit DNA gyrase and ESBL from E. coli which provided the potential against resistant bacteria and was promoted as lead compounds to be developed further.

Analysis of a novel class A β-lactamase OKP-B-6 of Klebsiella quasipneumoniae: structural characterisation and interaction with commercially available drugs

BACKGROUND

Gram-negative and Gram-positive bacteria produce beta-lactamase as factors to overcome beta-lactam antibiotics, causing their hydrolysis and impaired antimicrobial action. Class A beta-lactamase contains the chromosomal sulfhydryl reagent variable (SHV, point mutation variants of SHV-1), LEN (Klebsiella pneumoniae strain LEN-1), and other K. pneumoniae beta-lactamase (OKP) found mostly in Klebsiella’s phylogroups. The SHV known as extended-spectrum β-lactamase can inactivate most beta-lactam antibiotics. Class A also includes the worrisome plasmid-encoded Klebsiella pneumoniae carbapenemase (KPC-2), a carbapenemase that can inactivate most beta-lactam antibiotics, carbapenems, and some beta-lactamase inhibitors.

OBJECTIVES

So far, there is no 3D crystal structure for OKP-B, so our goal was to perform structural characterisation and molecular docking studies of this new enzyme.

METHODS

We applied a homology modelling method to build the OKP-B-6 structure, which was compared with SHV-1 and KPC-2 according to their electrostatic potentials at the active site. Using the DockThor-VS, we performed molecular docking of the SHV-1 inhibitors commercially available as sulbactam, tazobactam, and avibactam against the constructed model of OKP-B-6.

FINDINGS

From the point of view of enzyme inhibition, our results indicate that OKP-B-6 should be an extended-spectrum beta-lactamase (ESBL) susceptible to the same drugs as SHV-1.

MAIN CONCLUSIONS

This conclusion advantageously impacts the clinical control of the bacterial pathogens encoding OKP-B in their genome by using any effective, broad-spectrum, and multitarget inhibitor against SHV-containing bacteria.

A multiscale approach to predict the binding mode of metallo beta-lactamase inhibitors

Antibiotic resistance is a major threat to global public health. β-lactamases, which catalyze breakdown of β-lactam antibiotics, are a principal cause. Metallo β-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all β-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems.