Design of high-activity mutants of human butyrylcholinesterase against (-)-cocaine: structural and energetic factors affecting the catalytic efficiency

Long-lasting blocking of interoceptive effects of cocaine by a highly efficient cocaine hydrolase in rats

Cocaine dependence is a serious world-wide public health problem without an FDA-approved pharmacotherapy. We recently designed and discovered a highly efficient long-acting cocaine hydrolase CocH5-Fc(M6). The present study examined the effectiveness and duration of CocH5-Fc(M6) in blocking interoceptive effects of cocaine by performing cocaine discrimination tests in rats, demonstrating that the duration of CocH5-Fc(M6) in blocking cocaine discrimination was dependent on cocaine dose and CocH5-Fc(M6) plasma concentration. Particularly, a dose of 3 mg/kg CocH5-Fc(M6) effectively attenuated discriminative stimulus effects of 10 mg/kg cocaine, cumulative doses of 10 and 32 mg/kg cocaine, and cumulative doses of 10, 32 and 56 mg/kg cocaine by ≥ 20% for 41, 19, and 10 days, and completely blocked the discriminative stimulus effects for 30, 13, and 5 days with corresponding threshold plasma CocH5-Fc(M6) concentrations of 15.9, 72.2, and 221 nM, respectively, under which blood cocaine concentration was negligible. Additionally, based on the data obtained, cocaine discrimination model is more sensitive than the locomotor activity to reveal cocaine effects and that CocH5-Fc(M6) itself has no long-term toxicity regarding behavioral activities such as lever pressing and food consumption in rats, further demonstrating that CocH5-Fc(M6) has the desired properties as a promising therapeutic candidate for prevenance of cocaine dependence.

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.

Catalytic activities of a highly efficient cocaine hydrolase for hydrolysis of biologically active cocaine metabolites norcocaine and benzoylecgonine

Cocaine is a widely abused, hepatotoxic drug without an FDA-approved pharmacotherapy specific for cocaine addiction or overdose. It is recognized as a promising therapeutic strategy to accelerate cocaine metabolism which can convert cocaine to pharmacologically inactive metabolite(s) using an efficient cocaine-metabolizing enzyme. Our previous studies have successfully designed and discovered a highly efficient cocaine hydrolase, denoted as CocH5-Fc(M6), capable of rapidly hydrolyzing cocaine at the benzoyl ester moiety. In the present study, we determined the kinetic parameters of CocH5-Fc(M6) against norcocaine (k_cat = 9,210 min^-1, K_M = 20.9 μM, and k_cat/K_M = 1.87 × 10⁵ min^-1 M^-1) and benzoylecgonine (k_cat = 158 min^-1, K_M = 286 μM, and k_cat/K_M = 5.5 × 10⁵ min^-1 M^-1) for the first time. Further in vivo studies have demonstrated that CocH5-Fc(M6) can effectively accelerate clearance of not only cocaine, but also norcocaine (known as a cocaine metabolite which is more toxic than cocaine itself) and benzoylecgonine (known as an unfavorable long-lasting metabolite with some long-term toxicity concerns) in rats. Due to the desired high catalytic activity against norcocaine, CocH5-Fc(M6) is capable of quickly detoxifying both cocaine and its more toxic metabolite norcocaine after intraperitoneally administering lethal dose of 60 or 180 mg/kg cocaine. In addition, the ability of CocH5-Fc(M6) to accelerate clearance of benzoylecgonine should also be valuable for the use of CocH5-Fc(M6) in treatment of cocaine use disorder.

Gender differences in cocaine-induced hyperactivity and dopamine transporter trafficking to the plasma membrane

As well known, cocaine induces stimulant effects and dopamine transporter (DAT) trafficking to the plasma membrane of dopaminergic neurons. In the present study, we examined cocaine-induced hyperactivity along with cocaine-induced DAT trafficking and the recovery rate of the dopaminergic system in female rats in comparison with male rats, demonstrating interesting gender differences. Female rats are initially more sensitive to cocaine than male rats in terms of both the DAT trafficking and hyperactivity induced by cocaine. Particularly, intraperitoneal (i.p.) administration of 5 mg/kg cocaine induced significant hyperactivity and DAT trafficking in female rats but not in male rats. After repeated cocaine exposures (i.e., i.p. administration of 20 mg/kg cocaine every other day from Day 0 to Day 32), cocaine-induced hyperactivity in female rats gradually became a clear pattern of two phases, with the first phase of the hyperactivity lasting for only a few minutes and the second phase lasting for over an hour beginning at ~30 min, which is clearly different from that of male rats. It has also been demonstrated that the striatal DAT distribution of female rats may recover faster than that of male rats after multiple cocaine exposures. Nevertheless, despite the remarkable gender differences, our recently developed long-acting cocaine hydrolase, known as CocH5-Fc(M6), can similarly and effectively block cocaine-induced DAT trafficking and hyperactivity in both male and female rats.

