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

Diels-Alder Adducts of Morphinan-6,8-Dienes and Their Transformations

6,14-ethenomorphinans are semisynthetic opiate derivatives containing an ethylene bridge between positions 6 and 14 in ring-C of the morphine skeleton that imparts a rigid molecular structure. These compounds represent an important family of opioid receptor ligands in which the 6,14-etheno bridged structural motif originates from a [4 + 2] cycloaddition of morphinan-6,8-dienes with dienophiles. Certain 6,14-ethenomorphinans having extremely high affinity for opioid receptors are often non-selective for opioid receptor subtypes, but this view is now undergoing some revision. The agonist 20R-etorphine and 20R-dihydroetorphine are several thousand times more potent analgesics than morphine, whereas diprenorphine is a high-affinity non-selective antagonist. The partial agonist buprenorphine is used as an analgesic in the management of post-operative pain or in substitution therapy for opiate addiction, sometimes in combination with the non-selective antagonist naloxone. In the context of the current opioid crisis, we communicated a summary of several decades of work toward generating opioid analgesics with lesser side effects or abuse potential. Our summary placed a focus on Diels-Alder reactions of morphinan-6,8-dienes and subsequent transformations of the cycloadducts. We also summarized the pharmacological aspects of radiolabeled 6,14-ethenomorphinans used in molecular imaging of opioid receptors.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s^-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

Systematic Structure-Based Virtual Screening Approach to Antibody Selection and Design of a Humanized Antibody against Multiple Addictive Opioids without Affecting Treatment Agents Naloxone and Naltrexone

Opioid drug use, especially heroin, is known as a growing national crisis in America. Heroin itself is a prodrug and is converted to the most active metabolite 6-monoacetylmorphine (6-MAM) responsible for the acute toxicity of heroin and then to a relatively less-active metabolite morphine responsible for the long-term toxicity of heroin. Monoclonal antibodies (mAbs) are recognized as a potentially promising therapeutic approach in the treatment of opioid use disorders (OUDs). Due to the intrinsic challenges of discovering an mAb against multiple ligands, here we describe a general, systematic structure-based virtual screening and design approach which has been used to identify a known anti-morphine antibody 9B1 and a humanized antibody h9B1 capable of binding to multiple addictive opioids (including 6-MAM, morphine, heroin, and hydrocodone) without significant binding with currently available OUD treatment agents naloxone, naltrexone, and buprenorphine. The humanized antibody may serve as a promising candidate for the treatment of OUDs. The experimental binding affinities reasonably correlate with the computationally predicted binding free energies. The experimental activity data strongly support the computational predictions, suggesting that the systematic structure-based virtual screening and humanization design protocol is reliable. The general, systematic structure-based virtual screening and design approach will be useful for many other antibody selection and design efforts in the future.

Diagnoses of Pathological States Based on Acetylcholinesterase and Butyrylcholinesterase

Two cholinesterases exist: Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). While AChE plays a crucial role in neurotransmissions, BChE has no specific function apart from the detoxification of some drugs and secondary metabolites from plants. Thus, both AChE and BChE can serve as biochemical markers of various pathologies. Poisoning by nerve agents like sarin, soman, tabun, VX, novichok and overdosing by drugs used in some neurodegenerative disorders like Alzheimer´s disease and myasthenia gravis, as well as poisoning by organophosphorus pesticides are relevant to this issue. But it appears that changes in these enzymes take place in other processes including oxidative stress, inflammation, some types of cancer and genetically conditioned diseases. In this review, the cholinesterases are introduced, the mechanism of inhibitors action is explained and the relations between the cholinesterases and pathologies are explained.

Cebranopadol reduces cocaine self-administration in male rats: Dose, treatment and safety consideration

