Antibiotic recognition by binuclear metallo-beta-lactamases revealed by X-ray crystallography

Metallo-beta-lactamases are zinc-dependent enzymes responsible for resistance to beta-lactam antibiotics in a variety of host bacteria, usually Gram-negative species that act as opportunist pathogens. They hydrolyze all classes of beta-lactam antibiotics, including carbapenems, and escape the action of available beta-lactamase inhibitors. Efforts to develop effective inhibitors have been hampered by the lack of structural information regarding how these enzymes recognize and turn over beta-lactam substrates. We report here the crystal structure of the Stenotrophomonas maltophilia L1 enzyme in complex with the hydrolysis product of the 7alpha-methoxyoxacephem, moxalactam. The on-enzyme complex is a 3'-exo-methylene species generated by elimination of the 1-methyltetrazolyl-5-thiolate anion from the 3'-methyl group. Moxalactam binding to L1 involves direct interaction of the two active site zinc ions with the beta-lactam amide and C4 carboxylate, groups that are common to all beta-lactam substrates. The 7beta-[(4-hydroxyphenyl)malonyl]-amino substituent makes limited hydrophobic and hydrogen bonding contacts with the active site groove. The mode of binding provides strong evidence that a water molecule situated between the two metal ions is the most likely nucleophile in the hydrolytic reaction. These data suggest a reaction mechanism for metallo-beta-lactamases in which both metal ions contribute to catalysis by activating the bridging water/hydroxide nucleophile, polarizing the substrate amide bond for attack and stabilizing anionic nitrogen intermediates. The structure illustrates how a binuclear zinc site confers upon metallo-beta-lactamases the ability both to recognize and efficiently hydrolyze a wide variety of beta-lactam substrates.     

Hybrid QM/MM and DFT investigations of the catalytic mechanism and inhibition of the dinuclear zinc metallo-beta-lactamase CcrA from Bacteroides fragilis

Based on hybrid QM/MM molecular dynamics simulation and density functional theoretical (DFT) calculations, we investigate the mechanistic and energetic features of the catalytic action of dizinc metallo-beta-lactamase CcrA from Bacteroides fragilis. The 200 ps QM/MM simulation of the CcrA enzyme in complex with nitrocefin shows that the substrate beta-lactam moiety is directed toward the active site dizinc center through the interactions of aminocarbonyl and carboxylate groups with the two active site zinc ions and the two conserved residues, Lys167 and Asn176. From the determination of the potential energy profile of a relevant enzymatic reaction model, it is found that the nucleophilic displacement reaction step proceeds with a low-barrier height, leading to the formation of an energetically favored reaction intermediate. The results also show that the high catalytic activity of the CcrA enzyme stems from a simultaneous operation of three catalytic components: activation of the bridging hydroxide nucleophile by zinc-coordinated Asp86; polarization of the substrate aminocarbonyl group by the first zinc ion; stabilization of the negative charge developed on the departing amide nitrogen by the second zinc ion. Consistent with the previous experimental finding that the proton-transfer reaction step is rate-limiting, the activation energy of the second step is found to be 1.6 kcal/mol higher than that of the first step. Finally, through an examination of the structural and energetic features of binding of a thiazolidinecarboxylic acid inhibitor to the active site dizinc center, a two-step inhibition mechanism involving a protonation-induced ligand exchange reaction is proposed for the inhibitory action of a tight-binding inhibitor possessing a thiol group.     

New insights into the mechanism of the Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase from S. enterica: a theoretical study

The reaction pathway of Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase (DHQD) from S. enterica has been studied by performing molecular dynamics (MD) simulations and density functional theory (DFT) calculations and the corresponding potential energy profile has also been identified. On the basis of the results, the catalytic hydrolysis process for the wild-type enzyme consists of three major reaction steps, including nucleophilic attack on the carbon atom involved in the carbon-nitrogen double bond of the Schiff base intermediate by a water molecule, deprotonation of the His143 residue, and dissociation between the product and the Lys170 residue of the enzyme. The remarkable difference between this and the previously proposed reaction mechanism is that the second step here, absent in the previously proposed reaction mechanism, plays an important role in facilitating the reaction through a key proton transfer by the His143 residue, resulting in a lower energy barrier. Comparison with our recently reported results on the Schiff base formation and dehydration processes clearly shows that the Schiff base hydrolysis is rate-determining in the overall reaction catalyzed by type I DHQD, consistent with the experimental prediction, and the calculated energy barrier of ～16.0 kcal mol(-1) is in good agreement with the experimentally derived activation free energy of ～14.3 kcal mol(-1). When the imidazole group of His143 residue is missing, the Schiff base hydrolysis is initiated by a hydroxide ion in the solution, rather than a water molecule, and both the reaction mechanism and the kinetics of Schiff base hydrolysis have been remarkably changed, clearly elucidating the catalytic role of the His143 residue in the reaction. The new mechanistic insights obtained here will be valuable for the rational design of high-activity inhibitors of type I DHQD as non-toxic antimicrobials, anti-fungals, and herbicides.     

