Molecular dynamics simulation of cocaine binding with human butyrylcholinesterase and its mutants

Molecular dynamics (MD) simulations were carried out to study cocaine binding with wild-type human butyrylcholinesterase (BChE) and its mutants based on a recently reported X-ray crystal structure of human BChE. For each BChE-cocaine system, we simulated both the nonprereactive and prereactive complexes in water. Despite the significant difference found at the acyl binding pocket, the simulated structures confirm the fundamental structural and mechanistic insights obtained from earlier computational studies of wild-type BChE with cocaine based on a homology model, e.g. the rate-determining step for BChE-catalyzed hydrolysis of biologically active (-)-cocaine is the (-)-cocaine rotation in the active site from the nonprereactive BChE-(-)-cocaine complex to the prereactive complex. It has been demonstrated that the MD simulations on both the nonprereactive and prereactive BChE-cocaine complexes can clearly reveal whether specific mutations produce the desired BChE-(-)-cocaine binding structures in which the (-)-cocaine rotation is less hindered while the required prereactive BChE-(-)-cocaine binding is maintained. Based on the MD simulations, both A328W/Y332A and A328W/Y332G BChE's are expected to have catalytic activity for (-)-cocaine hydrolysis higher than that of wild-type BChE and the activity of A328W/Y332G BChE should be slightly higher than that of A328W/Y332A BChE due to the less-hindered (-)-cocaine rotation in the mutant BChE's. However, the less-hindered (-)-cocaine rotation is only a necessary condition for a higher activity mutant BChE. The (-)-cocaine rotation is also less hindered in A328W/Y332A/Y419S BChE, but (-)-cocaine binds with A328W/Y332A/Y419S BChE in a way that is not suitable for the catalysis. Thus, A328W/Y332A/Y419S BChE is expected to lose the catalytic activity. The computational predictions were confirmed by our experimental kinetic data, demonstrating that the MD simulation-based computational protocol used in this study is reliable in prediction of the catalytic activity of BChE mutants for (-)-cocaine hydrolysis.     

Modeling effects of oxyanion hole on the ester hydrolysis catalyzed by human cholinesterases

Molecular dynamics (MD) simulations and hydrogen bonding energy (HBE) calculations have been performed on the prereactive enzyme-substrate complexes (ES), transition states (TS1), and intermediates (INT1) for acetylcholinesterase (AChE)-catalyzed hydrolysis of acetylcholine (ACh), butyrylcholinesterase (BChE)-catalyzed hydrolysis of ACh, and BChE-catalyzed hydrolysis of (+)/(-)-cocaine to examine the protein environmental effects on the catalytic reactions. The hydrogen bonding of cocaine with the oxyanion hole of BChE is found to be remarkably different from that of ACh with AChE/BChE. Whereas G121/G116, G122/G117, and A204/A199 of AChE/BChE all can form hydrogen bonds with ACh to stabilize the transition state during the ACh hydrolysis, BChE only uses G117 and A199 to form hydrogen bonds with cocaine. The change of the estimated total HBE from ES to TS1 is ca. -5.4/-4.4 kcal/mol for AChE/BChE-catalyzed hydrolysis of ACh and ca. -1.7/-0.8 kcal/mol for BChE-catalyzed hydrolysis of (+)/(-)-cocaine. The remarkable difference of approximately 3 to 5 kcal/mol reveals that the oxyanion hole of AChE/BChE can lower the energy barrier of the ACh hydrolysis significantly more than that of BChE for the cocaine hydrolysis. These results help to understand why the catalytic activity of AChE against ACh is considerably higher than that of BChE against cocaine and provides valuable clues on how to improve the catalytic activity of BChE against cocaine.     

Rational design of an enzyme mutant for anti-cocaine therapeutics