As a novel first-in-class potent analgesic acting as an agonist of multiple opioid receptors, cebranopadol showed high efficacy and good tolerability in a broad range of preclinical models and clinical trials related to pain. In the present study, to evaluate the efficacy and safety of cebranopadol as a potential treatment of cocaine dependence, we tested the effects of cebranopadol with single and repeated doses (25, 50, 75, or 100 μg/kg, oral gavage) using rat models of cocaine fixed-ratio (FR) self-administration (SA), cocaine progressive-ratio (PR) SA, and sucrose pellet SA. In single-dosing treatment paradigm, cebranopadol significantly and dose-dependently reduced cocaine SA under FR and PR schedules and suppressed food intake under FR schedule without causing apparent side effects. In repeated-dosing treatment scheme, i.e. daily administration of 25, 50, 75, or 100 μg/kg cebranopadol for a week, the similar reduction in cocaine intake was detected, while non-negligible complications/side effects were observed at repeated high doses (75 and 100 μg/kg). The observed side effects were similar to the common toxic signs elicited by heroin at high doses, although cebranopadol did not fully substitute heroin's discriminative stimulant effects in our drug discriminative tests. These results demonstrated that the most appropriate oral dose of cebranopadol to balance the efficacy and safety is 50 μg/kg. Collectively, although cebranopadol may serve as a new treatment for cocaine dependence, more consideration, cautiousness, and a clear optimal dose window to dissociate its therapeutic effects from opioid side effects/complications in male and female subjects will be necessary to increase its practical clinical utility.

DARK Classics in Chemical Neuroscience: Heroin and Desomorphine

Opioids are arguably one of the most important pharmacologic classes, mainly due to their rich history, their useful and potent analgesic effects, and also, just as importantly, their "Dark Side", constituted by their reinforcing properties that have led countless of users to a spiral of addiction, biological dependence, tolerance, withdrawal syndromes, and death. Among the most significant abused and addictive known opioids are heroin and desomorphine, both synthetic derivatives of morphine that belong to the 4,5-epoxymorphinan structural chemical group of the opioid family drugs. These agents share not only structural, pharmacological, and epidemiological features but also a common geographical distribution. A drop in Afghan heroin production and its "exports" to Russia gave rise to widespread consumption of desomorphine in ex-Soviet republics during the first decade of the 21st century, representing an economical and accessible alternative for misusers to this sort of derivative. Herein we review the state of the art of history, chemistry and synthesis, pharmacology, and impact on society of these "cursed cousins".

Discovery of potent and selective butyrylcholinesterase inhibitors through the use of pharmacophore-based screening

Cholinesterase inhibitors have long been used in the treatment of Alzheimer's Disease (AD) via the protection of acetylcholine levels. However, recent research has shown that the specific inhibition of butyrylcholinesterase (BChE) could better ameliorate symptoms within patients. In addition, it has recently been shown that selective inhibition of BChE can also significantly attenuate the toxicity and physiological effects of heroin. Currently, there are no specific and potent inhibitors of BChE approved for use in AD or heroin abuse. Through a combined use of in silico and in vitro screening, we have found three compounds with sub-50 nM IC50 values that specifically target BChE. These newly discovered BChE inhibitors can act as the lead scaffolds for future development of the desirably potent and selective BChE inhibitors.

Cryo-EM structure of the native butyrylcholinesterase tetramer reveals a dimer of dimers stabilized by a superhelical assembly

The quaternary structures of the cholinesterases, acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), are essential for their localization and function. Of practical importance, BChE is a promising therapeutic candidate for intoxication by organophosphate nerve agents and insecticides, and for detoxification of addictive substances. Efficacy of the recombinant enzyme hinges on its having a long circulatory half-life; this, in turn, depends strongly on its ability to tetramerize. Here, we used cryoelectron microscopy (cryo-EM) to determine the structure of the highly glycosylated native BChE tetramer purified from human plasma at 5.7 Å. Our structure reveals that the BChE tetramer is organized as a staggered dimer of dimers. Tetramerization is mediated by assembly of the C-terminal tryptophan amphiphilic tetramerization (WAT) helices from each subunit as a superhelical assembly around a central lamellipodin-derived oligopeptide with a proline-rich attachment domain (PRAD) sequence that adopts a polyproline II helical conformation and runs antiparallel. The catalytic domains within a dimer are asymmetrically linked to the WAT/PRAD. In the resulting arrangement, the tetramerization domain is largely shielded by the catalytic domains, which may contribute to the stability of the human BChE (HuBChE) tetramer. Our cryo-EM structure reveals the basis for assembly of the native tetramers and has implications for the therapeutic applications of HuBChE. This mode of tetramerization is seen only in the cholinesterases but may provide a promising template for designing other proteins with improved circulatory residence times.

