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     

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.     

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

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

Oxidation of zinc-thiolate complexes of biological interest by hydrogen peroxide: a theoretical study

Zinc-thiolate complexes play a major structural and functional role in the living cell. Their stability is directly related to the thiolate reactivity toward reactive oxygen species naturally present in the cell. Oxidation of some zinc-thiolate complexes has a functional role, as is the case of zinc finger redox switches. Herein, we report a theoretical investigation on the oxidation of thiolate by hydrogen peroxide in zinc finger cores of CCCC, CCHC, and CCHH kinds containing either cysteine or histidine residues. In the case of the CCCC core, the calculated energy barrier for the oxidation to sulfenate of the complexed thiolate was found to be 16.0 kcal mol(-1), which is 2 kcal mol(-1) higher than that for the free thiolate. The energy barrier increases to 19.3 and 22.2 kcal mol(-1) for the monoprotonated and diprotonated CCCC cores, respectively. Substitution of cysteine by histidine also induces an increase in the magnitude of the reaction energy barrier: It becomes 20.0 and 20.9 kcal mol(-1) for the CCCH and CCHH cores, respectively. It is concluded that the energy barrier for the oxidation of zinc fingers is strictly dependent on the type of ligands coordinated to zinc and on the protonation state of the complex. These changes in the thiolate reactivity can be explained by the lowering of the nucleophilicity of complexed sulfur and by the internal reorganization of the complex (changes in the metal-ligand distances) upon oxidation. The next reaction steps subsequent to sulfenate formation are also considered. The oxidized thiolate (sulfenate) is predicted to dissociate very fast: For all complexes, the calculated dissociation energy barrier is lower than 3 kcal mol(-1). It is also shown that the dissociated sulfenic acid can interact with a free thiolate to form a sulfur-sulfur (SS) bridge in a reaction that is predicted to be quasi-diffusion limited. The interesting biological consequences of the modulation of thiolate reactivity by the chemical composition of the zinc finger cores are discussed.     

Mechanisms of antibiotic resistance: QM/MM modeling of the acylation reaction of a class A beta-lactamase with benzylpenicillin

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.     

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.     

Theoretical studies on the deacylation step of serine protease catalysis in the gas phase, in solution, and in elastase

The deacylation step of serine protease catalysis is studied using DFT and ab initio QM/MM calculations combined with MD/umbrella sampling calculations. Free energies of the entire reaction are calculated in the gas phase, in a continuum solvent, and in the enzyme elastase. The calculations show that a concerted mechanism in the gas phase is replaced by a stepwise mechanism when solvent effects or an acetate ion are added to the reference system, with the tetrahedral intermediate being a shallow minimum on the free energy surface. In the enzyme, the tetrahedral intermediate is a relatively stable species ( approximately 7 kcal/mol lower in energy than the transition state), mainly due to the electrostatic effects of the oxyanion hole and Asp102. It is formed in the first step of the reaction, as a result of a proton transfer from the nucleophilic water to His57 and of an attack of the remaining hydroxyl on the ester carbonyl. This is the rate-determining step of the reaction, which requires approximately 22 kcal/mol for activation, approximately 5 kcal/mol less than the reference reaction in water. In the second stage of the reaction, only small energy barriers are detected to facilitate the proton transfer from His57 to Ser195 and the breakdown of the tetrahedral intermediate. Those are attributed mainly to a movement of Ser195 and to a rotation of the His57 side chain. During the rotation, the imidazolium ion is stabilized by a strong H-bond with Asp102, and the C(epsilon)(1)-H...O H-bond with Ser214 is replaced by one with Thr213, suggesting that a "ring-flip mechanism" is not necessary as a driving force for the reaction. The movements of His57 and Ser195 are highly correlated with rearrangements of the binding site, suggesting that product release may be implicated in the deacylation process.     

Catalytic mechanism and product specificity of the histone lysine methyltransferase SET7/9: an ab initio QM/MM-FE study with multiple initial structures

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.     

Fixation of heterocumulenes. 1. A theoretical study on the irreversible reaction of CO(2) with a 2-lithio-1,3-dithiane

The mechanism of electrophilic addition of CO(2) to 2-lithio-2-phenyl-1, 3-dithiane-tetrahydrofuran-tetramethylethylendiamine, a Corey-Seebach "umpolungs" reagent for nucleophilic acylation, was investigated at the B3YLP/6-311+G//HF/6-31+G level of theory with specific solvation effects being included in the study. The overall reaction is exothermic by 18.8 kcal/mol and proceeds via precomplexation of the CO(2) (necessary for activation). The reaction barrier is calculated to be 11.9 kcal/mol (relative to the precomplex 2a) and represents a lower boundary for the activation energy. The reaction barrier is shown to originate from specific solvent effects. A detailed survey of the electronic effects governing the reaction is given. The ab initio results were then compared with semiempirical PM3 calculations, which were extended to include the heterocumulenes COS and CS(2). As expected, the reduced electrophilicity of these molecules leads to a higher activation barrier by retention of the same reaction mechanism.     

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.     

Influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine. A theoretical study

The influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine has been studied using density functional theory (DFT) methods. For the neutral system, two different pathways (with retention and inversion of configuration at the C1' anomeric carbon) have been found, both of them consisting of two steps and involving the formation of a dihydrofurane-like intermediate. The Gibbs free energy barrier for the first step is very high in both cases (53 and 46 kcal/mol for the process with inversion and with retention, respectively). However, the N7-protonated system shows a very different mechanism which consists of two steps. The first one leads to the formation of an oxacarbenium ion intermediate, with a Gibbs free energy barrier of 27 kcal/mol, and the second one corresponds to the nucleophilic attack of the water molecule to the oxacarbenium ion and takes place with a barrier of 1.3 kcal/mol. Thus, these results agree with a stepwise SN1 mechanism (DN*AN), with a discrete intermediate formed between the leaving group and the nucleophile approach, and show that N7 protonation strongly catalyzes the hydrolysis of the N-glycosidic bond, making the guanine a better leaving group. Finally, kinetic isotope effects have been calculated for the protonated system, and the results obtained are in very good agreement with experimental data for analogous systems.     

Theoretical study of 1,6-anhydrosugar formation from phenyl D-glucosides under basic condition: reasons for higher reactivity of β-anomer

Degradation of anomeric phenyl d-glucosides to levoglucosan under basic condition is theoretically studied. MP4(SDQ)//DFT(B3LYP)-computational results indicate that the degradation of phenyl α-glucoside (R(α)) occurs via the S(N)icB mechanism. In this mechanism, the oxyanion at the C6, which is formed through deprotonation of the OH group, directly attacks the anomeric carbon. On the other hand, the degradation of phenyl β-glucoside (R(β)) occurs via the S(N)icB(2) mechanism. In this mechanism, the oxyanion at the C2 attacks the anomeric carbon in a nucleophilic manner to afford 1,2-anhydride intermediate and then the oxyanion at the C6 attacks the anomeric carbon to afford levoglucosan. The activation barrier is much lower in the reaction of R(β) (ΔG(0++) = 25.6 kcal/mol and E(a) = 26.5 kcal/mol) than in the reaction of R(α) (ΔG(0++) = 38.1 kcal/mol and E(a) = 37.2 kcal/mol), which is consistent with the experimental observation that β-glucoside is generally much more reactive than the corresponding α-glucoside. The lower activation barrier of the reaction of R(β) arises from the stereoelectronic effect, which is induced by the charge transfer from the ring oxygen to the anomeric carbon, and the staggered conformation around the C1-C2 bond. When the stereoelectronic effect is absent, the degradation needs larger activation energy; for instance, the degradation of phenyl 5a-carba-β-d-glucoside (R(Cβ)) occurs with large ΔG(0++) and E(a) values like those of α-glucosides, because the methylene group of R(Cβ) does not contribute to the stereoelectronic effect. Also, the conformation around the C1-C2 bond is staggered in the transition state of the R(β) reaction but eclipsed in that of the R(α) reaction, which also leads to the larger reactivity of R(β).     

Fundamental reaction mechanism for cocaine hydrolysis in human butyrylcholinesterase

Butyrylcholinesterase (BChE)-cocaine binding and the fundamental pathway for BChE-catalyzed hydrolysis of cocaine have been studied by molecular modeling, molecular dynamics (MD) simulations, and ab initio calculations. Modeling and simulations indicate that the structures of the prereactive BChE/substrate complexes for (-)-cocaine and (+)-cocaine are all similar to that of the corresponding prereactive BChE/butyrylcholine (BCh) complex. The overall binding of BChE with (-)-cocaine and (+)-cocaine is also similar to that proposed with butyrylthiocholine and succinyldithiocholine, i.e., (-)- or (+)-cocaine first slides down the substrate-binding gorge to bind to Trp-82 and stands vertically in the gorge between Asp-70 and Trp-82 (nonprereactive complex) and then rotates to a position in the catalytic site within a favorable distance for nucleophilic attack and hydrolysis by Ser-198 (prereactive complex). In the prereactive complex, cocaine lies horizontally at the bottom of the gorge. The fundamental catalytic hydrolysis pathway, consisting of acylation and deacylation stages similar to those for ester hydrolysis by other serine hydrolases, was proposed on the basis of the simulated prereactive complex and confirmed theoretically by ab initio reaction coordinate calculations. Both the acylation and deacylation follow a double-proton-transfer mechanism. The calculated energetic results show that within the chemical reaction process the highest energy barrier and Gibbs free energy barrier are all associated with the first step of deacylation. The calculated ratio of the rate constant (k(cat)) for the catalytic hydrolysis to that (k(0)) for the spontaneous hydrolysis is approximately 9.0 x 10(7). The estimated k(cat)/k(0) value of approximately 9.0 x 10(7) is in excellent agreement with the experimentally derived k(cat)/k(0) value of approximately 7.2 x 10(7) for (+)-cocaine, whereas it is approximately 2000 times larger than the experimentally derived k(cat)/k(0) value of approximately 4.4 x 10(4) for (-)-cocaine. All of the results suggest that the rate-determining step of the BChE-catalyzed hydrolysis of (+)-cocaine is the first step of deacylation, whereas for (-)-cocaine the change from the nonprereactive complex to the prereactive complex is rate-determining and has a Gibbs free energy barrier higher than that for the first step of deacylation by approximately 4 kcal/mol. A further analysis of the structural changes from the nonprereactive complex to the prereactive complex reveals specific amino acid residues hindering the structural changes, providing initial clues for the rational design of BChE mutants with improved catalytic activity for (-)-cocaine.     

Understanding the participation of quadricyclane as nucleophile in polar [2sigma + 2sigma + 2pi] cycloadditions toward electrophilic pi molecules

The formal [2sigma + 2sigma + 2pi] cycloaddition of quadricyclane, 1, with dimethyl azodicarboxylate, 2, in water has been studied using DFT methods at the B3LYP/6-31G** and MPWB1K/6-31G** levels. In the gas phase, the reaction of 1 with 2 has a two-stage mechanism with a large polar character and an activation barrier of 23.2 kcal/mol. Inclusion of water through a combined discrete-continuum model changes the mechanism to a two-step model where the first nucleophilic attack of 1 to 2 is the rate-limiting step with an activation barrier of 14.7 kcal/mol. Analysis of the electronic structure of the transition state structures points out the large zwitterionic character of these species. A DFT analysis of the global electrophilicity and nucleophilicity of the reagents provides a sound explanation about the participation of 1 as a nucleophile in these cycloadditions. This behavior is reinforced by a further study of the reaction of 1 with 1,1-dicyanoethylene.     

On the origin of the stabilization of the zwitterionic resting state of cysteine proteases: a theoretical study

Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.     

A multilayered representation,quantum mechanical and molecular mechanics study of the CH3F + OH− reaction in water

The bimolecular nucleophilic substitution (S_N2) reaction of CH₃F + OH^? in aqueous solution was investigated using a combined quantum mechanical and molecular mechanics approach. Reactant complex, transition state, and product complex along the reaction pathway were analyzed in water. The potentials of mean force were calculated using a multilayered representation with the DFT and CCSD(T) level of theory for the reactive region. The obtained free energy activation barrier for this reaction at the CCSD(T)/MM representation is 18.3 kcal/mol which agrees well with the experimental value at ～21.6 kcal/mol. Both the solvation effect and solute polarization effect play key roles on raising the activation barrier height in aqueous solution, with the former raising the barrier height by 3.1 kcal/mol, the latter 1.5 kcal/mol. © 2013 Wiley Periodicals, Inc.     

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.     

Quantum mechanical and molecular mechanics approach with a multilayered-quantum representation study of solvent effects and potentials of mean force for the CH<Subscript>3</Subscript>CH<Subscript>2</Subscript>Cl + ClO<Superscript>−</Superscript> S<Subscript>N</Subscript>2 reaction in aqueous solution

The bimolecular nucleophilic substitution reaction of CH₃CH₂Cl + ClO^? in aqueous solution was investigated using a multilayered-quantum representation, quantum mechanical and molecular mechanics approach with an explicit water model. Ten configurations along the reaction pathway including reactant complex, transition state and product complex were analyzed in the presence of the aqueous solution. The obtained free energy activation barrier under the CCSD(T)/MM representation is 13.2 kcal/mol, while it is 11.7 kcal/mol under the DFT/MM representation which agrees very well with the DFT calculation, at 11.0 kcal/mol, with a polarizable continuum solvent model. The solvent effects including the solvation free energy contribution and the polarization effect raise the free activation barrier by 9.8 kcal/mol. The rate constant, at 298 K, is 5.27 × 10^?17 cm³/molecule/s which is about seven orders of magnitude smaller than that in the gas phase (1.10 × 10^?10 cm³/molecule/s). All in all, the aqueous solution plays an essential role in shaping the reaction pathway for this reaction in water.     

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.     

Effect of electrostatic interaction on the mechanism of dehalogenation catalyzed by haloalkane dehalogenase

Theoretical calculation has been carried out for the nucleophilic displacement reaction of 1,2‐dichloroethane catalyzed by haloalkane dehalogenase. The results indicate that different hydrogen bond patterns of the oxyanion hole and the halide‐stabilizing residues play an important role in the dehalogenation reaction. They cause concertedly an earlier transition state (TS) with the activation barrier of 16.60 kcal/mol. The stabilization effect of Trp125 and Trp175 on chlorine atom in the TS is larger than that of the reactant complex by 15.67 kcal/mol so that, they make contribution to the stabilization of the TS. Moreover, the reaction shows the enzymatic action can be attributed to a combination of reactant‐state destabilization and transition‐state electrostatic stabilization. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011