Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography

Drug Discovery in the Field of β-Lactams: An Academic Perspective

β-Lactams are the most widely prescribed class of antibiotics that inhibit penicillin-binding proteins (PBPs), particularly transpeptidases that function in peptidoglycan synthesis. A major mechanism of antibiotic resistance is the production of β-lactamase enzymes, which are capable of hydrolyzing β-lactam antibiotics. There have been many efforts to counter increasing bacterial resistance against β-lactams. These studies have mainly focused on three areas: discovering novel inhibitors against β-lactamases, developing new β-lactams less susceptible to existing resistance mechanisms, and identifying non-β-lactam inhibitors against cell wall transpeptidases. Drug discovery in the β-lactam field has afforded a range of research opportunities for academia. In this review, we summarize the recent new findings on both β-lactamases and cell wall transpeptidases because these two groups of enzymes are evolutionarily and functionally connected. Many efforts to develop new β-lactams have aimed to inhibit both transpeptidases and β-lactamases, while several promising novel β-lactamase inhibitors have shown the potential to be further developed into transpeptidase inhibitors. In addition, the drug discovery progress against each group of enzymes is presented in three aspects: understanding the targets, screening methodology, and new inhibitor chemotypes. This is to offer insights into not only the advancement in this field but also the challenges, opportunities, and resources for future research. In particular, cyclic boronate compounds are now capable of inhibiting all classes of β-lactamases, while the diazabicyclooctane (DBO) series of small molecules has led to not only new β-lactamase inhibitors but potentially a new class of antibiotics by directly targeting PBPs. With the cautiously optimistic successes of a number of new β-lactamase inhibitor chemotypes and many questions remaining to be answered about the structure and function of cell wall transpeptidases, non-β-lactam transpeptidase inhibitors may usher in the next exciting phase of drug discovery in this field.

Structure of a ribonucleotide reductase R2 protein radical

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Correction to "Bifurcated Hydrogen Bonds in a Peptide Crystal Unveiled by X-ray Diffraction and Polarized Raman Spectroscopy"

[This corrects the article DOI: 10.1021/acs.cgd.3c00302.].

Biochemical and Structural Characterization of CRH-1, a Carbapenemase from Chromobacterium haemolyticum Related to KPC β-Lactamases

KPC-2 is one of the most relevant serine-carbapenemases among the carbapenem-resistant Enterobacterales. We previously isolated from the environmental species Chromobacterium haemolyticum a class A CRH-1 β-lactamase displaying 69% amino acid sequence identity with KPC-2. The objective of this study was to analyze the kinetic behavior and crystallographic structure of this β-lactamase. Our results showed that CRH-1 can hydrolyze penicillins, cephalosporins (except ceftazidime), and carbapenems with similar efficacy compared to KPC-2. Inhibition kinetics showed that CRH-1 is not well inhibited by clavulanic acid, in contrast to efficient inhibition by avibactam (AVI). The high-resolution crystal of the apoenzyme showed that CRH-1 has a similar folding compared to other class A β-lactamases. The CRH-1/AVI complex showed that AVI adopts a chair conformation, stabilized by hydrogen bonds to Ser70, Ser237, Asn132, and Thr235. Our findings highlight the biochemical and structural similarities of CRH-1 and KPC-2 and the potential clinical impact of this carbapenemase in the event of recruitment by pathogenic bacterial species.

Mapping the determinants of catalysis and substrate specificity of the antibiotic resistance enzyme CTX-M β-lactamase

CTX-M β-lactamases are prevalent antibiotic resistance enzymes and are notable for their ability to rapidly hydrolyze the extended-spectrum cephalosporin, cefotaxime. We hypothesized that the active site sequence requirements of CTX-M-mediated hydrolysis differ between classes of β-lactam antibiotics. Accordingly, we use codon randomization, antibiotic selection, and deep sequencing to determine the CTX-M active-site residues required for hydrolysis of cefotaxime and the penicillin, ampicillin. The study reveals positions required for hydrolysis of all β-lactams, as well as residues controlling substrate specificity. Further, CTX-M enzymes poorly hydrolyze the extended-spectrum cephalosporin, ceftazidime. We further show that the sequence requirements for ceftazidime hydrolysis follow those of cefotaxime, with the exception that key active-site omega loop residues are not required, and may be detrimental, for ceftazidime hydrolysis. These results provide insights into cephalosporin hydrolysis and demonstrate that changes to the active-site omega loop are likely required for the evolution of CTX-M-mediated ceftazidime resistance.

Slow Protein Dynamics Elicits New Enzymatic Functions by Means of Epistatic Interactions

