Active Site Water Molecules Research Articles

Novel product chemistry from mechanistic analysis of ent-copalyl diphosphate synthases from plant hormone biosynthesis.

An active-site water molecule coordinated by conserved histidine and asparagine residues seems to serve as the catalytic base in all ent-copalyl diphosphate synthases (CPSs). When these residues are substituted by alanine, the mutant CPSs produce stereochemically novel ent-8-hydroxy-CPP. Given the requisite presence of CPSs in all land plants for gibberellin phytohormone biosynthesis, such plasticity presumably underlies the observed extensive diversification of the resulting labdane-related diterpenoids.

Arabinose 5-phosphate covalently inhibits transaldolase

Arabinose 5-phosphate (A5P) is the aldopentose version of the ketohexose fructose 6-phosphate (F6P), having identical stereochemistry but lacking atoms corresponding to the 1-carbon and 1-hydroxyl. Despite structural similarity and conservation of the reactive portion of F6P, F6P acts as a substrate whereas A5P is reported to be an inhibitor of transaldolase. To address the lack of A5P reactivity we determined a crystal structure of the Francisella tularensis transaldolase in complex with A5P. This structure reveals that like F6P, A5P forms a covalent Schiff base with active site Lys135. Unlike F6P, A5P binding fails to displace an ordered active site water molecule. Retaining this water necessitates conformational changes at the A5P-protein linkage that possibly hinder reactivity. The findings presented here show the basis of A5P inhibition and suggest an unusual mechanism of competitive, reversible-covalent transaldolase regulation.

Rapid decomposition and visualisation of protein-ligand binding free energies by residue and by water.

Recent advances in computational hardware, software and algorithms enable simulations of protein-ligand complexes to achieve timescales during which complete ligand binding and unbinding pathways can be observed. While observation of such events can promote understanding of binding and unbinding pathways, it does not alone provide information about the molecular drivers for protein-ligand association, nor guidance on how a ligand could be optimised to better bind to the protein. We have developed the waterswap (C. J. Woods et al., J. Chem. Phys., 2011, 134, 054114) absolute binding free energy method that calculates binding affinities by exchanging the ligand with an equivalent volume of water. A significant advantage of this method is that the binding free energy is calculated using a single reaction coordinate from a single simulation. This has enabled the development of new visualisations of binding affinities based on free energy decompositions to per-residue and per-water molecule components. These provide a clear picture of which protein-ligand interactions are strong, and which active site water molecules are stabilised or destabilised upon binding. Optimisation of the algorithms underlying the decomposition enables near-real-time visualisation, allowing these calculations to be used either to provide interactive feedback to a ligand designer, or to provide run-time analysis of protein-ligand molecular dynamics simulations.

Conformational Compression and Barrier Height Heterogeneity in theN-Acetylglutamate Kinase

The role of motions in the catalytic cycle of an enzyme is the subject of much debate. Crystallographic results for the enzyme N-acetyl-l-glutamate kinase (NAGK), which is a suitable target for antibacterial drugs, suggest that a conformational compression of the active site favors catalysis. We have used a QM/MM scheme to compute energy profiles of the phosphoryl transfer reaction for 20 conformations of NAGK, starting from four crystal structures that represent different stages of the catalytic process. All paths show a common associative mechanism but with a wide range of barrier heights. The position of several active site residues and water molecules are found to determine the energetic barrier of each conformation, as revealed by principal component and partial least-squares chemometric analyses. In particular, conformations in which the two substrates have a shorter distance separation and a more linear mutual orientation tend to have lower energetic barriers, thus supporting the putative role of conformational compressive motions in catalysis. Interestingly, these motions are the same that lead to opening of the active site, which molecular dynamics simulations indicate is a fast process when the enzyme is free of substrates. Despite the lack of extended sampling, the energy barrier we calculate for the chemical step lies significantly below the apparent energetic barrier derived from experiment. Although not conclusive, this result supports a previous hypothesis, also derived from experiment, that conformational motions, rather than the chemical step, are rate limiting.