Development of a Highly Efficient Long-Acting Cocaine Hydrolase Entity to Accelerate Cocaine Metabolism

It is particularly challenging to develop a truly effective pharmacotherapy for cocaine use disorder (CUD) treatment. Accelerating cocaine metabolism via hydrolysis at cocaine benzoyl ester using an efficient cocaine hydrolase (CocH) is known as a promising pharmacotherapeutic approach to CUD treatment. Preclinical and clinical studies on our first CocH (CocH1), in its human serum albumin-fused form known as TV-1380, have demonstrated the promise of a general concept of CocH-based pharmacotherapy for CUD treatment. However, the biological half-life of TV-1380 (t_1/2 = 8 h in rats, associated with t_1/2 = 43-77 h in humans) is not long enough for practical treatment of cocaine dependence, which requires enzyme injection for no more than once weekly. Through protein fusion of a human butyrylcholinesterase mutant (denoted as CocH5) with a mutant (denoted as Fc(M6)) of Fc from human IgG1, we have designed, prepared, and tested a new fusion protein (denoted as CocH5-Fc(M6)) for its pharmacokinetic profile and in vivo catalytic activity against (-)-cocaine. CocH5-Fc(M6) represents the currently most efficient long-acting cocaine hydrolase with both the highest catalytic activity against (-)-cocaine and the longest elimination half-life (t_1/2 = 229 ± 5 h) in rats. As a result, even at a single modest dose of 3 mg/kg, CocH5-Fc(M6) can significantly and effectively accelerate the metabolism of cocaine in rats for at least 60 days. In addition, ∼70 nM CocH5-Fc(M6) in plasma was able to completely block the toxicity and physiological effects induced by intraperitoneal injection of a lethal dose of cocaine (60 mg/kg).

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.

Clinical potential of a rationally engineered enzyme for treatment of cocaine dependence: Long-lasting blocking of the psychostimulant, discriminative stimulus, and reinforcing effects of cocaine

It is a grand challenge to develop a truly effective treatment of substance use disorder (SUD), particularly for cocaine and other drugs without an FDA-approved treatment available, because a truly effective therapy must effectively block the drug's physiological and reinforcing effects during the entire period of treatment in order to achieve the long-time abstinence required by the FDA. Whether a biologic, such as monoclonal antibody, vaccine, or therapeutic enzyme, can be truly effective for SUD treatment or not has been the subject of extensive debate. The main debate question is whether a biologic, particularly an exogenous enzyme, can effectively block the drug's reinforcing effect. In this report, we demonstrate that a modest dose of a recently redesigned long-acting cocaine hydrolase, CocH3-Fc(M6), can be used to effectively block the psychostimulant, discriminative stimulus, and reinforcing effects of cocaine for a sufficiently long period of time. For example, a dose of 3 mg/kg CocH3-Fc(M6) completely blocked the discriminative stimulus and reinforcing effects for 24/25 days and continued to significantly attenuate/decrease the cocaine effects for at least 29 days in rats. All the animal data consistently suggest that the long-acting cocaine hydrolase is a truly promising candidate of enzyme therapy for treatment of cocaine use disorder.

Structure-Based Design and Discovery of a Long-Acting Cocaine Hydrolase Mutant with Improved Binding Affinity to Neonatal Fc Receptor for Treatment of Cocaine Abuse

Despite decades of efforts to develop a pharmacotherapy for cocaine abuse treatment, there is still no FDA-approved treatment of diseases associated with this commonly abused drug. Our previously designed highly efficient cocaine hydrolases (CocHs) and the corresponding Fc-fusion proteins (e.g., CocH3-Fc) are recognized as potentially promising therapeutic enzyme candidates for cocaine abuse treatment, but all with limited biological half-lives. In order to prolong the biological half-life and, thus, decrease the required frequency of the enzyme administration for cocaine abuse treatment, we have modeled the Fc-fusion CocH binding with neonatal Fc receptor (FcRn) in the present study. This approach led to the design and testing of CocH3-Fc(M6), a CocH3-Fc mutant with nearly 100-fold increased binding affinity: from K_d = ~ 4 μM to K_d = 43 nM. As a result, CocH3-Fc(M6) indeed revealed a markedly prolonged biological half-life (t_1/2 = 206 ± 7 h or ~ 9 days) in rats, longer than other known Fc-fusion protein drugs such as abatacept and alefacept (for other therapeutic purposes) in the same species (rats). It has been demonstrated that a single dose of 3 mg/kg CocH3-Fc(M6) effectively blocked 20 mg/kg cocaine-induced hyperactivity on day 18 after CocH3-Fc(M6) administration. This is the first attempt to rationally design long-acting Fc-fusion enzyme mutant based on combined computational modeling and experimental measurement of the Fc-fusion CocH binding with FcRn. The similar structure-based design strategy may be used to prolong the biological half-lives of other Fc-fusion protein drugs.