Protein evolution depends on the adaptation of these molecules to different functional challenges. This occurs by tuning their biochemical, biophysical, and structural traits through the accumulation of mutations. While the role of protein dynamics in biochemistry is well recognized, there are limited examples providing experimental evidence of the optimization of protein dynamics during evolution. Here we report an NMR study of four variants of the CTX-M β-lactamases, in which the interplay of two mutations outside the active site enhances the activity against a cephalosporin substrate, ceftazidime. The crystal structures of these enzymes do not account for this activity enhancement. By using NMR, here we show that the combination of these two mutations increases the backbone dynamics in a slow timescale and the exposure to the solvent of an otherwise buried β-sheet. The two mutations located in this β-sheet trigger conformational changes in loops located at the opposite side of the active site. We postulate that the most active variant explores alternative conformations that enable binding of the more challenging substrate ceftazidime. The impact of the mutations in the dynamics is context-dependent, in line with the epistatic effect observed in the catalytic activity of the different variants. These results reveal the existence of a dynamic network in CTX-M β-lactamases that has been exploited in evolution to provide a net gain-of-function, highlighting the role of alternative conformations in protein evolution.

Small Molecules Promote Selective Denaturation and Degradation of Tubulin Heterodimers through a Low-Barrier Hydrogen Bond

Here, we report a novel mechanism to selectively degrade target proteins. 3-(3-Phenoxybenzyl)amino-β-carboline (PAC), a tubulin inhibitor, promotes selective degradation of αβ-tubulin heterodimers. Biochemical studies have revealed that PAC specifically denatures tubulin, making it prone to aggregation that predisposes it to ubiquitinylation and then degradation. The degradation is mediated by a single hydrogen bond formed between the pyridine nitrogen of PAC and βGlu198, which is identified as a low-barrier hydrogen bond (LBHB). In contrast, another two tubulin inhibitors that only form normal hydrogen bonds with βGlu198 exhibit no degradation effect. Thus, the LBHB accounts for the degradation. We then screened for compounds capable of forming an LBHB with βGlu198 and demonstrated that BML284, a Wnt signaling activator, also promotes tubulin heterodimer degradation through the LBHB. Our study provided a unique example of LBHB function and identified a novel approach to obtain tubulin degraders.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

Structural basis to repurpose boron-based proteasome inhibitors Bortezomib and Ixazomib as β-lactamase inhibitors

β-lactamases are a major cause of rapidly emerging and spreading antibiotic resistance. Currently β-lactamase inhibitors (BLIs) in clinical use act only on Ambler Class A, C and some class D lactamases. The urgent need to identify new BLIs recently lead to FDA approval of boron-based compounds BLIs, e.g. Vaborbactam. The boron-based proteasome inhibitors Bortezomib and Ixazomib are used in cancer therapy as multiple myeloma drugs but they also bind to Ser-/Thr- proteases. In this study we show the crystal structures of the β-lactamase CTX-M-14 with covalently bound Bortezomib and Ixazomib at high resolutions of 1.3 and 1.1 Å, respectively. Ixazomib is well defined in electron density whereas Bortezomib show some disorder which corresponds to weaker inhibition efficiency observed for Ixazomib. Both inhibitors mimic the deacylation transition state of β-lactam hydrolysis, because they replace the deacylating water molecule. We further investigate differences in binding of Bortezomib/Ixazomib to CTX-M-14 and its target proteases as well as known β-lactamase drugs. Our findings can help to use Bortezomib/Ixazomib as lead compounds for development of new BLIs.

Effective prediction of short hydrogen bonds in proteins via machine learning method

Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.

Mutation of the conserved Asp-Asp pair impairs the structure, function, and inhibition of CTX-M Class A β-lactamase

The Asp233-Asp246 pair is highly conserved in Class A β-lactamases, which hydrolyze β-lactam antibiotics. Here, we characterize its function using CTX-M-14 β-lactamase. The D233N mutant displayed decreased activity that is substrate-dependent, with reductions in k_cat /K_m ranging from 20% for nitrocefin to 6-fold for cefotaxime. In comparison, the mutation reduced the binding of a known reversible inhibitor by 10-fold. The mutant structures showed movement of the 213-219 loop and the loss of the Thr216-Thr235 hydrogen bond, which was restored by inhibitor binding. Mutagenesis of Thr216 further highlighted its contribution to CTX-M activity. These results demonstrate the importance of the aspartate pair to CTX-M hydrolysis of substrates with bulky side chains, while suggesting increased protein flexibility as a means to evolve drug resistance.

Protonation Equilibrium in the Active Site of the Photoactive Yellow Protein

The role and existence of low-barrier hydrogen bonds (LBHBs) in enzymatic and protein activity has been largely debated. An interesting case is that of the photoactive yellow protein (PYP). In this protein, two short HBs adjacent to the chromophore, p-coumaric acid (pCA), have been identified by X-ray and neutron diffraction experiments. However, there is a lack of agreement on the chemical nature of these H-bond interactions. Additionally, no consensus has been reached on the presence of LBHBs in the active site of the protein, despite various experimental and theoretical studies having been carried out to investigate this issue. In this work, we perform a computational study that combines classical and density functional theory (DFT)-based quantum mechanical/molecular mechanical (QM/MM) simulations to shed light onto this controversy. Furthermore, we aim to deepen our understanding of the chemical nature and dynamics of the protons involved in the two short hydrogen bonds that, in the dark state of PYP, connect pCA with the two binding pocket residues (E46 and Y42). Our results support the existence of a strong LBHB between pCA and E46, with the H fully delocalized and shared between both the carboxylic oxygen of E46 and the phenolic oxygen of pCA. Additionally, our findings suggest that the pCA interaction with Y42 can be suitably described as a typical short ionic H-bond of moderate strength that is fully localized on the phenolic oxygen of Y42.

