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.     

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.     

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.     

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.     

Ab initio and density functional theory studies of the catalytic mechanism for ester hydrolysis in serine hydrolases

We present results from ab initio and density functional theory studies of the mechanism for serine hydrolase catalyzed ester hydrolysis. A model system containing both the catalytic triad and the oxyanion hole was studied. The catalytic triad was represented by formate anion, imidazole, and methanol. The oxyanion hole was represented by two water molecules. Methyl formate was used as the substrate. In the acylation step, our computations show that the cooperation of the Asp group and oxyanion hydrogen bonds is capable of lowering the activation barrier by about 15 kcal/mol. The transition state leading to the first tetrahedral intermediate in the acylation step is rate limiting with an activation barrier (ΔE₀) of 13.4 kcal/mol. The activation barrier in the deacylation step is smaller. The double-proton-transfer mechanism is energetically unfavorable by about 2 kcal/mol. The bonds between the Asp group and the His group, and the hydrogen bonds in the oxyanion hole, increase in strength going from the Michaelis complex toward the transition state and the tetrahedral intermediate. In the acylation step, the tetrahedral intermediate is a very shallow minimum on the energy surface and is not viable when molecular vibrations are included. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 89–103, 1998     

Quantum chemical calculations on structural models of the catalytic site of chymotrypsin: comparison of calculated results with experimental data from NMR spectroscopy

Hybrid density functional quantum mechanical calculations were used to study the strength of the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in chymotrypsin and how it changes along the reaction coordinate. Comparison of experimental shifts with the results of chemical shift calculations on a variety of small molecules, including species containing very strong hydrogen bonds, has validated the overall approach and provided the means for calibrating and correcting the calculated values. Models of the active site of chymotrypsin in its resting state and tetrahedral intermediate state were derived from high-resolution X-ray structures. The distance between His(57) N(delta)(1) and Asp(102) O(delta)(1) in each model was varied between 2.77 A (weak hydrogen bond) and 2.50 A (extremely strong hydrogen bond), and the one-dimensional potential energy surface of the hydrogen-bonded proton (or deuteron/triton) was determined. The zero-point energy, probability distribution, and chemical shift were determined for each distance. Calculated values for NMR chemical shifts, NMR chemical shift differences between (1)H and (3)H, and (2)H/(1)H fractionation factors were compared with published experimental values. Energies provided by the calculations indicated that the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in the chymotrypsin active site increases in strength by 11 kcal mol(-)(1) in going from the resting state of the enzyme to the tetrahedral intermediate state. This result confirms the hypothesis that the strengthened hydrogen bond plays an important role in lowering the energy of the transition state and, hence, in the catalytic efficiency of the enzyme. Models of the transition state that best fit the experimental data are consistent with a "strong" hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) but apparently not a "low-barrier" or "very strong" hydrogen bond.     

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.     

The catalytic mechanism of mouse renin studied with QM/MM calculations

Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.     

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.     

Ab initio QM/MM dynamics simulation of the tetrahedral intermediate of serine proteases: insights into the active site hydrogen-bonding network

Ab initio QM/MM dynamics simulation is employed to examine the stability of the tetrahedral intermediate during the deacylation step in elastase-catalyzed hydrolysis of a simple peptide. An extended quantum region includes the catalytic triad, the tetrahedral structure, and the oxyanion hole. The calculations indicate that the tetrahedral intermediate of serine proteases is a stable species on the picosecond time scale. On the basis of geometrical and dynamical properties, and in agreement with many experimental and theoretical studies, it is suggested that the crucial hydrogen bonds involved in stabilizing this intermediate are between Asp-102 and His-57 and between the charged oxygen of the intermediate and the backbone N-H group of Gly-193 in the oxyanion hole. The mobility of the imidazolium ring between O(w) and O(gamma), two of the oxygens of the tetrahedral structure, shows how the intermediate could proceed toward the product state without a "ring-flip mechanism", proposed earlier on the basis of NMR data. In addition to the proposed C(epsilon)(1)-H.O hydrogen bond between the imidazolium ring and the backbone carbonyl of Ser-214, we observe an alternative C(epsilon)(1)-H.O hydrogen bond with the backbone carbonyl of Thr-213, that can stabilize the intermediate during the imidazolium movement. Proton hopping occurs between Asp-102 and His-57 during the simulation. The proton is, however, largely localized on the nitrogen, and hence it does not participate in a low-barrier hydrogen bond. The study also suggests factors that may be implicated in product release: breaking the hydrogen bond of the charged oxygen with the backbone of Ser-195 in the oxyanion hole and a loop opening between residues 216-225 that enables the breaking of a hydrogen bond in subsite S(3).     

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.     

A theoretical study on the substrate deacylation mechanism of class C beta-lactamase

The whole reaction of the deacylation of class C beta-lactamase was investigated by performing quantum chemical calculations under physiological conditions. In this study, the X-ray crystallographic structure of the inhibitor moxalactam-bound class C beta-lactamase (Patera et al. J. Am. Chem. Soc. 2000, 122, 10504-10512.) was utilized and moxalactam was changed into the substrate cefaclor. A model for quantum chemical calculations was constructed using an energy-minimized structure of the substrate-bound enzyme obtained by molecular mechanics calculation, in which the enzyme was soaked in thousands of TIP3P water molecules. It was found that the deacylation reaction consisted of two elementary processes. The first process was formation of a tetrahedral intermediate, which was initiated by the activation of catalytic water by Tyr150, and the second process was detachment of the hydroxylated substrate from the enzyme, which associated with proton transfer from the side chain of Lys67 to Ser64O(gamma). The first process is a rate-determining process, and the activation energy was estimated to be 30.47 kcal/mol from density functional theory calculations considering electron correlation (B3LYP/6-31G**). The side chain of Tyr150 was initially in a deprotonated state and was stably present in the active site of the acyl-enzyme complex, being held by Lys67 and Lys315 cooperatively.     

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.

     

