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

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.

Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates

A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.

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.

Coenzyme Coupling Boosts Charge Transport through Single Bioactive Enzyme Junctions

Oxidation of formate to CO₂ is catalyzed via the donation of electrons from formate dehydrogenase (FDH) to nicotinamide adenine dinucleotide (NAD⁺), and thus the charge transport characteristics of FDH become essential but remain unexplored. Here, we investigated the charge transport through single-enzyme junctions of FDH using the scanning tunneling microscope break junction technique (STM-BJ). We found that the coupling of NAD⁺ with FDH boosts the charge transport by ∼2,100%, and the single-enzyme conductance highly correlates with the enzyme activity. The combined flicker noise analysis demonstrated the switching of the coenzyme-mediated charge transport pathway and supported by the significantly reduced HOMO-LUMO gap from calculations. Site-specific mutagenesis analysis demonstrated that FDH-NAD⁺ stably combined own higher bioactivity and boosts charge transport, and the coupling has been optimized via the natural selection. Our work provides evidence of hydrogen bond coupling in bioactivity but also bridges the charge transport through single-enzyme junctions and enzyme activities.

•

Binding of NAD⁺ with FDH boosts the charge transport by more than 2,100%

•

Single-enzyme conductance highly correlates with the enzyme activity

•

Hydrogen bond bridges the charge transport and enzyme activities

•

Experiments combined with calculations probe switching of charge transport pathway

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.

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.

Study of the affinity between the protein kinase PKA and homoarginine-containing peptides derived from kemptide: Free energy perturbation (FEP) calculations

Protein kinases (PKs) discriminate between closely related sequences that contain serine, threonine, and/or tyrosine residues. Such specificity is defined by the amino acid sequence surrounding the phosphorylatable residue, so that it is possible to identify an optimal recognition motif (ORM) for each PK. The ORM for the protein kinase A (PKA), a well-known member of the PK family, is the sequence RRX(S/T)X, where arginines at the -3 and -2 positions play a key role with respect to the primed phosphorylation site. In this work, differential affinities of PKA for the peptide substrate Kemptide (LRRASLG) and mutants that substitute the arginine residues by the unnatural peptide homoarginine were evaluated through molecular dynamics (MD) and free energy perturbation (FEP) calculations. The FEP study for the homoarginine mutants required previous elaboration of a CHARMM "arginine to homoarginine" (R2B) hybrid topology file which is available in this manuscript as Supporting Information. Mutants substituting the arginine residues by alanine, lysine, and histidine were also considered in the comparison by using the same protocol. FEP calculations allowed estimating the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES ) when Kemptide structure was mutated. Both ΔΔG_S→ES values for homoarginine mutants were predicted with a difference below 1 kcal/mol. In addition, FEP correctly predicted that all the studied mutations decrease the catalytic efficiency of Kemptide for PKA. © 2018 Wiley Periodicals, Inc.

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.

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.

Energy Landscape Topography Reveals the Underlying Link Between Binding Specificity and Activity of Enzymes

Enzyme activity (often quantified by k_cat/K_m) is the main function of enzyme when it is active against the specific substrate. Higher or lower activities are highly desired for the design of novel enzyme and drug resistance. However, it is difficult to measure the activities of all possible variants and find the “hot-spot” within the limit of experimental time. In this study, we explore the underlying energy landscape of enzyme-substrate interactions and introduce the intrinsic specificity ratio (ISR), which reflects the landscape topography. By studying two concrete systems, we uncover the statistical correlation between the intrinsic specificity and the enzyme activity k_cat/K_m. This physics-based concept and method show that the energy landscape topography is valuable for understanding the relationship between enzyme specificity and activity. In addition, it can reveal the underlying mechanism of enzyme-substrate actions and has potential applications on enzyme design.

Understanding the non-catalytic behavior of human butyrylcholinesterase silent variants: Comparison of wild-type enzyme, catalytically active Ala328Cys mutant, and silent Ala328Asp variant