Catalytic activities of cocaine hydrolases against the most toxic cocaine metabolite norcocaethylene

A majority of cocaine users also consume alcohol. The concurrent use of cocaine and alcohol produces the pharmacologically active metabolites cocaethylene and norcocaethylene, in addition to norcocaine. Both cocaethylene and norcocaethylene are more toxic than cocaine itself. Hence, a truly valuable cocaine-metabolizing enzyme for cocaine abuse/overdose treatment should be effective for the hydrolysis of not only cocaine, but also its metabolites norcocaine, cocaethylene, and norcocaethylene. However, there has been no report on enzymes capable of hydrolyzing norcocaethylene (the most toxic metabolite of cocaine). The catalytic efficiency parameters (kcat and KM) of human butyrylcholinesterase (BChE) and two mutants (known as cocaine hydrolases E14-3 and E12-7) against norcocaethylene have been characterized in the present study for the first time, and they are compared with those against cocaine. According to the obtained kinetic data, wild-type human BChE showed a similar catalytic efficiency against norcocaethylene (kcat = 9.5 min-1, KM = 11.7 μM, and kcat/KM = 8.12 × 105 M-1 min-1) to that against (-)-cocaine (kcat = 4.1 min-1, KM = 4.5 μM, and kcat/KM = 9.1 × 105 M-1 min-1). E14-3 and E12-7 showed an improved catalytic activity against norcocaethylene compared to wild-type BChE. E12-7 showed a 39-fold improved catalytic efficiency against norcocaethylene (kcat = 210 min-1, KM = 6.6 μM, and kcat/KM = 3.18 × 107 M-1 min-1). It has been demonstrated that E12-7 as an exogenous enzyme can efficiently metabolize norcocaethylene in rats.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

Plant-expressed cocaine hydrolase variants of butyrylcholinesterase exhibit altered allosteric effects of cholinesterase activity and increased inhibitor sensitivity

Butyrylcholinesterase (BChE) is an enzyme with broad substrate and ligand specificities and may function as a generalized bioscavenger by binding and/or hydrolyzing various xenobiotic agents and toxicants, many of which target the central and peripheral nervous systems. Variants of BChE were rationally designed to increase the enzyme’s ability to hydrolyze the psychoactive enantiomer of cocaine. These variants were cloned, and then expressed using the magnICON transient expression system in plants and their enzymatic properties were investigated. In particular, we explored the effects that these site-directed mutations have over the enzyme kinetics with various substrates of BChE. We further compared the affinity of various anticholinesterases including organophosphorous nerve agents and pesticides toward these BChE variants relative to the wild type enzyme. In addition to serving as a therapy for cocaine addiction-related diseases, enhanced bioscavenging against other harmful agents could add to the practicality and versatility of the plant-derived recombinant enzyme as a multivalent therapeutic.

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.

Metabolic Enzymes of Cocaine Metabolite Benzoylecgonine

Cocaine is one of the most addictive drugs without a U.S. Food and Drug Administration (FDA)-approved medication. Enzyme therapy using an efficient cocaine-metabolizing enzyme is recognized as the most promising approach to cocaine overdose treatment. The actual enzyme, known as RBP-8000, under current clinical development for cocaine overdose treatment is our previously designed T172R/G173Q mutant of bacterial cocaine esterase (CocE). The T172R/G173Q mutant is effective in hydrolyzing cocaine but inactive against benzoylecgonine (a major, biologically active metabolite of cocaine). Unlike cocaine itself, benzoylecgonine has an unusually stable zwitterion structure resistant to further hydrolysis in the body and environment. In fact, benzoylecgonine can last in the body for a very long time (a few days) and, thus, is responsible for the long-term toxicity of cocaine and a commonly used marker for drug addiction diagnosis in pre-employment drug tests. Because CocE and its mutants are all active against cocaine and inactive against benzoylecgonine, one might simply assume that other enzymes that are active against cocaine are also inactive against benzoylecgonine. Here, through combined computational modeling and experimental studies, we demonstrate for the first time that human butyrylcholinesterase (BChE) is actually active against benzoylecgonine, and that a rationally designed BChE mutant can not only more efficiently accelerate cocaine hydrolysis but also significantly hydrolyze benzoylecgonine in vitro and in vivo. This sets the stage for advanced studies to design more efficient mutant enzymes valuable for the development of an ideal cocaine overdose enzyme therapy and for benzoylecgonine detoxification in the environment.

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.

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.

