Mechanisms of antibiotic resistance: QM/MM modeling of the acylation reaction of a class A beta-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.     

An ultrahigh resolution structure of TEM-1 beta-lactamase suggests a role for Glu166 as the general base in acylation

Although TEM-1 beta-lactamase is among the best studied enzymes, its acylation mechanism remains controversial. To investigate this problem, the structure of TEM-1 in complex with an acylation transition-state analogue was determined at ultrahigh resolution (0.85 A) by X-ray crystallography. The quality of the data was such as to allow for refinement to an R-factor of 9.1% and an R(free) of 11.2%. In the resulting structure, the electron density features were clear enough to differentiate between single and double bonds in carboxylate groups, to identify multiple conformations that are occupied by residues and loops, and to assign 70% of the protons in the protein. Unexpectedly, even at pH 8.0 where the protein was crystallized, the active site residue Glu166 is clearly protonated. This supports the hypothesis that Glu166 is the general base in the acylation half of the reaction cycle. This structure suggests that Glu166 acts through the catalytic water to activate Ser70 for nucleophilic attack on the beta-lactam ring of the substrate. The hydrolytic mechanism of class A beta-lactamases, such as TEM-1, appears to be symmetrical, as are the serine proteases. Apart from its mechanistic implications, this atomic resolution structure affords an unusually detailed view of the structure, dynamics, and hydrogen-bonding networks of TEM-1, which may be useful for the design of inhibitors against this key antibiotic resistance target.     

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.     

Ab initio QM/MM study of class A beta-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.     

Investigation of the mechanism of the cell wall DD-carboxypeptidase reaction of penicillin-binding protein 5 of Escherichia coli by quantum mechanics/molecular mechanics calculations

Penicillin-binding protein 5 (PBP 5) of Escherichia coli hydrolyzes the terminal D-Ala-D-Ala peptide bond of the stem peptides of the cell wall peptidoglycan. The mechanism of PBP 5 catalysis of amide bond hydrolysis is initial acylation of an active site serine by the peptide substrate, followed by hydrolytic deacylation of this acyl-enzyme intermediate to complete the turnover. The microscopic events of both the acylation and deacylation half-reactions have not been studied. This absence is addressed here by the use of explicit-solvent molecular dynamics simulations and ONIOM quantum mechanics/molecular mechanics (QM/MM) calculations. The potential-energy surface for the acylation reaction, based on MP2/6-31+G(d) calculations, reveals that Lys47 acts as the general base for proton abstraction from Ser44 in the serine acylation step. A discrete potential-energy minimum for the tetrahedral species is not found. The absence of such a minimum implies a conformational change in the transition state, concomitant with serine addition to the amide carbonyl, so as to enable the nitrogen atom of the scissile bond to accept the proton that is necessary for progression to the acyl-enzyme intermediate. Molecular dynamics simulations indicate that transiently protonated Lys47 is the proton donor in tetrahedral intermediate collapse to the acyl-enzyme species. Two pathways for this proton transfer are observed. One is the direct migration of a proton from Lys47. The second pathway is proton transfer via an intermediary water molecule. Although the energy barriers for the two pathways are similar, more conformers sample the latter pathway. The same water molecule that mediates the Lys47 proton transfer to the nitrogen of the departing D-Ala is well positioned, with respect to the Lys47 amine, to act as the hydrolytic water in the deacylation step. Deacylation occurs with the formation of a tetrahedral intermediate over a 24 kcal x mol(-1) barrier. This barrier is approximately 2 kcal x mol(-1) greater than the barrier (22 kcal x mol(-1)) for the formation of the tetrahedral species in acylation. The potential-energy surface for the collapse of the deacylation tetrahedral species gives a 24 kcal x mol(-1) higher energy species for the product, signifying that the complex would readily reorganize and pave the way for the expulsion of the product of the reaction from the active site and the regeneration of the catalyst. These computational data dovetail with the knowledge on the reaction from experimental approaches.     

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.     