Identification of catalytically important amino acid residues for enzymatic reduction of glyoxylate in plants

NADPH-dependent glyoxylate reductases from Arabidopsis thaliana (AtGLYR) convert both glyoxylate and succinic semialdehyde into their corresponding hydroxyacid equivalents. The primary sequence of cytosolic AtGLYR1 reveals several sequence elements that are consistent with the β-HAD (β-hydroxyacid dehydrogenase) protein family, whose members include 3-hydroxyisobutyrate dehydrogenase, tartronate semialdehyde reductase and 6-phosphogluconate dehydrogenase. Here, site-directed mutagenesis was utilized to identify catalytically important amino acid residues for glyoxylate reduction in AtGLYR1. Kinetic studies and binding assays established that Lys170 is essential for catalysis, Phe231, Asp239, Ser121 and Thr95 are more important in substrate binding than in catalysis, and Asn174 is more important in catalysis. The low activity of the mutant enzymes precluded kinetic studies with succinic semialdehyde. The crystal structure of AtGLYR1 in the absence of substrate was solved to 2.1Å by molecular replacement using a previously unrecognized member of the β-HAD family, cytokine-like nuclear factor, thereby enabling the 3-D structure of the protein to be modeled with substrate and co-factor. Structural alignment of AtGLYR1 with β-HAD family members provided support for the essentiality of Lys170, Phe173, Asp239, Ser121, Asn174 and Thr95 in the active site and preliminary support for an acid/base catalytic mechanism involving Lys170 as the general acid and a conserved active-site water molecule. This information established that AtGLYR1 is a member of the β-HAD protein family. Sequence and activity comparisons indicated that AtGLYR1 and the plastidial AtGLYR2 possess structural features that are absent in Arabidopsis hydroxypyruvate reductases and probably account for their stronger preference for glyoxylate over hydroxypyruvate.

Catalytic mechanism and origin of high activity of cellulase TmCel12A at high temperature: a quantum mechanical/molecular mechanical study

Understanding the factors that determine the catalytic efficiency of cellulases is of considerable importance in cellulosic ethanol production, especially at high temperature. The cellulase 12A from the hyperthermophile Thermotoga maritima (TmCel12A) is a possible candidate for accelerating the rate of hydrolysis via temperature elevation up to as high as 95 °C. However, the details of the catalytic mechanism and origin of the activity of TmCel12A at high temperature have not been well studied. Here, the enzyme-catalyzed reaction is explored using free energy simulations (potential of mean force) with umbrella sampling and quantum mechanical/molecular mechanical (SCC-DFTB/MM) potential at both relatively low (37 °C) and high (85 °C) temperatures. The free energy barriers for glycosylation and deglycosylation are calculated to be 22.5 ± 0.4 and 24.5 ± 0.7 kcal · mol−1 at 85 °C, respectively. The barrier for deglycosylation is found to decrease with increasing temperature or as a result of the Y61 → G mutation, consistent with experimental observations. The transition state for glycosylation and deglycosylation obtained from the simulations is in an oxocarbonium state with the −1 glucose ring having an E3 envelop (or 4H3 half-chair) conformation. A unique characteristic of the TmCel12A structure seems to be the existence of a stable moiety that may play a role in “holding” cellulose at the binding site with the correct orientation for the reaction even at 85 °C. This stable moiety (comprising hydrogen-bonded E116, E134, E227 and an active-site water molecule) may be one of the important factors for the relatively high activity of TmCel12A at high temperature.

Methylated Nω-Hydroxy-l-arginine Analogues as Mechanistic Probes for the Second Step of the Nitric Oxide Synthase-Catalyzed Reaction