Structural and Biochemical Characterization of the Novel CTX-M-151 Extended-Spectrum β-Lactamase and Its Inhibition by Avibactam

The diazabicyclooctane (DBO) inhibitor avibactam (AVI) reversibly inactivates most serine β-lactamases, including the CTX-M β-lactamases. Currently, more than 230 unique CTX-M members distributed in five clusters with less than 5% amino acid sequence divergence within each group have been described. Recently, a variant named CTX-M-151 was isolated from a Salmonella enterica subsp. enterica serovar Choleraesuis strain in Japan. This variant possesses a low degree of amino acid identity with the other CTX-Ms (63.2% to 69.7% with respect to the mature proteins), and thus it may represent a new subgroup within the family. CTX-M-151 hydrolyzes ceftriaxone better than ceftazidime (k _cat/K _m values 6,000-fold higher), as observed with CTX-Ms. CTX-M-151 is well inhibited by mechanism-based inhibitors like clavulanic acid (inactivation rate [k _inact]/inhibition constant [K_i ] = 0.15 μM^-1 · s^-1). For AVI, the apparent inhibition constant (K_{i app}), 0.4 μM, was comparable to that of KPC-2; the acylation rate (k₂/K) (37,000 M^-1 · s^-1) was lower than that for CTX-M-15, while the deacylation rate (k _off) (0.0015 s^-1) was 2- to 14-fold higher than those of other class A β-lactamases. The structure of the CTX-M-151/AVI complex (1.32 Å) reveals that AVI adopts a chair conformation with hydrogen bonds between the AVI carbamate and Ser70 and Ser237 at the oxyanion hole. Upon acylation, the side chain of Lys73 points toward Ser130, which is associated with the protonation of Glu166, supporting the role of Lys73 in the proton relay pathway and Glu166 as the general base in deacylation. To our knowledge, this is the first chromosomally encoded CTX-M in Salmonella Choleraesuis that shows similar hydrolytic preference toward cefotaxime (CTX) and ceftriaxone (CRO) to that toward ceftazidime (CAZ).

Lighting up solid states using a rubber

It is crucial and desirable to develop green and high-efficient strategies to regulate solid-state structures and their related material properties. However, relative to solution, it is more difficult to break and generate chemical bonds in solid states. In this work, a rubbing-induced photoluminescence on the solid states of ortho-pyridinil phenol family was achieved. This rubbing response relied on an accurately designed topochemical tautomerism, where a negative charge, exactly provided by the triboelectric effect of a rubber, can induce a proton transfer in a double H-bonded dimeric structure. This process instantaneously led to a bright-form tautomer that can be stabilized in the solid-state settings, leading to an up to over 450-fold increase of the fluorescent quantum yield of the materials. The property can be repeatedly used due to the reversibility of the tautomerism, enabling encrypted applications. Moreover, a further modification to the structure can be accomplished to achieve different properties, opening up more possibilities for the design of new-generation smart materials.

Changes in molecular properties due to stimuli response are critically important for the development of new materials. However, most processes are slow or inefficient in the solid state. Here the authors demonstrate property switching in the solid state using a rubbing-induced tautomerism in multiple hydrogen-bonded donor-acceptor couples.

Unraveling the origin of interactions of hydroxychloroquine with the receptor-binding domain of SARS-CoV-2 in aqueous medium

Interactions of hydroxychloroquin (HCQ) with the receptor binding domain (RBD) of SARS-CoV-2 were studied from atomistic simulation and ONIOM techniques. The key-residues of RBD responsible for the human transmission are recognized to be blocked in a heterogeneous manner with the favorable formation of key-residue:HCQ (1:1) complex. Such heterogeneity in binding was identified to be governed by the differential life-time of the hydrogen bonded water network anchoring HCQ and the key-residues. The intermolecular proton transfer facilitates the most favorable Lys417:HCQ complexation. The study demonstrates that off-target bindings of HCQ need to be minimized to efficiently prevent the transmission of SARS-CoV-2.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

A drug-resistant β-lactamase variant changes the conformation of its active-site proton shuttle to alter substrate specificity and inhibitor potency