Catalytic Reaction Mechanism for Drug Metabolism in Human Carboxylesterase-1: Cocaine Hydrolysis Pathway

Carboxylesterase-1 (CE-1) is a crucial enzyme responsible for metabolism/activation/inactivation of xenobiotics (therapeutic agents, prodrugs, abused drugs, and organophosphorus nerve agents etc.) and also involved in many other biological processes. In this study, we performed extensive computational modeling and simulations to understand the fundamental reaction mechanism of cocaine hydrolysis catalyzed by CE-1, revealing that CE-1-catalyzed cocaine hydrolysis follows a novel reaction pathway with only two reaction steps: a single-step acylation process and a single-step deacylation process. In the transition states of both single-step processes, the cocaine NH group joins the oxyanion hole to form an additional hydrogen bond with the negatively charged carbonyl oxygen atom of the cocaine. Thus, the transition states are stabilized by both intermolecular and intramolecular hydrogen bonds with the methyl ester of cocaine, specifically the carbonyl oxygen atom. The rate-limiting transition state is associated with the acylation process, and the activation free energy barrier was predicted to be 20.1 kcal/mol. Further, in vitro experimental kinetic analysis was performed for human CE-1-catalyzed cocaine hydrolysis. For CE-1-catalyzed cocaine hydrolysis, the computationally predicted free energy barrier (20.1 kcal/mol) is reasonably close to the experimentally derived turnover number ( k_cat = 0.058 min^-1), indicating the reasonability of the computational results. The obtained novel mechanistic insights are expected to benefit not only CE-1 related rational drug discovery but also future research on the catalytic mechanism of other esterases.

Kinetic characterization of cholinesterases and a therapeutically valuable cocaine hydrolase for their catalytic activities against heroin and its metabolite 6-monoacetylmorphine

As the most popularly abused one of opioids, heroin is actually a prodrug. In the body, heroin is hydrolyzed/activated to 6-monoacetylmorphine (6-MAM) first and then to morphine to produce its toxic and physiological effects. It has been known that heroin hydrolysis to 6-MAM and morphine is accelerated by cholinesterases, including acetylcholinesterase (AChE) and/or butyrylcholinesterase (BChE). However, there has been controversy over the specific catalytic activities and functional significance of the cholinesterases, which requires for the more careful kinetic characterization under the same experimental conditions. Here we report the kinetic characterization of AChE, BChE, and a therapeutically promising cocaine hydrolase (CocH1) for heroin and 6-MAM hydrolyses under the same experimental conditions. It has been demonstrated that AChE and BChE have similar k_cat values (2100 and 1840 min^-1, respectively) against heroin, but with a large difference in K_M (2170 and 120 μM, respectively). Both AChE and BChE can catalyze 6-MAM hydrolysis to morphine, with relatively lower catalytic efficiency compared to the heroin hydrolysis. CocH1 can also catalyze hydrolysis of heroin (k_cat = 2150 min^-1 and K_M = 245 μM) and 6-MAM (k_cat = 0.223 min^-1 and K_M = 292 μM), with relatively larger K_M values and lower catalytic efficiency compared to BChE. Notably, the K_M values of CocH1 against both heroin and 6-MAM are all much larger than previously reported maximum serum heroin and 6-MAM concentrations observed in heroin users, implying that the heroin use along with cocaine will not drastically affect the catalytic activity of CocH1 against cocaine in the CocH1-based enzyme therapy for cocaine abuse.

Insights into the isothiourea-catalyzed asymmetric [4 + 2] annulation of phenylacetic acid with alkylidene pyrazolone

A computational study on the isothiourea-catalyzed asymmetric [4 + 2] annulation of phenylacetic acid with alkylidene pyrazolone was performed using the DFT method.

Binding Space Concept: A New Approach To Enhance the Reliability of Docking Scores and Its Application to Predicting Butyrylcholinesterase Hydrolytic Activity