Role of zinc content on the catalytic efficiency of B1 metallo beta-lactamases

Metallo beta-lactamases (MbetaL) are enzymes naturally evolved by bacterial strains under the evolutionary pressure of beta-lactam antibiotic clinical use. They have a broad substrate spectrum and are resistant to all the clinically useful inhibitors, representing a potential risk of infection if massively disseminated. The MbetaL scaffold is designed to accommodate one or two zinc ions able to activate a nucleophilic hydroxide for the hydrolysis of the beta-lactam ring. The role of zinc content on the binding and reactive mechanism of action has been the subject of debate and still remains an open issue despite the large amount of data acquired. We report herein a study of the reaction pathway for binuclear CcrA from Bacteroides fragilis using density functional theory based quantum mechanics-molecular mechanics dynamical modeling. CcrA is the prototypical binuclear enzyme belonging to the B1 MbetaL family, which includes several harmful chromosomally encoded and transferable enzymes. The involvement of a second zinc ion in the catalytic mechanism lowers the energetic barrier for beta-lactam hydrolysis, preserving the essential binding features found in mononuclear B1 enzymes (BcII from Bacillus cereus) while providing a more efficient single-step mechanism. Overall, this study suggests that uptake of a second equivalent zinc ion is evolutionary favored.     

Theoretical Studies on the Binding Mode and Reaction Mechanism of TLP Hydrolase kpHIUH

In this work, we have investigated the binding conformations of the substrate in the active site of 5-HIU hydrolase kpHIUH and its catalytic hydrolysis mechanism. Docking calculations revealed that the substrate adopts a conformation in the active site with its molecular plane laying parallel to the binding interface of the protein dimer of kpHIUH, in which His7 and His92 are located adjacent to the hydrolysis site C6 and have hydrogen bond interactions with the lytic water. Based on this binding conformation, density functional theory calculations indicated that the optimal catalytic mechanism consists of two stages: (1) the lytic water molecule is deprotonated by His92 and carries out nucleophilic attack on C6=O of 5-HIU, resulting in an oxyanion intermediate; (2) by accepting a proton transferred from His92, C6–N5 bond is cleaved to completes the catalytic cycle. The roles of His7, His92, Ser108 and Arg49 in the catalytic reaction were revealed and discussed in detail.     

Antibiotic resistance: mono- and dinuclear zinc complexes as metallo-beta-lactamase mimics

Biomimetic systems containing one or two zinc(II) ions supported by phenolate ligands were developed as functional mimics of metallo-beta-lactamase. These complexes were shown to catalytically hydrolyze beta-lactam substrates, such as oxacillin and penicillin G. The dinuclear zinc complex 1, which has a coordinated water molecule, exhibits high beta-lactamase activity, whereas the dinuclear zinc complex 2, which has no water molecules, but labile chloride ligands, shows a much lower activity. The high beta-lactamase activity of complex 1 can be ascribed to the presence of a zinc-bound water molecule that is activated by being hydrogen bonded to acetate substituents. The kinetics of the hydrolysis of oxacillin by complex 1 and the effect of pH on the reaction rates are reported in detail. In addition, the kinetic parameters obtained for the synthetic analogues are compared with those of the natural metallo-beta-lactamase from Bacillus cereus (BcII). To understand the role of the second metal ion in hydrolysis, the syntheses and catalytic activities of two mononuclear complexes (3 and 4) that include coordinated water molecules are described. Interestingly, the mononuclear zinc complexes 3 and 4 also exhibit high activity, supporting the assumption that the second zinc ion is not crucial for the beta-lactamase activity.     

Computational modeling of the catalytic mechanism of human placental alkaline phosphatase (PLAP)