Nitric oxide synthase (NOS) catalyzes the conversion of L-arginine to L-citrulline through the intermediate N(ω)-hydroxy-L-arginine (NHA), producing nitric oxide, an important mammalian signaling molecule. Several disease states are associated with improper regulation of nitric oxide production, making NOS a therapeutic target. The first step of the NOS reaction has been well-characterized and is presumed to proceed through a compound I heme species, analogous to the cytochrome P450 mechanism. The second step, however, is enzymatically unprecedented and is thought to occur via a ferric peroxo heme species. To gain insight into the details of this unique second step, we report here the synthesis of NHA analogues bearing guanidinium methyl or ethyl substitutions and their investigation as either inhibitors of or alternate substrates for NOS. Radiolabeling studies reveal that N(ω)-methoxy-L-arginine, an alternative NOS substrate, produces citrulline, nitric oxide, and methanol. On the basis of these results, we propose a mechanism for the second step of NOS catalysis in which a methylated nitric oxide species is released and is further metabolized by NOS. Crystal structures of our NHA analogues bound to nNOS have been determined, revealing the presence of an active site water molecule only in the presence of singly methylated analogues. Bulkier analogues displace this active site water molecule; a different mechanism is proposed in the absence of the water molecule. Our results provide new insights into the steric and stereochemical tolerance of the NOS active site and substrate capabilities of NOS.

Docking challenge: protein sampling and molecular docking performance.

Computational tools are essential in the drug design process, especially in order to take advantage of the increasing numbers of solved X-ray and NMR protein-ligand structures. Nowadays, molecular docking methods are routinely used for prediction of protein-ligand interactions and to aid in selecting potent molecules as a part of virtual screening of large databases. The improvements and advances in computational capacity in the past decade have allowed for further developments in molecular docking algorithms to address more complicated aspects such as protein flexibility. The effects of incorporation of active site water molecules and implicit or explicit solvation of the binding site are other relevant issues to be addressed in the docking procedures. Using the right docking algorithm at the right stage of virtual screening is most important. We report a staged study to address the effects of various aspects of protein flexibility and inclusion of active site water molecules on docking effectiveness to retrieve (and to be able to predict) correct ligand poses and to rank docked ligands in relation to their biological activity for CHK1, ERK2, LpxC, and UPA. We generated multiple conformers for the ligand and compared different docking algorithms that use a variety of approaches to protein flexibility, including rigid receptor, soft receptor, flexible side chains, induced fit, and multiple structure algorithms. Docking accuracy varied from 1% to 84%, demonstrating that the choice of method is important.

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPDclosed ) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPDopen ) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPDclosed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPDopen conformation. When the active site water molecule network that is a part of the free WPDclosed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design.

Calculation of Vibrational Shifts of Nitrile Probes in the Active Site of Ketosteroid Isomerase upon Ligand Binding

The vibrational Stark effect provides insight into the roles of hydrogen bonding, electrostatics, and conformational motions in enzyme catalysis. In a recent application of this approach to the enzyme ketosteroid isomerase (KSI), thiocyanate probes were introduced in site-specific positions throughout the active site. This paper implements a quantum mechanical/molecular mechanical (QM/MM) approach for calculating the vibrational shifts of nitrile (CN) probes in proteins. This methodology is shown to reproduce the experimentally measured vibrational shifts upon binding of the intermediate analogue equilinen to KSI for two different nitrile probe positions. Analysis of the molecular dynamics simulations provides atomistic insight into the roles that key residues play in determining the electrostatic environment and hydrogen-bonding interactions experienced by the nitrile probe. For the M116C-CN probe, equilinen binding reorients an active-site water molecule that is directly hydrogen-bonded to the nitrile probe, resulting in a more linear C≡N--H angle and increasing the CN frequency upon binding. For the F86C-CN probe, equilinen binding orients the Asp103 residue, decreasing the hydrogen-bonding distance between the Asp103 backbone and the nitrile probe and slightly increasing the CN frequency. This QM/MM methodology is applicable to a wide range of biological systems and has the potential to assist in the elucidation of the fundamental principles underlying enzyme catalysis.

Screening of inhibitors for S130G inhibitor resistant mutants of TEM type beta-lactamase.