Modeling of the mechanism of hydrolysis of succinylcholine in the active site of native and modified (Asp70Gly) human butyrylcholinesterase

The quantum mechanics/molecular mechanics method was used to model the mechanism of hydrolysis of succinylcholine in the active site of native butyrylcholinesterase and the Asp70Gly single mutant enzyme. The energy diagram for the reaction pathway of the acylation was constructed, and the effect of the mutation on the reaction mechanism and energy characteristics of the process are discussed.     

Catalytic mechanism of class A beta-lactamase. I. The role of Glu166 and Serl30 in the deacylation reaction

The tetrahedral intermediate formation process, which is the first step in the deacylation reaction by class A beta-lactamase, was investigated by the ab initio molecular orbital method. In this study, benzyl penicillin was used as the substrate. From the results of our molecular dynamics study of the structure of beta-lactam antibiotics-beta-lactamase complex, the substrate, Ser70, Lys73, Ser130, Glu166 and a water molecule for the deacylation reaction were considered for construction of a model for calculation. The calculation results indicated that Glu166 plays a role in holding a water molecule, which is necessary for the deacylation reaction, and that the hydrogen bond network among Lys73Nzeta, Ser130Ogamma, and the carboxyl group of the beta-lactam antibiotics was formed by the uptake of beta-lactam antibiotics by beta-lactamase. The activation energy for this reaction was 33.3 kcal/mol, and it is very likely that the reaction occurred at body temperature. Subsequent calculation results obtained by using the model excluding Ser130 and the carboxyl group of the substrate indicated that the activation energy for this reaction was 40.8 kcal/mol, which is 7.5 kcal/mol higher than that of the previous reaction. It was found that the hydrogen bond network plays an important role in decreasing the activation energy for the tetrahedral intermediate formation reaction. Lys73Nzeta, which is located at the edge of the hydrogen bond network, played a role in forming a hydrogen bond with Glu166Oepsilon in order to help the deacylation reaction. The role of amino acid residues around the active site of class A beta-lactamase was also discussed.     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

On the possible reaction pathway for the acylation of AChE-catalyzed hydrolysis of ACh: Semiempirical quantum chemical study

The acylation process in the acetylcholinesterase (AChE)-catalyzed hydrolysis of the neurotransmitter, acetylcholine (ACh), has been determined with the semiempirical quantum chemical calculation method AM1 using the model molecules extracted from the X-ray crystal structure of Torpedo Californica AChE. For the sake of identifying the microfeatures of the mechanism of this reaction, two types of possible mechanisms, stepwise mechanism and cooperative mechanisms, were proposed and studied with AM1 methods. All the model molecules for the possible reactants, intermediates, transition states, and products in the reaction pathways of the two mechanisms were obtained. Energy profiles, the structural properties of the transition states, indicate that the acylation of AChE-catalyzed hydrolysis of ACh adopts the cooperative mechanism, i.e., the proton transfer from Ser220 of AChE to His440 occurs simultaneously with the nucleophilic attack of Ser200 to the carbonyl carbon atom of ACh. This result is in agreement with the kinetic data and the secondary isotope effects of AChE-catalyzed reactions. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 70: 515–525, 1998     

A quantum mechanics study on the reaction mechanism of chalcone formation from p-coumaroyl-CoA and malonyl-CoA catalyzed by chalcone synthase

Chalcone synthase catalyzes formation of phenylpropanoid chalcone from one p-coumaroyl-coenzyme A (CoA) and three malonyl-CoA molecules. In order to elucidate structural and energetic features of the reaction mechanism, we performed the quantum mechanics calculations and obtained the following results. In loading step, only a tetrahedral intermediate is located without transition state (TS). Our results indicate that His303 acts as a H₃₁ donor, but not a hydrogen bond donor, to stabilize the intermediate formation. In decarboxylation step, the reaction proceeds via a TS and is sensitive to the environment. In elongation step, a tetrahedral TS is located. All of the results above support the reaction mechanism and further complement the proposal of Noel JP et al.     