Conformational dynamics of wild-type human butyrylcholinesterase (BChE), two mutants of residue Ala328, the catalytically active Ala328Cys, and the catalytically inactive (silent) Ala328Asp, and their interactions with butyrylcholine were studied. The aim was to understand the molecular mechanisms by which point mutations may lead to silent BChE variant or alter catalytic activity. Importance of BChE natural variants is due to medical consequences, i.e. prolonged apnea, following administration of the myorelaxant esters, succinylcholine and mivacurium. Comparison of molecular dynamics (MD) simulations for the three model systems showed that: 1) the active mutant Ala328Cys mutant has some changes in configuration of catalytic residues, which do not prevent binding of butyrylcholine to the active site; 2) in the naturally-occurring silent variant Ala328Asp, the Asp328 carboxylate may either form a salt bridge with Lys339 or a H-bond with His438. In the first case, the Ω-loop swings off the gorge, disrupting the π-cation binding site and the catalytic triad. In the second case, binding of cationic substrates in the catalytic center is also impaired. MD simulations carried out in 0.15 M NaCl, close to physiological ionic strength conditions, favored the second situation. It was seen that Asp328 forms a H-bond with the catalytic triad His438, which in turn disrupts the catalytic machinery. Therefore, we concluded that the Ala328Asp variant is not catalytically active because of that dramatic event. Computational results, consistent with in vitro biochemical data and clinical observations, validate our MD approach.

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.

Accurate Estimates of Free Energy Changes in Charge Mutations

The ability to determine the effect of charge changes on the free energy is necessary for fundamental studies of the electrostatic contribution to protein binding and stability. Currently, calculations of differences in free energy for charge mutations carried out in systems with periodic boundary conditions must include an approximate self-energy correction that can be on the same order of magnitude as the calculated free energy change. Here, a new method for accurately calculating the free energy change associated with any alchemical mutation, regardless of charge, is presented. In this method, paired mutations of opposite charge exactly cancel the self-energy term because of its quadratic charge dependence. Since the self-energy term implicitly cancels within the method, a correction never needs to be applied, and the statistical uncertainty associated is thereby removed. An implementation procedure is described and applied to the free energy of ionic hydration and a charged amino acid mutation.

Molecular mechanism of HIV-1 Tat interacting with human dopamine transporter

Nearly 70% of HIV-1-infected individuals suffer from HIV-associated neurocognitive disorders (HAND). HIV-1 transactivator of transcription (Tat) protein is known to synergize with abused drugs and exacerbate the progression of central nervous system (CNS) pathology. Cumulative evidence suggest that the HIV-1 Tat protein exerts the neurotoxicity through interaction with human dopamine transporter (hDAT) in the CNS. Through computational modeling and molecular dynamics (MD) simulations, we develop a three-dimensional (3D) structural model for HIV-1 Tat binding with hDAT. The model provides novel mechanistic insights concerning how HIV-1 Tat interacts with hDAT and inhibits dopamine uptake by hDAT. In particular, according to the computational modeling, Tat binds most favorably with the outward-open state of hDAT. Residues Y88, K92, and Y470 of hDAT are predicted to be key residues involved in the interaction between hDAT and Tat. The roles of these hDAT residues in the interaction with Tat are validated by experimental tests through site-directed mutagensis and dopamine uptake assays. The agreement between the computational and experimental data suggests that the computationally predicted hDAT-Tat binding mode and mechanistic insights are reasonable and provide a new starting point to design further pharmacological studies on the molecular mechanism of HIV-1-associated neurocognitive disorders.

Kinetic characterization of a cocaine hydrolase engineered from mouse butyrylcholinesterase