Cocaine Hydrolase Gene Transfer Demonstrates Cardiac Safety and Efficacy against Cocaine-Induced QT Prolongation in Mice

Cocaine addiction is associated with devastating medical consequences, including cardiotoxicity and risk-conferring prolongation of the QT interval. Viral gene transfer of cocaine hydrolase engineered from butyrylcholinesterase offers therapeutic promise for treatment-seeking drug users. Although previous preclinical studies have demonstrated benefits of this strategy without signs of toxicity, the specific cardiac safety and efficacy of engineered butyrylcholinesterase viral delivery remains unknown. Here, telemetric recording of electrocardiograms from awake, unrestrained mice receiving a course of moderately large cocaine doses (30 mg/kg, twice daily for 3 weeks) revealed protection against a 2-fold prolongation of the QT interval conferred by pretreatment with cocaine hydrolase vector. By itself, this prophylactic treatment did not affect QT interval duration or cardiac structure, demonstrating that viral delivery of cocaine hydrolase has no intrinsic cardiac toxicity and, on the contrary, actively protects against cocaine-induced QT prolongation.

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.

Rational design, preparation, and characterization of a therapeutic enzyme mutant with improved stability and function for cocaine detoxification

Cocaine esterase (CocE) is known as the most efficient natural enzyme for cocaine hydrolysis. The major obstacle to the clinical application of wild-type CocE is the thermoinstability with a half-life of only ∼12 min at 37 °C. The previously designed T172R/G173Q mutant (denoted as enzyme E172-173) with an improved in vitro half-life of ∼6 h at 37 °C is currently in clinical trial Phase II for cocaine overdose treatment. Through molecular modeling and dynamics simulation, we designed and characterized a promising new mutant of E172-173 with extra L196C/I301C mutations (denoted as enzyme E196-301) to produce cross-subunit disulfide bonds that stabilize the dimer structure. The cross-subunit disulfide bonds were confirmed by X-ray diffraction. The designed L196C/I301C mutations have not only considerably extended the in vitro half-life at 37 °C to >100 days, but also significantly improved the catalytic efficiency against cocaine by ∼150%. In addition, the thermostable E196-301 can be PEGylated to significantly prolong the residence time in mice. The PEGylated E196-301 can fully protect mice from a lethal dose of cocaine (180 mg/kg, LD100) for at least 3 days, with an average protection time of ∼94h. This is the longest in vivo protection of mice from the lethal dose of cocaine demonstrated within all studies using an exogenous enzyme reported so far. Hence, E196-301 may be developed to become a more valuable therapeutic enzyme for cocaine abuse treatment, and it demonstrates that a general design strategy and protocol to simultaneously improve both the stability and function are feasible for rational protein drug design.

Kinetic characterization of human butyrylcholinesterase mutants for the hydrolysis of cocaethylene

It is known that the majority of cocaine users also consume alcohol. Alcohol can react with cocaine to produce a significantly more cytotoxic compound, cocaethylene. Hence a truly valuable cocaine-metabolizing enzyme as treatment for cocaine abuse/overdose should be efficient for not only cocaine itself, but also cocaethylene. The catalytic parameters (kcat and KM) of human BChE (butyrylcholinesterase) and two mutants (known as cocaine hydrolases E14-3 and E12-7) for cocaethylene are characterized in the present study, for the first time, in comparison with those for cocaine. On the basis of the obtained kinetic data, wild-type human BChE has a lower catalytic activity for cocaethylene (kcat=3.3 min(-1), KM=7.5 μM and kcat/KM=4.40 × 10(5) M(-1)·min(-1)) compared with its catalytic activity for (-)-cocaine. E14-3 and E12-7 have a considerably improved catalytic activity against cocaethylene compared with the wild-type BChE. E12-7 is identified as the most efficient enzyme for hydrolysing cocaethylene in addition to its high activity for (-)-cocaine. E12-7 has an 861-fold improved catalytic efficiency for cocaethylene (kcat=3600 min(-1), KM=9.5 μM and kcat/KM=3.79 × 10(8) M(-1)·min(-1)). It has been demonstrated that E12-7 as an exogenous enzyme can indeed rapidly metabolize cocaethylene in rats. Further kinetic modelling has suggested that E12-7 with an identical concentration as that of the endogenous BChE in human plasma can effectively eliminate (-)-cocaine, cocaethylene and norcocaine in simplified kinetic models of cocaine abuse and overdose associated with the concurrent use of cocaine and alcohol.

