Theoretical and experimental determination on two substrates turned over by 4-oxalocrotonate tautomerase

Quantum mechanical/molecular mechanical (QM/MM) calculations and experimental kinetic studies have been performed on 4-oxalocrotonate tautomerase (4OT) for two different substrates, 2-hydroxymuconate (2HM) and 2-oxo-4-hexenedioate (2o4hex). Potential (deltaE) and free energy (deltaG) paths for both steps of the reaction using both substrates were calculated to determine the free energy barriers and compared to the experimental values obtained from the kinetic studies via the transition state theory. In the first step, a proton from the hydroxyl oxygen on the second carbon of 2HM, or from the third carbon of 2o4hex, is abstracted by Pro-1. In the second step, the proton is transferred to the fifth carbon of the substrate to form the product, 2-oxo-3-hexenedioate (2o3hex). For both substrates we obtain a calculated deltaG of approximately 13 kcal/mol, in agreement with experimental determinations. The calculated free energy barrier difference deltaG2o4hex - deltaG2HM (deltadeltaG) is 0.87 kcal/mol. We obtained an experimental deltadeltaG of 0.85 kcal/mol. These results suggest that 2HM is turned over faster than 2o4hex by 4OT. However, these energy differences are so small that both 2HM and 2o4hex need to be taken into account in considering the mechanism of catalysis of 4OT.     

Reaction path determination for quantum mechanical/molecular mechanical modeling of enzyme reactions by combining first order and second order "chain-of-replicas" methods

A two-step procedure for the determination of reaction paths in enzyme systems is presented. This procedure combines two chain-of-states methods: a quantum mechanical/molecular mechanical (QM/MM) implementation of the nudged elastic band (NEB) method and a second order parallel path optimizer method both recently developed in our laboratory. In the first step, a reaction path determination is performed with the NEB method, along with a restrained minimization procedure for the MM environment to obtain a first approximation to the reaction path. In the second step, the calculated path is refined with the parallel path optimizer method. By combining these two methods the reaction paths are determined accurately, and in addition, the number of path optimization iterations are significantly reduced. This procedure is tested by calculating both steps of the isomerization of 2-oxo-4-hexenedioate by 4-oxalocrotonate tautomerase, which have been previously determined by our group. The calculated paths agree with the previously reported results and we obtain a reduction of 45%-55% in the number of path optimization cycles.     

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.     

Sequential quadratic programming method for determining the minimum energy path

A new method, referred to as the sequential quadratic programming method, is presented for determining minimum energy paths. The method is based on minimizing the points representing the path in the subspace perpendicular to the tangent of the path while using a penalty term to prevent kinks from forming. Rather than taking one full step, the minimization is divided into a number of sequential steps on an approximate quadratic surface. The resulting method can efficiently determine the reaction mechanism, from which transition state can be easily identified and refined with other methods. To improve the resolution of the path close to the transition state, points are clustered close to this region with a reparametrization scheme. The usefulness of the algorithm is demonstrated for the Muller-Brown potential, amide hydrolysis, and an 89 atom cluster taken from the active site of 4-oxalocrotonate tautomerase for the reaction which catalyzes 2-oxo-4-hexenedioate to the intermediate 2-hydroxy-2,4-hexadienedioate.     

