Electrostatic analysis of TEM1 beta-lactamase: effect of substrate binding, steep potential gradients and consequences of site-directed mutations

An in vivo selection system with tightly regulated gene expression enables directed evolution of highly efficient enzymes

In vivo selection systems are powerful tools for directed evolution of enzymes. The selection pressure of the systems can be tuned by regulating the expression levels of the catalysts. In this work, we engineered a selection system for laboratory evolution of highly active enzymes by incorporating a translationally suppressing cis repressor as well as an inducible promoter to impart stringent and tunable selection pressure. We demonstrated the utility of our selection system by performing directed evolution experiments using TEM β-lactamase as the model enzyme. Five evolutionary rounds afforded a highly active variant exhibiting 440-fold improvement in catalytic efficiency. We also showed that, without the cis repressor, the selection system cannot provide sufficient selection pressure required for evolving highly efficient TEM β-lactamase. The selection system should be applicable for the exploration of catalytic perfection of a wide range of enzymes.

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.

Tabtoxinine-β-lactam is a “stealth” β-lactam antibiotic that evades β-lactamase-mediated antibiotic resistance

Tabtoxinine-β-lactam (TβL) is a phytotoxin produced by plant pathogenic strains of Pseudomonas syringae.

Can molecular dynamics and QM/MM solve the penicillin binding protein protonation puzzle?

Benzylpenicillin, a member of the β-lactam antibiotic class, has been widely used to combat bacterial infections since 1947. The general mechanism is well-known: a serine protease enzyme (i.e., DD-peptidase) forms a long lasting intermediate with the lactam ring of the antibiotic known as acylation, effectively preventing biosynthesis of the bacterial cell wall. Despite this overall mechanistic understanding, many details of binding and catalysis are unclear. Specifically, there is ongoing debate about active site protonation states and the role of general acids/bases in the reaction. Herein, a unique combination of MD simulations, QM/MM minimizations, and QM/MM orbital analyses is combined with systematic variation of active site residue protonation states. Critical interactions that maximize the stability of the bound inhibitor are examined and used as metrics. This approach was validated by examining cefoxitin interactions in the CTX-M β-lactamase from E. coli and compared to an ultra high-resolution (0.88 Å) crystal structure. Upon confirming the approach used, an investigation of the preacylated Streptomyces R61 active site with bound benzylpenicillin was performed, varying the protonation states of His298 and Lys65. We concluded that protonated His298 and deprotonated Lys65 are most likely to exist in the R61 active site.

Can inhibitor-resistant substitutions in the Mycobacterium tuberculosis β-Lactamase BlaC lead to clavulanate resistance?: a biochemical rationale for the use of β-lactam-β-lactamase inhibitor combinations

The current emergence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis calls for novel treatment strategies. Recently, BlaC, the principal β-lactamase of Mycobacterium tuberculosis, was recognized as a potential therapeutic target. The combination of meropenem and clavulanic acid, which inhibits BlaC, was found to be effective against even extensively drug-resistant M. tuberculosis strains when tested in vitro. Yet there is significant concern that drug resistance against this combination will also emerge. To investigate the potential of BlaC to evolve variants resistant to clavulanic acid, we introduced substitutions at important amino acid residues of M. tuberculosis BlaC (R220, A244, S130, and T237). Whereas the substitutions clearly led to in vitro clavulanic acid resistance in enzymatic assays but at the expense of catalytic activity, transformation of variant BlaCs into an M. tuberculosis H37Rv background revealed that impaired inhibition of BlaC did not affect inhibition of growth in the presence of ampicillin and clavulanate. From these data we propose that resistance to β-lactam-β-lactamase inhibitor combinations will likely not arise from structural alteration of BlaC, therefore establishing confidence that this therapeutic modality can be part of a successful treatment regimen against M. tuberculosis.

Variations within class-A β-lactamase physiochemical properties reflect evolutionary and environmental patterns, but not antibiotic specificity

The bacterial enzyme β-lactamase hydrolyzes the β-lactam ring of penicillin and chemically related antibiotics, rendering them ineffective. Due to rampant antibiotic overuse, the enzyme is evolving new resistance activities at an alarming rate. Related, the enzyme's global physiochemical properties exhibit various amounts of conservation and variability across the family. To that end, we characterize the extent of property conservation within twelve different class-A β-lactamases, and conclusively establish that the systematic variations therein parallel their evolutionary history. Large and systematic differences within electrostatic potential maps and pairwise residue-to-residue couplings are observed across the protein, which robustly reflect phylogenetic outgroups. Other properties are more conserved (such as residue pKa values, electrostatic networks, and backbone flexibility), yet they also have systematic variations that parallel the phylogeny in a statistically significant way. Similarly, the above properties also parallel the environmental condition of the bacteria they are from in a statistically significant way. However, it is interesting and surprising that the only one of the global properties (protein charge) parallels the functional specificity patterns; meaning antibiotic resistance activities are not significantly constraining the global physiochemical properties. Rather, extended spectrum activities can emerge from the background of nearly any set of electrostatic and dynamic properties.

On 3LEZ, a deep-sea halophilic protein with in vitro class-a β-lactamase activity: molecular-dynamics, docking, and reactivity simulations

