Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step

Human Butyrylcholinesterase Mutants for (-)-Cocaine Hydrolysis: A Correlation Relationship between Catalytic Efficiency and Total Hydrogen Bonding Energy with an Oxyanion Hole

A combined computational and experimental study has been carried out to explore and test a quantitative correlation relationship between the relative catalytic efficiency (RCE) of human butyrylcholinesrase (BChE) mutant-catalyzed hydrolysis of substrate (-)-cocaine and the total hydrogen bonding energy (tHBE) of the carbonyl oxygen of the substrate with the oxyanion hole of the enzyme in the modeled transition-state structure (TS1), demonstrating a satisfactory linear correlation relationship between ln(RCE) and tHBE. The satisfactory correlation relationship has led us to computationally predict and experimentally confirm new human BChE mutants that have a further improved catalytic activity against (-)-cocaine, including the most active one (the A199S/F227S/S287G/A328W/Y332G mutant) with a 2790-fold improved catalytic efficiency (kcat/KM = 2.5 × 109 min-1 M-1) compared to the wild-type human BChE. Compared to the reference mutant (the A199S/S287G/A328W/Y332G mutant) tested in the reported clinical development of an enzyme therapy for cocaine dependence treatment, this new mutant (with a newly predicted additional F227S mutation) has an improved catalytic efficiency against (-)-cocaine by ∼2.6-fold. The good agreement between the computational and experimental ln(RCE) values suggests that the obtained correlation relationship is robust for computational prediction. A similar correlation relationship could also be explored in studying BChE or other serine hydrolases/esterases with an oxyanion hole stabilizing the carbonyl oxygen in the rate-determining reaction step of the enzymatic hydrolysis of other substrates.

Computational Design and Crystal Structure of a Highly Efficient Benzoylecgonine Hydrolase

Benzoylecgonine (BZE) is the major toxic metabolite of cocaine and is responsible for the long-term cocaine-induced toxicity owing to its long residence time in humans. BZE is also the main contaminant following cocaine consumption. Here, we identified the bacterial cocaine esterase (CocE) as a BZE-metabolizing enzyme (BZEase), which can degrade BZE into biological inactive metabolites (ecgonine and benzoic acid). CocE was redesigned by a reactant-state-based enzyme design theory. An encouraging mutant denoted as BZEase2, presented a >400-fold improved catalytic efficiency against BZE compared with wild-type (WT) CocE. In vivo, a single dose of BZEase2 (1 mg kg^-1 , IV) could eliminate nearly all BZE within only two minutes, suggesting the enzyme has the potential for cocaine overdose treatment and BZE elimination in the environment by accelerating BZE clearance. The crystal structure of a designed BZEase was also determined.

Catalytic Mechanisms for Cofactor-Free Oxidase-Catalyzed Reactions: Reaction Pathways of Uricase-Catalyzed Oxidation and Hydration of Uric Acid

First-principles quantum mechanical/molecular mechanical (QM/MM)-free energy calculations have been performed to uncover how uricase catalyzes metabolic reactions of uric acid (UA), demonstrating that the entire reaction process of UA in uricase consists of two stages-oxidation followed by hydration. The oxidation consists of four steps: (1) chemical transformation from 8-hydroxyxythine to an anionic radical via a proton transfer along with an electron transfer, which is different from the previously proposed electron-transfer mechanism that involves a dianion intermediate (UA^2-) during the catalytic reaction process; (2) proton transfer to the O₂^- anion (radical); (3) diradical recombination to form a peroxo intermediate; (4) dissociation of H₂O₂ to generate the dehydrourate. Hydration, for the most favorable pathway, is initiated by the nucleophilic attack of a water molecule on dehydrourate, along with a concerted proton transfer through residue Thr69 in the catalytic site. According to the calculated free energy profile, the hydration is the rate-determining step, and the corresponding free energy barrier of 16.2 kcal/mol is consistent with that derived from experimental kinetic data, suggesting that the computational insights into the catalytic mechanisms are reasonable. The mechanistic insights not only provide a mechanistic base for future rational design of uricase mutants with improved catalytic activity against uric acid as an improved enzyme therapy, but also are valuable for understanding a variety of other cofactor-free oxidase-catalyzed reactions involving an oxygen molecule.