(-)-Cocaine is a widely abused drug and there is no available anti-cocaine therapeutic. The disastrous medical and social consequences of cocaine addiction have made the development of an effective pharmacological treatment a high priority. An ideal anti-cocaine medication would be to accelerate (-)-cocaine metabolism producing biologically inactive metabolites. The main metabolic pathway of cocaine in body is the hydrolysis at its benzoyl ester group. Reviewed in this article is the state-of-the-art computational design of high-activity mutants of human butyrylcholinesterase (BChE) against (-)-cocaine. The computational design of BChE mutants have been based on not only the structure of the enzyme, but also the detailed catalytic mechanisms for BChE-catalyzed hydrolysis of (-)-cocaine and (+)-cocaine. Computational studies of the detailed catalytic mechanisms and the structure-and-mechanism-based computational design have been carried out through the combined use of a variety of state-of-the-art techniques of molecular modeling. By using the computational insights into the catalytic mechanisms, a recently developed unique computational design strategy based on the simulation of the rate-determining transition state has been employed to design high-activity mutants of human BChE for hydrolysis of (-)-cocaine, leading to the exciting discovery of BChE mutants with a considerably improved catalytic efficiency against (-)-cocaine. One of the discovered BChE mutants (i.e., A199S/S287G/A328W/Y332G) has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the computational design and discovery effort demonstrates that the unique computational design approach based on the transition-state simulation is promising for rational enzyme redesign and drug discovery.     

Structure-and-mechanism-based design and discovery of therapeutics for cocaine overdose and addiction

(-)-Cocaine is a widely abused drug and there is currently no available anti-cocaine therapeutic. Promising agents, such as anti-cocaine catalytic antibodies and high-activity mutants of human butyrylcholinesterase (BChE), for therapeutic treatment of cocaine overdose have been developed through structure-and-mechanism-based design and discovery. In particular, a unique computational design strategy based on the modeling and simulation of the rate-determining transition state has been developed and used to design and discover desirable high-activity mutants of BChE. One of the discovered high-activity mutants of BChE has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the structure-and-mechanism-based design and discovery effort demonstrates that the unique computational design approach based on transition state modeling and simulation is promising for rational enzyme redesign and drug discovery. The general approach of the structure-and-mechanism-based design and discovery may be used to design high-activity mutants of any enzyme or catalytic antibody.     

A promiscuous glutathione transferase transformed into a selective thiolester hydrolase

Human glutathione transferase A1-1 (hGST A1-1) can be reengineered by rational design into a catalyst for thiolester hydrolysis with a catalytic proficiency of 1.4 x 10(7) M(-1). The thiolester hydrolase, A216H that was obtained by the introduction of a single histidine residue at position 216 catalyzed the hydrolysis of a substrate termed GSB, a thiolester of glutathione and benzoic acid. Here we investigate the substrate requirements of this designed enzyme by screening a thiolester library. We found that only two thiolesters out of 18 were substrates for A216H. The A216H-catalyzed hydrolysis of GS-2 (thiolester of glutathione and naphthalenecarboxylic acid) exhibits a k(cat) of 0.0032 min(-1) and a KM of 41 microM. The previously reported catalysis of GSB has a k(cat) of 0.00078 min(-1) and KM of 5 microM. The k(cat) for A216H-catalyzed hydrolysis of GS-2 is thus 4.1 times higher than for GSB. The catalytic proficiency (k(cat)/KM)/k(uncat) for GS-2 is 3 x 10(6) M(-1). The promiscuous feature of the wt protein towards a range of different substrates has not been conserved in A216H but we have obtained a selective enzyme with high demands on the substrate.     

Formic acid catalyzed hydrolysis of SO3 in the gas phase: a barrierless mechanism for sulfuric acid production of potential atmospheric importance

Computational studies at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels are performed to explore the changes in reaction barrier height for the gas phase hydrolysis of SO(3) to form H(2)SO(4) in the presence of a single formic acid (FA) molecule. For comparison, we have also performed calculations for the reference reaction involving water assisted hydrolysis of SO(3) at the same level. Our results show that the FA assisted hydrolysis of SO(3) to form H(2)SO(4) is effectively a barrierless process. The barrier heights for the isomerization of the SO(3)···H(2)O···FA prereactive collision complex, which is the rate limiting step in the FA assisted hydrolysis, are found to be respectively 0.59 and 0.08 kcal/mol at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels. This is substantially lower than the ~7 kcal/mol barrier for the corresponding step in the hydrolysis of SO(3) by two water molecules--which is currently the accepted mechanism for atmospheric sulfuric acid production. Simple kinetic analysis of the relative rates suggests that the reduction in barrier height facilitated by FA, combined with the greater stability of the prereactive SO(3)···H(2)O···FA collision complex compared to SO(3)···H(2)O···H(2)O and the rather plentiful atmospheric abundance of FA, makes the formic acid mediated hydrolysis reaction a potentially important pathway for atmospheric sulfuric acid production.     