This work shows that a deep-sea protein, 3LEZ, with known in vitro β-lactamase activity, proved stable, substantially in the conformation detected by X-ray diffraction of the crystal, when subjected to molecular-dynamics (MD) simulations under conditions compatible with shallow seas. Docking simulations showed that the β-lactamase active site S85 of 3LEZ (S70 in Ambler numbering) is the preferential binding pocket for not only β-lactam antibiotics and inhibitors, but, surprisingly, also for a wide variety of other biologically active compounds in various chemical classes, including marine metabolites. In line with the in vitro β-lactamase activity, a) affinities on docking β-lactam antibiotics and inhibitors onto 3LEZ were found to roughly parallel published K(m) and K(i) values, obtained from MichaelisMenten kinetics under room conditions, and b) DFT calculations agreed with experiments that the irreversible reaction of the β-lactamase inhibitor clavulanic acid with the whole S85 catalytic center of 3LEZ is spontaneous. These observations must be viewed in the light that a) the compounds in other chemical classes showed comparable affinities, and, in some cases, even higher than β-lactams, for the S85 active site, b) K(m) and K(i) data are not available at the high hydrostatic pressure of the deep sea, where 3LEZ is believed to have evolved, c) an inverse order of affinities for the β-lactams, with respect to both experimentation and simulations at room conditions, was observed from comparative docking simulations with 3LEZ derived from MD under high hydrostatic pressure. Although MD requires a general assessment for high hydrostatic pressure before c) above is given the same weight as all other observations, this work questions the conclusion that the in vitro determined β-lactamase activity represents the ecological role of 3LEZ.

Three decades of beta-lactamase inhibitors

Since the introduction of penicillin, beta-lactam antibiotics have been the antimicrobial agents of choice. Unfortunately, the efficacy of these life-saving antibiotics is significantly threatened by bacterial beta-lactamases. beta-Lactamases are now responsible for resistance to penicillins, extended-spectrum cephalosporins, monobactams, and carbapenems. In order to overcome beta-lactamase-mediated resistance, beta-lactamase inhibitors (clavulanate, sulbactam, and tazobactam) were introduced into clinical practice. These inhibitors greatly enhance the efficacy of their partner beta-lactams (amoxicillin, ampicillin, piperacillin, and ticarcillin) in the treatment of serious Enterobacteriaceae and penicillin-resistant staphylococcal infections. However, selective pressure from excess antibiotic use accelerated the emergence of resistance to beta-lactam-beta-lactamase inhibitor combinations. Furthermore, the prevalence of clinically relevant beta-lactamases from other classes that are resistant to inhibition is rapidly increasing. There is an urgent need for effective inhibitors that can restore the activity of beta-lactams. Here, we review the catalytic mechanisms of each beta-lactamase class. We then discuss approaches for circumventing beta-lactamase-mediated resistance, including properties and characteristics of mechanism-based inactivators. We next highlight the mechanisms of action and salient clinical and microbiological features of beta-lactamase inhibitors. We also emphasize their therapeutic applications. We close by focusing on novel compounds and the chemical features of these agents that may contribute to a "second generation" of inhibitors. The goal for the next 3 decades will be to design inhibitors that will be effective for more than a single class of beta-lactamases.

Molecular mechanisms of antibiotic resistance: QM/MM modelling of deacylation in a class A beta-lactamase