Docking simulations are very popular approaches able to assess the capacity of a given ligand to interact with a target. Docking simulations are usually focused on a single best complex even though many studies showed that ligands retain a significant mobility within a binding pocket by assuming different binding modes all of which may contribute to the monitored ligand affinity. The present study describes an innovative concept, the binding space, which allows an exploration of the ligand mobility within the binding pocket by simultaneously considering several ligand poses as generated by docking simulations. The multiple poses and the relative docking scores can then be analyzed by taking advantage of the same concepts already used in the property space analysis; hence the binding space can be parametrized by (a) mean scores, (b) score ranges, and (c) score sensitivity values. The first parameter represents a very simple procedure to account for the contribution of the often neglected alternative binding modes, while the last two descriptors encode the degree of mobility which a given ligand retains within the binding cavity (score range) as well as the ease with which a ligand explores such a mobility (score sensitivity). Here, the binding space concept is applied to the prediction of the hydrolytic activity of BChE by synergistically considering multiple poses and multiple protein structures. The obtained results shed light on the remarkable potential of the binding space concept, whose parameters allow a significant increase of the predictive power of the docking results as revealed by extended correlative analyses. Mean scores are the parameters affording the largest statistical improvement, and all the here proposed docking-based descriptors show enhancing effects in developing predictive models. Finally, the study describes a new score function (Contacts score) simply based on the number of surrounding residues which appears to be particularly productive in the framework of the binding space.

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.

Site- and species-specific hydrolysis rates of heroin

The hydroxide-catalyzed non-enzymatic, simultaneous and consecutive hydrolyses of diacetylmorphine (DAM, heroin) are quantified in terms of 10 site- and species-specific rate constants in connection with also 10 site- and species-specific acid-base equilibrium constants, comprising all the 12 coexisting species in solution. This characterization involves the major and minor decomposition pathways via 6-acetylmorphine and 3-acetylmorphine, respectively, and morphine, the final product. Hydrolysis has been found to be 18-120 times faster at site 3 than at site 6, depending on the status of the amino group and the rest of the molecule. Nitrogen protonation accelerates the hydrolysis 5-6 times at site 3 and slightly less at site 6. Hydrolysis rate constants are interpreted in terms of intramolecular inductive effects and the concomitant local electron densities. Hydrolysis fraction, a new physico-chemical parameter is introduced and determined to quantify the contribution of the individual microspecies to the overall hydrolysis. Hydrolysis fractions are depicted as a function of pH.

Is it possible to reverse aged acetylcholinesterase inhibited by organophosphorus compounds? Insight from the theoretical study

The main treatment for organophosphorus (OP) compound poisoning in clinics is to restore the activity of acetylcholinesterase (AChE) through oxime-induced reactivation of the phosphorylated OP-AChE adduct. It suffers from a competitive and irreversible aging reaction of the phosphorylated OP-AChE adduct, resulting in permanent inactivity of AChE. However, it was recently reported that N-methyl-2-methoxypyridinium species can act as methylating agents to methylate the methyl methane-phosphonate monoanion, in which the reaction mimics the reverse of the aging reaction of the phosphorylated OP-AChE adduct. If the aging reaction could be really reversed, the efficiency for the OP detoxification should be significantly improved, bringing up the possibility to develop an agent to reverse the aging process of the phosphorylated OP-AChE adduct. However, such a reaction with the N-methyl-2-methoxypyridinium species in the enzyme is still not reported so far. It is of great interest to know whether or not this reaction is observable in the enzyme, and more importantly, if it turns out to be not observable in the enzyme, why such a reaction proceeds quickly in aqueous solution but not in the enzyme. In the present study, we performed DFT calculations and quantum mechanical/molecular mechanical (QM/MM) calculations to reveal the fundamental mechanism for the methylation of both the methyl methane-phosphonate monoanion and the aged sarin-AChE adduct by N-methyl-2-methoxypyridinium species, respectively. The obtained results support the SN2 reaction mechanism, not the stepwise mechanism, for the methylation of the methyl methane-phosphonate monoanion by 9 reported N-methyl-2-methoxypyridinium compounds. The calculated free energy barriers are in good agreement with the experimental data. The methylation of the aged sarin-AChE adduct by one N-methyl-2-methoxypyridinium compound (labeled as compound 2) also employs the SN2 reaction mechanism with an extremely high free energy barrier of 30.4 ± 3.5 (or 26.6) kcal mol(-1), implying that this reaction in the enzyme hardly occurs. Our results clearly show that compound 2 forms a strong π-π stacking interaction with the aromatic ring of the W86 residue of AChE, making itself unable to approach sarin for the reverse of the aging process. On the basis of the structure and mechanism, several possible strategies have been suggested for designing methylating agents with higher activity against the aged sarin-AChE adduct.

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.