Lys234 is one of the residues present in class A β-lactamases that is under selective pressure due to antibiotic use. Located adjacent to proton shuttle residue Ser130, it is suggested to play a role in proton transfer during catalysis of the antibiotics. The mechanism underpinning how substitutions in this position modulate inhibitor efficiency and substrate specificity leading to drug resistance is unclear. The K234R substitution identified in several inhibitor-resistant β-lactamase variants is associated with decreased potency of the inhibitor clavulanic acid, which is used in combination with amoxicillin to overcome β-lactamase-mediated antibiotic resistance. Here we show that for CTX-M-14 β-lactamase, whereas Lys234 is required for hydrolysis of cephalosporins such as cefotaxime, either lysine or arginine is sufficient for hydrolysis of ampicillin. Further, by determining the acylation and deacylation rates for cefotaxime hydrolysis, we show that both rates are fast, and neither is rate-limiting. The K234R substitution causes a 1500-fold decrease in the cefotaxime acylation rate but a 5-fold increase in kcat for ampicillin, suggesting that the K234R enzyme is a good penicillinase but a poor cephalosporinase due to slow acylation. Structural results suggest that the slow acylation by the K234R enzyme is due to a conformational change in Ser130, and this change also leads to decreased inhibition potency of clavulanic acid. Because other inhibitor resistance mutations also act through changes at Ser130 and such changes drastically reduce cephalosporin but not penicillin hydrolysis, we suggest that clavulanic acid paired with an oxyimino-cephalosporin rather than penicillin would impede the evolution of resistance.

Structural Basis and Binding Kinetics of Vaborbactam in Class A β-Lactamase Inhibition

Class A β-lactamases are a major cause of β-lactam resistance in Gram-negative bacteria. The recently FDA-approved cyclic boronate vaborbactam is a reversible covalent inhibitor of class A β-lactamases, including CTX-M extended-spectrum β-lactamase and KPC carbapenemase, both frequently observed in the clinic. Intriguingly, vaborbactam displayed different binding kinetics and cell-based activity for these two enzymes, despite their similarity. A 1.0-Å crystal structure of CTX-M-14 demonstrated that two catalytic residues, K73 and E166, are positively charged and neutral, respectively. Meanwhile, a 1.25-Å crystal structure of KPC-2 revealed a more compact binding mode of vaborbactam versus CTX-M-14, as well as alternative conformations of W105. Together with kinetic analysis of W105 mutants, the structures demonstrate the influence of this residue and the unusual conformation of the β3 strand on the inactivation rate, as well as the stability of the reversible covalent bond with S70. Furthermore, studies of KPC-2 S130G mutant shed light on the different impacts of S130 in the binding of vaborbactam versus avibactam, another recently approved β-lactamase inhibitor. Taken together, these new data provide valuable insights into the inhibition mechanism of vaborbactam and future development of cyclic boronate inhibitors.

Three-Step Spin State Transition and Hysteretic Proton Transfer in the Crystal of an Iron(II) Hydrazone Complex

A proton-electron coupling system, exhibiting unique bistability or multistability of the protonated state, is an attractive target for developing new switchable materials based on proton dynamics. Herein, we present an iron(II) hydrazone crystalline compound, which displays the stepwise transition and bistability of proton transfer at the crystal level. These phenomena are realized through the coupling with spin transition. Although the multi-step transition with hysteresis has been observed in various systems, the corresponding behavior of proton transfer has not been reported in crystalline systems; thus, the described iron(II) complex is the first example. Furthermore, because proton transfer occurs only in one of the two ligands and π electrons redistribute in it, the dipole moment of the iron(II) complexes changes with the proton transfer, wherein the total dipole moment in the crystal was canceled out owing to the antiferroelectric-like arrangement.

Discovery of an Orally Available Diazabicyclooctane Inhibitor (ETX0282) of Class A, C, and D Serine β-Lactamases

Multidrug resistant Gram-negative bacterial infections are an increasing public health threat due to rapidly rising resistance toward β-lactam antibiotics. The hydrolytic enzymes called β-lactamases are responsible for a large proportion of the resistance phenotype. β-Lactamase inhibitors (BLIs) can be administered in combination with β-lactam antibiotics to negate the action of the β-lactamases, thereby restoring activity of the β-lactam. Newly developed BLIs offer some advantage over older BLIs in terms of enzymatic spectrum but are limited to the intravenous route of administration. Reported here is a novel, orally bioavailable diazabicyclooctane (DBO) β-lactamase inhibitor. This new DBO, ETX1317, contains an endocyclic carbon–carbon double bond and a fluoroacetate activating group and exhibits broad spectrum activity against class A, C, and D serine β-lactamases. The ester prodrug of ETX1317, ETX0282, is orally bioavailable and, in combination with cefpodoxime proxetil, is currently in development as an oral therapy for multidrug resistant and carbapenem-resistant Enterobacterales infections.

The hydrolytic water molecule of Class A β-lactamase relies on the acyl-enzyme intermediate ES* for proper coordination and catalysis