Bacteria are remarkably adaptable organisms that acquire an almost limitless competence to survive under unpleasant conditions. The drastic emergence of antibiotic resistance among β-Lactamases is the most serious threat to hospitals and nosocomial settings. β-lactam inhibitors came into existence in order to overcome the problem of antibibiotic resistance in bacteria. The emergence of inhibitor resistant mutants has raised the alarms. In this study we have used structured based virtual screening approach and have screened out some inhibitors against S130G TEM mutant. All the compounds were tested in presence and absence of conserved active site water molecules. These compounds were found be showing much higher efficacy than known β-lactamase inhibitors. Amino acids G130, S70, N132, G130, Y105 and V216 were found crucial for the interaction of inhibitors within the active site.

The role of active-site Phe87 in modulating the organic co-solvent tolerance of cytochrome P450 BM3 monooxygenase.

Understanding the effects of organic co-solvents on protein structure and function is pivotal to engineering enzymes for biotransformation in non-aqueous solvents. The effects of DMSO on the catalytic activity of cytochrome P450 BM3 have previously been investigated and the importance of Phe87 in its organic co-solvent tolerance was identified. To probe the DMSO inactivation mechanism and the functional role of Phe87 in modulating the organic co-solvent tolerance of P450 BM3, the haem domain (Thr1-Leu455) of the F87A variant was cocrystallized in the presence of 14%(v/v) and 28%(v/v) DMSO. At both DMSO concentrations the protein retained the canonical structure of the P450 haem domain without any sign of partial or global unfolding. Interestingly, a DMSO molecule was found in the active site of both structures, with its O atom pointing towards the haem iron. The orientation of the DMSO molecule indicated a dynamic coordination process that was in competition with the active-site water molecule. The ability of the DMSO molecule to coordinate the haem iron is plausibly the main reason why P450 BM3 is inactivated at elevated DMSO concentrations. The data allowed an interesting comparison with the wild-type structures reported previously. A DMSO molecule was found when the wild-type protein was placed in 28%(v/v) DMSO, in which the DMSO molecule coordinated the haem iron directly via its S atom. Intriguingly, no DMSO molecule was observed at 14%(v/v) DMSO for the wild-type structure. These results suggested that the bulky phenyl side chain of Phe87 protects the haem from being accessed by the DMSO molecule and explains the higher tolerance of the wild-type enzyme towards organic co-solvents compared with its F87A variant.

The Catalytic Mechanism of HIV-1 Integrase for DNA 3′-End Processing Established by QM/MM Calculations

The development of HIV-1 integrase (INT) inhibitors has been hampered by incomplete structural and mechanistic information. Despite the efforts made to overcome these limitations, only one compound has been approved for clinical use so far. In this work, we have used all experimental information available for INT and similar enzymes, to build a model of the holo-integrase:DNA complex that includes an entire central core domain, a ssDNA GCAGT substrate, and two magnesium ions. Subsequently, we used a large array of computational techniques, which included molecular dynamics, thermodynamic integration, and high-level quantum mechanics/molecular mechanics (QM/MM) calculations to study the possible pathways for the mechanism of 3' end processing catalyzed by INT. We found that the only viable mechanism to hydrolyze the DNA substrate is a nucleophilic attack of an active site water molecule to the phosphorus atom of the scissile phosphoester bond, with the attacking water being simultaneously deprotonated by an Mg(2+)-bound hydroxide ion. The unstable leaving oxoanion is protonated by an Mg(2+)-bound water molecule within the same elementary reaction step. This reaction has an activation free energy of 15.4 kcal/mol, well within the limits imposed by the experimental turnover. This work significantly improves the fundamental knowledge on the integrase chemistry. It can also contribute to the discovery of leads against HIV-1 infection as it provides, for the first time, accurate transition states structures that can be successfully used as templates for high-throughput screening of new INT inhibitors.

Electrostatic Effects of Mutations of Ras Glutamine 61 Measured Using Vibrational Spectroscopy of a Thiocyanate Probe

Mutations of human oncoprotein p21(Ras) (hereafter Ras) at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.