Mouse butyrylcholinesterase (mBChE) and an mBChE-based cocaine hydrolase (mCocH, i.e. the A¹⁹⁹S/S²²⁷A/S²⁸⁷G/A³²⁸W/Y³³²G mutant) have been characterized for their catalytic activities against cocaine, i.e. naturally occurring (-)-cocaine, in comparison with the corresponding human BChE (hBChE) and an hBChE-based cocaine hydrolase (hCocH, i.e. the A¹⁹⁹S/F²²⁷A/S²⁸⁷G/A³²⁸W/Y³³²G mutant). It has been demonstrated that mCocH and hCocH have improved the catalytic efficiency of mBChE and hBChE against (-)-cocaine by ~8- and ~2000-fold respectively, although the catalytic efficiencies of mCocH and hCocH against other substrates, including acetylcholine (ACh) and butyrylthiocholine (BTC), are close to those of the corresponding wild-type enzymes mBChE and hBChE. According to the kinetic data, the catalytic efficiency (k(cat)/K(M)) of mBChE against (-)-cocaine is comparable with that of hBChE, but the catalytic efficiency of mCocH against (-)-cocaine is remarkably lower than that of hCocH by ~250-fold. The remarkable difference in the catalytic activity between mCocH and hCocH is consistent with the difference between the enzyme-(-)-cocaine binding modes obtained from molecular modelling. Further, both mBChE and hBChE demonstrated substrate activation for all of the examined substrates [(-)-cocaine, ACh and BTC] at high concentrations, whereas both mCocH and hCocH showed substrate inhibition for all three substrates at high concentrations. The amino-acid mutations have remarkably converted substrate activation of the enzymes into substrate inhibition, implying that the rate-determining step of the reaction in mCocH and hCocH might be different from that in mBChE and hBChE.

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.

Amino-acid mutations to extend the biological half-life of a therapeutically valuable mutant of human butyrylcholinesterase

Cocaine is a widely abused and addictive drug without an FDA-approved medication. Our recently designed and discovered cocaine hydrolase, particularly E12-7 engineered from human butyrylcholinesterase (BChE), has the promise of becoming a valuable cocaine abuse treatment. An ideal anti-cocaine therapeutic enzyme should have not only a high catalytic efficiency against cocaine, but also a sufficiently long biological half-life. However, recombinant human BChE and the known BChE mutants have a much shorter biological half-life compared to the native human BChE. The present study aimed to extend the biological half-life of the cocaine hydrolase without changing its high catalytic activity against cocaine. Our strategy was to design possible amino-acid mutations that can introduce cross-subunit disulfide bond(s) and, thus, change the distribution of the oligomeric forms and extend the biological half-life. Three new BChE mutants (E364-532, E377-516, and E535) were predicted to have a more stable dimer structure with the desirable cross-subunit disulfide bond(s) and, therefore, a different distribution of the oligomeric forms and a prolonged biological half-life. The rational design was followed by experimental tests in vitro and in vivo, confirming that the rationally designed new BChE mutants, i.e. E364-532, E377-516, and E535, indeed had a remarkably different distribution of the oligomeric forms and prolonged biological half-life in rats from ∼7 to ∼13h without significantly changing the catalytic activity against (-)-cocaine. This is the first demonstration that rationally designed amino-acid mutations can significantly prolong the biological half-life of a high-activity enzyme without significantly changing the catalytic activity.

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.

A model of glycosylated human butyrylcholinesterase

Human butyrylcholinesterase (BChE) and its mutants have shown great potential in treating cocaine overdose and addiction. In order to effectively suppress cocaine reward in the brain for a long period of time after an exogenous cocaine hydrolase administration, the therapeutic enzyme should have not only a high catalytic efficiency against cocaine, but also a sufficiently long circulation time. It has been known that PEGylation (covalent attachment of polyethylene glycol) modification of a therapeutic protein can prolong the biological half-life of the protein without affecting its biological function. However, the asparagine-linked glycans on the surface of glycosylated BChE may interfere with the PEGylation modification. In this study, we built a three-dimensional (3D) model of glycosylated human BChE to investigate the influence of glycans on the PEGylation modification. Glycans did not change the overall stability of the BChE structure, but could increase the flexibility of some local structures. For further evaluating the accessibility of the PEGylation reaction sites, particularly lysine residues, on the protein surface, we calculated the Solvent Accessible Surface Areas (SASAs) of these residues. The results indicate that some lysine residues show a significant decrease in SASA due to the direct or indirect influence of their surrounding glycans. The results also indicate that PEGylation reaction agents with smaller functional groups could have a better chance to react with lysine residues. This investigation provides a structural basis for rational engineering of human BChE and its mutants as therapeutic candidates.

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.