Quantum mechanical/molecular mechanical study of three stationary points along the deacylation step of the catalytic mechanism of elastase

A large amount of experimental as well as theoretical information is available about the mechanism of serine proteases, but many questions remain unanswered. Here we study the deacylation step of the reaction mechanism of elastase. The water molecule in the acyl-enzyme active site, the binding mode of the carbonyl oxygen in the oxyanion hole, the characteristics of the tetrahedral intermediate structure, and the mobility of the imidazole ring of His-57 were studied with quantum mechanical/molecular mechanical methods. The models are based on a recent high-resolution crystal structure of the acyl-enzyme intermediate. The nucleophilic water in the active site of the acyl-enzyme has been shown to have two minima that differ by only 2 kcalmol⁻¹ in energy. The carbonyl group of the acyl-enzyme is located in the oxyanion hole and is positioned for attack by the hydrolytic water. The tetrahedral intermediate is a weakly bonded system, which is electrostatically stabilized by short hydrogen bonds to the backbone NH groups of Gly-193 and Ser-195 in the oxyanion hole. The short distance between the N^ɛ2 of His-57 and the O^γ of Ser-195 in the tetrahedral intermediate indicates a small movement of the imidazole ring towards the product in the deacylation step. The carbonyl group of the enzyme-product complex is not held strongly in the oxyanion hole, which shows that the peptide is first released from the oxyanion hole before it leaves the active site to regenerate the native state of the enzyme. Received: 11 September 2000 / Accepted: 15 September 2000 / Published online: 21 March 2001     

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.     

Subangstrom crystallography reveals that short ionic hydrogen bonds, and not a His-Asp low-barrier hydrogen bond, stabilize the transition state in serine protease catalysis

To address questions regarding the mechanism of serine protease catalysis, we have solved two X-ray crystal structures of alpha-lytic protease (alphaLP) that mimic aspects of the transition states: alphaLP at pH 5 (0.82 A resolution) and alphaLP bound to the peptidyl boronic acid inhibitor, MeOSuc-Ala-Ala-Pro-boroVal (0.90 A resolution). Based on these structures, there is no evidence of, or requirement for, histidine-flipping during the acylation step of the reaction. Rather, our data suggests that upon protonation of His57, Ser195 undergoes a conformational change that destabilizes the His57-Ser195 hydrogen bond, preventing the back-reaction. In both structures the His57-Asp102 hydrogen bond in the catalytic triad is a normal ionic hydrogen bond, and not a low-barrier hydrogen bond (LBHB) as previously hypothesized. We propose that the enzyme has evolved a network of relatively short hydrogen bonds that collectively stabilize the transition states. In particular, a short ionic hydrogen bond (SIHB) between His57 Nepsilon2 and the substrate's leaving group may promote forward progression of the TI1-to-acylenzyme reaction. We provide experimental evidence that refutes use of either a short donor-acceptor distance or a downfield 1H chemical shift as sole indicators of a LBHB.     

Discussion of the catalytic pathway of cysteine proteases based on AM1 calculations

The deacylation reaction of the cysteine protease papain was examined by AM1 reaction-coordinate calculations. The transition-state (TS) structure was extracted from the reaction pathway as corresponding to the maximum point along this minimum-energy pathway. Consistent with experimental kinetic data revealing that deacylation is about 60 times faster for thioester (---C(O)---S---) than dithioester (---C(S)---S---) intermediates, calculated E_a values are about 10 kcal mol⁻¹ lower for the former than the latter. The calculated partial atomic charges indicate that the C=O carbon in the thioester is a good site for nucleophilic attack whereas the corresponding C=S carbon in the dithioester is a poor site. The present calculations reveal that the enzyme's oxyanion hole contributes about 9 kcal mol⁻¹ toward reducing E_a for the anionic tetrahedral intermediate and TS structure. On the other hand, the presence of Asn in the putative Asn-His-Cys catalytic triad contributes only about l kcal mol⁻¹ toward reducing their E_a value. The presence of this Asn, however, did appear to stabilize His in its protonated form (ImH⁺) over its unprotonated form (Im). Two novel mechanisms are introduced to explain the unusual effect of a remote X substituent on the deacylation kinetics of the substrate family under consideration. The first mechanism invokes a “field effect” while the second mechanism embodies the concepts of induction and homoconjugation. The unique feature of these two mechanisms is that, unlike other proposed models, they circumvent the requirement for a close N…S interaction which has stimulated controversy.     

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.     

A theoretical study on the catalytic mechanism of Mus musculus adenosine deaminase

The catalytic mechanism of Mus musculus adenosine deaminase (ADA) has been studied by quantum mechanics and two‐layered ONIOM calculations. Our calculations show that the previously proposed mechanism, involving His238 as the general base to activate the Zn‐bound water, has a high activation barrier of about 28 kcal/mol at the proposed rate‐determining nucleophilic addition step, and the corresponding calculated kinetic isotope effects are significantly different from the recent experimental observations. We propose a revised mechanism based on calculations, in which Glu217 serves as the general base to abstract the proton of the Zn‐bound water, and the protonated Glu217 then activates the substrate for the subsequent nucleophilic addition. The rate‐determining step is the proton transfer from Zn‐OH to 6‐NH₂ of the tetrahedral intermediate, in which His238 serves as a proton shuttle for the proton transfer. The calculated kinetic isotope effects agree well with the experimental data, and calculated activation energy is also consistent with the experimental reaction rate. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

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.