Roles of residues Cys69, Asn104, Phe160, Gly232, Ser237, and Asp240 in extended-spectrum beta-lactamase Toho-1

Mechanistic Insights into Enzyme Catalysis from Explaining Machine-Learned Quantum Mechanical and Molecular Mechanical Minimum Energy Pathways

With the increasing popularity of machine learning (ML) applications, the demand for explainable artificial intelligence techniques to explain ML models developed for computational chemistry has also emerged. In this study, we present the development of the Boltzmann-weighted cumulative integrated gradients (BCIG) approach for effective explanation of mechanistic insights into ML models trained on high-level quantum mechanical and molecular mechanical (QM/MM) minimum energy pathways. Using the acylation reactions of the Toho-1 β-lactamase and two antibiotics (ampicillin and cefalexin) as the model systems, we show that the BCIG approach could quantitatively attribute the energetic contribution in one system and the relative reactivity of individual steps across different systems to specific chemical processes such as the bond making/breaking and proton transfers. The proposed BCIG contribution attribution method quantifies chemistry-interpretable insights in terms of contributions from each elementary chemical process, which is in agreement with the validating QM/MM calculations and our intuitive mechanistic understandings of the model reactions.

QM/MM modeling of class A β-lactamases reveals distinct acylation pathways for ampicillin and cefalexin

Efficient mechanism-based design of antibiotics that are not susceptible to β-lactamases is hindered by the lack of comprehensive knowledge on the energetic landscapes for the hydrolysis of various β-lactams. Herein, we adopted efficient quantum mechanics/molecular mechanics simulations to explore the acylation reaction catalyzed by CTX-M-44 (Toho-1) β-lactamase. We show that the catalytic pathways for β-lactam hydrolysis are correlated to substrate scaffolds: using Glu166 as the only general base for acylation is viable for ampicillin but prohibitive for cefalexin. The present computational workflow provides quantitative insights to facilitate the optimization of future β-lactam antibiotics.

Non-catalytic-Region Mutations Conferring Transition of Class A β-Lactamases Into ESBLs

Despite class A ESBLs carrying substitutions outside catalytic regions, such as Cys69Tyr or Asn136Asp, have emerged as new clinical threats, the molecular mechanisms underlying their acquired antibiotics-hydrolytic activity remains unclear. We discovered that this non-catalytic-region (NCR) mutations induce significant dislocation of β3-β4 strands, conformational changes in critical residues associated with ligand binding to the lid domain, dynamic fluctuation of Ω-loop and β3-β4 elements. Such structural changes increase catalytic regions’ flexibility, enlarge active site, and thereby accommodate third-generation cephalosporin antibiotics, ceftazidime (CAZ). Notably, the electrostatic property around the oxyanion hole of Cys69Tyr ESBL is significantly changed, resulting in possible additional stabilization of the acyl-enzyme intermediate. Interestingly, the NCR mutations are as effective for antibiotic resistance by altering the structure and dynamics in regions mediating substrate recognition and binding as single amino-acid substitutions in the catalytic region of the canonical ESBLs. We believe that our findings are crucial in developing successful therapeutic strategies against diverse class A ESBLs, including the new NCR-ESBLs.

Dynamical Behavior of β-Lactamases and Penicillin- Binding Proteins in Different Functional States and Its Potential Role in Evolution

β-Lactamases are enzymes produced by bacteria to hydrolyze β-lactam-based antibiotics, and pose serious threat to public health through related antibiotic resistance. Class A β-lactamases are structurally and functionally related to penicillin-binding proteins (PBPs). Despite the extensive studies of the structures, catalytic mechanisms and dynamics of both β-lactamases and PBPs, the potentially different dynamical behaviors of these proteins in different functional states still remain elusive in general. In this study, four evolutionarily related proteins, including TEM-1 and TOHO-1 as class A β-lactamases, PBP-A and DD-transpeptidase as two PBPs, are subjected to molecular dynamics simulations and various analyses to characterize their dynamical behaviors in different functional states. Penicillin G and its ring opening product serve as common ligands for these four proteins of interest. The dynamic analyses of overall structures, the active sites with penicillin G, and three catalytically important residues commonly shared by all four proteins reveal unexpected cross similarities between Class A β-lactamases and PBPs. These findings shed light on both the hidden relations among dynamical behaviors of these proteins and the functional and evolutionary relations among class A β-lactamases and PBPs.

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