Acylation is rate-limiting in glycosylasparaginase-catalyzed hydrolysis of N4-(4'-substituted phenyl)-L-asparagines

Glycosylasparaginase catalyzes the hydrolysis of the N-glycosylic bond between N-acetyl-D-glucosamine and L-asparagine in the catabolism of glycoproteins. The mechanism has been proposed to resemble that of serine proteases involving an acylation step where a nucleophilic attack by a catalytic Thr residue on the carbonyl carbon of the N-glycosylic bond gives rise to a covalent beta-aspartyl-enzyme intermediate, and a deacylation step to give the final products. The question posed in this study was: Is the acylation step the rate-limiting step in the hydrolysis reaction as in serine proteases? To answer this question a series of mostly new substituted anilides was synthesized and characterized, and their hydrolysis reactions catalyzed by glycosylasparaginase from human amniotic fluid were studied. Five N4-(4'-substituted phenyl)-L-asparagine compounds were synthesized and characterized: 4'-hydrogen, 4'-ethyl, 4'-bromo, 4'-nitro, and 4'-methoxy. Each of these anilides was a substrate for the enzyme. Hammett plots of the kinetic parameters showed that acylation is the rate-limiting step in the reaction and that upon binding the electron distribution of the substrate is perturbed toward the transition state. This is the first direct evidence that acylation is the rate-limiting step in the enzyme-catalyzed reaction. A Br?nsted plot indicates a small, negative charge (-0.25) on the nitrogen atom of the leaving group anilines containing electron-withdrawing groups, and a small, positive charge (0.43) on the nitrogen atom of the leaving group anilines containing electron-donating groups. The free energy (incremental) change of binding (delta deltaGb) in the enzyme-substrate transition state complexes shows that substitution of a substituted phenyl group for the pyranosyl group in the natural substrate results in an overall loss of binding energy equivalent to a weak hydrogen bond, the magnitude of which is dependent on the substituent group. The data are consistent with a mechanism for glycosylasparaginase involving rapid formation of a tetrahedral structure upon substrate binding, and a rate-limiting breakdown of the tetrahedral structure to a covalent beta-aspartyl-enzyme intermediate that is dependent on the electronic properties of the substituent group and on the degree of protonation of the leaving group in the transition state by a general acid.     

磷酰化丝氨酸形成六配位磷中间体的理论研究

用MNDO方法对磷酰化丝氨酸仿生化反应机理中六配位磷中间体的形成过程进行了研究.磷酰化丝氨酸(1)形成分子内磷酸-羧酸分子内混酐的五配位磷中间体(2)后,其酸性质子解离,分子经过具有氢桥键结构的过渡态,使氨基酸侧链羟基上的氢通过氢键作用向磷上的O₁进行转移,然后再经过构型由三角双锥向八面体的转变,形成六配位磷中间体(3).氢桥键的存在使反应过渡态能量降低,其相对能量为148.5kJ/mol.理论计算较成功的解释了六配位磷中间体的形成机理以及磷酰化丝氨酸仿生化反应中羧基和侧链羟基共同参与的实验结果.     

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.     

Hybrid QM/MM study of thio effects in transphosphorylation reactions

Results of a series of hybrid quantum mechanical/molecular mechanical (QM/MM) activated dynamics simulations of thio effects in the transphosphorylation (methanolysis) of a 2'-ribose, 5'-methyl phosphate-diester under basic conditions are presented. Single and double substitutions in the nonbridging oxygen positions exhibit thio effects in accord with experimental data and show the existence of a stable intermediate. Thio substitution at the 2' and 5' positions resulted in reactions having a single transition state with increased and decreased free energy barriers, respectively, relative to the unsubstituted reaction. In all of the reactions except for the 5' substitution, the rate-limiting step corresponds to exocyclic cleavage. In the 5' substitution reaction, the rate-limiting step corresponds to endocyclic cleavage and shows a considerable reverse thio effect, in accord with experimental observations of phosphates with enhanced leaving groups. Thio substitution at the 3' position results in a mild reverse thio effect that arises from electronic stabilization of the dianionic transition state. The results presented here provide an important step toward the development and application of new hybrid QM/MM methods that, combined with experiment, may provide a detailed picture of the molecular mechanisms of RNA catalysis.     

Computational determination of fundamental pathway and activation barriers for acetohydroxyacid synthase‐catalyzed condensation reactions of α‐keto acids

Acetohydroxyacid synthase (AHAS) is the first common enzyme in the biosynthetic pathway leading to the production of various branched‐chain amino acids. AHAS is recognized as a promising target for new antituberculosis drugs, antibacterial drugs, and herbicides. Extensive first‐principles quantum mechanical (QM) and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have enabled us, in this study, to uncover the fundamental reaction pathway, determine the activation barriers, and obtain valuable insights concerning the specific roles of key amino acid residues for the common steps of AHAS‐catalyzed condensation reactions of α‐keto acids. The computational results reveal that the rate‐determining step of the AHAS‐catalyzed reactions is the second reaction step and that the most important amino acid residues involved in the catalysis include Glu144′, Gln207′, Gly121′, and Gly511 that form favorable hydrogen bonds with the reaction center (consisting of atoms from the substrate and cofactor) during the reaction process. In addition, Glu144′ also accepts a proton from cofactor thiamin diphosphate (ThDP) through hydrogen bonding during the catalytic reaction. The favorable interactions between the reaction center and protein environment remarkably stabilize the transition state and, thus, lower the activation barrier for the rate‐determining reaction step by ～20 kcal/mol. The activation barrier calculated for the rate‐determining step is in good agreement with the experimental activation barrier. The detailed structural and mechanistic insights should be valuable for rational design of novel, potent AHAS inhibitors that may be used as promising new anti‐tuberculosis drugs, antibacterial drugs, and/or herbicides to overcome drug resistance problem. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

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

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

Cycloreversion of cyclobutanes

Semiempirical MO calculations with the method SINDO1 were performed to study the potential energy surface of cyclobutane and several substituted cyclobutanes with substituents F, OCH₃ and CN. The reaction pathway with the lowest activation energy leading to two ethylenic fragments is nonconcerted. One carbon bond is broken after symmetric opening of two adjacent bond angles and twisting of the carbon framework. The first transition state is asymmetric and diradicaloid. The reaction proceeds to a diradicaloid, non-zwitterionic intermediate. The second transition state is characterized by bond breaking of the inner carbon-carbon bond. For the unsubstituted case, the barrier for free rotation of the outer methylenic groups was also calculated. In comparison, the unsubstituted reaction is characterized by transition states of almost equal energy, whereas in the substituted reactions the barriers for the second bond breaking are much higher than for the first bond breaking step.     

The 4-oxalocrotonate tautomerase- and YwhB-catalyzed hydration of 3E-haloacrylates: implications for the evolution of new enzymatic activities

4-Oxalocrotonate tautomerase (4-OT) catalyzes the conversion of 2-oxo-4E-hexenedioate to 2-oxo-3E-hexenedioate through the intermediate, 2-hydroxy-2,4E-hexadienedioate. 4-OT and a homologue found in Bacillus subtilis (designated YwhB) share sequence identity and two key catalytic groups, Pro-1 and Arg-11, with the two subunits comprising trans-3-chloroacrylic acid dehalogenase (CaaD). 4-OT and YwhB have now been found to display a low-level hydratase activity, resulting in the dehalogenation of 3E-haloacrylates. The enzymes are highly selective for the (E)-isomer, and Pro-1 is critical for the activity while an arginine is likely required. Two mechanisms are proposed in which Pro-1 functions as a general base or a general acid catalyst and, along with the arginine, facilitates the Michael addition of water. Both mechanisms suggest an intriguing route for the evolution of the CaaD activity. One or more mutations could decrease the hydrophobic environment of the active site, which would make it more favorable for a hydrolytic reaction, thereby raising the pKa of Pro-1 and increasing the concentration of enzyme in the reactive form.     

AM1 study of acid-catalyzed hydrolysis of maleamic (4-amino-4-oxo-2-butenoic) acids

The acid-catalyzed hydrolysis mechanisms of maleamic (4-amino-4-oxo-2-butenoic) acids were studied using AM1 method. The reaction proceeded mainly in two steps: (1) nucleophilic attack of the undissociated carboxyl group on the adjacent aminocarbonyl carbon via a zwitterionic intermediate; and (2) the rate-determining proton transfer to form the zwitterionic tetrahedral intermediate. In each step, the hydration of water and hydronium ion molecules was important in stabilizing the polarized intermediates. The substituent effects at the amide moiety and the 2,3-positions of the maleamic acids were qualitatively estimated for each step.     

Molecular details from computational reaction dynamics for the cellobiohydrolase I glycosylation reaction

Glycosylation of cellobiose hydrolase I (CBHI), is a key step in the processing and degradation of cellulose. Here the pathways and barriers of the reaction are explored using the free energy from adaptive reaction coordinate forces (FEARCF) reaction dynamics method coupled with SCC-DFTB/MM. In many respects CBHI follows the expected general GH7 family mechanism that involves the Glu-X-Asp-X-X-Glu motif. However, critical electronic and conformational details, previously not known, were discovered through our computations. The central feature that ensures the success of the glycosylation reaction are the Glu212 nucleophile's hydrogen bond to the hydroxyl on C2, of the glucose in the -1 position of the cellulosic strand. This Glu212 function restricts the C2 hydroxyl in such a way as to favor the formation of the (4)E ring pucker of the -1 position glucose. A frontier molecular orbital analysis of the structures along the reaction surface proves the existence of an oxocarbenium ion, which has both transition state and intermediate character. The transition state structure is able to descend down the glycosylation pathway through the critical combination of Asp214 (HOMO), ring oxygen (LUMO), and Glu212 (HOMO), anomeric carbon (LUMO) interactions. Using the fully converged FEARCF SCC-DFTB/MM reaction surface, we find a barrier of 17.48 kcal/mol separating bound cellulose chain from the glycosylated CBHI. Taking recrossing into account gives k(cat) = 0.415 s(-1) for cellobiohydrolase glycosylation.     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

Hydrolysis of cyclic phosphates by ribonuclease A: a computational study using a simplified ab initio quantum model

The second step in the enzyme-catalyzed hydrolysis of phosphate esters by ribonuclease A (RNase A) was studied using an ab initio quantum-based model of the active site including constrained parts of three critical residues, His-12, His-119, and Lys-41, and a small substrate. The competition between release of the cyclic phosphate intermediate and subsequent hydrolysis following transphosphorylation was explored to determine the electronic factors that contribute to preferential intermediate product release observed experimentally. The structural and energetic results obtained at both the RHF and MP2 levels reveal several contributing factors consistent with experimental observation. Although the intrinsic electronic effects tend to favor hydrolysis slightly with an overall activation free energy of approximately 70 kJ mol(-1), entropic and environmental effects favor release of the cyclic phosphate intermediate over hydrolysis. Exploration of the second, hydrolysis step also revealed interesting similarity with the transphosphorylation step, including the observation of autocatalysis by the substrate. Moreover, both steps of the overall RNase A reaction reveal multiple pathways involving proton transfers to sites of similar proton affinities. The anionic phosphate in both steps can act as a stable proton binding site as protons are moved around the active site throughout the progress of the reaction. These results suggest autocatalysis may be representative of more general behavior in enzymes containing highly charged substrates, especially phosphates.     

Fundamental reaction pathway and free energy profile for inhibition of proteasome by Epoxomicin

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to provide the first detailed computational study on the possible mechanisms for reaction of proteasome with a representative peptide inhibitor, Epoxomicin (EPX). The calculated results reveal that the most favorable reaction pathway consists of five steps. The first is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is nucleophilic attack on the carbonyl carbon of EPX by the negatively charged Thr1-O(γ) atom, followed by a proton transfer from Thr1-N(z) to the carbonyl oxygen of EPX (third step). Then, Thr1-N(z) attacks on the carbon of the epoxide group of EPX, accompanied by the epoxide ring-opening (S(N)2 nucleophilic substitution) such that a zwitterionic morpholino ring is formed between residue Thr1 and EPX. Finally, the product of morpholino ring is generated via another proton transfer. Noteworthy, Thr1-O(γ) can be activated directly by Thr1-N(z) to form the zwitterionic intermediate (with a free energy barrier of only 9.9 kcal/mol), and water cannot assist the rate-determining step, which is remarkably different from the previous perception that a water molecule should mediate the activation process. The fourth reaction step has the highest free energy barrier (23.6 kcal/mol) which is reasonably close to the activation free energy (～21-22 kcal/mol) derived from experimental kinetic data. The obtained novel mechanistic insights should be valuable for not only future rational design of more efficient proteasome inhibitors but also understanding the general reaction mechanism of proteasome with a peptide or protein.     

Investigation of a general base mechanism for ester hydrolysis in C-C hydrolase enzymes of the alpha/beta-hydrolase superfamily: a novel mechanism for the serine catalytic triad

Previous mechanistic and crystallographic studies on two C-C hydrolase enzymes, Escherichia coli MhpC and Burkholderia xenovorans BphD, support a general base mechanism for C-C hydrolytic cleavage, rather than the nucleophilic mechanism expected for a serine hydrolase. The role of the active site serine residue could be to form a hydrogen bond with a gem-diolate intermediate, or to protonate such an intermediate. Hydrolase BphD is able to catalyse the hydrolysis of p-nitrophenyl benzoate ester substrates, which has enabled an investigation of these mechanisms using a Hammett analysis, and comparative studies upon five serine esterases and lipases from the alpha/beta-hydrolase family. A reaction parameter (rho) value of +0.98 was measured for BphD-catalysed ester hydrolysis, implying a build-up of negative charge in the transition state, consistent with a general base mechanism. Values of +0.31-0.61 were measured for other serine esterases and lipases, for the same series of esterase substrates. Pre-steady state kinetic studies of ester hydrolysis, using p-nitrophenyl acetate as the substrate, revealed a single step kinetic mechanism for BphD-catalysed ester hydrolysis, with no burst kinetics. A general base mechanism for BphD-catalysed ester hydrolysis is proposed, in which Ser-112 stabilises an oxyanion intermediate through hydrogen bonding, and assists the rotation of this oxyanion intermediate via proton transfer, a novel reaction mechanism for the serine catalytic triad.     

Substrate orientation in 4-oxalocrotonate tautomerase and its effect on QM/MM energy profiles

The tautomerization of 2-oxo-4E-hexendioate by 4-oxalocrotonate tautomerase has been studied by quantum mechanical/molecular mechanical (QM/MM) methods using three models, A-C, with different substrate orientations. The computed QM/MM energy profiles are rather different. Various energy partitioning analyses indicate the origin of these differences and the role of the active site residues for different substrate orientations. The proposed new model C is preferred over the previously used models A and B because it combines favorable substrate binding geometries with reasonable barriers and is consistent with the experimental evidence from mutation studies concerning the catalytic ability of specific residues in the binding site, especially R11'.     

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.