Accumulation of tetrahedral intermediates in cholinesterase catalysis: a secondary isotope effect study

In a previous communication, kinetic β-deuterium secondary isotope effects were reported that support a mechanism for substrate-activated turnover of acetylthiocholine by human butyrylcholinesterase (BuChE) wherein the accumulating reactant state is a tetrahedral intermediate ( Tormos , J. R. ; et al. J. Am. Chem. Soc. 2005 , 127 , 14538 - 14539 ). In this contribution additional isotope effect experiments are described with acetyl-labeled acetylthiocholines (CL(3)COSCH(2)CH(2)N(+)Me(3); L = H or D) that also support accumulation of the tetrahedral intermediate in Drosophila melanogaster acetylcholinesterase (DmAChE) catalysis. In contrast to the aforementioned BuChE-catalyzed reaction, for this reaction the dependence of initial rates on substrate concentration is marked by pronounced substrate inhibition at high substrate concentrations. Moreover, kinetic β-deuterium secondary isotope effects for turnover of acetylthiocholine depended on substrate concentration, and gave the following: (D3)k(cat)/K(m) = 0.95 ± 0.03, (D3)k(cat) = 1.12 ± 0.02 and (D3)βk(cat) = 0.97 ± 0.04. The inverse isotope effect on k(cat)/K(m) is consistent with conversion of the sp(2)-hybridized substrate carbonyl in the E + A reactant state into a quasi-tetrahedral transition state in the acylation stage of catalysis, whereas the markedly normal isotope effect on k(cat) is consistent with hybridization change from sp(3) toward sp(2) as the reactant state for deacylation is converted into the subsequent transition state. Transition states for Drosophila melanogaster AChE-catalyzed hydrolysis of acetylthiocholine were further characterized by measuring solvent isotope effects and determining proton inventories. These experiments indicated that the transition state for rate-determining decomposition of the tetrahedral intermediate is stabilized by multiple protonic interactions. Finally, a simple model is proposed for the contribution that tetrahedral intermediate stabilization provides to the catalytic power of acetylcholinesterase.     

Beta-D-xylosidase from Selenomonas ruminantium: thermodynamics of enzyme-catalyzed and noncatalyzed reactions

Beta-D-Xylosidase/alpha-L-arabinofuranosidase from Selenomonas ruminantium is the most active enzyme known for catalyzing hydrolysis of 1,4-beta-D: -xylooligosaccharides to D-xylose. Temperature dependence for hydrolysis of 4-nitrophenyl-beta-D-xylopyranoside (4NPX), 4-nitrophenyl-alpha-L-arabinofuranoside (4NPA), and 1,4-beta-D-xylobiose (X2) was determined on and off (k (non)) the enzyme at pH 5.3, which lies in the pH-independent region for k (cat) and k (non). Rate enhancements (k (cat)/k (non)) for 4NPX, 4NPA, and X2 are 4.3 x 10(11), 2.4 x 10(9), and 3.7 x 10(12), respectively, at 25 degrees C and increase with decreasing temperature. Relative parameters k (cat) (4NPX)/k (cat) (4NPA), k (cat) (4NPX)/k (cat) (X2), and (k (cat)/K (m))(4NPX)/(k (cat)/K (m))(X2) increase and (k (cat)/K (m))(4NPX)/(k (cat)/K (m))(4NPA), (1/K (m))(4NPX)/(1/K (m))(4NPA), and (1/K (m))(4NPX)/(1/K (m))(X2) decrease with increasing temperature.     

Purification and characterization of a glycoside hydrolase family 43 beta-xylosidase from Geobacillus thermoleovorans IT-08