β-lactams are one of the most common antimicrobials used to treat bacterial infections. However, bacterial resistance has compromised their efficacy, mainly due to the β-lactamase enzyme production. To overcome this resistance, β-lactamase inhibitors can be used in association with these antimicrobials. Herein, we analyzed the structural characteristics of β-lactamases and their interactions with classical inhibitors, such as clavulanic acid (CA), sulbactam (SB) and tazobactam (TZ) to gain insights into resistance. The homology models of five class A β-lactamases, namely CARB-3, IMI-1, SFO-1, SHV-5 and TEM-10, were constructed and validated and revealed an overall 3D structural conservation, but with significant differences in the electrostatic potential maps, especially at important regions in the catalytic site. Molecular dockings of CA, SB and TZ with these enzymes revealed a covalent bond with the S70 in all complexes, except Carb-3 which is in agreement with experimental data reported so far. This is likely related to the less voluminous active site of Carb-3 model. Although few specific contacts were observed in the β-lactamase-inhibitor complexes, all compounds interacted with the residues in positions 73, 130, 132, 236 and 237. Therefore, this study provides new perspectives for the design of innovative compounds with broad-spectrum inhibitory profiles against β-lactamases.

Antibacterial Spectrum of a Tetrazole-Based Reversible Inhibitor of Serine β-Lactamases

CTX-M is the most prevalent family of extended-spectrum β-lactamases. We recently developed a tetrazole-derived noncovalent inhibitor of CTX-M-9. Here, we present the biochemical and microbiological activity of this inhibitor across a representative panel of serine β-lactamases and Gram-negative bacteria. The compound displayed significant activity against all major subgroups of CTX-M, including CTX-M-15, while it exhibited some low-level inhibition of other serine β-lactamases. Complex crystal structures with the CTX-M-14 S237A mutant and CTX-M-27 illustrate the binding contribution of specific active-site residues on the β3 strand. In vitro pharmacokinetic studies revealed drug-like properties and positive prospects for further optimization. These studies suggest that tetrazole-based compounds can provide novel chemotypes for future serine β-lactamase inhibitor discovery.

Defining Substrate Specificity in the CTX-M Family: the Role of Asp240 in Ceftazidime Hydrolysis

The natural diversification of CTX-M β-lactamases led to the emergence of Asp240Gly variants in the clinic that confer reduced susceptibility to ceftazidime (CAZ). In this study, we compared the impact of this substitution on CAZ and ceftazidime-avibactam (CZA) MICs against isogenic Escherichia coli strains with different porin deficiencies. Our results show a noticeable increase in CAZ resistance in clones expressing Asp240Gly-harboring CTX-M when combined with OmpF porin deficiency. Kinetic analysis revealed that the k_cat/K_m for CAZ was 5- to 15-fold higher for all Asp240Gly variants but remained 200- to 725-fold lower than that for cefotaxime (CTX). In vitro selection of CAZ-resistant clones yielded nonsusceptible CTX-M producers (MIC of >16 μg/ml) only after overnight incubation; the addition of avibactam (AVI) decreased MICs to a susceptible range against these variants. In contrast, the use of CZA as a selective agent did not yield resistant clones. AVI inactivated both CTX-M-12 and CTX-M-96, with an apparent inhibition constant comparable to that of SHV-2 and 1,000-fold greater than that of PER-2 and CMY-2, and k₂/K for CTX-M-12 was 24- and 35-fold higher than that for CTX-M-96 and CTX-M-15, respectively. Molecular modeling suggests that AVI interacts similarly with CTX-M-96 and CTX-M-15. We conclude that the impact of Asp240Gly in resistance may arise when other mechanisms are also present (i.e., OmpF deficiency). Additionally, CAZ selection could favor the emergence of CAZ-resistant subpopulations. These results define the role of Asp240 and the impact of the -Gly substitution and allow us to hypothesize that the use of CZA is an effective preventive strategy to delay the development of resistance in this family of extended-spectrum β-lactamases.

A Structure-Based Classification of Class A β-Lactamases, a Broadly Diverse Family of Enzymes

For medical biologists, sequencing has become a commonplace technique to support diagnosis. Rapid changes in this field have led to the generation of large amounts of data, which are not always correctly listed in databases. This is particularly true for data concerning class A β-lactamases, a group of key antibiotic resistance enzymes produced by bacteria. Many genomes have been reported to contain putative β-lactamase genes, which can be compared with representative types. We analyzed several hundred amino acid sequences of class A β-lactamase enzymes for phylogenic relationships, the presence of specific residues, and cluster patterns. A clear distinction was first made between dd-peptidases and class A enzymes based on a small number of residues (S70, K73, P107, 130SDN132, G144, E166, 234K/R, 235T/S, and 236G [Ambler numbering]). Other residues clearly separated two main branches, which we named subclasses A1 and A2. Various clusters were identified on the major branch (subclass A1) on the basis of signature residues associated with catalytic properties (e.g., limited-spectrum β-lactamases, extended-spectrum β-lactamases, and carbapenemases). For subclass A2 enzymes (e.g., CfxA, CIA-1, CME-1, PER-1, and VEB-1), 43 conserved residues were characterized, and several significant insertions were detected. This diversity in the amino acid sequences of β-lactamases must be taken into account to ensure that new enzymes are accurately identified. However, with the exception of PER types, this diversity is poorly represented in existing X-ray crystallographic data.