Physiologic and metabolic safety of butyrylcholinesterase gene therapy in mice

In continuing efforts to develop gene transfer of human butyrylcholinesterase (BChE) as therapy for cocaine addiction, we conducted wide-ranging studies of physiological and metabolic safety. For that purpose, mice were given injections of adeno-associated virus (AAV) vector or helper-dependent adenoviral (hdAD) vector encoding human or mouse BChE mutated for optimal cocaine hydrolysis. Age-matched controls received saline or AAV-luciferase control vector. At times when transduced BChE was abundant, physiologic and metabolic parameters in conscious animals were evaluated by non-invasive Echo-MRI and an automated "Comprehensive Laboratory Animal Monitoring System" (CLAMS). Despite high vector doses (up to 10(13) particles per mouse) and high levels of transgene protein in the plasma (∼1500-fold above baseline), the CLAMS apparatus revealed no adverse physiologic or metabolic effects. Likewise, body composition determined by Echo-MRI, and glucose tolerance remained normal. A CLAMS study of vector-treated mice given 40 mg/kg cocaine showed none of the physiologic and metabolic fluctuations exhibited in controls. We conclude that neither the tested vectors nor great excesses of circulating BChE affect general physiology directly, while they protect mice from disturbance by cocaine. Hence, viral gene transfer of BChE appears benign and worth exploring as a therapy for cocaine abuse and possibly other disorders as well.

A highly efficient cocaine-detoxifying enzyme obtained by computational design

Compared to naturally occurring enzymes, computationally designed enzymes are usually much less efficient, with their catalytic activities being more than six orders of magnitude below the diffusion limit. Here we use a two-step computational design approach, combined with experimental work, to design a highly efficient cocaine hydrolising enzyme. We engineer E30-6 from human butyrylcholinesterase (BChE), which is specific for cocaine hydrolysis, and obtain a much higher catalytic efficiency for cocaine conversion than for conversion of the natural BChE substrate, acetylcholine (ACh). The catalytic efficiency of E30-6 for cocaine hydrolysis is comparable to that of the most efficient known naturally-occurring hydrolytic enzyme, acetylcholinesterase, the catalytic activity of which approaches the diffusion limit. We further show that E30-6 can protect mice from a subsequently administered lethal dose of cocaine, suggesting the enzyme may have therapeutic potential in the setting of cocaine detoxification or cocaine abuse.

Kinetic characterization of high-activity mutants of human butyrylcholinesterase for the cocaine metabolite norcocaine

It has been known that cocaine produces its toxic and physiological effects through not only cocaine itself, but also norcocaine formed from cocaine oxidation catalysed by microsomal CYP (cytochrome P450) 3A4 in the human liver. The catalytic parameters (kcat and Km) of human BChE (butyrylcholinesterase) and its three mutants (i.e. A199S/S287G/A328W/Y332G, A199S/F227A/S287G/A328W/E441D and A199S/F227A/S287G/A328W/Y332G) for norcocaine have been characterized in the present study for the first time and compared with those for cocaine. On the basis of the obtained kinetic data, wild-type human BChE has a significantly lower catalytic activity for norcocaine (kcat=2.8 min(-1), Km=15 μM and kcat/Km=1.87 × 10(5) M(-1)·min(-1)) compared with its catalytic activity for (-)-cocaine. The BChE mutants examined in the present study have considerably improved catalytic activities against both cocaine and norcocaine compared with the wild-type enzyme. Within the enzymes examined in the present study, the A199S/F227A/S287G/A328W/Y332G mutant (CocH3) is identified as the most efficient enzyme for hydrolysing both cocaine and norcocaine. CocH3 has a 1080-fold improved catalytic efficiency for norcocaine (kcat=2610 min(-1), Km=13 μM and kcat/Km=2.01 × 10(8) M(-1)·min(-1)) and a 2020-fold improved catalytic efficiency for cocaine. It has been demonstrated that CocH3 as an exogenous enzyme can rapidly metabolize norcocaine, in addition to cocaine, in rats. Further kinetic modelling has suggested that CocH3 with an identical concentration with that of the endogenous BChE in human plasma can effectively eliminate both cocaine and norcocaine in a simplified kinetic model of cocaine abuse.

Substrate selectivity of high-activity mutants of human butyrylcholinesterase