Serine-based β-lactamases of Class A, C and D all rely on a key water molecule to hydrolyze and inactivate β-lactam antibiotics. This process involves two conserved catalytic steps. In the first acylation step, the β-lactam antibiotic forms an acyl-enzyme intermediate (ES*) with the catalytic serine residue. In the second deacylation step, an activated water molecule serves as nucleophile (WAT_Nu) to attack ES* and release the inactivated β-lactam. The coordination and activation of WAT_Nu is not fully understood. Using time-resolved x-ray crystallography and QM/MM simulations, we analyzed three intermediate structures of Class A β-lactamase PenP as it slowly hydrolyzed cephaloridine. WAT_Nu is centrally located in the apo structure but becomes slightly displaced away by ES* in the post-acylation structure. In the deacylation structure, WAT_Nu moves back and is positioned along the Bürgi–Dunitz trajectory with favorable energetic profile to attack ES*. Unexpectedly, WAT_Nu is also found to adopt a catalytically incompetent conformation in the deacylation structure forming a hydrogen bond with ES*. Our results reveal that ES* plays a significant role in coordinating and activating WAT_Nu through subtle yet distinct interactions at different stages of the catalytic process. These interactions may serve as potential targets to circumvent β-lactamase-mediated antibiotic resistance.

Extreme Catalytic Power of Ketosteroid Isomerase Related to the Reversal of Proton Dislocations in Hydrogen-Bond Network

Dynamic electrostatic catalytic field (DECF) vectors derived from transition state and reactant wavefunctions for the two-step reaction occurring within ketosteroid isomerase (KSI) have been calculated using MP2/aug-cc-pVTZ and lower theory levels to determine the magnitude of the catalytic effect and the optimal directions of proton transfers in the KSI hydrogen-bond network. The most surprising and meaningful finding is that the KSI catalytic activity is enhanced by proton dislocations proceeding in opposite directions for each of the two consecutive reaction steps in the same hydrogen network. Such a mechanism allows an ultrafast switching of the catalytic proton wire environment, possibly related to the exceptionally high KSI catalytic power.

Mechanism of proton transfer in class A β-lactamase catalysis and inhibition by avibactam

Gram-negative bacteria expressing class A β-lactamases pose a serious health threat due to their ability to inactivate all β-lactam antibiotics. The acyl-enzyme intermediate is a central milestone in the hydrolysis reaction catalyzed by these enzymes. However, the protonation states of the catalytic residues in this complex have never been fully analyzed experimentally due to inherent difficulties. To help unravel the ambiguity surrounding class A β-lactamase catalysis, we have used ultrahigh-resolution X-ray crystallography and the recently approved β-lactamase inhibitor avibactam to trap the acyl-enzyme complex of class A β-lactamase CTX-M-14 at varying pHs. A 0.83-Å-resolution CTX-M-14 complex structure at pH 7.9 revealed a neutral state for both Lys73 and Glu166. Furthermore, the avibactam hydroxylamine-O-sulfonate group conformation varied according to pH, and this conformational switch appeared to correspond to a change in the Lys73 protonation state at low pH. In conjunction with computational analyses, our structures suggest that Lys73 has a perturbed acid dissociation constant (pK_a) compared with acyl-enzyme complexes with β-lactams, hindering its function to deprotonate Glu166 and the initiation of the deacylation reaction. Further NMR analysis demonstrated Lys73 pK_a to be ∼5.2 to 5.6. Together with previous ultrahigh-resolution crystal structures, these findings enable us to follow the proton transfer process of the entire acylation reaction and reveal the critical role of Lys73. They also shed light on the stability and reversibility of the avibactam carbamoyl acyl-enzyme complex, highlighting the effect of substrate functional groups in influencing the protonation states of catalytic residues and subsequently the progression of the reaction.

β-Lactamase of Mycobacterium tuberculosis Shows Dynamics in the Active Site That Increase upon Inhibitor Binding

The Mycobacterium tuberculosis β-lactamase BlaC is a broad-spectrum β-lactamase that can convert a range of β-lactam antibiotics. Enzymes with low specificity are expected to exhibit active-site flexibility. To probe the motions in BlaC, we studied the dynamic behavior in solution using nuclear magnetic resonance (NMR) spectroscopy. ¹⁵N relaxation experiments show that BlaC is mostly rigid on the pico- to nanosecond timescale. Saturation transfer experiments indicate that also on the high-millisecond timescale BlaC is not dynamic. Using relaxation dispersion experiments, clear evidence was obtained for dynamics in the low-millisecond range, with an exchange rate of ca. 860 s^-1 The dynamic amide groups are localized in the active site. Upon formation of an adduct with the inhibitor avibactam, extensive line broadening occurs, indicating an increase in magnitude of the active-site dynamics. Furthermore, the rate of the motions increases significantly. Upon reaction with the inhibitor clavulanic acid, similar line broadening is accompanied by duplication of NMR signals, indicative of at least one additional, slower exchange process (exchange rate, k _ex, of <100 s^-1), while for this inhibitor also loss of pico- to nanosecond timescale rigidity is observed for some amides in the α domain. Possible sources of the observed dynamics, such as motions in the omega loop and rearrangements of active-site residues, are discussed. The increase in dynamics upon ligand binding argues against a model of inhibitor binding through conformational selection. Rather, the induced dynamics may serve to maximize the likelihood of sampling the optimal conformation for hydrolysis of the bound ligand.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Structural Basis for Substrate Specificity and Carbapenemase Activity of OXA-48 Class D β-Lactamase