A theoretical approach for the acylation/deacylation mechanisms of avibactam in the reversible inhibition of KPC-2

Klebsiella pneumoniae carbapenemase (KPC-2) is the most commonly encountered class A β-lactamase variant worldwide, which confer high-level resistance to most available antibiotics. In this article we address the issue by a combined approach involving molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations. The study contributes to improve the understanding, at molecular level, of the acylation and deacylation stages of avibactam involved in the inhibition of KPC-2. The results show that both mechanisms, acylation and deacylation, the reaction occur via the formation of a tetrahedral intermediate. The formation of this intermediate corresponds to the rate limiting stage. The activation barriers are 19.5 kcal/mol and 23.0 kcal/mol for the acylation and deacylation stages, respectively. The associated rate constants calculated, using the Eyring equation, are 1.2?×?10^?1 and 3.9?×?10^?4 (s^?1). These values allow estimating a value of 3.3?×?10^?3 for the inhibition constant, in good agreement with the experimental value.

     

Ab initio model study on acetylcholinesterase catalysis: potential energy surfaces of the proton transfer reactions

Ab initio molecular orbital (MO) and hybrid density functional theory (DFT) calculations have been applied to the initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE), which is the nucleophiric addition of Ser200 in catalytic triads to a neurotransmitter acetylcholine (ACh). We focus our attention mainly on the effects of oxyanion hole and Glu327 on the potential energy surfaces (PESs) for the proton transfer reactions in the catalytic triad Ser200-His440-Glu327. The activation barrier for the addition reaction of Ser200 to ACh was calculated to be 23.4 kcal/mol at the B3LYP/6-31G(d)//HF/3-21G(d) level of theory. The barrier height under the existence of oxyanion hole, namely, Ser200-His440-Glu327-ACh-(oxyanion hole) system, decreased significantly to 14.2 kcal/mol, which is in reasonable agreement with recent experimental value (12.0 kcal/mol). Removal of Glu327 from the catalytic triad caused destabilization of both energy of transition state for the reaction and tetrahedral intermediate (product). PESs calculated for the proton transfer reactions showed that the first proton transfer process is the most important in the stabilization of tetrahedral intermediate complex. The mechanism of addition reaction of ACh was discussed on the basis of theoretical results.     

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.     

Regiodivergent C−H Acylation of Arenes by Switching from Ionic- to Radical-Type Chemistry Using NHC Catalysis

The Friedel–Crafts acylation reaction, which belongs to the class of electrophilic aromatic substitutions is a highly valuable and versatile reaction in synthesis. Regioselectivity is predictable and determined by electronic as well as steric factors of the (hetero)arene substrate. Herein, a radical approach for the acylation of arenes and heteroarenes is presented. C−H acylation is achieved through mild cooperative photoredox/NHC radical catalysis with the cross-coupling of an arene radical cation with an NHC-bound ketyl radical as a key step. As compared to the classical Friedel–Crafts acylation, a regiodivergent outcome is observed upon switching from the ionic to the radical mode. In these divergent reactions, aroyl fluorides act as the acylation reagents in both the ionic as well as the radical process.     

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α -carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.     

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

Kinetics and modeling of (R,S)‐1‐phenylethanol acylation over lipase