Cocaine is one of the most addictive drugs, and there is still no FDA (Food and Drug Administration)-approved medication specific for cocaine abuse. A promising therapeutic strategy is to accelerate cocaine metabolism, producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e. cocaine hydrolysis catalyzed by butyrylcholinesterase (BChE) in plasma. However, the native BChE has a low catalytic efficiency against the abused cocaine, i.e. (-)-cocaine. Our recently designed and discovered A199S/F227A/S287G/A328W/Y332G mutant and other mutants of human BChE have a considerably improved catalytic efficiency against (-)-cocaine. In the present study, we carried out both computational modeling and experimental kinetic analysis on the catalytic activities of these promising new BChE mutants against other known substrates, including neurotransmitter acetylcholine (ACh), acetylthiocholine (ATC), butyrylthiocholine (BTC), and (+)-cocaine, in comparison with the corresponding catalytic activity against (-)-cocaine. Both the computational modeling and kinetic analysis have consistently revealed that all the examined amino acid mutations only considerably improve the catalytic efficiency of human BChE against (-)-cocaine, without significantly improving the catalytic efficiency of the enzyme against any of the other substrates examined. In particular, all the examined BChE mutants have a slightly lower catalytic efficiency against neurotransmitter ACh compared to the wild-type BChE. This observation gives us confidence in developing an anti-cocaine enzyme therapy by using one of these BChE mutants, particularly the A199S/F227A/S287G/A328W/Y332G mutant.

Modeling in vitro inhibition of butyrylcholinesterase using molecular docking, multi-linear regression and artificial neural network approaches

Butyrylcholinesterase (BChE) has been an important protein used for development of anti-cocaine medication. Through computational design, BChE mutants with ∼2000-fold improved catalytic efficiency against cocaine have been discovered in our lab. To study drug-enzyme interaction it is important to build mathematical model to predict molecular inhibitory activity against BChE. This report presents a neural network (NN) QSAR study, compared with multi-linear regression (MLR) and molecular docking, on a set of 93 small molecules that act as inhibitors of BChE by use of the inhibitory activities (pIC₅₀ values) of the molecules as target values. The statistical results for the linear model built from docking generated energy descriptors were: r(2)=0.67, rmsd=0.87, q(2)=0.65 and loormsd=0.90; the statistical results for the ligand-based MLR model were: r(2)=0.89, rmsd=0.51, q(2)=0.85 and loormsd=0.58; the statistical results for the ligand-based NN model were the best: r(2)=0.95, rmsd=0.33, q(2)=0.90 and loormsd=0.48, demonstrating that the NN is powerful in analysis of a set of complicated data. As BChE is also an established drug target to develop new treatment for Alzheimer's disease (AD). The developed QSAR models provide tools for rationalizing identification of potential BChE inhibitors or selection of compounds for synthesis in the discovery of novel effective inhibitors of BChE in the future.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Plants as a source of butyrylcholinesterase variants designed for enhanced cocaine hydrolase activity

Cocaine addiction affects millions of people with disastrous personal and social consequences. Cocaine is one of the most reinforcing of all drugs of abuse, and even those who undergo rehabilitation and experience long periods of abstinence have more than 80% chance of relapse. Yet there is no FDA-approved treatment to decrease the likelihood of relapse in rehabilitated addicts. Recent studies, however, have demonstrated a promising potential treatment option with the help of the serum enzyme butyrylcholinesterase (BChE), which is capable of breaking down naturally occurring (-)-cocaine before the drug can influence the reward centers of the brain or affect other areas of the body. This activity of wild-type (WT) BChE, however, is relatively low. This prompted the design of variants of BChE which exhibit significantly improved catalytic activity against (-)-cocaine. Plants are a promising means to produce large amounts of these cocaine hydrolase variants of BChE, cheaply, safely with no concerns regarding human pathogens and functionally equivalent to enzymes derived from other sources. Here, in expressing cocaine-hydrolyzing mutants of BChE in Nicotiana benthamiana using the MagnICON virus-assisted transient expression system, and in reporting their initial biochemical analysis, we provide proof-of-principle that plants can express engineered BChE proteins with desired properties.

Catalytic activities of a cocaine hydrolase engineered from human butyrylcholinesterase against (+)- and (-)-cocaine