Kinetic and crystallographic studies of extended-spectrum GES-11, GES-12, and GES-14 β-lactamases

GES-1 is a class A extended-spectrum β-lactamase conferring resistance to penicillins, narrow- and expanded-spectrum cephalosporins, and ceftazidime. However, GES-1 poorly hydrolyzes aztreonam and cephamycins and exhibits very low k(cat) values for carbapenems. Twenty-two GES variants have been discovered thus far, differing from each other by 1 to 3 amino acid substitutions that affect substrate specificity. GES-11 possesses a Gly243Ala substitution which seems to confer to this variant an increased activity against aztreonam and ceftazidime. GES-12 differs from GES-11 by a single Thr237Ala substitution, while GES-14 differs from GES-11 by the Gly170Ser mutation, which is known to confer increased carbapenemase activity. GES-11 and GES-12 were kinetically characterized and compared to GES-1 and GES-14. Purified GES-11 and GES-12 showed strong activities against most tested β-lactams, with the exception of temocillin, cefoxitin, and carbapenems. Both variants showed a significantly increased rate of hydrolysis of cefotaxime, ceftazidime, and aztreonam. On the other hand, GES-11 and GES-12 (and GES-14) variants all containing Ala243 exhibited increased susceptibility to classical inhibitors. The crystallographic structures of the GES-11 and GES-14 β-lactamases were solved. The overall structures of GES-11 and GES-14 are similar to that of GES-1. The Gly243Ala substitution caused only subtle local rearrangements, notably in the typical carbapenemase disulfide bond. The active sites of GES-14 and GES-11 are very similar, with the Gly170Ser substitution leading only to the formation of additional hydrogen bonds of the Ser residue with hydrolytic water and the Glu166 residue.

Epidemiological expansion, structural studies, and clinical challenges of new β-lactamases from gram-negative bacteria

β-Lactamase evolution presents to the infectious disease community a major challenge in the treatment of infections caused by multidrug-resistant gram-negative bacteria. Because over 1,000 of these naturally occurring β-lactamases exist, attempts to correlate structure and function have become daunting. Although new enzymes in the extended-spectrum β-lactamase (ESBL) families are frequently identified, the older CTX-M-14 and CTX-M-15 enzymes have become the most prevalent ESBLs in global surveillance. Carbapenemases with either serine-based or zinc-facilitated hydrolysis mechanisms are posing some of the most critical problems. Most geographical regions now report KPC serine carbapenemases and the metallo-β-lactamases VIM, IMP, and NDM-1, even though NDM-1 was only recently identified. The rapid emergence of these newer enzymes, with multiple β-lactamases appearing in a single organism, makes the design of new β-lactamase inactivators or β-lactamase-stable β-lactams all the more difficult. Combination therapy will likely be required to counteract the continuing evolution of these insidious enzymes in multidrug-resistant pathogens.

Distant and new mutations in CTX-M-1 beta-lactamase affect cefotaxime hydrolysis

The CTX-M β-lactamases are an increasingly prevalent group of extended-spectrum β-lactamases (ESBL). Point mutations in CTX-M β-lactamases are considered critical for enhanced hydrolysis of cefotaxime. In order to clarify the structural determinants of the activity against cefotaxime in CTX-M β-lactamases, screening for random mutations was carried out to search for decreased activity against cefotaxime, with the CTX-M-1 gene as a model. Thirteen single mutants with a considerable reduction in cefotaxime MICs were selected for biochemical and stability studies. The 13 mutated genes of the CTX-M-1 β-lactamase were expressed, and the proteins were purified for kinetic studies against cephalothin and cefotaxime (as the main antibiotics). Some of the positions, such as Val103Asp, Asn104Asp, Asn106Lys, and Pro107Ser, are located in the (103)VNYN(106) loop, which had been described as important in cefotaxime hydrolysis, although this has not been experimentally confirmed. There are four mutations located close to catalytic residues-Thr71Ile, Met135Ile, Arg164His, and Asn244Asp-that may affect the positioning of these residues. We show here that some distant mutations, such as Ala219Val, are critical for cefotaxime hydrolysis and highlight the role of this loop at the top of the active site. Other distant substitutions, such as Val80Ala, Arg191, Ala247Ser, and Val260Leu, are in hydrophobic cores and may affect the dynamics and flexibility of the enzyme. We describe here, in conclusion, new residues involved in cefotaxime hydrolysis in CTX-M β-lactamases, five of which are in positions distant from the catalytic center.