N–O bond cleavage mechanism(s) in nitrous oxide reductase

Quantum chemical calculations of active-site models of nitrous oxide reductase (N(2)OR) have been undertaken to elucidate the mechanism of N-O bond cleavage mediated by the supported tetranuclear Cu(4)S core (Cu(Z)) found in the enzymatic active site. Using either a minimal model previously employed by Gorelsky et al. (J. Am. Chem. Soc. 128:278-290, 2006) or a more extended model including key residue side chains in the active-site second shell, we found two distinct mechanisms. In the first model, N(2)O binds to the fully reduced Cu(Z) in a bent μ-(1,3)-O,N bridging fashion between the Cu(I) and Cu(IV) centers and subsequently extrudes N(2) while generating the corresponding bridged μ-oxo species. In the second model, substrate N(2)O binds loosely to one of the coppers of Cu(Z) in a terminal fashion, i.e., using only the oxygen atom; loss of N(2) generates the same μ-oxo copper core. The free energies of activation predicted for these two alternative pathways are sufficiently close to one another that theory does not provide decisive support for one over the other, posing an interesting problem with respect to experiments that might be designed to distinguish between the two. Effects of nearby residues and active-site water molecules are also explored.

Monoamine oxidase inhibition by selected anilide derivatives

A series of anilide derivatives were synthesized and evaluated as inhibitors of recombinant human monoamine oxidase (MAO) A and B. The most potent inhibitors among the derivatives that were initially evaluated were (2 E)-N-(3-chlorophenyl)-3-phenylprop-2-enamide ( 2c) and (2 E)-N-(3-bromophenyl)-3-phenylprop-2-enamide ( 2d) with IC 50 values of 0.53 μM and 0.45 μM, respectively. These derivatives exhibited reversible and selective inhibition of MAO-B with binding affinities 37 fold higher for MAO-B than for MAO-A. Analysis of the possible binding interactions of these inhibitors with active site models of human MAO-A and –B led to the design of phenolic and benzonitrile derivatives of 2c and 2d. Among these were (2E)-N-(3-chlorophenyl)-3-(4-hydroxyphenyl)prop-2-enamide ( 7c) and (2E)-N-(3-bromophenyl)-3-(4-hydroxyphenyl)prop-2-enamide ( 7d) which inhibited MAO-B selectively and reversibly with IC 50 values of 0.032 μM and 0.026 μM, respectively. These inhibitors were at least 14 fold more potent than 2c and 2d. This study concludes that N,3-diphenylprop-2-enamide is a suitable scaffold for the design of selective MAO-B inhibitors and structural modifications to enhance the binding affinities of the inhibitors for the MAO-B active site include substitution with halogens on the N-phenyl ring and substitution with hydroxyl and nitrile functional groups on the para and meta positions, respectively, of the C3 phenyl ring. Possible binding modes of these structures within the MAO-B active site are proposed with the emphasis on the interactions of the inhibitor halogens and the hydroxyl and nitrile functional groups with active site residues and water molecules.

Water in the Active Site of Ketosteroid Isomerase

Classical molecular dynamics simulations were utilized to investigate the structural and dynamical properties of water in the active site of ketosteroid isomerase (KSI) to provide insight into the role of these water molecules in the enzyme-catalyzed reaction. This reaction is thought to proceed via a dienolate intermediate that is stabilized by hydrogen bonding with residues Tyr16 and Asp103. A comparative study was performed for the wild-type (WT) KSI and the Y16F, Y16S, and Y16F/Y32F/Y57F (FFF) mutants. These systems were studied with three different bound ligands: equilenin, which is an intermediate analog, and the intermediate states of two steroid substrates. Several distinct water occupation sites were identified in the active site of KSI for the WT and mutant systems. Three additional sites were identified in the Y16S mutant that were not occupied in WT KSI or the other mutants studied. The number of water molecules directly hydrogen bonded to the ligand oxygen was approximately two in the Y16S mutant and one in the Y16F and FFF mutants, with intermittent hydrogen bonding of one water molecule in WT KSI. The molecular dynamics trajectories of the Y16F and FFF mutants reproduced the small conformational changes of residue 16 observed in the crystal structures of these two mutants. Quantum mechanical/molecular mechanical calculations of (1)H NMR chemical shifts of the protons in the active site hydrogen-bonding network suggest that the presence of water in the active site does not prevent the formation of short hydrogen bonds with far-downfield chemical shifts. The molecular dynamics simulations indicate that the active site water molecules exchange much more frequently for WT KSI and the FFF mutant than for the Y16F and Y16S mutants. This difference is most likely due to the hydrogen-bonding interaction between Tyr57 and an active site water molecule that is persistent in the Y16F and Y16S mutants but absent in the FFF mutant and significantly less probable in WT KSI.