Free energy profiles of cocaine esterase-cocaine binding process by molecular dynamics and potential of mean force simulations

The combined molecular dynamics (MD) and potential of mean force (PMF) simulations have been performed to determine the free energy profile of the CocE)-(+)-cocaine binding process in comparison with that of the corresponding CocE-(-)-cocaine binding process. According to the MD simulations, the equilibrium CocE-(+)-cocaine binding mode is similar to the CocE-(-)-cocaine binding mode. However, based on the simulated free energy profiles, a significant free energy barrier (∼5 kcal/mol) exists in the CocE-(+)-cocaine binding process whereas no obvious free energy barrier exists in the CocE-(-)-cocaine binding process, although the free energy barrier of ∼5 kcal/mol is not high enough to really slow down the CocE-(+)-cocaine binding process. In addition, the obtained free energy profiles also demonstrate that (+)-cocaine and (-)-cocaine have very close binding free energies with CocE, with a negligible difference (∼0.2 kcal/mol), which is qualitatively consistent with the nearly same experimental K_M values of the CocE enzyme for (+)-cocaine and (-)-cocaine. The consistency between the computational results and available experimental data suggests that the mechanistic insights obtained from this study are reasonable.

Fundamental reaction pathway and free energy profile of proteasome inhibition by syringolin A (SylA)

In this study, molecular dynamics (MD) simulations and first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations have been performed to uncover the fundamental reaction pathway of proteasome with a representative inhibitor syringolin A (SylA). The calculated results reveal that the reaction process consists of three steps. The first step is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is a nucleophilic attack on the olefin carbon of SylA by the negatively charged Thr1-O(γ) atom. The last step is a proton transfer from Thr1-N(z) to another olefin carbon of SylA to complete the inhibition reaction process. The calculated free energy profile demonstrates that the second step should be the rate-determining step and has the highest free energy barrier of 24.6 kcal mol(-1), which is reasonably close to the activation free energy (∼22.4-23.0 kcal mol(-1)) derived from the available experimental kinetic data. In addition, our computational results indicate that no water molecule can assist the rate-determining step, since the second step is not involved in a proton transfer process. The obtained mechanistic insights should be valuable for understanding the inhibition process of proteasome by SylA and structurally related inhibitors at a molecular level, and thus provide a solid mechanistic base and valuable clues for future rational design of novel, more potent inhibitors of proteasome.

Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine

A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.

Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson-Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.

Fundamental reaction pathway for peptide metabolism by proteasome: insights from first-principles quantum mechanical/molecular mechanical free energy calculations

Proteasome is the major component of the crucial non-lysosomal protein degradation pathway in the cells, but the detailed reaction pathway is unclear. In this study, first-principles quantum mechanical/molecular mechanical free energy calculations have been performed to explore, for the first time, possible reaction pathways for proteasomal proteolysis/hydrolysis of a representative peptide, succinyl-leucyl-leucyl-valyl-tyrosyl-7-amino-4-methylcoumarin (Suc-LLVY-AMC). The computational results reveal that the most favorable reaction pathway consists of six steps. The first is a water-assisted proton transfer within proteasome, activating Thr1-O(γ). The second is a nucleophilic attack on the carbonyl carbon of a Tyr residue of substrate by the negatively charged Thr1-O(γ), followed by the dissociation of the amine AMC (third step). The fourth step is a nucleophilic attack on the carbonyl carbon of the Tyr residue of substrate by a water molecule, accompanied by a proton transfer from the water molecule to Thr1-N(z). Then, Suc-LLVY is dissociated (fifth step), and Thr1 is regenerated via a direct proton transfer from Thr1-N(z) to Thr1-O(γ). According to the calculated energetic results, the overall reaction energy barrier of the proteasomal hydrolysis is associated with the transition state (TS3(b)) for the third step involving a water-assisted proton transfer. The determined most favorable reaction pathway and the rate-determining step have provided a reasonable interpretation of the reported experimental observations concerning the substituent and isotopic effects on the kinetics. The calculated overall free energy barrier of 18.2 kcal/mol is close to the experimentally derived activation free energy of ∼18.3-19.4 kcal/mol, suggesting that the computational results are reasonable.

Fundamental reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of heroin