Carbapenem-hydrolyzing class D β-lactamases (CHDLs) are a diverse family of enzymes that are rapidly becoming the predominant cause of bacterial resistance against β-lactam antibiotics in many regions of the world. OXA-48, an atypical member of CHDLs, is one of the most frequently observed in the clinic and exhibits a unique substrate profile. We applied X-ray crystallography to OXA-48 complexes with multiple β-lactam antibiotics to elucidate this enzyme's carbapenemase activity and its preference of imipenem over meropenem and other substrates such as cefotaxime. In particular, we obtained acyl-enzyme complexes of OXA-48 with imipenem, meropenem, faropenem, cefotaxime, and cefoxitin, and a product complex with imipenem. Importantly, the product complex captures a key reaction milestone with the newly generated carboxylate group still in the oxyanion hole, and represents the first such complex with a wild-type serine β-lactamase. A potential hydrogen bond is observed between the two carboxylate groups from the product and the carbamylated Lys73, representing the stage immediately after the breakage of the acyl-enzyme bond where the product carboxylate would be neutral. The placement of the product carboxylate also illustrates the approximate transient location of the deacylation water that has long eluded structural characterization in class D β-lactamases. Additionally, comparing the product complex with the acyl-enzyme intermediates provides new insights into the various mechanisms by which specific side chain groups hinder the access of the deacylation water to the acyl-enzyme linkage, especially in meropenem. Taken together, these data offer valuable information on the substrate specificity of OXA-48 and the catalytic mechanism of CHDLs.

Probing the role of the conserved residue Glu166 in a class A β-lactamase using neutron and X-ray protein crystallography

The amino-acid sequence of the Toho-1 β-lactamase contains several conserved residues in the active site, including Ser70, Lys73, Ser130 and Glu166, some of which coordinate a catalytic water molecule. This catalytic water molecule is essential in the acylation and deacylation parts of the reaction mechanism through which Toho-1 inactivates specific antibiotics and provides resistance to its expressing bacterial strains. To investigate the function of Glu166 in the acylation part of the catalytic mechanism, neutron and X-ray crystallographic studies were performed on a Glu166Gln mutant. The structure of this class A β-lactamase mutant provides several insights into its previously reported reduced drug-binding kinetic rates. A joint refinement of both X-ray and neutron diffraction data was used to study the effects of the Glu166Gln mutation on the active site of Toho-1. This structure reveals that while the Glu166Gln mutation has a somewhat limited impact on the positions of the conserved amino acids within the active site, it displaces the catalytic water molecule from the active site. These subtle changes offer a structural explanation for the previously observed decreases in the binding of non-β-lactam inhibitors such as the recently developed diazobicyclooctane inhibitor avibactam.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

Cysteine SH and Glutamate COOH Contributions to [NiFe] Hydrogenase Proton Transfer Revealed by Highly Sensitive FTIR Spectroscopy

A [NiFe] hydrogenase (H₂ ase) is a proton-coupled electron transfer enzyme that catalyses reversible H₂ oxidation; however, its fundamental proton transfer pathway remains unknown. Herein, we observed the protonation of Cys546-SH and Glu34-COOH near the Ni-Fe site with high-sensitivity infrared difference spectra by utilizing Ni-C-to-Ni-L and Ni-C-to-Ni-SI_a photoconversions. Protonated Cys546-SH in the Ni-L state was verified by the observed SH stretching frequency (2505 cm^-1 ), whereas Cys546 was deprotonated in the Ni-C and Ni-SI_a states. Glu34-COOH was double H-bonded in the Ni-L state, as determined by the COOH stretching frequency (1700 cm^-1 ), and single H-bonded in the Ni-C and Ni-SI_a states. Additionally, a stretching mode of an ordered water molecule was observed in the Ni-L and Ni-C states. These results elucidate the organized proton transfer pathway during the catalytic reaction of a [NiFe] H₂ ase, which is regulated by the H-bond network of Cys546, Glu34, and an ordered water molecule.

Underestimated Noncovalent Interactions in Protein Data Bank

Noncovalent interactions (NCIs) play essential roles in the structure and function of biomacromolecules. There are various NCIs, e.g., hydrogen bonds (HBs), cation-π and π-π interactions, and ionic bonds, among which HBs are the most widespread and well-studied. By utilizing the ratio of the observed HBs over pseudo HBs (1.0 Å longer than the HB distance criteria without angle constraints), we demonstrated that HBs in both protein-ligand and protein-protein interfaces are overlooked in structures deposited in PDB. After the QM/MM optimization of 12 protein-ligand complexes, we showed that the overlooked HBs could be recovered. With a systematic search in the PDB, we found that the HB number per residue (N_HB/R) in proteins decreases as structural resolution becomes lower, implying that HBs are overlooked even today, regardless of the type of refinement approach used. Similarly, cation-π, π-π, and ionic interactions were found to be significantly lost, manifesting the universal underestimation of various NCIs. Considering the vital role of NCIs, it is important to recover the NCIs to facilitate drug design, to explore protein-protein interaction, and to study protein structure and function.

Unraveling the structural and chemical features of biological short hydrogen bonds

Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.

The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.