Preparation and in vivo characterization of a cocaine hydrolase engineered from human butyrylcholinesterase for metabolizing cocaine

Cocaine is a widely abused drug without an FDA (Food and Drug Administration)-approved medication. It has been recognized that an ideal anti-cocaine medication would accelerate cocaine metabolism producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e. human BChE (butyrylcholinesterase)-catalysed hydrolysis. However, the native human BChE has a low catalytic activity against cocaine. We recently designed and discovered a BChE mutant (A199S/F227A/S287G/A328W/Y332G) with a high catalytic activity (kcat=5700 min-1, Km=3.1 μM) specifically for cocaine, and the mutant was proven effective in protecting mice from acute cocaine toxicity of a lethal dose of cocaine (180 mg/kg of body weight, LD100). Further characterization in animal models requires establishment of a high-efficiency stable cell line for the BChE mutant production at a relatively larger scale. It has been extremely challenging to develop a high-efficiency stable cell line expressing BChE or its mutant. In the present study, we successfully developed a stable cell line efficiently expressing the BChE mutant by using a lentivirus-based repeated-transduction method. The scaled-up protein production enabled us to determine for the first time the in vivo catalytic activity and the biological half-life of this high-activity mutant of human BChE in accelerating cocaine clearance. In particular, it has been demonstrated that the BChE mutant (administered to mice 1 min prior to cocaine) can quickly metabolize cocaine and completely eliminate cocaine-induced hyperactivity in rodents, implying that the BChE mutant may be developed as a promising therapeutic agent for cocaine abuse treatment.

Accurate calculation of mutational effects on the thermodynamics of inhibitor binding to p38α MAP kinase: a combined computational and experimental study

A major current challenge for drug design efforts focused on protein kinases is the development of drug resistance caused by spontaneous mutations in the kinase catalytic domain. The ubiquity of this problem means that it would be advantageous to develop fast, effective computational methods that could be used to determine the effects of potential resistance-causing mutations before they arise in a clinical setting. With this long-term goal in mind, we have conducted a combined experimental and computational study of the thermodynamic effects of active-site mutations on a well-characterized and high-affinity interaction between a protein kinase and a small-molecule inhibitor. Specifically, we developed a fluorescence-based assay to measure the binding free energy of the small-molecule inhibitor, SB203580, to the p38α MAP kinase and used it measure the inhibitor's affinity for five different kinase mutants involving two residues (Val38 and Ala51) that contact the inhibitor in the crystal structure of the inhibitor-kinase complex. We then conducted long, explicit-solvent thermodynamic integration (TI) simulations in an attempt to reproduce the experimental relative binding affinities of the inhibitor for the five mutants; in total, a combined simulation time of 18.5 μs was obtained. Two widely used force fields - OPLS-AA/L and Amber ff99SB-ILDN - were tested in the TI simulations. Both force fields produced excellent agreement with experiment for three of the five mutants; simulations performed with the OPLS-AA/L force field, however, produced qualitatively incorrect results for the constructs that contained an A51V mutation. Interestingly, the discrepancies with the OPLS-AA/L force field could be rectified by the imposition of position restraints on the atoms of the protein backbone and the inhibitor without destroying the agreement for other mutations; the ability to reproduce experiment depended, however, upon the strength of the restraints' force constant. Imposition of position restraints in corresponding simulations that used the Amber ff99SB-ILDN force field had little effect on their ability to match experiment. Overall, the study shows that both force fields can work well for predicting the effects of active-site mutations on small molecule binding affinities and demonstrates how a direct combination of experiment and computation can be a powerful strategy for developing an understanding of protein-inhibitor interactions.

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.

Interception of cocaine by enzyme or antibody delivered with viral gene transfer: a novel strategy for preventing relapse in recovering drug users