It can be argued that an ideal anti-cocaine medication would be one that accelerates cocaine metabolism producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., hydrolysis catalyzed by butyrylcholinesterase (BChE) in plasma. However, wild-type BChE has a low catalytic efficiency against naturally occurring (-)-cocaine. Interestingly, wild-type BChE has a much higher catalytic activity against unnatural (+)-cocaine. According to available positron emission tomography (PET) imaging analysis using [(11)C](-)-cocaine and [(11)C](+)-cocaine tracers in human subjects, only [(11)C](-)-cocaine was observed in the brain, whereas no significant [(11)C](+)-cocaine signal was observed in the brain. The available PET data imply that an effective therapeutic enzyme for treatment of cocaine abuse could be an exogenous cocaine-metabolizing enzyme with a catalytic activity against (-)-cocaine comparable to that of wild-type BChE against (+)-cocaine. Our recently designed A199S/F227A/S287G/A328 W/Y332G mutant of human BChE has a considerably improved catalytic efficiency against (-)-cocaine and has been proven active in vivo. In the present study, we have characterized the catalytic activities of wild-type BChE and the A199S/F227A/S287G/A328 W/Y332G mutant against both (+)- and (-)-cocaine at the same time under the same experimental conditions. Based on the obtained kinetic data, the A199S/F227A/S287G/A328 W/Y332G mutant has a similarly high catalytic efficiency (kcat/KM) against (+)- and (-)-cocaine, and indeed has a catalytic efficiency (k(cat/)K(M) = 1.84 × 10(9) M(-1) min(-1)) against (-)-cocaine comparable to that (k(cat)/K(M) = 1.37 × 10(9) M(-1) min(-1)) of wild-type BChE against (+)-cocaine. Thus, the mutant may be used to effectively prevent (-)-cocaine from entering brain and producing physiological effects in the enzyme-based treatment of cocaine abuse.

Modeling of pharmacokinetics of cocaine in human reveals the feasibility for development of enzyme therapies for drugs of abuse

A promising strategy for drug abuse treatment is to accelerate the drug metabolism by administration of a drug-metabolizing enzyme. The question is how effectively an enzyme can actually prevent the drug from entering brain and producing physiological effects. In the present study, we have developed a pharmacokinetic model through a combined use of in vitro kinetic parameters and positron emission tomography data in human to examine the effects of a cocaine-metabolizing enzyme in plasma on the time course of cocaine in plasma and brain of human. Without an exogenous enzyme, cocaine half-lives in both brain and plasma are almost linearly dependent on the initial cocaine concentration in plasma. The threshold concentration of cocaine in brain required to produce physiological effects has been estimated to be 0.22±0.07 µM, and the threshold area under the cocaine concentration versus time curve (AUC) value in brain (denoted by AUC2_∞) required to produce physiological effects has been estimated to be 7.9±2.7 µM·min. It has been demonstrated that administration of a cocaine hydrolase/esterase (CocH/CocE) can considerably decrease the cocaine half-lives in both brain and plasma, the peak cocaine concentration in brain, and the AUC2_∞. The estimated maximum cocaine plasma concentration which a given concentration of drug-metabolizing enzyme can effectively prevent from entering brain and producing physiological effects can be used to guide future preclinical/clinical studies on cocaine-metabolizing enzymes. Understanding of drug-metabolizing enzymes is key to the science of pharmacokinetics. The general insights into the effects of a drug-metabolizing enzyme on drug kinetics in human should be valuable also in future development of enzyme therapies for other drugs of abuse.

In this computational study, we have examined, for the first time, the potential effects of a drug-metabolizing enzyme on drug pharmacokinetics in human, showing that a high-activity drug-metabolizing enzyme can completely/effectively prevent the drug of abuse from entering brain to produce physiological effects. Based on this encouraging insight, it is feasible to develop enzyme therapies for drugs of abuse. Through pharmacokinetic modeling, we have demonstrated that, without an exogenous enzyme, the drug half-lives in both brain and plasma are almost linearly dependent on the initial drug concentration in plasma. This finding indicates that one may not simply say the half-life of a drug without clearly indicating the actual dose condition. We have also demonstrated for the first time how a high-activity drug-metabolizing enzyme can considerably decrease the peak concentration of drug in brain and drug half-lives in both brain and plasma. In addition, we have calculated the minimum (threshold) concentration of cocaine in brain required to produce physiological effects. The predicted threshold concentration, along with all of the general insights obtained in this study, will provide a rational base for future design of further experimental studies required for the enzyme therapy development.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

His-tag truncated butyrylcholinesterase as a useful construct for in vitro characterization of wild-type and variant butyrylcholinesterases

