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

Abstract

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.

Key words: class A β-lactamase; OKP-B-6; multitarget antibiotics; Klebsiella quasipneumoniae; structure-based drug design; molecular docking

Beta-lactamases are enzymes produced by Gram-negative and Gram-positive bacteria as an evolutionary response against naturally occurring beta-lactams and constitute a significant resistance factor against beta-lactam antibiotics to control bacterial infections. The beta-lactam antibiotics target the bacterial cell wall, leading to cell lysis and death. However, in the presence of beta-lactamases, the beta-lactam drug is hydrolysed and inactivated, thereby preventing it from reaching the cell wall target. A thriving pharmaceutical industry strategy has been the beta-lactamase inhibitors (BLIs) designed to overcome this beta-lactamase-mediated inactivation of beta-lactam drugs. Thus, while the BLI binds to the beta-lactamase and interrupts its function, it enables the bactericidal effect of the beta-lactam antibiotics against pathogenic bacteria.¹

Structurally, beta-lactamases are similar to penicillin-binding proteins (PBPs).^2,3 Indeed, both enzymes acquired deacylation machinery to inactivate or hydrolyse beta-lactam antibiotics and have an α-helix domain and a mixed α/β domain, being the active site localised near the amino terminus of the central α-helix.⁴ The Ambler classification recognises four classes of beta-lactamases, namely A-D.² While classes A, C, and D are active-site serine beta-lactamases, class B are structurally unrelated metallo beta-lactamases that use one or two Zn⁺² catalytically functional atoms for beta-lactam hydrolysis.^2,5 Mainly, class A includes an alarming group of extended-spectrum beta-lactamases (ESBLs) that can inactivate most beta-lactam antibiotics and carbapenemases that inactivate even carbapenems and some BLIs yielding a variant enzyme designated “inhibitor-resistant”,⁶ making infections caused by ESBL-producing and carbapenemase-producing bacteria tough to treat. The most worrying carbapenemase is the carbapenem-hydrolysing enzyme Klebsiella pneumoniae carbapenemase (KPC-2) found in Klebsiella phylogroups and other enteric bacteria.^7,8,9

To overcome the activity of the ESBLs and carbapenemases, pharmaceutical challenges have focused on either the new generation of beta-lactams less susceptible to cleavage (i.e., fifth-generation cephalosporins and ceftolozane, fourth-generation carbapenems) or BLIs with better inhibitory activity (i.e., avibactam and vaborbactam).¹⁰ In the context of small molecules inhibiting beta-lactamases, structure-based computational methods have been applied to provide a rationale for the design of selective and broad-spectrum BLIs.^11,12 For example, all “first generation” BLIs (e.g., clavulanic acid, sulbactam, and tazobactam) function as mechanism-based compounds that utilise various strategies to inactivate the beta-lactamase. The “first generation” BLIs can promote semi-stable intermediates to inhibit serine beta-lactamases by preventing the deacylation machinery of this enzyme, or these BLIs yield covalent modifications on either Ser130 or the catalytic Ser70 residue on the active site of beta-lactamases.^{11,12,13,14,15} In the class of “second generation” BLIs, i.e., avibactam (NXL104), approved by the FDA in 2015 as a ceftazidime/avibactam formulation, it can form an acyl-enzyme complex with a serine residue at the beta-lactamase binding site. Interestingly, this BLI can be removed from the beta-lactamase via deacylation, yielding an intact liberated inhibitor molecule.¹⁶ However, the carbapenemase KPC-2, carrying some amino acid substitutions in the Ω loop of the enzyme, is even capable of slowly desulfating avibactam, triggering the production of inactive BLI fragments.¹⁷ Currently, for treating strains producing the KPC-2 variant conferring resistance to ceftazidime-avibactam, the meropenem-vaborbactam combination becomes a valuable agent.¹⁸ It is important to note that clavulanic acid and penicillanic acid sulfones, sulbactam, and tazobactam exhibit a similar spectrum of inhibition. These inhibitors are recognised to inhibit most class A beta-lactamases, including ESBLs, and serine carbapenemases to a lesser extent.¹⁹ On the other hand, avibactam has a broader spectrum of activity than clavulanic acid and the sulfone inhibitors, potently inhibiting class-A penicillinases, ESBLs, serine carbapenemases, class-C cephalosporinases and some class-D oxacillinases.²⁰