The gene encoding a glycoside hydrolase family 43 beta-xylosidase (GbtXyl43A) from the thermophilic bacterium Geobacillus thermoleovorans strain IT-08 was synthesized and cloned with a C-terminal His-tag into a pET29b expression vector. The recombinant gene product termed GbtXyl43A was expressed in Escherichia coli and purified to apparent homogeneity. Michaelis-Menten kinetic parameters were obtained for the artificial substrates p-nitrophenyl-beta-D: -xylopyranose (4NPX) and p-nitrophenyl-alpha-L: -arabinofuranose (4NPA), and it was found that the ratio k (cat)/K (m) 4NPA/k (cat)/K (m) 4NPX was approximately 7, indicting greater catalytic efficiency for 4NP hydrolysis from the arabinofuranose aglycon moiety. Substrate inhibition was observed for the substrates 4-methylumbelliferyl xylopyranoside (muX) and the arabinofuranoside cogener (muA), and the ratio k (cat)/K (m) muA/k (cat)/K (m) muX was approximately 5. The enzyme was competitively inhibited by monosaccharides, with an arabinose K (i) of 6.8 +/- 0.62 mM and xylose K (i) of 76 +/- 8.5 mM. The pH maxima was 5.0, and the enzyme was not thermally stable above 54 degrees C, with a t (1/2) of 35 min at 57.5 degrees C. GbtXyl43A showed a broad substrate specificity for hydrolysis of xylooligosaccharides up to the highest degree of polymerization tested (xylopentaose), and also released xylose from birch and beechwood arabinoxylan.     

Proton inventory studies of alpha-thrombin-catalyzed reactions of substrates with selected P and P' sites

Deuterium kinetic solvent isotope effects for the human alpha-thrombin-catalyzed hydrolysis of (1) substrates with selected P(1)-P(3) sites, Z-Pro-Arg-7-amido-4-methylcoumarin (7-AMC), N-t-Boc-Val-Pro-Arg-7-AMC, Bz-Phe-Val-Arg-4-nitroanilide (pNA), and H-D-Phe-L-Pip-Arg-pNA, are (DOD)k(cat) = (2.8-3.3) +/- 0.1 and (DOD)(k(cat)/K(m)) = (0.8-2.1) +/- 0.1 and (2) internally fluorescence-quenched substrates (a) (AB)Val-Phe-Pro-Arg-Ser-Phe-Arg-Leu-Lys(DNP)-Asp-OH, an optimal sequence, and (b) (AB)Val-Ser-Pro-Arg-Ser-Phe-Gln-Lys(DNP)-Asp-OH, recognition sequence for factor VIII, are (DOD)k(cat) = 2.2 +/- 0.2 and (DOD)(k(cat)/K(m)) = (0.8-0.9) +/- 0.1, at the pL (L = H, D) maximum, 8.4-9.0, and (25.0-26.0) +/- 0.1 degrees C. The most plausible models fitting the partial isotope effect (proton inventory) data have been selected on the basis of lowest values of the reduced chi squared and consistency of fractionation factors at all substrate concentrations, assuming rate-determining acylation. The data for Z-Pro-Arg-7-AMC are consistent with a single-proton bridge at the transition state phi(TS) = 0.39 +/- 0.05 and components for solvent reorganization phi(S) = 0.8 +/- 0.1 and phi(S) = 1.22 for k(cat) and k(cat)/K(m), respectively. The data for tripeptide amides fit bowl-shaped curves; an example is N-t-Boc-Val-Pro-Arg-7-AMC: phi(TS)(1) = phi(TS)(2) = 0.57 +/- 0.01 and phi(S) = 1 for k(cat) and 1.6 +/- 0.1 for k(cat)/K(m). Proton inventories for the nonapeptide (2b) are linear. The data for k(cat) for H-D-Phe-L-Pip-Arg-pNA and the decapeptide (2a) are most consistent with two identical fractionation factors for catalytic proton bridging, phi(TS)(1) = phi(TS)(2) = 0.68 +/- 0.02 and a large inverse component (phi(S) = 3.1 +/- 0.5) for the latter, indicative of substantial solvent reorganization upon leaving group departure. Proton inventory curves for k(cat)/K(m) for nearly all substrates are dome-shaped with an inverse isotope effect component (phi(S) = 1.2-2.4) originating from solvent reorganization during association of thrombin with substrate. These large contributions from medium effects are in full accord with the conformational adjustments required for the fulfillment of the dual, hemostatic and thrombolytic, functions of thrombin.     