Human butyrylcholinesterase (BChE) can scavenge and thereby provide protection against various toxic esters, including organophosphate-based chemical warfare agents and the recreational drug cocaine. It is currently being used in molecular evolution studies to generate novel enzymes with improved ability to hydrolyze toxic ester compounds. Currently, the most commonly used purification strategies for recombinant BChE enzymes involve using affinity resins based on small molecule interactions with the enzyme's substrate binding site. However, as BChE variants are discovered and developed, a generic purification protocol that is insensitive to amino acid substitutions is necessary. In the current manuscript, an expression vector encoding a C-terminal truncation and a His₆-tag was designed for BChE and used to express recombinant "wild-type" enzyme and two variants (i.e., G117H BChE and G117H/E197Q BChE). All the three His₆-tagged enzymes were successfully purified via metal-affinity columns using similar procedures with good recovery. Steady-state kinetic parameters were determined for each enzyme, and values were compared to those obtained with the corresponding non-truncated non-His₆-tagged enzymes. Rates of inhibition by echothiophate, a model compound for organophosphate-based pesticides, and rates of oxime-mediated reactivation after inhibition with a nerve agent model compound were also determined for selected enzymes. Rates of spontaneous reactivation from ETP inhibition were determined for the G117H variants. In all instances examined, truncation of the C-terminus of BChE and introduction of a His₆-tag had no significant effects on the observed kinetic parameters, making this a highly useful construct for in vitro characterization of wild-type and variant BChEs.

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

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to examine the catalytic mechanism for cocaine esterase (CocE)-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. It has been shown that the acylation of (+)-cocaine consists of nucleophilic attack of the hydroxyl group of Ser117 on the carbonyl carbon of (+)-cocaine benzoyl ester and the dissociation of (+)-cocaine benzoyl ester. The first reaction step of deacylation of (+)-cocaine, which is identical to that of (-)-cocaine, is rate-determining, indicating that CocE-catalyzed hydrolyses of (+)- and (-)-cocaine have a common rate-determining step. The computational results predict that the catalytic rate constant of CocE against (+)-cocaine should be the same as that of CocE against (-)-cocaine, in contrast with the remarkable difference between human butyrylcholinesterase-catalyzed hydrolyses of (+)- and (-)-cocaine. The prediction has been confirmed by experimental kinetic analysis on CocE-catalyzed hydrolysis of (+)-cocaine in comparison with CocE-catalyzed hydrolysis of (-)-cocaine. The determined common rate-determining step indicates that rational design of a high-activity mutant of CocE should be focused on the first reaction step of the deacylation. Furthermore, the obtained mechanistic insights into the detailed differences in the acylation between the (+)- and (-)-cocaine hydrolyses provide indirect clues for rational design of amino acid mutations that could more favorably stabilize the rate-determining transition state in the deacylation and, thus, improve the catalytic activity of CocE. This study provides a valuable mechanistic base for rational design of an improved esterase for therapeutic treatment of cocaine abuse.

Design, preparation, and characterization of high-activity mutants of human butyrylcholinesterase specific for detoxification of cocaine

Cocaine is a widely abused drug without a U.S. Food and Drug Administration-approved medication. There is a recognized, promising anticocaine medication to accelerate cocaine metabolism, producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway [i.e., cocaine hydrolysis catalyzed by butyrylcholinesterase (BChE) in plasma]. An ideal, therapeutically valuable mutant of human BChE should have not only a significantly improved catalytic activity against (-)-cocaine but also certain selectivity for (-)-cocaine over neurotransmitter acetylcholine (ACh), such that one would not expect systemic administration of the BChE mutant to interrupt cholinergic transmission. The present study accounting for the mutation-caused changes of the catalytic activities of BChE against both (-)-cocaine and ACh by means of molecular modeling and site-directed mutagenesis has led to identification of three BChE mutants that have not only a considerably improved catalytic efficiency against (-)-cocaine but also the desirable selectivity for (-)-cocaine over ACh. Two representative BChE mutants have been confirmed to be potent in actual protection of mice from acute toxicity (convulsion and lethality) of a lethal dose of cocaine (180 mg/kg). Pretreatment with the BChE mutant (i.e., 1 min before cocaine administration) dose-dependently protected mice against cocaine-induced convulsions and lethality. In particular, all mice pretreated with the mutant (e.g., 0.02 mg or more of A199S/F227A/S287G/A328W/E441D BChE) survived. The in vivo data reveal the primary factor (i.e., the relative catalytic efficiency), determining the efficacy in practical protection of mice from the acute cocaine toxicity and future direction for further improving the efficacy of the enzyme in the cocaine overdose treatment.

Enzyme-therapy approaches for the treatment of drug overdose and addiction

“Recent progress in the study of cocaine-metabolizing enzymes demonstrates that enzyme-therapy approaches using appropriately designed enzymes show promise for the treatment of drug overdose and addiction.”