Unexpected Reaction Pathway for butyrylcholinesterase-catalyzed inactivation of "hunger hormone" ghrelin

Extensive computational modeling and simulations have been carried out, in the present study, to uncover the fundamental reaction pathway for butyrylcholinesterase (BChE)-catalyzed hydrolysis of ghrelin, demonstrating that the acylation process of BChE-catalyzed hydrolysis of ghrelin follows an unprecedented single-step reaction pathway and the single-step acylation process is rate-determining. The free energy barrier (18.8 kcal/mol) calculated for the rate-determining step is reasonably close to the experimentally-derived free energy barrier (~19.4 kcal/mol), suggesting that the obtained mechanistic insights are reasonable. The single-step reaction pathway for the acylation is remarkably different from the well-known two-step acylation reaction pathway for numerous ester hydrolysis reactions catalyzed by a serine esterase. This is the first time demonstrating that a single-step reaction pathway is possible for an ester hydrolysis reaction catalyzed by a serine esterase and, therefore, one no longer can simply assume that the acylation process must follow the well-known two-step reaction pathway.

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.

Mechanistic insights into cobalt(ii/iii)-catalyzed C-H oxidation: a combined theoretical and experimental study

Cobalt-mediated C-H functionalization has been the subject of extensive interest in synthetic chemistry, but the mechanisms of many of these reactions (such as the cobalt-catalyzed C-H oxidation) are poorly understood. In this paper, possible mechanisms including single electron transfer (SET) and the concerted metalation-deprotonation (CMD) pathways of the CoII/CoIII-catalyzed alkoxylation of C(sp2)-H bonds have been investigated for the first time using the DFT method. CoII(OAc)2 has been employed as an efficient catalyst in our previous experimental study, but the calculated results unexpectedly indicated that the intermolecular SET pathway with CoIII as the actual catalyst might be the most favorable pathway. To support this theoretical prediction, we have explored a series of Cp*CoIII(CO)I2 catalyzed C(sp2)-H bond alkoxylations, extending the application of cobalt-catalyzed functionalization of C-H bonds. Furthermore, kinetic isotope effect (KIE) data, electron paramagnetic resonance (EPR) data, and TEMPO inhibition experiments also support the SET mechanism in both the Co-catalyzed alkoxylation reactions. Thus, this work should support an understanding of the possible mechanisms of the CoII/CoIII-catalyzed C(sp2)-H functionalization, and also provide an example of the rational design of novel catalytic reactions guided by theoretical calculations.

Reaction pathways and free energy profiles for cholinesterase-catalyzed hydrolysis of 6-monoacetylmorphine

As the most active metabolite of heroin, 6-monoacetylmorphine (6-MAM) can penetrate into the brain for the rapid onset of heroin effects. The primary enzymes responsible for the metabolism of 6-MAM to the less potent morphine in humans are acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The detailed reaction pathways for AChE- and BChE-catalyzed hydrolysis of 6-MAM to morphine have been explored, for the first time, in the present study by performing first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the two enzymatic reaction processes follow similar catalytic reaction mechanisms, and the whole catalytic reaction pathway for each enzyme consists of four reaction steps. According to the calculated results, the second reaction step associated with the transition state TS2(a)/TS2(b) should be rate-determining for the AChE/BChE-catalyzed hydrolysis, and the free energy barrier calculated for the AChE-catalyzed hydrolysis (18.3 kcal mol(-1)) is 2.5 kcal mol(-1) lower than that for the BChE-catalyzed hydrolysis (20.8 kcal mol(-1)). The free energy barriers calculated for the AChE- and BChE-catalyzed reactions are in good agreement with the experimentally derived activation free energies (17.5 and 20.7 kcal mol(-1) for the AChE- and BChE-catalyzed reactions, respectively). Further structural analysis reveals that the aromatic residues Phe295 and Phe297 in the acyl pocket of AChE (corresponding to Leu286 and Val288 in BChE) contribute to the lower energy of TS2(a) relative to TS2(b). The obtained structural and mechanistic insights could be valuable for use in future rational design of a novel therapeutic treatment of heroin abuse.