Hydrolysis of N-alkyl sulfamates and the catalytic efficiency of an S-N cleaving sulfamidase

The final step in the degradation of heparin sulfate involves the enzymatic hydrolysis of its 2-sulfamido groups. To evaluate the power of the corresponding sulfamidases as catalysts, we examined the reaction of N-neopentyl sulfamate at elevated temperatures and found it to undergo specific acid catalyzed hydrolysis even at alkaline pH. A rate constant of 10(-16) s(-1) was calculated using the Eyring equation for water attack on the N-protonated species at pH 7, 25 °C. As a model for the pH neutral reaction, a rate constant for hydroxide attack on (CH(3))(3)CCH(2)N(+)H(2)SO(3)(-) at pH 7, 25 °C was calculated to be 10(-19) s(-1). The corresponding rate enhancement (k(cat)/k(non)) produced by the N-sulfamidase of F. heparinum is approximately 10(16)-fold, which is somewhat larger than those generated by most hydrolytic enzymes but considerably smaller than those generated by S-O cleaving sulfatases.     

Design of a dinuclear nickel(II) bioinspired hydrolase to bind covalently to silica surfaces: synthesis, magnetism, and reactivity studies

Presented herein is the design of a dinuclear Ni(II) synthetic hydrolase [Ni(2)(HBPPAMFF)(μ-OAc)(2)(H(2)O)]BPh(4) (1) (H(2)BPPAMFF = 2-[(N-benzyl-N-2-pyridylmethylamine)]-4-methyl-6-[N-(2-pyridylmethyl)aminomethyl)])-4-methyl-6-formylphenol) to be covalently attached to silica surfaces, while maintaining its catalytic activity. An aldehyde-containing ligand (H(2)BPPAMFF) provides a reactive functional group that can serve as a cross-linking group to bind the complex to an organoalkoxysilane and later to the silica surfaces or directly to amino-modified surfaces. The dinuclear Ni(II) complex covalently attached to the silica surfaces was fully characterized by different techniques. The catalytic turnover number (k(cat)) of the immobilized Ni(II)Ni(II) catalyst in the hydrolysis of 2,4-bis(dinitrophenyl)phosphate is comparable to the homogeneous reaction; however, the catalyst interaction with the support enhanced the substrate to complex association constant, and consequently, the catalytic efficiency (E = k(cat)/K(M)) and the supported catalyst can be reused for subsequent diester hydrolysis reactions.     

Free energy perturbation (FEP) simulation on the transition states of cocaine hydrolysis catalyzed by human butyrylcholinesterase and its mutants

A novel computational protocol based on free energy perturbation (FEP) simulations on both the free enzyme and transition state structures has been developed and tested to predict the mutation-caused shift of the free energy change from the free enzyme to the rate-determining transition state for human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The calculated shift, denoted by DeltaDeltaG(1 --> 2), of such kind of free energy change determines the catalytic efficiency (kcat/KM) change caused by the simulated mutation transforming enzyme 1 to enzyme 2. By using the FEP-based computational protocol, the DeltaDeltaG(1 --> 2) values for the mutations A328W/Y332A --> A328W/Y332G and A328W/Y332G --> A328W/Y332G/A199S were calculated to be -0.22 and -1.94 kcal/mol, respectively. The calculated DeltaDeltaG(1 --> 2) values predict that the change from the A328W/Y332A mutant to the A328W/Y332G mutant should slightly improve the catalytic efficiency and that the change from the A328W/Y332G mutant to the A328W/Y332G/A199S mutant should significantly improve the catalytic efficiency of the enzyme for the (-)-cocaine hydrolysis. The predicted catalytic efficiency increases are supported by the experimental data showing that kcat/KM = 8.5 x 10(6), 1.4 x 10(7), and 7.2 x 10(7) min(-1) M(-1) for the A328W/Y332A, A328W/Y332G, and A328W/Y332G/A199S mutants, respectively. The qualitative agreement between the computational and experimental data suggests that the FEP simulations may provide a promising protocol for rational design of high-activity mutants of an enzyme. The general computational strategy of the FEP simulation on a transition state can be used to study the effects of a mutation on the activation free energy for any enzymatic reaction.     

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.     