Recent progress in enzyme engineering has led to versions of human butyrylcholinesterase (BChE) that hydrolyze cocaine efficiently in plasma, reduce concentrations reaching reward neurocircuity in the brain, and weaken behavioral responses to this drug. Along with enzyme advances, increasingly avid anti-cocaine antibodies and potent anti-cocaine vaccines have also been developed. Here we review these developments and consider the potential advantages along with the risks of delivering drug-intercepting proteins via gene transfer approaches to treat cocaine addiction.

Computational approaches for rational design of proteins with novel functionalities

Proteins are the most multifaceted macromolecules in living systems and have various important functions, including structural, catalytic, sensory, and regulatory functions. Rational design of enzymes is a great challenge to our understanding of protein structure and physical chemistry and has numerous potential applications. Protein design algorithms have been applied to design or engineer proteins that fold, fold faster, catalyze, catalyze faster, signal, and adopt preferred conformational states. The field of de novo protein design, although only a few decades old, is beginning to produce exciting results. Developments in this field are already having a significant impact on biotechnology and chemical biology. The application of powerful computational methods for functional protein designing has recently succeeded at engineering target activities. Here, we review recently reported de novo functional proteins that were developed using various protein design approaches, including rational design, computational optimization, and selection from combinatorial libraries, highlighting recent advances and successes.

Combined cocaine hydrolase gene transfer and anti-cocaine vaccine synergistically block cocaine-induced locomotion

Mice and rats were tested for reduced sensitivity to cocaine-induced hyper-locomotion after pretreatment with anti-cocaine antibody or cocaine hydrolase (CocH) derived from human butyrylcholinesterase (BChE). In Balb/c mice, direct i.p. injection of CocH protein (1 mg/kg) had no effect on spontaneous locomotion, but it suppressed responses to i.p. cocaine up to 80 mg/kg. When CocH was injected i.p. along with a murine cocaine antiserum that also did not affect spontaneous locomotion, there was no response to any cocaine dose. This suppression of locomotor activity required active enzyme, as it was lost after pretreatment with iso-OMPA, a selective BChE inhibitor. Comparable results were obtained in rats that developed high levels of CocH by gene transfer with helper-dependent adenoviral vector, and/or high levels of anti-cocaine antibody by vaccination with norcocaine hapten conjugated to keyhole limpet hemocyanin (KLH). After these treatments, rats were subjected to a locomotor sensitization paradigm involving a “training phase" with an initial i.p. saline injection on day 1 followed by 8 days of repeated cocaine injections (10 mg/kg, i.p.). A 15-day rest period then ensued, followed by a final “challenge" cocaine injection. As in mice, the individual treatment interventions reduced cocaine-stimulated hyperactivity to a modest extent, while combined treatment produced a greater reduction during all phases of testing compared to control rats (with only saline pretreatment). Overall, the present results strongly support the view that anti-cocaine vaccine and cocaine hydrolase vector treatments together provide enhanced protection against the stimulatory actions of cocaine in rodents. A similar combination therapy in human cocaine users might provide a robust therapy to help maintain abstinence.

Cocaine hydrolase gene therapy for cocaine abuse

Rapid progress in the past decade with re-engineering of human plasma butyrylcholinesterase has led to enzymes that destroy cocaine so efficiently that they prevent or interrupt drug actions in the CNS even though confined to the blood stream. Over the same time window, improved gene-transfer technology has made it possible to deliver such enzymes by endogenous gene transduction at high levels for periods of a year or longer after a single treatment. This article reviews recent advances in this field and considers prospects for development of a robust therapy aimed at aiding recovering drug users avoid addiction relapse.

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.

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.

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).

Uricases as therapeutic agents to treat refractory gout: Current states and future directions