Influence of substrates and inhibitors on the structure of Klebsiella pneumoniae carbapenemase-2

The hydrolysis of last resort carbapenem antibiotics by Klebsiella pneumoniae carbapenemase-2 (KPC-2) presents a significant danger to global health. Combined with horizontal gene transfer, the emergence KPC-2 threatens to quickly expand carbapenemase activity to ever increasing numbers of pathogens. Our understanding of KPC-2 has greatly increased over the past decade thanks, in great part, to 20 crystal structures solved by groups around the world. These include apo KPC-2 structures, along with structures featuring a library of 10 different inhibitors representing diverse structural and functional classes. Herein we focus on cataloging the available KPC-2 structures and presenting a discussion of key aspects of each structure and important relationships between structures. Although the available structures do not provide information on dynamic motions with KPC-2, and the family of structures indicates small conformational changes across a wide array of bound inhibitors, substrates, and products, the structures provide a strong foundation for additional studies in the coming years to discover new KPC-2 inhibitors.

Impact statement

The work herein is important to the field as it provides a clear and succinct accounting of available KPC-2 structures. The work advances the field by collecting and analyzing differences and similarities across the available structures. This work features new analyses and interpretations of the existing structures which will impact the field in a positive way by making structural insights more widely available among the beta-lactamase community.

An Empirical Study of Amide-Heteroarene π-Stacking Interactions Using Reversible Inhibitors of a Bacterial Serine Hydrolase

Compared to aryl-aryl π-stacking interactions, the analogous stacking of heteroarenes on amide π systems is less well understood and vastly underutilized in structure-based drug design. Recent theoretical studies have delineated the important geometric coordinates of the interaction, some of which have been confirmed with synthetic model systems based on Rebek imides. Unfortunately, a broadly useful and tractable protein-ligand model system of this interaction has remained elusive. Here we employed a known inhibitor scaffold to study π-stacking of diverse heteroarene substituents on the amide face of Gly238 in the cephalosporinases CTX-M-14 and CTX-M-27. Biochemical inhibition constants (K _i) and biophysical binding constants (K _d) were determined for nineteen new analogues against both enzymes, while multiple high-resolution co-crystal structures revealed remarkably consistent placement of the probe heteroarene on Gly238. The data presented support the predicted importance of opposing dipoles in amide-heteroarene interactions and should be useful for evaluating other theoretical predictions concerning these interactions.

β-Lactamases and β-Lactamase Inhibitors in the 21st Century

The β-lactams retain a central place in the antibacterial armamentarium. In Gram-negative bacteria, β-lactamase enzymes that hydrolyze the amide bond of the four-membered β-lactam ring are the primary resistance mechanism, with multiple enzymes disseminating on mobile genetic elements across opportunistic pathogens such as Enterobacteriaceae (e.g., Escherichia coli) and non-fermenting organisms (e.g., Pseudomonas aeruginosa). β-Lactamases divide into four classes; the active-site serine β-lactamases (classes A, C and D) and the zinc-dependent or metallo-β-lactamases (MBLs; class B). Here we review recent advances in mechanistic understanding of each class, focusing upon how growing numbers of crystal structures, in particular for β-lactam complexes, and methods such as neutron diffraction and molecular simulations, have improved understanding of the biochemistry of β-lactam breakdown. A second focus is β-lactamase interactions with carbapenems, as carbapenem-resistant bacteria are of grave clinical concern and carbapenem-hydrolyzing enzymes such as KPC (class A) NDM (class B) and OXA-48 (class D) are proliferating worldwide. An overview is provided of the changing landscape of β-lactamase inhibitors, exemplified by the introduction to the clinic of combinations of β-lactams with diazabicyclooctanone and cyclic boronate serine β-lactamase inhibitors, and of progress and strategies toward clinically useful MBL inhibitors. Despite the long history of β-lactamase research, we contend that issues including continuing unresolved questions around mechanism; opportunities afforded by new technologies such as serial femtosecond crystallography; the need for new inhibitors, particularly for MBLs; the likely impact of new β-lactam:inhibitor combinations and the continuing clinical importance of β-lactams mean that this remains a rewarding research area.

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.

Nuclear quantum effect and H/D isotope effect on Cl· + (H2O)n → HCl + OH·(H2O)n−1 (n = 1–3) reactions

Cl· + (H₂O)_n → HCl + OH·(H₂O)_n−1 (n = 1–3) reactions are fundamental and important ones in atmospheric chemistry. In this study, we focused on the nuclear quantum effect (NQE) of the hydrogen nucleus on these reactions with the aid of the multicomponent quantum mechanics (MC_QM) method, which can directly take account of NQE of light nuclei. Our study reveals that the NQE of the hydrogen nucleus lowers the activation barriers of the reactions and enhances the catalytic effects of second and third water molecules. In particular, we find that (i) the NQE of the proton removes the activation barrier of the reverse reaction of HCl + OH· → Cl· + H₂O, and (ii) the catalytic effect of the third water molecule appears in only our MC_QM calculation. We also analyze the H/D isotope effects on these reactions by using the MC_QM method.

