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.     

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.     

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.     

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.     

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.     

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.     

Role of the catalytic triad and oxyanion hole in acetylcholinesterase catalysis: an ab initio QM/MM study

The initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE) has been studied by a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach. The reaction proceeds through the nucleophilic addition of the Ser203 O to the carbonyl C of acetylcholine, and the reaction is facilitated by simultaneous proton transfer from Ser203 to His447. The calculated potential energy barrier at the MP2(6-31+G) QM/MM level is 10.5 kcal/mol, consistent with the experimental reaction rate. The third residue of the catalytic triad, Glu334, is found to be essential in stabilizing the transition state through electrostatic interactions. The oxyanion hole, formed by peptidic NH groups from Gly121, Gly122, and Ala204, is also found to play an important role in catalysis. Our calculations indicate that, in the AChE-ACh Michaelis complex, only two hydrogen bonds are formed between the carbonyl oxygen of ACh and the peptidic NH groups of Gly121 and Gly122. As the reaction proceeds, the distance between the carbonyl oxygen of ACh and NH group of Ala204 becomes smaller, and the third hydrogen bond is formed both in the transition state and in the tetrahedral intermediate.     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

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.     

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.     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

Force-field parameters for the simulation of tetrahedral intermediates of serine hydrolases

CHARMM force-field parameters are reported for the tetrahedral intermediate of serine hydrolases. The fitting follows the standard protocol proposed for CHARMM22. The reference data include ab initio (RHF/6-31G*) interaction energies of complexes between water and the model compound 1,1-dimethoxyethoxide, torsional profiles of related model compounds from correlated ab initio (MP2/6-311+G*//B3LYP/6-31+G*) calculations, as well as molecular geometries and vibrational frequencies from density functional theory (B3LYP/6-31+G*). The optimized parameters reproduce the target data well. Their utility is demonstrated by a QM/MM study of the tetrahedral intermediate in Bacillus subtilis lipase A, and by classical molecular modeling of enantioselectivity in Pseudomonas aeruginosa lipase and its mutants.     

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.     

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.     

Toward accurate barriers for enzymatic reactions: QM/MM case study on p-hydroxybenzoate hydroxylase

The hydroxylation reaction catalyzed by p-hydroxybenzoate hydroxylase has been investigated by quantum mechanical/molecular mechanical (QM/MM) calculations at different levels of QM theory. The solvated enzyme was modeled (approximately 23,000 atoms in total, 49 QM atoms). The geometries of reactant and transition state were optimized for ten representative pathways using semiempirical (AM1) and density functional (B3LYP) methods as QM components. Single-point calculations at B3LYP/MM optimized geometries were performed with local correlation methods [LMP2, LCCSD(T0)] and augmented triple-zeta basis sets. A careful validation of the latter approach with regard to all computational parameters indicates convergence of the QM contribution to the computed barriers to within approximately 1 kcal mol(-1). Comparison with the available experimental data supports this assessment.     

Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydrates

An extensive quantum mechanical study of a water dimer suggests that the introduction of a diffuse function into the basis set, which significantly reduces the basis set superposition error (BSSE) in the hydrogen bonding energy calculation, is the key to better calculations of the potential energy surfaces of carbohydrates. This article examines the potential energy surfaces of selected d -aldo- and d -ketohexoses (a total of 82 conformers) by quantum mechanics (QM) and molecular mechanics (MM) methods. In contrast to the results with a smaller basis set (B3LYP/6-31G** 5d), we found at the higher level calculation (B3LYP/6-311++G(2d,2p)//B3LYP/6-31G** 5d) that, in most cases, the furanose forms are less stable than the pyranose forms. These discrepancies are mainly due to the fact that intramolecular hydrogen bonding energies are overestimated in the lower level calculations. The higher level QM calculations of the potential energy surfaces of d -aldo- and d -ketohexoses now are more comparable to the MM3 results. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1593–1603, 1999     

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.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

A theoretical study of the aminolysis reaction of lysine 199 of human serum albumin with benzylpenicillin: consequences for immunochemistry of penicillins

Herein, we present results of a computational study on benzylpenicillin attachment to Lys199 of human serum albumin via an aminolysis reaction. The internal geometry of the reactive part of the system was taken from previous work at the B3LYP/6-31+G* level on the water-assisted aminolysis reaction of a penicillin model compound (Díaz, N.; Suárez, D.; Sordo, T. L. J. Am. Chem. Soc. 2000, 122, 6710--6719). The protein environment around Lys199, the 6-acylamino side chain, and the 2-methyl groups of benzylpenicillin were relaxed by carrying out geometry optimizations with a hybrid QM/MM method (PM3/AMBER). Two different mechanistic routes were explored: a one-step water-assisted process and a carboxylate and water-assisted route in which the beta-lactam carboxylate and the ancillary water molecule mediate the proton transfer from the epsilon-amino group of Lys199 to the beta-lactam leaving N atom. The corresponding energy profiles in the protein combine the B3LYP/6-31+G* and PM3 energies of the reactive subsystem (benzylpenicillin + Lys199 side chain + the ancillary water molecule) and semiempirical PM3 energies of the entire system evaluated with a "divide and conquer" linear-scaling method. It is observed that penicillin haptenation to HSA can proceed through the water-assisted concerted mechanism which is calculated to have a high energy barrier of approximately 38 kcal/mol, in agreement with the experimentally observed slow reaction kinetics.