Treatment of refractory gout remains a challenge on drug development. While pegloticase, a recombinant mammalian uricase modified with monomethoxyl-poly(ethylene glycol) (mPEG) is effective in treating refractory gout, after continued treatment for three months biweekly at a therapeutic dose of 0.14 mg/kg body weight, it elicits an immune response against mPEG in nearly 20% of patients. For continued treatment of refractory gout PEGylated uricases at monthly therapeutic doses below 4 μg/kg body weight have promise. To formulate uricases to achieve monthly therapeutic regimens requires pharmacodynamics simulation and experimentation including: (a) molecular engineering of uricases based on rational design and evolution biotechnology in combination to improve their inherent catalytic efficiency, thermostability and selectivity for urate over xanthine and; (b) optimization of the number and distribution of accessible reactive amino acid residues in native uricases for site-specific PEGylation with PEG derivatives with lower of immunogenicity than mPEG to retain activity, minimize immunogenicity and enhance the pharmacokinetics of the PEGylated uricase. These issues are briefly reviewed as a means to stimulate the development of safer uricase formulations for continued treatment of refractory gout.

Recent progress in protein drug design and discovery with a focus on novel approaches to the development of anti-cocaine medications

Cocaine is highly addictive and no anti-cocaine medication is currently available. Accelerating cocaine metabolism, producing biologically inactive metabolites, is recognized as an ideal anti-cocaine medication strategy, especially for the treatment of acute cocaine toxicity. However, currently known wild-type enzymes have either too low a catalytic efficiency against the abused cocaine, in other words (-)-cocaine, or the in vivo half-life is too short. Novel computational strategies and design approaches have been developed recently to design and discover thermostable or high-activity mutants of enzymes based on detailed structures and catalytic/inactivation mechanisms. The structure- and mechanism-based computational design efforts have led to the discovery of high-activity mutants of butyrylcholinesterase and thermostable mutants of cocaine esterase as promising anti-cocaine therapeutics. The structure- and mechanism-based computational strategies and design approaches may be used to design high-activity and/or thermostable mutants of many other proteins that have clear therapeutic potentials and to design completely new therapeutic enzymes.

The evolution of catalytic function in the HIV-1 protease

The evolution of species is a complex phenomenon based on the optimization of a multidimensional function referred to as fitness. At the level of biomolecular evolution, the fitness function can be reduced to include physiochemical properties relevant to the biological function of a particular molecule. In this work, questions involving the physical-chemical mechanisms underlying the evolution of HIV-1 protease are addressed through molecular simulation and subsequent analysis of thermodynamic properties related to the activity of the enzyme. Specifically, the impact of 40 single amino acid mutations on the binding affinity toward the matrix/capsid (MA/CA) substrate and corresponding transition state intermediate has been characterized using a molecular mechanics Poisson-Boltzmann surface area approach. We demonstrate that this approach is capable of extracting statistically significant information relevant to experimentally determined catalytic activity. Further, no correlation was observed between the effect of mutations on substrate and transition state binding, suggesting independent evolutionary pathways toward optimizing substrate specificity and catalytic activity. In addition, a detailed analysis of calculated binding affinity data suggests that ground-state destabilization (reduced binding affinity for the substrate) could be a contributing factor in the evolutionary optimization of HIV-1 protease. A numerical model is developed to demonstrate that ground-state destabilization is a valid mechanism for activity optimization given the high concentrations of substrate experienced by the functional enzyme in vivo.

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.

Reaction pathway and free energy profile for prechemical reaction step of human butyrylcholinesterase-catalyzed hydrolysis of (-)-cocaine by combined targeted molecular dynamics and potential of mean force simulations

Combined targeted molecular dynamics (TMD) and potential of mean force (PMF) simulations have been carried out to uncover the detailed pathway and determine the corresponding free energy profile for the structural transformation from the nonprereactive butyrylcholinesterase (BChE)-(-)-cocaine binding to the prereactive BChE-(-)-cocaine binding associated with the (-)-cocaine rotation in the binding pocket of BChE. It has been shown that the structural transformation involves two transition states (TS1(rot) and TS2(rot)). TS1(rot) is mainly associated with the deformation of the nonprereactive complex, whereas TS2(rot) is mainly associated with the formation of the prereactive complex. It has also been demonstrated that the A328W/Y332G mutation significantly reduces the steric hindrance for (-)-cocaine rotation in the binding pocket of BChE and, thus, decreases the free energy barrier for the structural transformation from the nonprereactive binding to the prereactive binding. The calculated relative free energy barriers are all consistent with available experimental kinetic data. The new mechanistic insights obtained and the novel computational protocol tested in this study should be valuable for future computational design of high-activity mutants of BChE. The general computational strategy and approach based on the combined TMD and PMF simulations may be also valuable in computational studies of detailed pathways and free energy profiles for other similar mechanistic problems involving ligand rotation or another type of structural transformation in the binding pocket of a protein.

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.”

Reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of acetylcholine

A catalytic mechanism for the butyrylcholinesterase (BChE)-catalyzed hydrolysis of acetylcholine (ACh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy calculations on both acylation and deacylation of BChE. It has been shown that the acylation with ACh includes two reaction steps, including nucleophilic attack on the carbonyl carbon of ACh and dissociation of choline ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the hydroxyl oxygen of the Ser198 side chain. Notably, despite the fact that acetylcholinesterase (AChE) and BChE are very similar enzymes, the acylation of BChE with ACh is rate-determining, which is remarkably different from the AChE-catalyzed hydrolysis of ACh, in which the deacylation is rate-determining. The computational prediction is consistent with available experimental kinetic data. The overall free energy barrier calculated for BChE-catalyzed hydrolysis of ACh is 13.8 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 13.3 kcal/mol.

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

The present study was aimed to explore the correlation between the protein structure and catalytic efficiency of butyrylcholinesterase (BChE) mutants against (-)-cocaine by modeling the rate-determining transition state (TS1), i.e., the transition state for the first step of chemical reaction process, of (-)-cocaine hydrolysis catalyzed by various mutants of human BChE in comparison with the wild type. Molecular modeling of the TS1 structures revealed that mutations on certain nonactive site residues can indirectly affect the catalytic efficiency of the enzyme against (-)-cocaine through enhancing or weakening the overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. Computational insights and predictions were supported by the catalytic activity data obtained from wet experimental tests on the mutants of human BChE, including five new mutants reported for the first time. The BChE mutants with at least ∼1000-fold improved catalytic efficiency against (-)-cocaine compared to the wild-type BChE are all associated with the TS1 structures having stronger overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. The combined computational and experimental data demonstrate a reasonable correlation relationship between the hydrogen-bonding distances in the TS1 structure and the catalytic efficiency of the enzyme against (-)-cocaine.

Characterization of a high-activity mutant of human butyrylcholinesterase against (-)-cocaine

Cocaine addiction and overdose are a well-known public health problem. There is no approved medication available for cocaine abuse treatment. Our recently designed and discovered high-activity mutant (A199S/S287G/A328W/Y332G) of human butyrylcholinesterase (BChE) has been recognized to be worth exploring for clinical application in humans as a potential anti-cocaine medication. The catalytic rate constant (k(cat)) and Michaelis-Menten constant (K(M)) for (-)-cocaine hydrolysis catalyzed by A199S/S287G/A328W/Y332G BChE (without fusion with any other peptide) have been determined to be 3,060 min(-1) and 3.1 microM, respectively, in the present study. The determined kinetic parameters reveal that the un-fused A199S/S287G/A328W/Y332G mutant has a approximately 1,080-fold improved catalytic efficiency (k(cat)/K(M)) against (-)-cocaine compared to the wild-type BChE. The approximately 1,080-fold improvement in the catalytic efficiency of the un-fused A199S/S287G/A328W/Y332G mutant is very close to the previously reported the approximately 1,000-fold improvement in the catalytic efficiency of the A199S/S287G/A328W/Y332G mutant fused with human serum albumin. These results suggest that the albumin fusion did not significantly change the catalytic efficiency of the BChE mutant while extending the plasma half-life. In addition, we have also examined the catalytic activities of the A199S/S287G/A328W/Y332G mutant against two other substrates, acetylthiocholine (ATC) and butyrylthiocholine (BTC). It has been shown that the A199S/S287G/A328W/Y332G mutations actually decreased the catalytic efficiencies of BChE against ATC and BTC, while considerably improving the catalytic efficiency of BChE against (-)-cocaine.