Cl·+ (H₂O)_n → HCl + OH(H₂O)_n−1 (n = 1–3) reactions have been investigated using multicomponent quantum mechanics method, which can take account of the nuclear quantum effect of proton and deuteron.

Mechanisms of proton relay and product release by Class A β-lactamase at ultrahigh resolution

The β-lactam antibiotics inhibit penicillin-binding proteins (PBPs) by forming a stable, covalent, acyl-enzyme complex. During the evolution from PBPs to Class A β-lactamases, the β-lactamases acquired Glu166 to activate a catalytic water and cleave the acyl-enzyme bond. Here we present three product complex crystal structures of CTX-M-14 Class A β-lactamase with a ruthenocene-conjugated penicillin-a 0.85 Å resolution structure of E166A mutant complexed with the penilloate product, a 1.30 Å resolution complex structure of the same mutant with the penicilloate product, and a 1.18 Å resolution complex structure of S70G mutant with a penicilloate product epimer-shedding light on the catalytic mechanisms and product inhibition of PBPs and Class A β-lactamases. The E166A-penilloate complex captured the hydrogen bonding network following the protonation of the leaving group and, for the first time, unambiguously show that the ring nitrogen donates a proton to Ser130, which in turn donates a proton to Lys73. These observations indicate that in the absence of Glu166, the equivalent lysine would be neutral in PBPs and therefore capable of serving as the general base to activate the catalytic serine. Together with previous results, this structure suggests a common proton relay network shared by Class A β-lactamases and PBPs, from the catalytic serine to the lysine, and ultimately to the ring nitrogen. Additionally, the E166A-penicilloate complex reveals previously unseen conformational changes of key catalytic residues during the release of the product, and is the first structure to capture the hydrolyzed product in the presence of an unmutated catalytic serine.

DATABASE

Structural data are available in the PDB database under the accession numbers 5TOP, 5TOY, and 5VLE.

High Pressure Spectroscopic Investigation on Proton Transfer in Squaric Acid and 4,4'-Bipyridine Co-crystal

In attempt to the obtain detailed geometric information of proton transfer compound (subsequently denote as SQBP) formed between squaric acid (SQ)and 4,4′-bipyridine(BP), and to investigate the mechanisms of pressure-induced double proton transfer and related structural phase transition, we carried out in-situ high pressure Raman spectroscopy of SQBP up to 20 GPa. A solid-solid phase transition together with double proton transfer phenomenon was confirmed by Raman spectroscopy at about 1.5 GPa, and the activation of C = O stretching mode in Raman spectra indicates a square-ring structure of SQ with four symmetric C = O bond formation. These results are further supported by first-principals calculations and in-situ high pressure infrared absorption spectroscopy. Additionally, Raman intensity analysis suggests that a higher-order phase transition with planar BP molecular structure occurred in the pressure range of 3~6 GPa. As a result, the π electron delocalization effect in BP dominated the intensity enhancement of C = O stretching mode in SQ. To the best of our knowledge, this is the first time observation of the intensity enhancement of proton donor’s normal modes induced by proton acceptor’s π electron delocalization.

A QM/MM study on the enzymatic inactivation of cefotaxime

The reaction between the antibiotic cefotaxime and the CTX-M-14 class A serine hydrolase is addressed from a theoretical point of view, by means of hybrid quantum mechanics/molecular mechanical (QM/MM) calculations, adopting a new approach that postulates that the residue Ser70 itself should play the role of the acid-base species required for the cefotaxime acylation. The proposed mechanism differs from earlier proposals existing in literature for other class A β-lactamases. The results confirm the hypothesis, and show that the reaction should occur via a concerted mechanism in which the acylation of the lactam carbonyl carbon, protonation of the N₇ lactam atom, and opening of the β-lactam ring occurs simultaneously. Exploration of the potential energy surface shows three critical points, associated with reactants, transition state and product. The transition state is characterized by frequency, intrinsic reaction coordinate, atomic charge, and bond orders calculations. The calculated activation barrier is 20 kcal mol^-1, and the reaction appears to be slightly endothermic by about 12 kcal mol^-1. We conclude that this approach is feasible, and should be considered as an alternative mechanism or may exist in competition with others already published in the literature. This information should be useful for the design of novel antibiotics and β-lactamase inhibitors. Graphical abstract Three-dimensional view of the potential energy surface of cefotaxime.

Molecular Basis of Substrate Recognition and Product Release by the Klebsiella pneumoniae Carbapenemase (KPC-2)

Carbapenem-resistant Enterobacteriaceae are resistant to most β-lactam antibiotics due to the production of the Klebsiella pneumoniae carbapenemase (KPC-2) class A β-lactamase. Here, we present the first product complex crystal structures of KPC-2 with β-lactam antibiotics containing hydrolyzed cefotaxime and faropenem. They provide experimental insights into substrate recognition by KPC-2 and its unique cephalosporinase/carbapenemase activity. These structures also represent the first product complexes for a wild-type serine β-lactamase, elucidating the product release mechanism of these enzymes in general.