Alkaline phosphatases (APs) catalyze the hydrolysis and transphosphorylation of phosphate monoesters. Quantum mechanical, molecular dynamics, and molecular docking techniques were applied to computationally model the catalytic mechanism of human placental AP (PLAP). Kinetic and thermodynamic evaluations were performed for each reaction step. The functional significances of the more important residues within the active site were analyzed. The role of the metal ion at the metal binding site M3 was also examined. The calculated activation and reaction energy and free energy values obtained suggested the nucleophilic attack of the Ser92 alkoxide on the phosphorus atom of the substrate would be the rate-limiting step of the catalytic hydrolysis of alkyl phosphate monoesters by PLAP. The reactivities of the wild-type M3-Mg enzyme and the M3-Zn protein were compared, and the main difference observed was a change in the coordination number of the M3 metal for the M3-Zn enzyme. This modification in the active site structure lowered the free energy profile for the second chemical step of the catalytic mechanism (hydrolysis of the covalent phosphoserine intermediate). Consequently, a greater stabilization of the phosphoseryl moiety resulted in a small increment in the activation free energy of the phosphoserine hydrolysis reaction. These computational results suggest that the activation of APs by magnesium at the M3 site is caused by the preference of Mg(2+) for octahedral coordination, which structurally stabilizes the active site into a catalytically most active conformation. The present theoretical results are in good agreement with previously reported experimental studies.     

Antibiotic deactivation by a dizinc beta-lactamase: mechanistic insights from QM/MM and DFT studies

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods and density functional theory (DFT) were used to investigate the initial ring-opening step in the hydrolysis of moxalactam catalyzed by the dizinc L1 beta-lactamase from Stenotrophomonas maltophilia. Anchored at the enzyme active site via direct metal binding as suggested by a recent X-ray structure of an enzyme-product complex (Spencer, J.; et al. J. Am. Chem. Soc. 2005, 127, 14439), the substrate is well aligned with the nucleophilic hydroxide that bridges the two zinc ions. Both QM/MM and DFT results indicate that the addition of the hydroxide nucleophile to the carbonyl carbon in the substrate lactam ring leads to a metastable intermediate via a dominant nucleophilic addition barrier. The potential of mean force obtained by SCC-DFTB/MM simulations and corrected by DFT/MM calculations yields a reaction free energy barrier of 23.5 kcal/mol, in reasonable agreement with the experimental value of 18.5 kcal/mol derived from kcat of 0.15 s(-1). It is further shown that zinc-bound Asp120 plays an important role in aligning the nucleophile, but accepts the hydroxide proton only after the nucleophilic addition. The two zinc ions are found to participate intimately in the catalysis, consistent with the proposed mechanism. In particular, the Zn(1) ion is likely to serve as an "oxyanion hole" in stabilizing the carbonyl oxygen, while the Zn(2) ion acts as an electrophilic catalyst to stabilize the anionic nitrogen leaving group.     

Insights into the reaction mechanism of soluble epoxide hydrolase from theoretical active site mutants

Density functional theory calculations of active site mutants are used to gain insights into the reaction mechanism of the soluble epoxide hydrolases (sEHs). The quantum chemical model is based on the X-ray crystal structure of the human soluble epoxide hydrolase. The role of two conserved active site tyrosines is explored through in silico single and double mutations to phenylalanine. Full potential energy curves for hydrolysis of (1S,2S)-beta-methylstyrene oxide are presented. The results indicate that the two active site tyrosines act in concert to lower the activation barrier for the alkylation step. For the wild-type and three different tyrosine mutant models, the regioselectivity of epoxide opening is compared for the substrates (1S,2S)-beta-methylstyrene oxide and (S)-styrene oxide. An additional part of our study focuses on the importance of the catalytic histidine for the alkylation half-reaction. Different models are presented to explore the protonation state of the catalytic histidine in the alkylation step and to evaluate the possibility of an interaction between the nucleophilic aspartate and the catalytic histidine.     

Promiscuity in alkaline phosphatase superfamily. Unraveling evolution through molecular simulations