Although most beta-lactamases described to date are plasmid-encoded enzymes associated with hospital outbreaks, like the KPCs that increasingly spread among clinical isolates,^21,22,23 there are specific chromosomally encoded beta-lactamase genes also identified from clinical isolates.^9,24 The major variants of a chromosomal beta-lactamase gene, like bla _SHV, bla _LEN, and bla _OKP (OKP for other K. pneumoniae beta-lactamase)²⁵ evolved in parallel with the diversification of three K. pneumoniae phylogenetic groups, namely KpI (K. pneumoniae), KpII (K. quasipneumoniae), and KpIII (K. variicola), but without concomitant phenotypic changes in beta-lactam susceptibility.^24,26 Usually, the beta-lactamases variants could be a distinctive feature of these Klebsiella species since sulfhydryl reagent variable (SHV), point mutation variants of SHV-1 is predominantly found in K. pneumoniae, OKP in K. quasipneumoniae, and LEN (K. pneumoniae strain LEN-1) in K. variicola.⁹ Concerning the K. quasipneumoniae, the first isolates were described in 2014, belonging to subspecies K. quasipneumoniae subsp. quasipneumoniae (KpII-A) and K. quasipneumoniae subsp. similipneumoniae (KpII-B), carrying the divergent genes bla _OKP-A and bla _OKP-B in their chromosomes.^25,26,27 A few studies have reported isolates belonging to this new subspecies describing their antimicrobial resistance genes repertoires.^28,29,30,31 In Brazil, our group described at the genomic level a novel multidrug-resistant (MDR) clinical isolate of K. quasipneumoniae subsp. similipneumoniae, carrying both the OKP-B-6 variant enzyme and the plasmid-encoded beta-lactamase KPC-2. The new isolate called Kqps142 was associated with a hospital outbreak in southeast Brazil, representing a significant concern in public health.³²