The Mechanism of Peptidyl Transfer Catalysis by the Ribosome

The ribosome catalyzes two fundamental biological reactions: peptidyl transfer, the formation of a peptide bond during protein synthesis, and peptidyl hydrolysis, the release of the complete protein from the peptidyl tRNA upon completion of translation. The ribosome is able to utilize and distinguish the two different nucleophiles for each reaction, the α-amine of the incoming aminoacyl tRNA versus the water molecule. The correct binding of substrates induces structural rearrangements of ribosomal active-site residues and the substrates themselves, resulting in an orientation suitable for catalysis. In addition, active-site residues appear to provide further assistance by ordering active-site water molecules and providing an electrostatic environment via a hydrogen network that stabilizes the reaction intermediates and possibly shuttles protons. Major questions remain concerning the timing, components, and mechanism of the proton transfer steps. This review summarizes the recent progress in structural, biochemical, and computational advances and presents the current mechanistic models for these two reactions.

SET7/9 Catalytic Mutants Reveal the Role of Active Site Water Molecules in Lysine Multiple Methylation

SET domain lysine methyltransferases (KMTs) methylate specific lysine residues in histone and non-histone substrates. These enzymes also display product specificity by catalyzing distinct degrees of methylation of the lysine ε-amino group. To elucidate the molecular mechanism underlying this specificity, we have characterized the Y245A and Y305F mutants of the human KMT SET7/9 (also known as KMT7) that alter its product specificity from a monomethyltransferase to a di- and a trimethyltransferase, respectively. Crystal structures of these mutants in complex with peptides bearing unmodified, mono-, di-, and trimethylated lysines illustrate the roles of active site water molecules in aligning the lysine ε-amino group for methyl transfer with S-adenosylmethionine. Displacement or dissociation of these solvent molecules enlarges the diameter of the active site, accommodating the increasing size of the methylated ε-amino group during successive methyl transfer reactions. Together, these results furnish new insights into the roles of active site water molecules in modulating lysine multiple methylation by SET domain KMTs and provide the first molecular snapshots of the mono-, di-, and trimethyl transfer reactions catalyzed by these enzymes.

Catalysis of Carboxypeptidase A: Promoted-Water versus Nucleophilic Pathways

The catalytic mechanism of carboxypeptidase A (CPA) for the hydrolysis of ester substrates is investigated using hybrid quantum mechanical/molecular mechanical (QM/MM) methods and high-level density functional theory. The prevailing mechanism was found to utilize an active-site water molecule assisted by Glu270, and this so-called promoted-water pathway is similar to that in the CPA catalyzed proteolytic reaction (D. Xu and H. Guo, J. Am. Chem. Soc. 2009, 131, 9780). On the other hand, our simulations indicated the existence of an alternative pathway due to direct nucleophilic attack of Glu270 on the scissile carbonyl carbon. This so-called nucleophilic pathway, which is not viable in proteolytic reactions, leads to a stable acyl-enzyme complex. However, the nucleophilic pathway is nonproductive as it is blocked by a high barrier in the deacylation step. On the basis of results reported here and in our earlier publication, a unified model is proposed to account for nearly all experimental observations concerning the catalysis of CPA.