We here present a theoretical study of the alkaline hydrolysis of a phosphodiester (methyl p-nitrophenyl phosphate or MpNPP) in the active site of Escherichia coli alkaline phosphatase (AP), a monoesterase that also presents promiscuous activity as a diesterase. The analysis of our simulations, carried out by means of molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials, shows that the reaction takes place through a D(N)A(N) or dissociative mechanism, the same mechanism employed by AP in the hydrolysis of monoesters. The promiscuous activity observed in this superfamily can be then explained on the basis of a conserved reaction mechanism. According to our simulations the specialization in the hydrolysis of phosphomonoesters or phosphodiesters, developed in different members of the superfamily, is a consequence of the interactions established between the protein and the oxygen atoms of the phosphate group and, in particular, with the oxygen atom that bears the additional alkyl group when the substrate is a diester. A water molecule, belonging to the coordination shell of the Mg(2+) ion, and residue Lys328 seem to play decisive roles stabilizing a phosphomonoester substrate, but the latter contributes to increase the energy barrier for the hydrolysis of phosphodiesters. Then, mutations affecting the nature or positioning of Lys328 lead to an increased diesterase activity in AP. Finally, the capacity of this enzymatic family to catalyze the reaction of phosphoesters having different leaving groups, or substrate promiscuity, is explained by the ability of the enzyme to stabilize different charge distributions in the leaving group using different interactions involving either one of the zinc centers or residues placed on the outer side of the catalytic site.     

Catalytic mechanism of 6-phosphogluconate dehydrogenase: a theoretical investigation

Density functional calculations are employed to theoretically explore the mechanism of all elementary reaction steps involved in the catalytic reaction of 6-phosphogluconate dehydrogenase (6PGDH). The model systems we choose for the enzyme contain the essential parts of the cofactor (NADP+), the substrate 6-phosphogluconate (6PG), and some key residues (Lys183 and Glu190) in the active site of sheep liver 6PGDH. The effect of the apoenzyme electrostatic environment on the studied reaction is treated by the self-consistent reaction-field method. Our calculations demonstrate that the first step of the catalytic reaction is the formation of a 3-keto 6PG intermediate, which proceeds through a concerted transition state involving a hydride transfer from 6PG to NADP+, and a proton transfer from 6PG to Lys183. The second step is the elimination of a CO2 molecule from 6-PG, concomitant with a proton transfer from Lys183 to 6-PG. In the final step, a concerted double proton transfer (one from Glu190 to the substrate, another from the substrate to Lys183) results in the final product, the keto form of ribulose 5-phosphate (Ru5P). The rate-limiting step is the formation of a 3-keto 6PG intermediate, with a free energy barrier of 22.7 kcal/mol at room temperature in the protein environment, and all three steps are calculated to be thermodynamically favorable. These results are in good agreement with the general acid/general base mechanism suggested from previous experiments for the 6PGDH reaction.     

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.     

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.     

Computational investigation of irreversible inactivation of the zinc-dependent protease carboxypeptidase A

Zinc proteases are ubiquitous and the zinc ion plays a central function in the catalytic mechanism of these enzymes. A novel class of mechanism-based inhibitors takes advantage of the zinc ion chemistry in carboxypeptidase A (CPA) to promote covalent attachment of an inhibitor to the carboxylate of Glu-270, resulting in irreversible inhibition of the enzyme. The effect of the active site zinc ion on irreversible inactivation of CPA was probed by molecular orbital (MO) calculations on a series of active site models and the Cl(-) + CH(3)Cl S(N)2 reaction fragment. Point charge models representing the active site reproduced energetics from full MO calculations at 12.0 A separation between the zinc and the central carbon of the S(N)2 reaction, but at 5.0 A polarization played an important role in moderating barrier suppression. ONIOM MO/MO calculations that included the residues within 10 A of the active site zinc suggest that about 75% of the barrier suppression arises from the zinc ion and its ligands. A model of the pre-reactive complex of the 2-benzyl-3-iodopropanoate inactivator with CPA was constructed from the X-ray structure of l-phenyl lactate bound in the active site of the enzyme. The model was fully solvated and minimized by using the AMBER force field to generate the starting structure for the ONIOM QM/MM calculations. Optimization of this structure led to the barrierless S(N)2 displacement of the iodide of the inhibitor by Glu-270, assisted by interaction of the zinc ion with the leaving group. The resulting product is in good agreement with the X-ray structure of the covalently modified enzyme obtained by irreversible inhibition of CPA by 2-benzyl-3-iodopropanoate.     

Computational Insights on the Hydride and Proton Transfer Mechanisms of D-Arginine Dehydrogenase

D-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is an amine oxidase which catalyzes the conversion of D-arginine into iminoarginine. It contains a non-covalent FAD cofactor that is involved in the oxidation mechanism. Based on substrate, solvent, and multiple kinetic isotope effects studies, a stepwise hydride transfer mechanism is proposed. It was shown that D-arginine binds to the active site of enzyme as α-amino group protonated, and it is deprotonated before a hydride ion is transferred from its α-C to FAD. Based on a mutagenesis study, it was concluded that a water molecule is the most likely catalytic base responsible from the deprotonation of α-amino group. In this study, we formulated computational models based on ONIOM method to elucidate the oxidation mechanism of D-arginine into iminoarginine using the crystal structure of enzyme complexed with iminoarginine. The calculations showed that Arg222, Arg305, Tyr249, Glu87, His 48, and two active site water molecules play key roles in binding and catalysis. Model systems showed that the deprotonation step occurs prior to hydride transfer step, and active site water molecule(s) may have participated in the deprotonation process.     