Modelling of the first step of the deacylation reaction of benzylpenicillin in the E. coli TEM1 beta-lactamase (with B3LYP/6-31G + (d)//AM1-CHARMM22 quantum mechanics/molecular mechanics methods) shows that a mechanism in which Glu166 acts as the base to deprotonate a conserved water molecule is both energetically and structurally consistent with experimental data; the results may assist the design of new antibiotics and beta-lactamase inhibitors.

Dynamical aspects of TEM-1 beta-lactamase probed by molecular dynamics

The dynamical aspects of the fully hydrated TEM-1 beta-lactamase have been determined by a 5 ns Molecular Dynamics simulation. Starting from the crystallographic coordinates, the protein shows a relaxation in water with an overall root mean square deviation from the crystal structure increasing up to 0.17 nm, within the first nanosecond. Then a plateau is reached and the molecule fluctuates around an equilibrium conformation. The results obtained in the first nanosecond are in agreement with those of a previous simulation (Diaz et al., J. Am. Chem. Soc., (2003) 125, 672-684). The successive equilibrium conformation in solution shows an increased mobility characterized by the following aspects. A flap-like translational motion anchors the omega-loop to the body of the enzyme. A relevant part of the backbone dynamics implies a rotational motion of one domain relative to the other. The water molecules in the active site can exchange with different residence times. The H-bonding networks formed by the catalytic residues are frequently interrupted by water molecules that could favour proton transfer reactions. An additional simulation, where the aspartyl dyad D214-D233 was considered fully deprotonated, shows that the active site is destabilized.

Ab Initio QM/MM Study of Class A β-Lactamase Acylation: Dual Participation of Glu166 and Lys73 in a Concerted Base Promotion of Ser70

Beta-lactamase acquisition is the most prevalent basis for Gram-negative bacteria resistance to the beta-lactam antibiotics. The mechanism used by the most common class A Gram-negative beta-lactamases is serine acylation followed by hydrolytic deacylation, destroying the beta-lactam. The ab initio quantum mechanical/molecular mechanical (QM/MM) calculations, augmented by extensive molecular dynamics simulations reported herein, describe the serine acylation mechanism for the class A TEM-1 beta-lactamase with penicillanic acid as substrate. Potential energy surfaces (based on approximately 350 MP2/6-31+G calculations) reveal the proton movements that govern Ser70 tetrahedral formation and then collapse to the acyl-enzyme. A remarkable duality of mechanism for tetrahedral formation is implicated. Following substrate binding, the pathway initiates by a low energy barrier (5 kcal mol(-1)) and an energetically favorable transfer of a proton from Lys73 to Glu166, through the catalytic water molecule and Ser70. This gives unprotonated Lys73 and protonated Glu166. Tetrahedral formation ensues in a concerted general base process, with Lys73 promoting Ser70 addition to the beta-lactam carbonyl. Moreover, the three-dimensional potential energy surface also shows that the previously proposed pathway, involving Glu166 as the general base promoting Ser70 through a conserved water molecule, exists in competition with the Lys73 process. The existence of two routes to the tetrahedral species is fully consistent with experimental data for mutant variants of the TEM beta-lactamase.

Mechanisms of Antibiotic Resistance: QM/MM Modeling of the Acylation Reaction of a Class A β-Lactamase with Benzylpenicillin

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.

Crystal Structures of the Apo and Penicillin-acylated Forms of the BlaR1 β-Lactam Sensor of Staphylococcus aureus

Staphylococcus aureus is among the most prevalent and antibiotic-resistant of pathogenic bacteria. The resistance of S. aureus to prototypal beta-lactam antibiotics is conferred by two mechanisms: (i) secretion of hydrolytic beta-lactamase enzymes and (ii) production of beta-lactam-insensitive penicillin-binding proteins (PBP2a). Despite their distinct modes of resistance, expression of these proteins is controlled by similar regulation systems, including a repressor (BlaI/MecI) and a multidomain transmembrane receptor (BlaR1/MecR1). Resistance is triggered in response to a covalent binding event between a beta-lactam antibiotic and the extracellular sensor domain of BlaR1/MecR1 by transduction of the binding signal to an intracellular protease domain capable of repressor inactivation. This study describes the first crystal structures of the sensor domain of BlaR1 (BlaRS) from S. aureus in both the apo and penicillin-acylated forms. The structures show that the sensor domain resembles the beta-lactam-hydrolyzing class D beta-lactamases, but is rendered a penicillin-binding protein due to the formation of a very stable acyl-enzyme. Surprisingly, conformational changes upon penicillin binding were not observed in our structures, supporting the hypothesis that transduction of the antibiotic-binding signal into the cytosol is mediated by additional intramolecular interactions of the sensor domain with an adjacent extracellular loop in BlaR1.

The importance of a critical protonation state and the fate of the catalytic steps in class A beta-lactamases and penicillin-binding proteins

Beta-lactamases and penicillin-binding proteins are bacterial enzymes involved in antibiotic resistance to beta-lactam antibiotics and biosynthetic assembly of cell wall, respectively. Members of these large families of enzymes all experience acylation by their respective substrates at an active site serine as the first step in their catalytic activities. A Ser-X-X-Lys sequence motif is seen in all these proteins, and crystal structures demonstrate that the side-chain functions of the serine and lysine are in contact with one another. Three independent methods were used in this report to address the question of the protonation state of this important lysine (Lys-73) in the TEM-1 beta-lactamase from Escherichia coli. These techniques included perturbation of the pK(a) of Lys-73 by the study of the gamma-thialysine-73 variant and the attendant kinetic analyses, investigation of the protonation state by titration of specifically labeled proteins by nuclear magnetic resonance, and by computational treatment using the thermodynamic integration method. All three methods indicated that the pK(a) of Lys-73 of this enzyme is attenuated to 8.0-8.5. It is argued herein that the unique ground-state ion pair of Glu-166 and Lys-73 of class A beta-lactamases has actually raised the pK(a) of the active site lysine to 8.0-8.5 from that of the parental penicillin-binding protein. Whereas we cannot rule out that Glu-166 might activate the active site water, which in turn promotes Ser-70 for the acylation event, such as proposed earlier, we would like to propose as a plausible alternative for the acylation step the possibility that the ion pair would reconfigure to the protonated Glu-166 and unprotonated Lys-73. As such, unprotonated Lys-73 could promote serine for acylation, a process that should be shared among all active-site serine beta-lactamases and penicillin-binding proteins.

Understanding the acylation mechanisms of active-site serine penicillin-recognizing proteins: a molecular dynamics simulation study

Beta-lactam antibiotics inhibit enzymes involved in the last step of peptidoglycan synthesis. These enzymes, also identified as penicillin-binding proteins (PBPs), form a long-lived acyl-enzyme complex with beta-lactams. Antibiotic resistance is mainly due to the production of beta-lactamases, which are enzymes that hydrolyze the antibiotics and so prevent them reaching and inactivating their targets, and to mutations of the PBPs that decrease their affinity for the antibiotics. In this study, we present a theoretical study of several penicillin-recognizing proteins complexed with various beta-lactam antibiotics. Hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations have been performed to understand the role of several residues, and pK(a) calculations have also been done to determine their protonation state. We analyze the differences between the beta-lactamase TEM-1, the membrane-bound PBP2x of Streptococcus pneumoniae, and the soluble DD-transpeptidase of Streptomyces K15.

A new TEM-derived extended-spectrum beta-lactamase (TEM-91) with an R164C substitution at the omega-loop confers ceftazidime resistance

A new plasmid-mediated TEM-derived extended-spectrum beta-lactamase, TEM-91, was identified in a ceftazidime-resistant (MIC, >128 microg per ml) Escherichia coli strain isolated in 1996 in Japan. TEM-91 has three amino acid substitutions, R164C, M184T, and E240K, compared with TEM-1 penicillinase. The isoelectric point (pI), K(m), and k(cat) of TEM-91 for ceftazidime were 5.7, 179 microM, and 29.0 s(-1), respectively. The K(i) of clavulanic acid for ceftazidime hydrolysis was 30.3 nM.

Ultrahigh resolution structure of a class A beta-lactamase: on the mechanism and specificity of the extended-spectrum SHV-2 enzyme

Bacterial beta-lactamases hydrolyze beta-lactam antibiotics such as penicillins and cephalosporins. The TEM-type class A beta-lactamase SHV-2 is a natural variant that exhibits activity against third-generation cephalosporins normally resistant to hydrolysis by class A enzymes. SHV-2 contains a single Gly238Ser change relative to the wild-type enzyme SHV-1. Crystallographic refinement of a model including hydrogen atoms gave R and R(free) of 12.4% and 15.0% for data to 0.91 A resolution. The hydrogen atom on the O(gamma) atom of the reactive Ser70 is clearly seen for the first time, bridging to the water molecule activated by Glu166. Though hydrogen atoms on the nearby Lys73 are not seen, this observation of the Ser70 hydrogen atom and the hydrogen bonding pattern around Lys73 indicate that Lys73 is protonated. These findings support a role for the Glu166-water couple, rather than Lys73, as the general base in the deprotonation of Ser70 in the acylation process of class A beta-lactamases. Overlay of SHV-2 with SHV-1 shows a significant 1-3 A displacement in the 238-242 beta-strand-turn segment, making the beta-lactam binding site more open to newer cephalosporins with large C7 substituents and thereby expanding the substrate spectrum of the variant enzyme. The OH group of the buried Ser238 side-chain hydrogen bonds to the main-chain CO of Asn170 on the Omega loop, that is unaltered in position relative to SHV-1. This structural role for Ser238 in protein-protein binding makes less likely its hydrogen bonding to oximino cephalosporins such as cefotaxime or ceftazidime.

pKa, MM, and QM studies of mechanisms of beta-lactamases and penicillin-binding proteins: acylation step

The acylation step of the catalytic mechanism of beta-lactamases and penicillin-binding proteins (PBPs) has been studied with various approaches. The methods applied range from molecular dynamics (MD) simulations to multiple titration calculations using the Poisson-Boltzmann approach to quantum mechanical (QM) methods. The mechanism of class A beta-lactamases was investigated in the greatest detail. Most approaches support the critical role of Glu-166 and hydrolytic water in the acylation step of the enzymatic catalysis in class A beta-lactamases. The details of the catalytic mechanism have been revealed by the QM approach, which clearly pointed out the critical role of Glu-166 acting as a general base in the acylation step with preferred substrates. Lys-73 shuffles a proton abstracted by Glu-166 O(epsilon ) to the beta-lactam nitrogen through Ser-130 hydroxyl. This proton is transferred from O(gamma) of the catalytic Ser-70 through the bridging hydrolytic water to Glu-166 O(epsilon ). Then the hydrogen is simultaneously passed through S(N)2 inversion mechanism at Lys-73 N(zeta) to Ser-130 O(gamma), which loses its proton to the beta-lactam nitrogen. The protonation of beta-lactam nitrogen proceeds with an immediate ring opening and collapse of the first tetrahedral species into an acyl-enzyme intermediate. However, the studies that considered the effect of solvation lower the barrier for the pathway, which utilizes Lys-73 as a general base, thus creating a possibility of multiple mechanisms for the acylation step in the class A beta-lactamases. These findings help explain the exceptional efficiency of these enzymes. They emphasize an important role of Glu-166, Lys-73, and Ser-130 for enzymatic catalysis and shed light on details of the acylation step of class A beta-lactamase mechanism. The acylation step for class C beta-lactamases and six classes of PBPs were also considered with continuum solvent models and MD simulations.

EstB from Burkholderia gladioli: a novel esterase with a beta-lactamase fold reveals steric factors to discriminate between esterolytic and beta-lactam cleaving activity

Esterases form a diverse class of enzymes of largely unknown physiological role. Because many drugs and pesticides carry ester functions, the hydrolysis of such compounds forms at least one potential biological function. Carboxylesterases catalyze the hydrolysis of short chain aliphatic and aromatic carboxylic ester compounds. Esterases, D-alanyl-D-alanine-peptidases (DD-peptidases) and beta-lactamases can be grouped into two distinct classes of hydrolases with different folds and topologically unrelated catalytic residues, the one class comprising of esterases, the other one of beta-lactamases and DD-peptidases. The chemical reactivities of esters and beta-lactams towards hydrolysis are quite similar, which raises the question of which factors prevent esterases from displaying beta-lactamase activity and vice versa. Here we describe the crystal structure of EstB, an esterase isolated from Burkholderia gladioli. It shows the protein to belong to a novel class of esterases with homology to Penicillin binding proteins, notably DD-peptidase and class C beta-lactamases. Site-directed mutagenesis and the crystal structure of the complex with diisopropyl-fluorophosphate suggest Ser75 within the "beta-lactamase" Ser-x-x-Lys motif to act as catalytic nucleophile. Despite its structural homology to beta-lactamases, EstB shows no beta-lactamase activity. Although the nature and arrangement of active-site residues is very similar between EstB and homologous beta-lactamases, there are considerable differences in the shape of the active site tunnel. Modeling studies suggest steric factors to account for the enzyme's selectivity for ester hydrolysis versus beta-lactam cleavage.

Role of protein flexibility in enzymatic catalysis: quantum mechanical-molecular mechanical study of the deacylation reaction in class A beta-lactamases

We present a theoretical study of a mechanism for the hydrolysis of the acyl-enzyme complex formed by a class A beta-lactamase (TEM1) and an antibiotic (penicillanate), as a part of the process of antibiotic's inactivation by this type of enzymes. In the presented mechanism the carboxylate group of a particular residue (Glu166) activates a water molecule, accepting one of its protons, and afterward transfers this proton directly to the acylated serine residue (Ser70). In our study we employed a quantum mechanics (AM1)-molecular mechanics partition scheme (QM/MM) where all the atoms of the system were allowed to relax. For this purpose we used the GRACE procedure in which part of the system is used to define the Hessian matrix while the rest is relaxed at each step of the stationary structures search. By use of this computational scheme, the hydrolysis of the acyl-enzyme is described as a three-step process: The first step corresponds to the proton transfer from the hydrolytic water molecule to the carboxylate group of Glu166 and the subsequent formation of a tetrahedral adduct as a consequence of the attack of this activated water molecule to the carbonyl carbon atom of the beta-lactam. In the second step, the acyl-enzyme bond is broken, obtaining a negatively charged Ser70. In the last step this residue is protonated by means of a direct proton transfer from Glu166. The large mobility of Glu166, a residue that is placed in a Ohms-loop, is essential to facilitate this mechanism. The geometry of the acyl-enzyme complex shows a large distance between Glu166 and Ser70 and thus, if protein coordinates were kept frozen during the reaction path, it would be difficult to get a direct proton transfer between these two residues. This computational study shows how a flexible treatment suggests the feasibility of a mechanism that could have been discounted on the basis of crystallographic positions.

Insights into the molecular basis for the carbenicillinase activity of PSE-4 beta-lactamase from crystallographic and kinetic studies

PSE-4 is a class A beta-lactamase produced by strains of Pseudomonas aeruginosa and is highly active for the penicillin derivative carbenicillin. The crystal structure of the wild-type PSE-4 carbenicillinase has been determined to 1.95 A resolution by molecular replacement and represents the first structure of a carbenicillinase published to date. A superposition of the PSE-4 structure with that of TEM-1 shows a rms deviation of 1.3 A for 263 Calpha atoms. Most carbenicillinases are unique among class A beta-lactamases in that residue 234 is an arginine (ABL standard numbering scheme), while in all other class A enzymes this residue is a lysine. Kinetic characterization of a R234K PSE-4 mutant reveals a 50-fold reduction in k(cat)/K(m) and confirms the importance of Arg 234 for carbenicillinase activity. A comparison of the structure of the R234K mutant refined to 1.75 A resolution with the wild-type structure shows that Arg 234 stabilizes an alternate conformation of the Ser 130 side chain, not seen in other class A beta-lactamase structures. Our molecular modeling studies suggest that the position of a bound carbenicillin would be shifted relative to that of a bound benzylpenicillin in order to avoid a steric clash between the carbenicillin alpha-carboxylate group and the conserved side chain of Asn 170. The alternate conformation of the catalytic Ser 130 in wild-type PSE-4 may be involved in accommodating this shift in the bound substrate position.

Insights into class D beta-lactamases are revealed by the crystal structure of the OXA10 enzyme from Pseudomonas aeruginosa

BACKGROUND

beta-lactam antibiotic therapies are commonly challenged by the hydrolytic activities of beta-lactamases in bacteria. These enzymes have been grouped into four classes: A, B, C, and D. Class B beta-lactamases are zinc dependent, and enzymes of classes A, C, and D are transiently acylated on a serine residue in the course of the turnover chemistry. While class A and C beta-lactamases have been extensively characterized by biochemical and structural methods, class D enzymes remain the least studied despite their increasing importance in the clinic.

RESULTS

The crystal structure of the OXA10 class D beta-lactamase has been solved to 1.66 A resolution from a gold derivative and MAD phasing. This structure reveals that beta-lactamases from classes D and A, despite very poor sequence similarity, share a similar overall fold. An additional beta strand in OXA10 mediates the association into dimers characterized by analytical ultracentrifugation. Major differences are found when comparing the molecular details of the active site of this class D enzyme to the corresponding regions in class A and C beta-lactamases. In the native structure of the OXA10 enzyme solved to 1.8 A, Lys-70 is carbamylated.

CONCLUSIONS

Several features were revealed by this study: the dimeric structure of the OXA10 beta-lactamase, an extension of the substrate binding site which suggests that class D enzymes may bind other substrates beside beta-lactams, and carbamylation of the active site Lys-70 residue. The CO2-dependent activity of the OXA10 enzyme and the kinetic properties of the natural OXA17 mutant protein suggest possible relationships between carbamylation, inhibition of the enzyme by anions, and biphasic behavior of the enzyme.

The high resolution crystal structure for class A beta-lactamase PER-1 reveals the bases for its increase in breadth of activity

The treatment of infectious diseases by beta-lactam antibiotics is continuously challenged by the emergence and dissemination of new beta-lactamases. In most cases, the cephalosporinase activity of class A enzymes results from a few mutations in the TEM and SHV penicillinases. The PER-1 beta-lactamase was characterized as a class A enzyme displaying a cephalosporinase activity. This activity was, however, insensitive to the mutations of residues known to be critical for providing extended substrate profiles to TEM and SHV. The x-ray structure of the protein, solved at 1.9-A resolution, reveals that two of the most conserved features in class A beta-lactamases are not present in this enzyme: the fold of the Omega-loop and the cis conformation of the peptide bond between residues 166 and 167. The new fold of the Omega-loop and the insertion of four residues at the edge of strand S3 generate a broad cavity that may easily accommodate the bulky substituents of cephalosporin substrates. The trans conformation of the 166-167 bond is related to the presence of an aspartic acid at position 136. Selection of class A enzymes based on the occurrence of both Asp(136) and Asn(179) identifies a subgroup of enzymes with high sequence homology.

Protonation of the -lactam nitrogen is the trigger event in the catalytic action of class A -lactamases

The pH dependence of the pK(a) values of all ionizable groups and of the electrostatic potential at grid points corresponding to catalytically important atoms in the active site of TEM-1 beta-lactamase has been calculated by a mean-field approach for reaction intermediates modeled on the basis of energy minimized x-ray crystallographic coordinates. By estimating electrostatic contributions to the free energy changes accompanying the conversion of the free enzyme into the acylenzyme reaction intermediate, we found that acid-catalyzed protonation of the beta-lactam nitrogen is energetically favored as the initiating event, followed by base-catalyzed nucleophilic attack on the carbonyl carbon of the beta-lactam group. N-protonation is catalyzed through a hydrogen-bonded cluster involving the 2-carboxylate group of the substrate, the side chains of S130 and K234, and a solvent molecule. Nucleophilic attack on the carbonyl carbon is carried out by the side chain of S70 with proton abstraction catalyzed by a water molecule hydrogen-bonded to the side chain of E166. Stabilization of ion pairs in the active site through interactions with distant clusters of charged residues in the enzyme was concluded to be an important driving force of the catalytic mechanism.

Protonation of the beta-lactam nitrogen is the trigger event in the catalytic action of class A beta-lactamases

The pH dependence of the pK(a) values of all ionizable groups and of the electrostatic potential at grid points corresponding to catalytically important atoms in the active site of TEM-1 beta-lactamase has been calculated by a mean-field approach for reaction intermediates modeled on the basis of energy minimized x-ray crystallographic coordinates. By estimating electrostatic contributions to the free energy changes accompanying the conversion of the free enzyme into the acylenzyme reaction intermediate, we found that acid-catalyzed protonation of the beta-lactam nitrogen is energetically favored as the initiating event, followed by base-catalyzed nucleophilic attack on the carbonyl carbon of the beta-lactam group. N-protonation is catalyzed through a hydrogen-bonded cluster involving the 2-carboxylate group of the substrate, the side chains of S130 and K234, and a solvent molecule. Nucleophilic attack on the carbonyl carbon is carried out by the side chain of S70 with proton abstraction catalyzed by a water molecule hydrogen-bonded to the side chain of E166. Stabilization of ion pairs in the active site through interactions with distant clusters of charged residues in the enzyme was concluded to be an important driving force of the catalytic mechanism.

Selection of beta-lactamases and penicillin binding mutants from a library of phage displayed TEM-1 beta-lactamase randomly mutated in the active site omega-loop

A combinatorial library of mutants of the phage displayed TEM-1 lactamase was generated in the region encompassing residues 163 to 171 of the active site Omega-loop. Two in vitro selection protocols were designed to extract from the library phage-enzymes characterised by a fast acylation by benzyl-penicillin (PenG) to yield either stable or very unstable acyl-enzymes. The critical step of the selections was the kinetically controlled labelling of the phages by reaction with either a biotinylated penicillin derivative or a biotinylated penicillin sulfone, i.e. a beta-lactamase suicide substrate; the biotinylated phages were recovered by panning on immobilised streptavidin. As labelling with biotinylated suicide substrates tends to select enzymes that do not turnover, a counter-selection against penicillin binding mutants was introduced to extract the beta-lactamases. The selected phage-enzymes were characterised by sequencing to identify conserved residues and by kinetic analysis of the reaction with benzyl-penicillin. Several penicillin binding mutants, in which the essential Glu166 is replaced by Asn, were shown to be acylated very fast by PenG, the acylation being characterised by biphasic kinetics. These data are interpreted by a kinetic scheme in which the enzymes exist in two interconvertible conformations. The rate constant of the conformational change suggests that it involves an isomerisation of the peptide bond between residues 166 and 167 and controls a conformation of the Omega-loop compatible with fast acylation of the active site serine residue.

X-ray structure of the Asn276Asp variant of the Escherichia coli TEM-1 beta-lactamase: direct observation of electrostatic modulation in resistance to inactivation by clavulanic acid

The clinical use of beta-lactam antibiotics combined with beta-lactamase inactivators, such as clavulanate, has resulted in selection of beta-lactamases that are insensitive to inactivation by these molecules. Therefore, therapeutic combinations of an enzyme inactivator and a penicillin are harmless for bacteria harboring such an enzyme. The TEM beta-lactamase variants are the most frequently encountered enzymes of this type, and presently, 20 variants are designated as inhibitor-resistant TEM ("IRT") enzymes. Three mutations appear to account for the phenotype of the majority of IRT enzymes, one of them being the Asn276Asp substitution. In this study, we have characterized the kinetic properties of the inhibition process of the wild-type TEM-1 beta-lactamase and of its Asn276Asp variant with the three clinically used inactivators, clavulanic acid (clavulanate), sulbactam, and tazobactam, and we report the X-ray structure for the mutant variant at 2.3 A resolution. The changes in kinetic parameters for the interactions of the inhibitors with the wild-type and the mutant enzymes were more pronounced for clavulanate, and relatively inconsequential for sulbactam and tazobactam. The structure of the Asn276Asp mutant enzyme revealed a significant movement of Asp276 and the formation of a salt bridge of its side chain with the guanidinium group of Arg244, the counterion of the inhibitor carboxylate. A water molecule critical for the inactivation chemistry by clavulanate, which is observed in the wild-type enzyme structure, is not present in the crystal structure of the mutant variant. Such structural changes favor the turnover process over the inactivation chemistry for clavulanate, with profound phenotypic consequences. The report herein represents the best studied example of inhibitor-resistant beta-lactamases.

pKa calculations for class A beta-lactamases: influence of substrate binding

Beta-Lactamases are responsible for bacterial resistance to beta-lactams and are thus of major clinical importance. However, the identity of the general base involved in their mechanism of action is still unclear. Two candidate residues, Glu166 and Lys73, have been proposed to fulfill this role. Previous studies support the proposal that Glu166 acts during the deacylation, but there is no consensus on the possible role of this residue in the acylation step. Recent experimental data and theoretical considerations indicate that Lys73 is protonated in the free beta-lactamases, showing that this residue is unlikely to act as a proton abstractor. On the other hand, it has been proposed that the pKa of Lys73 would be dramatically reduced upon substrate binding and would thus be able to act as a base. To check this hypothesis, we performed continuum electrostatic calculations for five wild-type and three beta-lactamase mutants to estimate the pKa of Lys73 in the presence of substrates, both in the Henri-Michaelis complex and in the tetrahedral intermediate. In all cases, the pKa of Lys73 was computed to be above 10, showing that it is unlikely to act as a proton abstractor, even when a beta-lactam substrate is bound in the enzyme active site. The pKa of Lys234 is also raised in the tetrahedral intermediate, thus confirming a probable role of this residue in the stabilization of the tetrahedral intermediate. The influence of the beta-lactam carboxylate on the pKa values of the active-site lysines is also discussed.

Catalytic properties of class A beta-lactamases: efficiency and diversity

beta-Lactamases are the main cause of bacterial resistance to penicillins, cephalosporins and related beta-lactam compounds. These enzymes inactivate the antibiotics by hydrolysing the amide bond of the beta-lactam ring. Class A beta-lactamases are the most widespread enzymes and are responsible for numerous failures in the treatment of infectious diseases. The introduction of new beta-lactam compounds, which are meant to be 'beta-lactamase-stable' or beta-lactamase inhibitors, is thus continuously challenged either by point mutations in the ubiquitous TEM and SHV plasmid-borne beta-lactamase genes or by the acquisition of new genes coding for beta-lactamases with different catalytic properties. On the basis of the X-ray crystallography structures of several class A beta-lactamases, including that of the clinically relevant TEM-1 enzyme, it has become possible to analyse how particular structural changes in the enzyme structures might modify their catalytic properties. However, despite the many available kinetic, structural and mutagenesis data, the factors explaining the diversity of the specificity profiles of class A beta-lactamases and their amazing catalytic efficiency have not been thoroughly elucidated. The detailed understanding of these phenomena constitutes the cornerstone for the design of future generations of antibiotics.

Crystal structure of an acylation transition-state analog of the TEM-1 beta-lactamase. Mechanistic implications for class A beta-lactamases

The crystal structure of a phosphonate complex of the class A TEM-1 beta-lactamase has been determined to a resolution of 2.0 A. The phosphonate appears stoichiometrically at the active site, bound covalently to Ser70Ogamma, with one phosphonyl oxygen in the oxyanion hole. Although the overall structure is very similar to that of the native enzyme (rms difference 0.37 A for all heavy atoms), changes have occurred in the position of active site functional groups. The active site is also not in the conformation observed in the complex of another class A beta-lactamase, that of Staphylococcus aureus PC1, with the same phosphonate [Chen, C. C. H., et al. (1993) J. Mol. Biol. 234,165-178]. Both phosphonate structures, however, can be seen to represent models of acylation transition-states since in each the deacylating water molecule appears firmly bound to the Glu166 carboxylate group. The major difference between the structures lies in the positioning of Lys73Nzeta and Ser130Ogamma. In the S. aureus structure, the closest interaction of these functional groups is between Lys73Nzeta and Ser70Ogamma (2.8 A), while in the TEM-1 structure it is between Ser130Ogamma and the second phosphonyl oxygen of the bound inhibitor (2.8 A). The former structure therefore may resemble a transition state for formation of the tetrahedral species in acylation by nucleophilic attack on the substrate, where Lys73Nzeta presumably catalyzes the reaction as a general base. The TEM-1 structure can then be seen as an analogue of the transition state for breakdown of the tetrahedral species, where Ser130Ogamma is acting as a general acid, assisting the departure of the leaving group. The class A beta-lactamase crystal structures now available lead to a self-consistent proposal for a mechanism of catalysis by these enzymes.

pKa calculations for class A beta-lactamases: methodological and mechanistic implications

Beta-lactamases are responsible for resistance to penicillins and related beta-lactam compounds. Despite numerous studies, the identity of the general base involved in the acylation step is still unclear. It has been proposed, on the basis of a previous pKa calculation and analysis of structural data, that the unprotonated Lys73 in the active site could act as the general base. Using a continuum electrostatic model with an improved treatment of the multiple titration site problem, we calculated the pKa values of all titratable residues in the substrate-free TEM-1 and Bacillus licheniformis class A beta-lactamases. The pKa of Lys73 in both enzymes was computed to be above 10, in good agreement with recent experimental data on the TEM-1 beta-lactamase, but inconsistent with the proposal that Lys73 acts as the general base. Even when the closest titratable residue, Glu166, is mutated to a neutral residue, the predicted downward shift of the pKa of Lys73 shows that it is unlikely to act as a proton abstractor in either enzyme. These results support a mechanism in which the proton of the active Ser70 is transferred to the carboxylate group of Glu166.

Structure-based design of a potent transition state analogue for TEM-1 beta-lactamase

The structure of the plasmid-mediated beta-lactamase TEM-1 has been solved in complex with a designed boronic acid inhibitor (1R)-1-acetamido-2-(3-carboxyphenyl)ethane boronic acid at 1.7 A resolution. The boronate inhibitor was designed based on the crystallographic coordinates of the acyl-enzyme intermediate of TEM-1 bound to the substrate penicillin G. The boronate-TEM-1 complex is highly ordered and defines a novel transition state analogue of the deacylation step in the beta-lactamase reaction pathway. The design principles of this highly effective inhibitor (Ki = 110 nM) and the resulting structural and mechanistic implications are presented.

Structural basis of extended spectrum TEM beta-lactamases. Crystallographic, kinetic, and mass spectrometric investigations of enzyme mutants

The E166Y and the E166Y/R164S TEM-1 beta-lactamase mutant enzymes display extended spectrum substrate specificities. Electrospray mass spectrometry demonstrates that, with penicillin G as substrate, the rate-limiting step in catalysis is the hydrolysis of the E166Y acyl-enzyme complex. Comparison of the 1.8-A resolution x-ray structures of the wild-type and of the E166Y mutant enzymes shows that the binding of cephalosporin substrates is improved, in the mutant enzyme, by the enlargement of the substrate binding site. This enlargement is due to the rigid body displacement of 60 residues driven by the movement of the omega-loop. These structural observations strongly suggest that the link between the position of the omega-loop and that of helix H5, plays a central role in the structural events leading to extended spectrum TEM-related enzymes. The increased omega-loop flexibility caused by the R164S mutation, which is found in several natural mutant TEM enzymes, may lead to similar structural effects. Comparisons of the kinetic data of the E166Y, E166Y/R164S, and R164S mutant enzymes supports this hypothesis.

The catalytic mechanism of beta-lactamases: NMR titration of an active-site lysine residue of the TEM-1 enzyme

Beta-Lactamases are widespread in the bacterial world, where they are responsible for resistance to penicillins, cephalosporins, and related compounds, currently the most widely used antibacterial agents. Detailed structural and mechanistic understanding of these enzymes can be expected to guide the design of new antibacterial compounds resistant to their action. A number of high-resolution structures are available for class A beta-lactamases, whose catalytic mechanism involves the acylation of a serine residue at the active site. The identity of the general base which participates in the activation of this serine residue during catalysis has been the subject of controversy, both a lysine residue and a glutamic acid residue having been proposed as candidates for this role. We have used the pH dependence of chemical modification of epsilon-amino groups by 2,4,6,-trinitrobenzenesulfonate and the pH dependence of the epsilon-methylene 1H and 13C chemical shifts (in enzyme selectively labeled with [epsilon-13C]lysine) to estimate the pKa of the relevant lysine residue, lysine-73, of TEM-1 beta-lactamase. Both methods show that the pKa of this residue is > 10, making it very unlikely that this residue could act as a proton acceptor in catalysis. An alternative mechanism in which this role is performed by glutamate-166 through an intervening water molecule is described.