RNA-cleaving DNA enzymes with altered regio- or enantioselectivity

In vitro evolution methods were used to obtain DNA enzymes that cleave either a 2',5'-phosphodiester following a D-ribonucleotide or a 3',5'-phosphodiester following an L-ribonucleotide. Both enzymes can operate in an intermolecular reaction format with multiple turnover. The DNA enzyme that cleaves a 2',5'-phosphodiester exhibits a k(cat) of approximately 0.01 min(-1) and catalytic efficiency, k(cat)/K(m), of approximately 10(8) M(-1) min(-1). The enzyme that cleaves an L-ribonucleotide is about 10-fold slower and has a catalytic efficiency of approximately 4 x 10(5) M(-1) min(-1). Both enzymes require a divalent metal cation for their activity and have optimal catalytic rate at pH 7-8 and 35-50 degrees C. In a comparison of each enzyme's activity with either its corresponding substrate that contains an unnatural ribonucleotide or a substrate that instead contains a standard ribonucleotide, the 2',5'-phosphodiester-cleaving DNA enzyme exhibited a regioselectivity of 6000-fold, while the L-ribonucleotide-cleaving DNA enzyme exhibited an enantioselectivity of 40-fold. These molecules demonstrate how in vitro evolution can be used to obtain regio- and enantioselective catalysts that exhibit specificities for nonnatural analogues of biological compounds.     

The "neutral" hydrolysis of simple carboxylic esters in water and the rate enhancements produced by acetylcholinesterase and other carboxylic acid esterases

Experiments at elevated temperatures permit the determination of rate constant and thermodynamic activation parameters for the neutral hydrolysis of the neurotransmitter acetylcholine in water. At 25 °C, the extrapolated rate constant for the uncatalyzed (or neutral) hydrolysis of acetylcholine is 3.9 × 10(-7) s(-1) at 25 °C (ΔH(?) = 20.0 kcal/mol; TΔS(?) = -6.1 kcal/mol). Acetylcholine is more susceptible to neutral and base-catalyzed hydrolysis than ethyl acetate but less susceptible to acid-catalyzed hydrolysis. For acetylcholinesterase from the electric eel, the catalytic proficiency [(k(cat)/K(m))/k(neutral)] is 2 × 10(16) M(-1), comparable in magnitude with the catalytic proficiencies of aminohydrolases that act on peptides and nucleosides.     

Supramolecular catalysis of orthoformate hydrolysis in basic solution: an enzyme-like mechanism

A water-soluble self-assembled supramolecular host molecule catalyzes the hydrolysis of orthoformates in basic solution. Comparison of the rate constants of the catalyzed and uncatalyzed reactions for hydrolysis displays rate accelerations of up to 3900 for tri- n-propyl orthoformate. Kinetic analysis shows that the mechanism of hydrolysis with the supramolecular host obeys the Michaelis-Menten model. Mechanistic studies, including (13)C-labeling experiments, revealed that the resting state of the catalytic system is the neutral substrate encapsulated in the host. Activation parameters for the k cat step of the reaction revealed that upon substrate encapsulation in the assembly, the entropy of activation becomes more negative in contrast to the uncatalyzed reaction. Furthermore, solvent isotope effects reveal a normal k(H 2O)/ k(D 2O) = 1.6, confirming an A-S E2 mechanism in which proton transfer occurs in the rate-limiting step. This is in contrast with the A1 mechanism of the uncatalyzed reaction in which decomposition of the protonated substrate is rate-limiting.     

Decay of the peroxide intermediate in laccase: reductive cleavage of the O-O bond.