Considering the medical importance of beta-lactamases, the structures of a considerable number of these enzymes were elucidated, even various in complex with approved BLIs (more than 5,370 beta-lactamase structures from Bacteria in protein data bank (PDB), http://www.rcsb.org/pdb/ accessed on March 02, 2022). Most of the beta-lactamase PDB entries contain structural data from the ESBL class A, such as TEM [named after the patient (Temoneira) infected with the first isolate expressing TEM-1], cefotaxime-hydrolysing beta-lactamase isolated in Munich (CTX-M), SHV, and the carbapenemase KPC-2. However, no structure is yet available in PDB for OKP variant enzymes. Regarding the emergence of a novel variant beta-lactamase, such as OKP-B-6 of K. quasipneumoniae subsp. similipneumoniae Kqps142, diverse approaches to gain structural insight into these new enzymes are critically required. Since no 3D crystal structure has been determined for OKP-B so far, our goal was to perform structural characterisation and molecular docking studies of this new variant enzyme through its interaction with commercial BLIs.

MATERIALS AND METHODS

Sequence availability - We retrieved the OKP-B-6 sequence from the genome of K. quasipneumoniae subsp. similipneumoniae Kqps142 (NCBI Bioproject number PRJNA383559), annotated as an OKP-B-6 beta-lactamase corresponding to GenBank accession no. AVR37864.1 and locus tag KPC142_02249 (https://www.ncbi.nlm.nih.gov/protein/AVR37864.1). We confirmed the final annotation as an OKP-B-6 beta-lactamase by the best match reference hit found in K. pneumoniae by “The Comprehensive Antibiotic Resistance Database” (CARD) with E-value = 0.0, 99% identity, and 100% coverage, which is also in agreement with Fevre et al.²⁷ In the text, the number of the amino acid residues from the template sequence will be mentioned in parenthesis when discussing about the OKP-B-6 sequence and structure.

Signal peptide cleavage - The β-lactamases are secreted proteins (preproteins) containing an N-terminal signal peptide (SP) essential for targeting the cytoplasmic membrane and the translocation of the preprotein across the membrane. Then, upon translocating the secretory protein, the SP is cleaved off.³³ According to the post-translational modification (PTM), we removed the SP from the original OKP-B-6 amino acid sequence using the SignalP 4.1 server (www.cbs.dtu.dk/services/SignalP) to proceed with the molecular modelling experiments.

Comparative modelling - We used the Swiss-Model server to construct the molecular 3D model of OKP-B-6 (https://swissmodel.expasy.org/).³⁴ The selection of the template structure was based on the amino acid sequence identity and the presence of ligands in the active site. After the generation of the tridimensional model, we verified the stereochemical qualities by MOLPROBITY^35,36 (Available at http://molprobity.biochem.duke.edu), PROCHECK,^37,38 and ERRAT³⁹ (available at https://saves.mbi.ucla.edu/). We compared the catalytic site of OKP-B-6 with both SHV-1 (PDB code 1VM1)⁴⁰ and KPC-2 (PDB code 4ZBE)¹⁶ structures according to their electrostatic potentials calculated using the PDB2PQR⁴¹ and the Adaptive Poisson Boltzmann Solver server (APBS) available at http://server.poissonboltzmann.org/. The potential surface of the structures was generated with the Pymol v.2.0.7⁴² using the CHARMM force field considering a probe radius of 1.4 Å. The sulbactam binding cavity was characterised by DogSiteScorer, an automated pocket detection and analysis tool that used a grid-based method and support vector machine (SVM) for druggability assessment based on the 3D structure of the target protein.⁴³

Preparation of the OKP-B-6 model and the structures of inhibitors - We prepared the modeled structure of OKP-B-6 using the Protein Preparation Wizard tool from the Maestro suite (Schrödinger, LLC, New York, NY, 2020) considering two conserved water molecules in the binding site (Wat1 and Wat2, respectively numbered 636 and 656 in the 1VM1 template structure).⁴⁰ These water molecules were considered since they are bind tightly at the binding site interacting with at least two amino acid residues, mediate hydrogen bonds between protein and ligand, and are conserved in other beta-lactamase structures [e.g., PDB codes 4FH2 (SHV-1 S70C mutant complexed with sulbactam), 4ZAM (SHV-1 complexed with avibactam) and 4ZBE (KPC-2 complexed with avibactam)]. To predict the protonation states of the amino acid residues, we used PROPKA at pH = 7 and optimised the hydrogen bond network to adjust the orientation of the polar hydrogen atoms.⁴⁴ According to the proposed mechanism provided by a recent study based on ultrahigh-resolution X-ray crystallography of the CTX-M-14 covalently bound to avibactam at different pH values, the apo form of class A β-lactamases has Lys48(73) positively charged and Glu141 (166) deprotonated,⁴⁵ in agreement with previous studies.^46,47 We also adopted this protonation state for Lys48 (73) and Glu141 (166) in OKP-B-6, which were predicted by PROPKA with pKa values of 9.48 and 4.95, respectively. The three-dimensional structures of the SHV-1 inhibitors tazobactam, sulbactam, and avibactam were retrieved from PubChem⁴⁸ and prepared using Maestro to predict the protonation states and tautomers with Epik at pH 7.2 ± 0.2⁴⁹ (Fig. 1).

Fig. 1:
3D and 2D representation of the docked inhibitors with their experimentally measured binding affinities.^50,60 The predicted pKa values of ionizable groups of sulbactam (CID 130313), tazobactam (CID 123630), and avibactam (CID 9835049) are labeled in red and blue.

Docking methodology - Molecular docking is a structure-based drug design strategy that aims to predict potential binding modes and the associated affinity of protein-ligand complexes, being extensively used to better understand the interaction profile between small molecules and therapeutic targets from many pathogens of clinical interest.^{11,51,52,53,54} In this work, we predicted the most favorable binding modes of the beta-lactamases inhibitors sulbactam, tazobactam, and avibactam against the constructed OKP-B-6 model through molecular docking with the DockThor program^55,56 (freely available at www.dockthor.lncc.br). We compared the binding mode obtained for the OKP-B-sulbactam complex with the pre-acylation state of this inhibitor complexed with the Ser70Cys variant of SHV-1 (PDB code 4FH2)⁵⁷ to validate the model constructed for OKP-B-6. The docking experiments were submitted using the standard configuration of the search algorithm: 1,000,000 evaluations, population size of 750, and 24 independent runs. The grid was centered on the sulbactam complexed with SHV1 in the pre-acylation state aligned to the OKP-B-6 model (X = 61.90, Y = 64.15, Z = 38.48) with a total size of 22 Å on each dimension and discretisation of 0.25 Å. We clustered the binding modes found using the dtstatistic tool from the DockThor platform. We ranked the top-energy binding modes according to the Total Energy (E_total, the default DockThor scoring function for pose prediction composed of the torsional, van der Waals, and electrostatic terms from the MMFF94S force field).^58,59

RESULTS AND DISCUSSION

Comparative modelling - Firstly, the SP sequence was removed, which was predicted at the N-terminal region from residues 1 to 21. Then, for the model building step, we considered the amino acid sequence from residues Ser22 to Arg286 as the mature chain of the OKP-B-6 enzyme mapped in the interval from 1 to 265 amino acids (Fig. 2A). The OKP-B-6 model was constructed using as template the crystal structure of SHV-1 beta-lactamase from K. pneumoniae inhibited by tazobactam (PDB code 1VM1, resolution of 2.02 Å and R-value work of 0.17),⁴⁰ sharing 90.19% sequence identity. The Ramachandran plots confirmed good stereochemical quality, with 91.3% of residues in the most favored regions. The normalised QMEAN score is in the range [0.5,1.0], and the Z-score indicates values in typical ranges for similar proteins in the PDB with values of 0.22 ± 0.52. The ERRAT showed a quality factor of 94.16, comparable to other models obtained for beta-lactamases, where the value of ERRAT was 80.92.³⁹ Due to the high identity of sequence with the SHV-1 template, the structure of OKP-B-6 conserved the secondary structures, composed of 13 α-helices (H1 to H13) and three β-sheets (A, B and C) (Fig. 2). The predicted protein structure is a monomer of 265 amino acids with one disulfide bond between residues Cys52 and Cys98, differing from the SHV-1 protein in only 25 residues that fall outside the binding pockets of interest.

Fig. 2:
comparison of the 2D topology and 3D structure of other Klebsiella pneumoniae beta-lactamase (OKP-B-6) (coloured deep teal) andsulfhydryl reagent variable, point mutation variants of SHV-1 (coloured salmon). In the top is the alignment of sequences of the template SHV-1 and the OKP-B-6. Numbering in SHV-1 is according to the Ambler consensus for class A beta-lactamases, where the catalytic serine is set as Ser70. Non-identical residues are shaded in purple. The α-helices are depicted as H1 to H13. Turns are represented as β, and gamma turns are represented as γ. É highlights beta-hairpin. Strands are labeled as A, B, and C. The number 1 in yellow shows the disulfide bonds. Black spheres indicate the position at each 10 amino acid residues. Below is the structural alignment of the 3D structures of the template SHV-1 protein data bank (PDB) code 1VM1, coloured salmon and the constructed model of OKP-B-6 (coloured teal). The non-identical residues are represented in sticks. The structure of the pre-acylation state of sulbactam (PDB code 4FH2, aligned to the OKP-B-6 model) is represented as transparent spheres and sticks to highlight the binding site location.

As expected from the sequence alignment and structural similarity, the residues at the binding site are highly conserved between SHV-1 and OKP-B-6 (Figs 2-3), including the catalytically important residues Lys48, Ser105, and Glu141, suggesting that OKP-B-6 is also susceptible to the same substrate and inhibitors. Furthermore, the residues involved with beta-lactamases resistance are also conserved in OKP-B-6, mainly Met44, Ser105, Lys206, Arg218, and Asn249.⁶⁰

Fig. 3:
three-dimensional model of other Klebsiella pneumoniae beta-lactamase (OKP-B-6) (left) and comparison between OKP-B-6 and sulfhydryl reagent variable, point mutation variants of SHV-1 binding sites (right). Structural waters present are represented as red spheres. The sticks of main chains are hidden for clarity purposes, except for Ala212 (SHV-1 Ala237). The structure of the pre-acylation state of sulbactam protein data bank (PDB) code 4FH2, aligned to the OKP-B-6 model, highlights the binding site location.

We also analysed the differences between the catalytic sites of OKP-B-6, SVH-1, and KPC-2 through the comparison of their electrostatic profiles and volumes (Fig. 4). All enzymes exhibited a positively charged environment in the carboxylate binding pocket, with OKP-B-6 and SVH-1 exhibiting a larger pocket with an electrostatic potential more pronounced than that of KPC-2 in this region. The positive charge on both catalytic sites is known to be critical for guiding the negatively charged substrates. In Fig. 5, the sulbactam binding cavity, as predicted with DogSiteScorer, indicates that this pocket is slightly smaller in OKP-B-6 than in the SHV-1 structure, with a volume of 255.87 ų (292.54 ų in SHV-1) and a surface of 261.18 Ų (280.16 Ų in SHV-1), and similar druggability scores of DrugScore = 0.53 (0.62 in SHV-1) and SimpleScore = 0.05 (0.09 in SHV-1).

Fig. 4:
electrostatic potential of (A) other Klebsiella pneumoniae beta-lactamase (OKP-B-6), (B) sulfhydryl reagent variable, point mutation variants of SHV-1 and (C) K. pneumoniae carbapenemase (KPC-2) focusing on the sulbactam binding site. The electrostatic potential was calculated in Pymol/Adaptive Poisson-Boltzmann Solver (APBS) at pH 7, ranging from blue (+3 kT/e, positively charged) to red (-3 kT/e, negatively charged). The structure of the pre-acylation state of sulbactam protein data bank (PDB) code 4FH2, aligned to the OKP-B-6 model, is present to highlight the binding site location, but it was not present during the potential calculation.

Fig. 5:
sulbactam binding cavity predicted with DogSiteScorer for other Klebsiella pneumoniae beta-lactamase (OKP-B-6) and sulfhydryl reagent variable, point mutation variants of SHV-1. DrugScore and SimpleScore range from 0 to 1, so the higher the score, the more druggable the pocket is estimated to be. The binding site location is highlighted by the structure of the pre-acylation state of sulbactam protein data bank (PDB) code code 4FH2,⁵² aligned to the OKP-B-6 model.

Docking of the inhibitors against OKP-B-6 - Sulbactam and tazobactam have similar mechanisms of action, inhibiting beta-lactamases through a serine-acylation mechanism.^15,57 Avibactam, an inhibitor that has the beta-lactam core replaced by a 5-membered ring with a carbonyl bond adjacent to ring nitrogen, has been proposed to inhibit the enzyme through a distinct mechanism of action, despite interacting with the same key residues within the binding site.^12,16

Our docking protocol against OKP-B-6 successfully predicted the experimental binding mode of sulbactam observed in the pre-acylation state of SHV-1-sulbactam (PDB code 4FH2) with an RMSD of 0.97 Å (E_total = -21.88 kcal/mol), which is also similar to the binding mode proposed for the recognition of beta-lactam substrates (PDB code 1KVL)^57,61 (Fig. 6 and Table). After the acylation step, the ligand undergoes a significantly conformational change,⁶¹ as observed in many experimentally solved structures of the acyl-enzyme intermediates for different inhibitors. Since sulbactam and tazobactam have the same mechanism of action, it is expected that their binding modes are similar before the covalent bond formation. Despite no experimental structure of the pre-acylation state of tazobactam is available, its top-ranked predicted pose (E_total = -31.03 kcal/mol) perfectly fits the experimental binding mode of sulbactam, being characterised by the key interactions with the binding site, namely: (i) hydrogen bonds from the carbonyl oxygen from the beta-lactam ring with the oxyanion hole formed by Ala212 (237) and the catalytic Ser45 (70), (ii) the stabilisation of the negatively charged carboxyl moiety by strong hydrogen bonds with Thr210 (235), Wat2 and the catalytic Ser105 (130), and through salt bridges with Arg218 (244) and Lys209 (234) (Fig. 6). The top-ranked docked pose of avibactam (E_total = -22.91 kcal/mol, Fig. 6 and Table) is characterised by the positioning of the sulfate moiety at the carboxyl binding pocket, interacting with Ser45 (70), Ser105 (130), Thr210 (235) and Arg218 (244) through hydrogen bonds, whereas the carbonyl oxygen from the terminal amide group is oriented to the side chains of Asn107 (132) and Asn145 (170). However, in this binding mode, the carbonyl carbon atom of the 5-membered ring, expected to suffer the nucleophilic attack, is exposed to the solvent, opposed to the Ser45 (70) (4.6 Å). Exploring the different binding modes predicted by DockThor, we found that the top-2 best-docked conformation of avibactam (E_total = -19.10 kcal/mol) is more likely to be the pre-acylation state, with the carbonyl oxygen atom located at the oxyanion hole interacting with the side chain of Ser45 (70) at 2.6 Å and the main chain of Ala212 (237), as observed in the experimental structures of avibactam complexed with different beta-lactamases.^{16,46,62,63,64} The docking results obtained in this work for the evaluated SHV-1 inhibitors against the OKP-B-6 binding site, associated with the highly conserved residues when compared to the SHV-1 target, suggest that OKP-B-6 is susceptible to the same drugs.

Fig. 6:
docking results for the inhibitors sulbactam, tazobactam, and avibactam against the other Klebsiella pneumoniae beta-lactamase (OKP-B-6) model. The pre-acylation state of sulbactam observed experimentally protein data bank (PDB) code 4FH2 is represented as cyan sticks. Conserved waters are represented as red spheres. Hydrogen bonds are represented as yellow dashes, and the distance between the carbonyl oxygen and Ser45 is highlighted. Salmon dashed line represents the distance between the carbonyl carbon atom of avibactam (top-1 pose) and Ser45 hydroxyl oxygen.

TABLE
Docking results of the inhibitors sulbactam, tazobactam and avibactam against the other Klebsiella pneumoniae beta-lactamase (OKP-B-6) model. The two top-ranked binding modes of avibactam were selected according to the total energy, which is the default scoring function for pose prediction in the DockThor program

Compound	Total energya	Elec. energya	vdW energya
Sulbactam	-21.88 (0.97b )	-39.70	1.88
Tazobactam	-31.03	-39.00	-0.40
Avibactam (top-1)	-22.91	-38.16	-4.32
Avibactam (top-2)	-19.10	-26.53	-11.39

• a: given in kcal/mol; b: root-mean-square deviation from the pre-acylation state of sulbactam protein data bank (PDB) code 4FH2, given in Å.

In conclusion - Our results indicate that OKP-B-6 should be an ESBL like SHV-1, with both beta-lactamases having the amino acid residues within the binding cavities highly conserved. We identified that the predicted binding modes with the docking of either sulbactam, tazobactam, or avibactam with OKP-B-6 are similar to those observed experimentally with SHV-1, which is expected since these inhibitors have broad activity against equivalent ESBLs. The docking results obtained in this work validated the constructed OKP-B-6 model and provided useful insights for its inhibition. These results also have important implications for clinical microbiologists in the control of pathogens that carry one of these types of class A beta-lactamases, SHV, or OKP-B. Particularly, the formulation prescribed with one of the inhibitors tested here would also be effective in controlling K. quasipneumoniae carrying OKP-B.

ACKNOWLEDGEMENTS

To the support of the Brazilian Sistema Nacional de Processamento de Alto Desempenho (SINAPAD) and the availability of the computational resources provided by the Supercomputer SDumont (LNCC/MCTI). We also acknowledge the Instituto Nacional de Ciência e Tecnologia em Fármacos e Medicamentos (INCT-INOFAR).

REFERENCES

• ¹ Docquier JD, Mangani S. An update on ß-lactamase inhibitor discovery and development. Drug Resist Updat. 2018; 36: 13-29.
• ² Ambler RP. The structure of beta-lactamases. Philos Trans R Soc Lond B Biol Sci. 1980; 289(1036): 321-31.
• ³ Ghuysen JM. Serine beta-lactamases and penicillin-binding proteins. Annu Rev Microbiol. 1991; 45: 37-67.
• ⁴ Frère JM, Page MG. Penicillin-binding proteins: evergreen drug targets. Curr Opin Pharmacol. 2014; 18: 112-9.
• ⁵ Lamotte-Brasseur J, Knox J, Kelly JA, Charlier P, Fonzé E, Dideberg O, et al. The structures and catalytic mechanisms of active-site serine beta-lactamases. Biotechnol Genet Eng Rev. 1994; 12: 189-230.
• ⁶ Bush K, Jacoby GA. Updated functional classification of beta-lactamases. Antimicrob Agents Chemother. 2010; 54(3): 969-76.
• ⁷ Ke W, Bethel CR, Thomson JM, Bonomo RA, van den Akker F. Crystal structure of KPC-2: insights into carbapenemase activity in class A beta-lactamases. Biochemistry. 2007; 46(19): 5732-40.
• ⁸ Nordmann P, Cuzon G, Naas T. The real threat of Klebsiella pneumoniae carbapenemase-producing bacteria. Lancet Infect Dis. 2009; 9(4): 228-36.
• ⁹ Long SW, Linson SE, Saavedra MO, Cantu C, Davis JJ, Brettin T, et al. Whole-genome sequencing of human clinical Klebsiella pneumoniae isolates reveals misidentification and misunderstandings of Klebsiella pneumoniae, Klebsiella variicola, and Klebsiella quasipneumoniae. mSphere. 2017; 2(4): e00290-17.
• ¹⁰ Tehrani KHME, Martin NI. ß-lactam/ß-lactamase inhibitor combinations: an update. MedChemComm. 2018; 9(9): 1439-56.
• ¹¹ Pereira R, Rabelo VWH, Sibajev A, Abreu PA, Castro HC. Class A ß-lactamases and inhibitors: In silico analysis of the binding mode and the relationship with resistance. J Biotechnol. 2018; 279: 37-46.
• ¹² van den Akker F, Bonomo RA. Exploring additional dimensions of complexity in inhibitor design for serine ß-lactamases: mechanistic and intra- and inter-molecular chemistry approaches. Front Microbiol. 2018; 9: 622.
• ¹³ Padayatti PS, Helfand MS, Totir MA, Carey MP, Hujer AM, Carey PR, et al. Tazobactam forms a stoichiometric trans-enamine intermediate in the E166A variant of SHV-1 beta-lactamase: 1.63 A crystal structure. Biochemistry. 2004; 43(4): 843-8.
• ¹⁴ Sun T, Bethel CR, Bonomo RA, Knox JR. Inhibitor-resistant class A beta-lactamases: consequences of the Ser130-to-Gly mutation seen in Apo and tazobactam structures of the SHV-1 variant. Biochemistry. 2004; 43(44): 14111-7.
• ¹⁵ Totir MA, Padayatti PS, Helfand MS, Carey MP, Bonomo RA, Carey PR, et al. Effect of the inhibitor-resistant M69V substitution on the structures and populations of trans-enamine ß-lactamase intermediates. Biochemistry. 2006; 45(39): 11895-904.
• ¹⁶ Krishnan NP, Nguyen NQ, Papp-Wallace KM, Bonomo RA, van den Akker F. Inhibition of Klebsiella ß-lactamases (SHV-1 and KPC-2) by avibactam: a dtructural study. PLoS One. 2015; 10(9): e0136813.
• ¹⁷ Ehmann DE, Jahic H, Ross PL, Gu RF, Hu J, Kern G, et al. Avibactam is a covalent, reversible, non- -lactam -lactamase inhibitor. Proc Natl Acad Sci. 2012; 109(29): 11663-8.
• ¹⁸ Tsivkovski R, Lomovskaya O. Potency of vaborbactam is less affected than that of avibactam in strains producing KPC-2 mutations that confer resistance to ceftazidime-avibactam. Antimicrob Agents Chemother. 2020; 64(4): e01936-19.
• ¹⁹ Bush K, Bradford PA. ß-lactams and ß-lactamase inhibitors: an overview. Cold Spring Harb Perspect Med. 2016; 6(8): a025247.
• ²⁰ Ehmann DE, Jahic H, Ross PL, Gu RF, Hu J, Durand-Réville TF, et al. Kinetics of avibactam inhibition against class A, C, and D ß-lactamases. J Biol Chem. 2013; 288(39): 27960-71.
• ²¹ Ramos PIP, Picão RC, de Almeida LGP, Lima NCB, Girardello R, Vivan ACP, et al. Comparative analysis of the complete genome of KPC-2-producing Klebsiella pneumoniae Kp13 reveals remarkable genome plasticity and a wide repertoire of virulence and resistance mechanisms. BMC Genomics. 2014; 15(1): 54.
• ²² Ambrose SJ, Harmer CJ, Hall RM. Evolution and typing of IncC plasmids contributing to antibiotic resistance in Gram-negative bacteria. Plasmid. 2018; 99: 40-55.
• ²³ van Dorp L, Wang Q, Shaw LP, Acman M, Brynildsrud OB, Eldholm V, et al. Rapid phenotypic evolution in multidrug-resistant Klebsiella pneumoniae hospital outbreak strains. Microb Genom. 2019; 5(4): e000263.
• ²⁴ Haeggman S, Löfdahl S, Paauw A, Verhoef J, Brisse S. Diversity and evolution of the class A chromosomal beta-lactamase gene in Klebsiella pneumoniae. Antimicrob Agents Chemother. 2004; 48(7): 2400-8.
• ²⁵ Siebor E, Péchinot A, Duez JM, Neuwirth C. One new LEN enzyme and two new OKP enzymes in Klebsiella pneumoniae clinical isolates and proposed nomenclature for chromosomal ß-lactamases of this species. Antimicrob Agents Chemother. 2005; 49(7): 3097-8.
• ²⁶ Brisse S, Passet V, Grimont PAD. Description of Klebsiella quasipneumoniae sp. nov., isolated from human infections, with two subspecies, Klebsiella quasipneumoniae subsp. quasipneumoniae subsp. nov. and Klebsiella quasipneumoniae subsp. similipneumoniae subsp. nov., and demonstration that Klebsiella singaporensis is a junior heterotypic synonym of Klebsiella variicola. Int J Syst Evol Microbiol. 2014; 64(Pt 9): 3146-52.
• ²⁷ Fevre C, Passet V, Weill FX, Grimont PAD, Brisse S. Variants of the Klebsiella pneumoniae OKP chromosomal beta-lactamase are divided into two main groups, OKP-A and OKP-B. Antimicrob Agents Chemother. 2005; 49(12): 5149-52.
• ²⁸ Martínez-Romero E, Rodríguez-Medina N, Beltrán-Rojel M, Toribio-Jiménez J, Garza-Ramos U. Klebsiella variicola and Klebsiella quasipneumoniae with capacity to adapt to clinical and plant settings. Salud Publica Mex. 2018; 60(1): 29-40.
• ²⁹ Arabaghian H, Salloum T, Alousi S, Panossian B, Araj GF, Tokajian S. Molecular characterization of carbapenem resistant Klebsiella pneumoniae and Klebsiella quasipneumoniae isolated from Lebanon. Sci Rep. 2019; 9(1): 531.
• ³⁰ Mathers AJ, Crook D, Vaughan A, Barry KE, Vegesana K, Stoesser N, et al. Klebsiella quasipneumoniae provides a window into carbapenemase gene transfer, plasmid rearrangements, and patient interactions with the hospital environment. Antimicrob Agents Chemother. 2019; 63(6): e02513-18.
• ³¹ Octavia S, Kalisvar M, Venkatachalam I, Ng OT, Xu W, Sridatta PSR, et al. Klebsiella pneumoniae and Klebsiella quasipneumoniae define the population structure of blaKPC-2Klebsiella: a 5 year retrospective genomic study in Singapore. J Antimicrob Chemother. 2019; 74(11): 3205-10.
• ³² Nicolás MF, Ramos PIP, de Carvalho FM, Camargo DRA, Alves CFM, de Morais GL, et al. Comparative genomic analysis of a clinical isolate of Klebsiella quasipneumoniae subsp. similipneumoniae, a KPC-2 and OKP-B-6 beta-lactamases producer harboring two drug-resistance plasmids from southeast Brazil. Front Microbiol. 2018; 9: 220.
• ³³ Paetzel M. Bacterial signal peptidases. Subcell Biochem. 2019; 92: 187-219.
• ³⁴ Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018; 46(W1): W296-303.
• ³⁵ Williams CJ, Headd JJ, Moriarty NW, Prisant MG, Videau LL, Deis LN, et al. MolProbity: more and better reference data for improved all-atom structure validation. Protein Sci. 2018; 27(1): 293-315.
• ³⁶ Chen VB, Arendall WB 3rd, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, et al. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr D Biol Crystallogr. 2010; 66(Pt 1): 12-21.
• ³⁷ Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr. 1993; 26(2): 283-91.
• ³⁸ Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996; 8(4): 477-86.
• ³⁹ Colovos C, Yeates TO. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993; 2(9): 1511-9.
• ⁴⁰ Kuzin AP, Nukaga M, Nukaga Y, Hujer A, Bonomo RA, Knox JR. Inhibition of the SHV-1 ß-lactamase by sulfones: crystallographic observation of two reaction intermediates with tazobactam. Biochemistry. 2001; 40(6): 1861-6.
• ⁴¹ Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004; 32: W665-7.
• ⁴² The PyMOL Molecular Graphics System. Version 2.5.0. Schrödinger, LLC. Available from: https://pymol.org/
  » https://pymol.org/
• ⁴³ Volkamer A, Kuhn D, Rippmann F, Rarey M. DoGSiteScorer: a web server for automatic binding site prediction, analysis and druggability assessment. Bioinformatics. 2012; 28(15): 2074-5.
• ⁴⁴ Olsson MHM, Søndergaard CR, Rostkowski M, Jensen JH. PROPKA3: consistent treatment of internal and surface residues in empirical pKa predictions. J Chem Theory Comput. 2011; 7(2): 525-37.
• ⁴⁵ Pemberton OA, Noor RE, Kumar MVV, Sanishvili R, Kemp MT, Kearns FL, et al. Mechanism of proton transfer in class A ß-lactamase catalysis and inhibition by avibactam. Proc Natl Acad Sci. 2020; 117(11): 5818-25.
• ⁴⁶ Lahiri SD, Mangani S, Durand-Reville T, Benvenuti M, De Luca F, Sanyal G, et al. Structural insight into potent broad-spectrum inhibition with reversible recyclization mechanism: avibactam in complex with CTX-M-15 and Pseudomonas aeruginosa AmpC ß-lactamases. Antimicrob Agents Chemother. 2013; 57(6): 2496-505.
• ⁴⁷ Sgrignani J, Grazioso G, De Amici M, Colombo G. Inactivation of TEM-1 by avibactam (NXL-104): insights from quantum mechanics/molecular mechanics metadynamics simulations. Biochemistry. 2014; 53(31): 5174-85.
• ⁴⁸ Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, et al. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021; 49(D1): D1388-95.
• ⁴⁹ Shelley JC, Cholleti A, Frye LL, Greenwood JR, Timlin MR, Uchimaya M. Epik: a software program for pKa prediction and protonation state generation for drug-like molecules. J Comput Aided Mol Des. 2007; 21(12): 681-91.
• ⁵⁰ Che T, Rodkey EA, Bethel CR, Shanmugam S, Ding Z, Pusztai-Carey M, et al. Detecting a quasi-stable imine species on the reaction pathway of SHV-1 ß-lactamase and 6ß-(hydroxymethyl)penicillanic acid sulfone. Biochemistry. 2015; 54(3): 734-43.
• ⁵¹ Guedes IA, Pereira FSS, Dardenne LE. Empirical scoring functions for structure-based virtual screening: applications, critical aspects, and challenges. Front Pharmacol. 2018; 9: 1-18.
• ⁵² Guedes IA, de Magalhães CS, Dardenne LE. Receptor-ligand molecular docking. Biophys Rev. 2014; 6(1): 75-87.
• ⁵³ Capriles PVSZ, Baptista LPR, Guedes IA, Guimarães ACR, Custódio FL, Alves-Ferreira M, et al. Structural modeling and docking studies of ribose 5-phosphate isomerase from Leishmania major and Homo sapiens: a comparative analysis for Leishmaniasis treatment. J Mol Graph Model. 2015; 55: 134-47.
• ⁵⁴ Guedes IA, Costa LSC, dos Santos KB, Karl ALM, Rocha GK, Teixeira IM, et al. Drug design and repurposing with DockThor-VS web server focusing on SARS-CoV-2 therapeutic targets and their non-synonym variants. Sci Rep. 2021; 11(1): 5543.
• ⁵⁵ dos Santos KB, Guedes IA, Karl ALM, Dardenne L. Highly flexible ligand docking: benchmarking of the DockThor Program on the LEADS-PEP Protein-peptide Dataset. J Chem Inf Model. 2020; 60(2): 667-83.
• ⁵⁶ de Magalhães CS, Almeida DM, Barbosa HJC, Dardenne LE. A dynamic niching genetic algorithm strategy for docking highly flexible ligands. Inf Sci. 2014; 289: 206-24.
• ⁵⁷ Rodkey EA, Drawz SM, Sampson JM, Bethel CR, Bonomo RA, van den Akker F. Crystal structure of a preacylation complex of the ß-lactamase inhibitor sulbactam bound to a sulfenamide bond-containing thiol-ß-lactamase. J Am Chem Soc. 2012; 134(40): 16798-804.
• ⁵⁸ Halgren TA. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J Comput Chem. 1996; 17(5-6): 490-519.
• ⁵⁹ Guedes IA, Barreto AMS, Marinho D, Krempser E, Kuenemann MA, Sperandio O, et al. New machine learning and physics-based scoring functions for drug discovery. Sci Rep. 2021; 11(1): 3198.
• ⁶⁰ Winkler ML, Papp-Wallace KM, Taracila MA, Bonomo RA. Avibactam and inhibitor-resistant SHV ß-lactamases. Antimicrob Agents Chemother. 2015; 59(7): 3700-9.
• ⁶¹ Beadle BM, Trehan I, Focia PJ, Shoichet BK. Structural milestones in the reaction pathway of an amide hydrolase: substrate, acyl, and product complexes of cephalothin with AmpC beta-lactamase. Structure. 2002; 10(3): 413-24.
• ⁶² Xu H, Hazra S, Blanchard JS. NXL104 irreversibly inhibits the ß-lactamase from Mycobacterium tuberculosis. Biochemistry. 2012; 51(22): 4551-7.
• ⁶³ King DT, King AM, Lal SM, Wright GD, Strynadka NCJ. Molecular mechanism of avibactam-mediated ß-lactamase inhibition. ACS Infect Dis. 2015; 1(4): 175-84.
• ⁶⁴ Stewart NK, Toth M, Stasyuk A, Vakulenko SB, Smith CA. In crystallo time-resolved interaction of the clostridioides difficile CDD-1 enzyme with avibactam provides new insights into the catalytic mechanism of class D ß-lactamases. ACS Infect Dis. 2021; 7(6): 1765-76.

• Financial support: CNPq (Grant no. 306894/2019-0, 307634/2019-1, 307145/2021-2), FAPERJ (Grant no. E-26/010.001229/2015, E-26/202.826/2018, E-26/210.935/2019, E-26/211.357/2021, E-26/200.608/2022 and E-26/210.372/2022), CAPES (Grant no. 88887.368759/2019-00), PCI-LNCC/CNPq (Grant no. 302201/2021-1).

Publication Dates

• Publication in this collection
  23 Sept 2022
• Date of issue
  2022

History

• Received
  30 Apr 2022
• Accepted
  22 Aug 2022