The kinetics of the acylation of (R,S)‐1‐phenylethanol was investigated using lipase as a catalyst. The main parameters were temperature, reaction atmosphere, different acyl donors, and different amounts of acyl donor as well as the presence of some additives in the reaction mixture. The initial reaction rate increased with increasing temperature and with a decreasing amount of an acyl donor. The activated esters, such as isopropenyl‐ and vinyl acetate, exhibited very high acylation rates for R‐1‐phenylethanol, whereas low rates were obtained with ethyl acetate and 2‐methoxyethyl acetate. The addition of water and acetophenone decreased the acylation rate. A kinetic model was developed based on a sequential step mechanism, in which enzyme was reacting in the first step with an acyl donor followed by the reaction of a modified enzyme complex with the reactant, R‐1‐phenylethanol. Comparison with experimental data obtained at different temperatures allowed simplification of this model, leading to a kinetic equation with just one apparent parameter. The influence of the amount of acyl donor, ethyl acetate, could be quantitatively described by taking into account the competitive inhibition of the ethanol produced. The rate constants and apparent activation energy for experiments performed under different temperatures and the amounts of acylation agent were determined. The apparent activation energy was 24.5 kJ/mol. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 629–639, 2010     

脂肪酰胺水解酶催化三联体膦酰化失活的机理：一个模型体系的理论研究

甲基花生四烯基氟代膦酸酯(MAFP)是脂肪酰胺水解酶(FAAH)的一个抑制剂. FAAH的丝氨酸241(Ser241)-丝氨酸217(Ser217)-赖氨酸142(Lys142)催化三联体被MAFP膦酰化后将导致FAAH失活. 本文采用B3LYP/6-311G(d,p)和MP2/6-311G(d,p)方法及一个简化的计算模型体系对这个膦酰化抑制反应进行理论研究. 考虑了两种反应途径. Path A涉及FAAH的催化三联体的所有残基, 是一个分步的加成-消除过程, 形成两性离子的三角双锥中间体, 其中第一步反应是决速步骤. 在这个反应途径中, Ser217和Lys142对亲核试剂Ser241起到碱催化活化的作用, 而Ser217充作Lys142和Ser241之间的桥梁. 此外, 溶剂中的一个水分子作为Lys142和MAFP间的“氢桥”具有关键的作用, 通过给出和接收质子促进了长距离的质子转移. Path B是催化三联体中的残基Lys142被突变为丙氨酸以后的膦酰化反应, 也是一个分步过程. 水的本体溶剂效应通过极化连续介质模型(PCM)估算. 计算结果显示膦酰化反应的Path A是优势途径, 在水溶液中其决速步骤的活化能垒为64.9 kJ·mol-1. FAAH催化三联体中残基Lys142的变异会降低膦酰化反应的速率, 这与实验结果相一致.     

Structure, function, and inhibition along the reaction coordinate of CTX-M beta-lactamases

CTX-M enzymes are an emerging group of extended spectrum beta-lactamases (ESBLs) that hydrolyze not only the penicillins but also the first-, second-, and third-generation cephalosporins. Although they have become the most frequently observed ESBLs in certain areas, there are few effective inhibitors and relatively little is known about their detailed mechanism. Here we describe the X-ray crystal structures of CTX-M enzymes in complex with different transition-state analogues and beta-lactam inhibitors, representing the enzyme as it progresses from its acylation transition state to its acyl enzyme complex to the deacylation transition state. As the enzyme moves along this reaction coordinate, two key catalytic residues, Lys73 and Glu166, change conformations, tracking the state of the reaction. Unexpectedly, the acyl enzyme complex with the beta-lactam inhibitor cefoxitin still has the catalytic water bound; this water had been predicted to be displaced by the unusual 7alpha-methoxy of the inhibitor. Instead, the 7alpha-group appears to inhibit by preventing the formation of the deacylation transition state through steric hindrance. From an inhibitor design standpoint, we note that the best of the reversible inhibitors, a ceftazidime-like boronic acid compound, binds to CTX-M-16 with a K(i) value of 4 nM. When used together in cell culture, this inhibitor reversed cefotaxime resistance in CTX-M-producing bacteria. The structure of its complex with CTX-M enzyme and the structural view of the reaction coordinate described here provide templates for inhibitor design and intervention to combat this family of antibiotic resistance enzymes.