Laccase is a multicopper oxidase that contains four Cu ions, one type 1, one type 2, and a coupled binuclear type 3 Cu pair. The type 2 and type 3 centers form a trinuclear Cu cluster that is the active site for O(2) reduction to H(2)O. To examine the reaction between the type 2/type 3 trinuclear cluster and dioxygen, the type 1 Cu was removed and replaced with Hg(2+), producing the T1Hg derivative. When reduced T1Hg laccase is reacted with dioxygen, a peroxide intermediate (P) is formed. The present study examines the kinetics and mechanism of formation and decay of P in T1HgLc. The formation of P was found to be independent of pH and did not involve a kinetic solvent isotope effect, indicating that no proton is involved in the rate-determining step of formation of P. Alternatively, pH and isotope studies on the decay of P revealed that a proton enhances the rate of decay by 10-fold at low pH. This process shows an inverse k(H)/k(D) kinetic solvent isotope effect and involves protonation of a nearby residue that assists in catalysis, rather than direct protonation of the peroxide. Decay of P also involves a significant oxygen isotope effect (k(16)O(2)/k(18)O(2)) of 1.11 +/- 0.05, indicating that reductive cleavage of the O-O bond is the rate-determining step in the decay of P. The activation energy for this process was found to be approximately 9.0 kcal/mol. The exceptionally slow rate of decay of P is explained by the fact that this process involves a 1e(-) reductive cleavage of the O-O bond and there is a large Franck-Condon barrier associated with this process. Alternatively, the 2e(-) reductive cleavage of the O-O bond has a much larger driving force which minimizes this barrier and accelerates the rate of this reaction by approximately 10(7) in the native enzyme. This large difference in rate for the 2e(-) versus 1e(-) process supports a molecular mechanism for multicopper oxidases in which O(2) is reduced to H(2)O in two 2e(-) steps.     

Aristolochene synthase: mechanistic analysis of active site residues by site-directed mutagenesis

Incubation of farnesyl diphosphate (1) with Penicillium roqueforti aristolochene synthase yielded (+)-aristolochene (4), accompanied by minor quantities of the proposed intermediate (S)-(-)germacrene A (2) and the side-product (-)-valencene (5) in a 94:4:2 ratio. By contrast, the closely related aristolochene synthase from Aspergillus terreus cyclized farnesyl diphosphate only to (+)-aristolochene (4). Site-directed mutagenesis of amino acid residues in two highly conserved Mg(2+)-binding domains led in most cases to reductions in both k(cat) and k(cat)/K(m) as well as increases in the proportion of (S)-(-)germacrene A (2), with the E252Q mutant of the P. roqueforti aristolochene synthase producing only (-)-2. The P. roqueforti D115N, N244L, and S248A/E252D mutants were inactive, as was the A. terreus mutant E227Q. The P. roqueforti mutant Y92F displayed a 100-fold reduction in k(cat) that was offset by a 50-fold decrease in K(m), resulting in a relatively minor 2-fold decrease in catalytic efficiency, k(cat)/K(m). The finding that Y92F produced (+)-aristolochene (4) as 81% of the product, accompanied by 7% 5 and 12% 2, rules out Tyr-92 as the active site Lewis acid that is responsible for protonation of the germacrene A intermediate in the formation of aristolochene (4).     

Most efficient cocaine hydrolase designed by virtual screening of transition states

Cocaine is recognized as the most reinforcing of all drugs of abuse. There is no anticocaine medication available. The disastrous medical and social consequences of cocaine addiction have made the development of an anticocaine medication a high priority. It has been recognized that an ideal anticocaine medication is one that accelerates cocaine metabolism producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., cocaine hydrolysis catalyzed by plasma enzyme butyrylcholinesterase (BChE). However, wild-type BChE has a low catalytic efficiency against the abused cocaine. Design of a high-activity enzyme mutant is extremely challenging, particularly when the chemical reaction process is rate-determining for the enzymatic reaction. Here we report the design and discovery of a high-activity mutant of human BChE by using a novel, systematic computational design approach based on transition-state simulations and activation energy calculations. The novel computational design approach has led to discovery of the most efficient cocaine hydrolase, i.e., a human BChE mutant with an approximately 2000-fold improved catalytic efficiency, promising for therapeutic treatment of cocaine overdose and addiction as an exogenous enzyme in human. The encouraging discovery resulted from the computational design not only provides a promising anticocaine medication but also demonstrates that the novel, generally applicable computational design approach is promising for rational enzyme redesign and drug discovery.