Catalytic mechanism of limonene epoxide hydrolase, a theoretical study

The catalytic mechanism of limonene epoxide hydrolase (LEH) was investigated theoretically using the density functional theory method B3LYP. LEH is part of a novel limonene degradation pathway found in Rhodococcus erythropolis DCL14, where it catalyzes the hydrolysis of limonene-1,2-epoxide to give limonene-1,2-diol. The recent crystal structure of LEH was used to build a model of the LEH active site composed of five amino acids and a crystallographically observed water molecule. With this model, hydrolysis of different substrates was investigated. It is concluded that LEH employs a concerted general acid/general base-catalyzed reaction mechanism involving protonation of the substrate by Asp101, nucleophilic attack by water on the epoxide, and abstraction of a proton from water by Asp132. Furthermore, we provide an explanation for the experimentally observed regioselective hydrolysis of the four stereoisomers of limonene-1,2-epoxide.     

Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study

The MMP-2 reaction mechanism is investigated by using different computational methodologies. First, quantum mechanical (QM) calculations are carried out on a cluster model of the active site bound to an Ace-Gly approximately Ile-Nme peptide. Along the QM reaction path, a Zn-bound water molecule attacks the Gly carbonyl group to give a tetrahedral intermediate. The breaking of the C-N bond is completed thanks to the Glu 404 residue that shuttles a proton from the water molecule to Ile-N atom. The gas-phase QM energy barrier is quite low ( approximately 14 kcal/mol), thus suggesting that the essential catalytic machinery is included in the cluster model. A similar reaction path occurs in the MMP-2 catalytic domain bound to an octapeptide substrate according to hybrid QM and molecular mechanical (QM/MM) geometry optimizations. However, the rupture of the Gly( P 1) approximately Ile( P 1') amide bond is destabilized in the static QM/MM calculations, owing to the positioning of the Ile( P 1') side chain inside the MMP-2 S 1' pocket and to the inability of simple energy miminization methodologies to properly relax complex systems. Molecular dynamics simulations show that these steric limitations are overcome easily through structural fluctuations. The energetic effect of structural fluctuations is taken into account by combining QM energies with average MM Poisson-Boltzmann free energies, resulting in a total free energy barrier of 14.8 kcal/mol in good agreement with experimental data. The rate-determining event in the MMP-2 mechanism corresponds to a H-bond rearrangement involving the Glu 404 residue and/or the Glu 404-COOH --> N-Ile( P 1') proton transfer. Overall, the present computational results and previous experimental data complement each other well in order to provide a detailed view of the MMPs catalytic mechanism.     

Computational study on hydrolysis of cefotaxime in gas phase and in aqueous solution

We are presenting a theoretical study of the hydrolysis of a β‐lactam antibiotic in gas phase and in aqueous solution by means of hybrid quantum mechanics/molecular mechanics potentials. After exploring the potential energy surfaces at semiempirical and density functional theory (DFT) level, potentials of mean force have been computed for the reaction in solution with hybrid PM3/TIP3P calculations and corrections with the B3LYP and M06‐2X functionals. Inclusion of the full molecule of the antibiotic, Cefotaxime, in the gas phase molecular model has been demonstrated to be crucial since its carboxylate group can activate a nucleophilic water molecule. Moreover, the flexibility of the substrate implies the existence of a huge number of possible conformers, some of them implying formation of intramolecular hydrogen bond interaction that can determine the energetics of the conformers defining the different states along the reaction profile. The results show PM3 provides results that are in qualitative agreement with DFT calculations. The free energy profiles show a step‐wise mechanism that is kinetically determined by the nucleophilic attack of a water molecule activated by the proton transfer to the carboxylate group of the substrate (the first step). However, since the main role of the β‐lactamase would be reducing the free energy barrier of the first step, and keeping in mind the barrier obtained from second intermediate to products, population of this second intermediate could be significant and consequently experimentally detected in β‐lactamases, as shown in the literature. © 2012 Wiley Periodicals, Inc.     

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.     

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.