The pharmacological function of heroin requires an activation process that transforms heroin into 6-monoacetylmorphine (6-MAM), which is the most active form. The primary enzyme responsible for this activation process in human plasma is butyrylcholinesterase (BChE). The detailed reaction pathway of the activation process via BChE-catalyzed hydrolysis has been explored computationally, for the first time, in this study via molecular dynamics simulation and first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the whole reaction process includes acylation and deacylation stages. The acylation consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the 3-acetyl group of heroin by the hydroxyl oxygen of the Ser198 side chain and the dissociation of 6-MAM. The deacylation also consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the acyl-enzyme intermediate by a water molecule and the dissociation of the acetic acid from Ser198. The calculated free energy profile reveals that the second transition state (TS2) should be rate-determining. The structural analysis reveals that the oxyanion hole of BChE plays an important role in the stabilization of rate-determining TS2. The free energy barrier (15.9 ± 0.2 or 16.1 ± 0.2 kcal/mol) calculated for the rate-determining step is in good agreement with the experimentally derived activation free energy (~16.2 kcal/mol), suggesting that the mechanistic insights obtained from this computational study are reliable. The obtained structural and mechanistic insights could be valuable for use in the future rational design of a novel therapeutic treatment of heroin abuse.

Catalytic mechanism of cytochrome P450 for N-methylhydroxylation of nicotine: reaction pathways and regioselectivity of the enzymatic nicotine oxidation

The fundamental reaction mechanism of cytochrome P450 2A6 (CYP2A6)-catalyzed N-methylhydroxylation of (S)-(-)-nicotine and the free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical (QM/MM) reaction-coordinate calculations. In the CYP2A6-(S)-(-)-nicotine binding structures that allow for 5'-hydroxylation, the N-methyl group is also sufficiently close to the oxygen of Cpd I for the N-methylhydroxylation reaction to occur. It has been demonstrated that the CYP2A6-catalyzed N-methylhydroxylation reaction is a concerted process involving a hydrogen-transfer transition state on both the quartet and the doublet states. The N-methylhydroxylation reaction proceeds mainly in the doublet state, since the free energy barriers on the doublet state are lower than the corresponding ones on the quartet state. The calculated free energy barriers indicate that (S)-(-)-nicotine oxidation catalyzed by CYP2A6 proceeds with a high regioselective abstraction of the hydrogen at the 5'-position, rather than the hydrogen at the N-methyl group. The predicted regioselectivity of 93% is in agreement with the most recent experimentally reported regioselectivity of 95%. The binding mode of (S)-(-)-nicotine in the active site of CYP2A6 is an important determinant for the stereoselectivity of nicotine (S)-(-)-oxidation, whereas the regioselectivity of (S)-(-)-nicotine oxidation is determined mainly by the free energy barrier difference between the 5'-hydroxylation and N-methylhydroxylation reactions.

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.

Cocaine esterase-cocaine binding process and the free energy profiles by molecular dynamics and potential of mean force simulations

The combined molecular dynamics (MD) and potential of mean force (PMF) simulations have been performed to determine the free energy profiles for the binding process of (-)-cocaine interacting with wild-type cocaine esterase (CocE) and its mutants (T172R/G173Q and L119A/L169K/G173Q). According to the MD simulations, the general protein-(-)-cocaine binding mode is not affected by the mutations; e.g.. the benzoyl group of (-)-cocaine is always bound in a subsite composed of aromatic residues W151, W166, F261, and F408 and hydrophobic residue L407, while the carbonyl oxygen on the benzoyl group of (-)-cocaine is hydrogen-bonded with the oxyanion-hole residues Y44 and Y118. According to the PMF-calculated free energy profiles for the binding process, the binding free energies for (-)-cocaine with the wild-type, T172R/G173Q, and L119A/L169K/G173Q CocEs are predicted to be -6.4, -6.2, and -5.0 kcal/mol, respectively. The computational predictions are supported by experimental kinetic data, as the calculated binding free energies are in good agreement with the experimentally derived binding free energies, i.e., -7.2, -6.7, and -4.8 kcal/mol for the wild-type, T172R/G173Q, and L119A/L169K/G173Q, respectively. The reasonable agreement between the computational and experimental data suggests that the PMF simulations may be used as a valuable tool in new CocE mutant design that aims to decrease the Michaelis-Menten constant of the enzyme for (-)-cocaine.