Computational redesign of human butyrylcholinesterase for anticocaine medication

Stereoselectivity of phosphotriesterase with paraoxon derivatives: a computational study

The bacterial enzyme phosphotriesterase (PTE) exhibits stereoselectivity toward hydrolysis of chiral substrates with a preference for the Sp enantiomer. In this work, docking analysis and two explicit-solvent molecular dynamics (MD) simulations were performed to characterize and differentiate the structural dynamics of PTE bound to the Sp and Rp paraoxon derivative enantiomers (Rp-1 and Sp-1) hydrolyzed with distinct catalytic efficiencies. Comparative analysis of the molecular trajectories for PTE bound to Rp-1 and Sp-1 suggested that substrate binding induced conformational changes in the loops near the active site. After 100 ns of MD simulation, the Zn β(2+) metal ion formed hexacoordinated- and tetracoordinated geometries in the Sp-1-PTE and Rp-1-PTE ensembles, respectively. Simulation results further showed that the hydrogen bond between Asp301 and His254 occurred with a higher probability after Sp-1 binding to PTE (47.5%) than that after Rp-1 binding (22.2%). These results provide a qualitative and molecular-level explanation for the 10 orders of magnitude increase in the catalytic efficiency of PTE toward the Sp enantiomer of paraoxon.

Metabolism of β-valine via a CoA-dependent ammonia lyase pathway

Pseudomonas species strain SBV1 can rapidly grow on medium containing β-valine as a sole nitrogen source. The tertiary amine feature of β-valine prevents direct deamination reactions catalyzed by aminotransferases, amino acid dehydrogenases, and amino acid oxidases. However, lyase- or aminomutase-mediated conversions would be possible. To identify enzymes involved in the degradation of β-valine, a PsSBV1 gene library was prepared and used to complement the β-valine growth deficiency of a closely related Pseudomonas strain. This resulted in the identification of a gene encoding β-valinyl-coenzyme A ligase (BvaA) and two genes encoding β-valinyl-CoA ammonia lyases (BvaB1 and BvaB2). The BvaA protein demonstrated high sequence identity to several known phenylacetate CoA ligases. Purified BvaA enzyme did not convert phenyl acetic acid but was able to activate β-valine in an adenosine triphosphate (ATP)- and CoA-dependent manner. The substrate range of the enzyme appears to be narrow, converting only β-valine and to a lesser extent, 3-aminobutyrate and β-alanine. Characterization of BvaB1 and BvaB2 revealed that both enzymes were able to deaminate β-valinyl-CoA to produce 3-methylcrotonyl-CoA, a common intermediate in the leucine degradation pathway. Interestingly, BvaB1 and BvaB2 demonstrated no significant sequence identity to known CoA-dependent ammonia lyases, suggesting they belong to a new family of enzymes. BLAST searches revealed that BvaB1 and BvaB2 show high sequence identity to each other and to several enoyl-CoA hydratases, a class of enzymes that catalyze a similar reaction with water instead of amine as the leaving group.

Safety, pharmacokinetics, and pharmacodynamics of TV-1380, a novel mutated butyrylcholinesterase treatment for cocaine addiction, after single and multiple intramuscular injections in healthy subjects

Human plasma butyrylcholinesterase (BChE) contributes to cocaine metabolism and has been considered for use in treating cocaine addiction and cocaine overdose. TV-1380 is a recombinant protein composed of the mature form of human serum albumin fused at its amino terminus to the carboxy-terminus of a truncated and mutated BChE. In preclinical studies, TV-1380 has been shown to rapidly eliminate cocaine in the plasma thus forestalling entry of cocaine into the brain and heart. Two randomized, blinded phase I studies were conducted to evaluate the safety, pharmacokinetics, and pharmacodynamics of TV-1380, following single and multiple administration in healthy subjects. TV-1380 was found to be safe and well tolerated with a long half-life (43-77 hours) and showed a dose-proportional increase in systemic exposure. Consistent with preclinical results, the ex vivo cocaine hydrolysis, TV-1380 activity clearly increased upon treatment in a dose-dependent manner. In addition, there was a direct relationship between ex vivo cocaine hydrolysis (kel ) and TV-1380 serum concentrations. There was no evidence that TV-1380 affected heart rate, the uncorrected QT interval, or the heart-rate-corrected QTcF interval. TV-1380, therefore, offers a safe once-weekly therapy to increase cocaine hydrolysis.

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.

Computationally efficient and accurate enantioselectivity modeling by clusters of molecular dynamics simulations

Computational approaches could decrease the need for the laborious high-throughput experimental screening that is often required to improve enzymes by mutagenesis. Here, we report that using multiple short molecular dynamics (MD) simulations makes it possible to accurately model enantioselectivity for large numbers of enzyme-substrate combinations at low computational costs. We chose four different haloalkane dehalogenases as model systems because of the availability of a large set of experimental data on the enantioselective conversion of 45 different substrates. To model the enantioselectivity, we quantified the frequency of occurrence of catalytically productive conformations (near attack conformations) for pairs of enantiomers during MD simulations. We found that the angle of nucleophilic attack that leads to carbon-halogen bond cleavage was a critical variable that limited the occurrence of productive conformations; enantiomers for which this angle reached values close to 180° were preferentially converted. A cluster of 20-40 very short (10 ps) MD simulations allowed adequate conformational sampling and resulted in much better agreement to experimental enantioselectivities than single long MD simulations (22 ns), while the computational costs were 50-100 fold lower. With single long MD simulations, the dynamics of enzyme-substrate complexes remained confined to a conformational subspace that rarely changed significantly, whereas with multiple short MD simulations a larger diversity of conformations of enzyme-substrate complexes was observed.

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.

Long-term reduction of cocaine self-administration in rats treated with adenoviral vector-delivered cocaine hydrolase: evidence for enzymatic activity

A new pharmacokinetic approach treating cocaine addiction involves rapidly metabolizing cocaine before it reaches brain reward centers using mutated human butyrylcholinesterase (BChE) or cocaine hydrolase (CocH). Recent work has shown that helper-dependent adenoviral (hdAD) vector-mediated plasma CocH reduced the locomotor-activating effects of cocaine and prevented reinstatement of cocaine-seeking behavior up to 6 months in rats. The present study investigated whether hdAD-CocH could decrease ongoing intravenous cocaine (0.4 mg/kg) self-administration. The hdAD-CocH vector was injected into self-administering rats, and after accumulation of plasma CocH, there was a dramatic reduction in cocaine infusions earned under a fixed ratio 1 schedule of reinforcement that lasted for the length of the study (>2 months). Pretreatment with the selective BChE and CocH inhibitor iso-OMPA (1.5 mg/kg) restored cocaine intake; therefore, the decline in self-administration was likely due to rapid CocH-mediated cocaine metabolism. Direct measurements of cocaine levels in plasma and brain samples taken after the conclusion of behavioral studies provided strong support for this conclusion. Further, rats injected with hdAD-CocH did not experience a deficit in operant responding for drug reinforcement and self-administered methamphetamine (0.05 mg/kg) at control levels. Overall, these outcomes suggest that viral gene transfer can yield plasma CocH levels that effectively diminish long-term cocaine intake and may have potential treatment implications for cocaine-dependent individuals seeking to become and remain abstinent.

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.

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.

Pharmacokinetic approaches to treatment of drug addiction

The pharmacokinetic approach to treatment targets the drug molecules themselves, aiming to reduce their concentration at the site of action, thereby reducing or preventing any pharmacodynamic effect. This approach might be useful in the treatment of acute drug toxicity/overdose and in the long-term treatment of addiction. Early clinical trials with anticocaine and antinicotine vaccines have shown reduced drug use and good tolerability. Also showing promise in animal studies are monoclonal antibodies against cocaine, methamphetamine and phencyclidine, as well as the enhancment of cocaine metabolism with genetic variants of human butyrylcholinesterase, using a bacterial esterase or catalytic monoclonal antibodies. Pharmacokinetic treatments offer potential advantages in terms of patient compliance, absence of medication interactions and benefit for patients who cannot take standard medications.

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.

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.

Prospects, promise and problems on the road to effective vaccines and related therapies for substance abuse

This review addresses potential new treatments for stimulant drugs of abuse, especially cocaine. Clinical trials of vaccines against cocaine and nicotine have been completed with the generally encouraging result that subjects showing high titers of antidrug antibody experience a reduction in drug reward, which may aid in cessation. New vaccine technologies, including gene transfer of highly optimized monoclonal antibodies, are likely to improve such outcomes further. In the special case of cocaine abuse, a metabolic enzyme is emerging as an alternative or added therapeutic intervention, which would also involve gene transfer. Such approaches still require extensive studies of safety and efficacy, but they may eventually contribute to a robust form of in vivo drug interception that greatly reduces the risks of addiction relapse.

Computational and experimental insights into the mechanism of substrate recognition and feedback inhibition of protoporphyrinogen oxidase

Protoporphyrinogen IX oxidase (PPO; EC 1.3.3.4) is an essential enzyme catalyzing the last common step in the pathway leading to heme and chlorophyll biosynthesis. Great interest in PPO inhibitors arises from both its significance to agriculture and medicine. However, the discovery of PPO inhibitors with ultrahigh potency and selectivity is hampered due to lack of structural and mechanistic understanding about the substrate recognition, which remains a longstanding question central in porphyrin biology. To understand the mechanism, a novel binding model of protogen (protoporphyrinogen IX, the substrate) was developed through extensive computational simulations. Subsequently, amino acid residues that are critical for protogen binding identified by computational simulations were substituted by mutagenesis. Kinetic analyses of these mutants indicated that these residues were critical for protogen binding. In addition, the calculated free energies of protogen binding with these mutants correlated well with the experimental data, indicating the reasonability of the binding model. On the basis of this novel model, the fundamental mechanism of substrate recognition was investigated by performing potential of mean force (PMF) calculations, which provided an atomic level description of conformational changes and pathway intermediates. The free energy profile revealed a feedback inhibition mechanism of proto (protoporphyrin IX, the product), which was also in agreement with experimental evidence. The novel mechanistic insights obtained from this study present a new starting point for future rational design of more efficient PPO inhibitors based on the product-bound PPO structure.

Structure-based engineering of streptavidin monomer with a reduced biotin dissociation rate

We recently reported the engineering of monomeric streptavidin, mSA, corresponding to one subunit of wild type (wt) streptavidin tetramer. The monomer was designed by homology modeling, in which the streptavidin and rhizavidin sequences were combined to engineer a high affinity binding pocket containing residues from a single subunit only. Although mSA is stable and binds biotin with nanomolar affinity, its fast off rate (koff ) creates practical challenges during applications. We obtained a 1.9 Å crystal structure of mSA bound to biotin to understand their interaction in detail, and used the structure to introduce targeted mutations to improve its binding kinetics. To this end, we compared mSA to shwanavidin, which contains a hydrophobic lid containing F43 in the binding pocket and binds biotin tightly. However, the T48F mutation in mSA, which introduces a comparable hydrophobic lid, only resulted in a modest 20-40% improvement in the measured koff . On the other hand, introducing the S25H mutation near the bicyclic ring of bound biotin increased the dissociation half life (t½ ) from 11 to 83 min at 20°C. Molecular dynamics (MD) simulations suggest that H25 stabilizes the binding loop L3,4 by interacting with A47, and protects key intermolecular hydrogen bonds by limiting solvent entry into the binding pocket. Concurrent T48F or T48W mutation clashes with H25 and partially abrogates the beneficial effects of H25. Taken together, this study suggests that stabilization of the binding loop and solvation of the binding pocket are important determinants of the dissociation kinetics in mSA.

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.

Effects of anti-cocaine vaccine and viral gene transfer of cocaine hydrolase in mice on cocaine toxicity including motor strength and liver damage

In developing an vivo drug-interception therapy to treat cocaine abuse and hinder relapse into drug seeking provoked by re-encounter with cocaine, two promising agents are: (1) a cocaine hydrolase enzyme (CocH) derived from human butyrylcholinesterase and delivered by gene transfer; (2) an anti-cocaine antibody elicited by vaccination. Recent behavioral experiments showed that antibody and enzyme work in a complementary fashion to reduce cocaine-stimulated locomotor activity in rats and mice. Our present goal was to test protection against liver damage and muscle weakness in mice challenged with massive doses of cocaine at or near the LD50 level (100-120 mg/kg, i.p.). We found that, when the interceptor proteins were combined at doses that were only modestly protective in isolation (enzyme, 1mg/kg; antibody, 8 mg/kg), they provided complete protection of liver tissue and motor function. When the enzyme levels were ~400-fold higher, after in vivo transduction by adeno-associated viral vector, similar protection was observed from CocH alone.

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.

Modification of pharmacokinetic and abuse-related effects of cocaine by human-derived cocaine hydrolase in monkeys

Although substantial research effort has focused on developing pharmacological treatments for cocaine abuse, no effective medications have been developed. Recent studies show that enzymes that metabolize cocaine in the periphery, forestalling its entry into the brain, can prevent cocaine toxicity and its behavioral effects in rodents. Here we report on effects of one such enzyme (Albu-CocH) on the pharmacokinetic and behavioral effects of cocaine in squirrel monkeys. Albu-CocH was developed from successive mutations of human butyrylcholinesterase (BChE) and has 1000-fold greater catalytic activity against cocaine than naturally occurring BChE. Pharmacokinetic studies showed that Albu-CocH (5 mg/kg) had a half-life of 56.6 hours in squirrel monkeys. In these studies, plasma levels of cocaine following i.v. 1 mg/kg cocaine were reduced 2 hours after administration of Albu-CocH, whereas plasma levels of the cocaine metabolite ecgonine methyl ester were increased. These effects were still evident 72 hours following Albu-CocH administration. In behavioral experiments in monkeys, pre-treatment with 5 mg/kg Albu-CocH dramatically decreased self-administration of a reinforcing dose of i.v. cocaine (30 µg/kg/injection) for over 24 hours. Pre-treatment with 5 mg/kg Albu-CocH also attenuated the reinstatement of extinguished cocaine self-administration by an i.v. priming injection of cocaine (0.1 or 0.3 mg/kg) and, in separate studies, attenuated the discriminative-stimulus effects of cocaine. The ability of Albu-CocH to attenuate the abuse-related effects of cocaine in squirrel monkeys indicates that further investigation of BChE mutants as potential treatment for cocaine abuse and toxicity is warranted.

Anti-cocaine antibody and butyrylcholinesterase-derived cocaine hydrolase exert cooperative effects on cocaine pharmacokinetics and cocaine-induced locomotor activity in mice

We are investigating treatments for cocaine abuse based on viral gene transfer of a cocaine hydrolase (CocH) derived from human butyrylcholinesterase, which can reduce cocaine-stimulated locomotion and cocaine-primed reinstatement of drug-seeking behavior in rats for many months. Here, in mice, we explored the possibility that anti-cocaine antibodies can complement the actions of CocH to reduce cocaine uptake in brain and block centrally-evoked locomotor stimulation. Direct injections of test proteins showed that CocH (0.3 or 1mg/kg) was effective by itself in reducing drug levels in plasma and brain of mice given cocaine (10mg/kg, s.c., or 20mg/kg, i.p). Administration of cocaine antibody per se at a low dose (8 mg/kg, i.p.) exerted little effect on cocaine distribution. However, a higher dose of antibody (12 mg/kg) caused peripheral trapping (increased plasma drug levels), which led to increased cocaine metabolism by CocH, as evidenced by a 6-fold rise in plasma benzoic acid. Behavioral tests with small doses of CocH and antibody (1 and 8 mg/kg, respectively) showed that neither agent alone reduced mouse locomotor activity triggered by a very large cocaine dose (100mg/kg, i.p.). However, dual treatment completely suppressed the locomotor stimulation. Altogether, we found cooperative and possibly synergistic actions that warrant further exploration of dual therapies for treatment of cocaine abuse.

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.

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.

Antibody-antigen pair probed by combinatorial approach and rational design: bringing together structural insights, directed evolution, and novel functionality

The unique hypervariability of the immunoglobulin (Ig) superfamily provides a means to create both binding and catalytic antibodies with almost any desired specificity and activity. The diversity of antigens and concept of adaptive response suggest that it is possible to find an antigen pair to any raised Ig. In the current review we discuss combinatorial approaches, which makes it possible to obtain an antibody with predefined properties, followed by 3D structure-based rational design to enhance or dramatically change its characteristics. A similar strategy, but applied to the second partner of the antibody-antigen pair, may result in selection of complementary substrates to the chosen Ig. Finally, 2D screening may be performed solving the "Chicken and Egg" problem when neither antibody nor antigen is known.

Repeated administration of a mutant cocaine esterase: effects on plasma cocaine levels, cocaine-induced cardiovascular activity, and immune responses in rhesus monkeys

Previous studies have demonstrated the capacity of a long-acting mutant form of a naturally occurring bacterial double mutant cocaine esterase (DM CocE) to antagonize the reinforcing, discriminative, convulsant, and lethal effects of cocaine in rodents and reverse the increases in mean arterial pressure (MAP) and heart rate (HR) produced by cocaine in rhesus monkeys. This study was aimed at characterizing the immunologic responses to repeated dosing with DM CocE and determining whether the development of anti-CocE antibodies altered the capacity of DM CocE to reduce plasma cocaine levels and ameliorate the cardiovascular effects of cocaine in rhesus monkeys. Under control conditions, intravenous administration of cocaine (3 mg/kg) resulted in a rapid increase in the plasma concentration of cocaine (n = 2) and long-lasting increases in MAP and HR (n = 3). Administration of DM CocE (0.32 mg/kg i.v.) 10 min after cocaine resulted in a rapid hydrolysis of cocaine with plasma levels below detection limits within 5 to 8 min. Elevations in MAP and HR were significantly reduced within 25 and 50 min of DM CocE administration, respectively. Although slight (10-fold) increases in anti-CocE antibodies were observed after the fourth administration of DM CocE, these antibodies did not alter the capacity of DM CocE to reduce plasma cocaine levels or ameliorate cocaine's cardiovascular effects. Anti-CocE titers were transient and generally dissipated within 8 weeks. Together, these results suggest that highly efficient cocaine esterases, such as DM CocE, may provide a novel and effective therapeutic for the treatment of acute cocaine intoxication in humans.

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.

Long-lasting effects of a PEGylated mutant cocaine esterase (CocE) on the reinforcing and discriminative stimulus effects of cocaine in rats

Recent mutagenesis studies have identified a mutant G4C/S10C/T172R/G173Q cocaine esterase (CCRQ CocE) with an in vitro duration of action of >40 days. Although the in vivo duration of CCRQ CocE's action was <24 h, modification of this enzyme with polyethylene glycol (PEG) polymers resulted in a CocE (PEG-CCRQ CocE) capable of preventing cocaine-induced lethality for up to 72 h. The current studies were aimed at providing a detailed characterization of the effectiveness, selectivity, and duration of PEG-CCRQ CocE's actions in cocaine self-administration and discrimination assays in rats. Pretreatment with PEG-CCRQ CocE produced dose-dependent rightward shifts in the dose-response curves for cocaine self-administration and discrimination, with the highest dose of PEG-CCRQ CocE capable of producing an initial shift of cocaine's reinforcing and interoceptive effects of >30-fold to the right, with significant inhibition of these effects observed for up to 72 h. Although PEG-CCRQ CocE also produced slight reductions in the rates of methylphenidate- and food-reinforced responding, these effects were short-lived, lasting <24 h. Finally, when taken together with the finding that PEG-CCRQ CocE failed to alter the cocaine-like interoceptive effects of either methylphenidate or d-amphetamine, these results suggest that PEG-CCRQ CocE possesses a high degree of pharmacologic specificity for cocaine and a prolonged in vivo duration of action. In conclusion, these studies provide strong evidence to support the further development of long-lasting, highly efficient CocEs, such as PEG-CCRQ CocE, as a potential therapeutic option for the treatment of cocaine abuse in humans.

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

Are pharmacokinetic approaches feasible for treatment of cocaine addiction and overdose?

“… we discuss the main challenges for developing therapeutic treatment of cocaine addiction and explain why pharmacokinetic approaches, particularly those based on our recently developed efficient cocaine-metabolizing enzymes, are feasible for treatment of cocaine addiction and overdose.”

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.

Evaluation of the hydrolytic activity of a long-acting mutant bacterial cocaine in the presence of commonly co-administered drugs

BACKGROUND

Cocaine toxicity is a prevalent problem in the Unites States for which there is currently no FDA-approved pharmacotherapy. We have developed a bacterial cocaine esterase (CocE) towards this indication. A thermostabilized mutant of CocE (DM-CocE) retains the hydrolytic activity of the wild-type esterase, rapidly hydrolyzing cocaine into the inactive metabolites ecgonine methyl ester and benzoic acid, and can prevent cocaine toxicities in rodent and non-human primate models. To advance DM-CocE towards clinical use, we examine here how the hydrolytic activity of DM-CocE is altered by some drugs commonly co-administered with cocaine.

METHODS

We employed a spectrophotometric cocaine hydrolysis assay to evaluate whether pharmacologically relevant doses of alcohol, nicotine, morphine, phencyclidine, ketamine, methamphetamine, naltrexone, naloxone, or midazolam would alter the Michaelis-Menten kinetics of DM-CocE for cocaine. Mass spectrometry was used to evaluate interaction with diazepam as this drug interferes with the absorbance spectra of cocaine critical for the spectrophotometric assay.

RESULTS

Alcohol, nicotine, morphine, phencyclidine, ketamine, methamphetamine, naltrexone, naloxone, and midazolam did not alter cocaine hydrolysis by DM-CocE. However, diazepam significantly slowed DM-CocE cocaine hydrolysis at very high concentrations, most likely through interaction of the phenyl ring of the benzodiazepine with the active site of DM-CocE.

CONCLUSIONS

DM-CocE does not display significant drug interactions, with the exception of diazepam, which may warrant further study as DM-CocE progresses towards a clinically used pharmacotherapy for cocaine toxicity. Alternate benzodiazepines, e.g., midazolam could be used to avoid this potential interaction.

Clair C, Nichols J, Narasimhan D, Ko MC, Sunahara RK, Woods JH. Effects of a long-acting mutant bacterial cocaine esterase on acute cocaine toxicity in rats

BACKGROUND

A longer acting, double mutant bacterial cocaine esterase (CocE T172R/G173Q; DM CocE) has been shown to protect mice from cocaine-induced lethality, inhibit the reinforcing effects of cocaine in rats, and reverse cocaine's cardiovascular effects in rhesus monkeys. The current studies evaluated the effectiveness of DM CocE to protect against, and reverse cocaine's cardiovascular, convulsant, and lethal effects in male and female rats.

METHODS

Pretreatment studies were used to determine the effectiveness and in vivo duration of action for DM CocE to protect rats against the occurrence of cardiovascular changes, convulsion and lethality associated with acute cocaine toxicity. Posttreatment studies were used to evaluate the capacity of DM CocE to rescue rats from the cardiovascular and lethal effects of large doses of cocaine. In addition, male and female rats were studied to determine if there were any potential effects of sex on the capacity of DM CocE to protect against, or reverse acute cocaine toxicity in rats.

RESULTS

Pretreatment with DM CocE dose-dependently protected rats against cocaine-induced cardiovascular changes, convulsion and lethality, with higher doses active for up to 4h, and shifting cocaine-induced lethality at least 10-fold to the right. In addition to dose-dependently recovering rats from an otherwise lethal dose of cocaine, post-treatment with DM CocE also reversed the cardiovascular effects of cocaine. There were no sex-related differences in the effectiveness of DM CocE to protect against, or reverse acute cocaine toxicity.

CONCLUSIONS

Together, these results support the development of DM CocE for the treatment of acute cocaine toxicity.

The ability of bacterial cocaine esterase to hydrolyze cocaine metabolites and their simultaneous quantification using high-performance liquid chromatography-tandem mass spectrometry

Cocaine toxicity is a widespread problem in the United States, responsible for more than 500,000 emergency department visits a year. There is currently no U.S. Food and Drug Administration-approved pharmacotherapy to directly treat cocaine toxicity. To this end, we have developed a mutant bacterial cocaine esterase (DM-CocE), which has been previously shown to rapidly hydrolyze cocaine into inert metabolites, preventing and reversing toxicity with limited immunogenic potential. Herein we describe the ability of DM-CocE to hydrolyze the active cocaine metabolites norcocaine and cocaethylene and its inability to hydrolyze benzoylecgonine. DM-CocE hydrolyzes norcocaine and cocaethylene with 58 and 45% of its catalytic efficiency for cocaine in vitro as measured by a spectrophotometric assay. We have developed a mass spectrometry method to simultaneously detect cocaine, benzoylecgonine, norcocaine, and ecgonine methyl ester to quantify the effect of DM-CocE on normal cocaine metabolism in vivo. DM-CocE administered to rats 10 min after a convulsant dose of cocaine alters the normal metabolism of cocaine, rapidly decreasing circulating levels of cocaine and norcocaine while increasing ecgonine methyl ester formation. Benzoylecgonine was not hydrolyzed in vivo, but circulating concentrations were reduced, suggesting that DM-CocE may bind and sequester this metabolite. These findings suggest that DM-CocE may reduce cocaine toxicity by eliminating active and toxic metabolites along with the parent cocaine molecule.

Active site gating and substrate specificity of butyrylcholinesterase and acetylcholinesterase: insights from molecular dynamics simulations

Butyrylcholinesterase (BChE) and acetylcholinesterase (AChE) are highly homologous proteins with distinct substrate preferences. In this study we compared the active sites of monomers and tetramers of human BChE and human AChE after performing molecular dynamics (MD) simulations in water-solvated systems. By comparing the conformational dynamics of gating residues of AChE and BChE, we found that the gating mechanisms of the main door of AChE and BChE are responsible for their different substrate specificities. Our simulation of the tetramers of AChE and BChE indicates that both enzymes could have two dysfunctional active sites due to their restricted accessibility to substrates. The further study on catalytic mechanisms of multiple forms of AChE and BChE would benefit from our comparison of the active sites of the monomers and tetramers of both enzymes.

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.

Computational design of a thermostable mutant of cocaine esterase via molecular dynamics simulations

Cocaine esterase (CocE) has been known as the most efficient native enzyme for metabolizing naturally occurring cocaine. A major obstacle to the clinical application of CocE is the thermoinstability of native CocE with a half-life of only ∼11 min at physiological temperature (37 °C). It is highly desirable to develop a thermostable mutant of CocE for therapeutic treatment of cocaine overdose and addiction. To establish a structure-thermostability relationship, we carried out molecular dynamics (MD) simulations at 400 K on wild-type CocE and previously known thermostable mutants, demonstrating that the thermostability of the active form of the enzyme correlates with the fluctuation (characterized as the root-mean square deviation and root-mean square fluctuation of atomic positions) of the catalytic residues (Y44, S117, Y118, H287, and D259) in the simulated enzyme. In light of the structure-thermostability correlation, further computational modelling including MD simulations at 400 K predicted that the active site structure of the L169K mutant should be more thermostable. The prediction has been confirmed by wet experimental tests showing that the active form of the L169K mutant had a half-life of 570 min at 37 °C, which is significantly longer than those of the wild-type and previously known thermostable mutants. The encouraging outcome suggests that the high-temperature MD simulations and the structure-thermostability relationship may be considered as a valuable tool for the computational design of thermostable mutants of an enzyme.

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.

Amelioration of the cardiovascular effects of cocaine in rhesus monkeys by a long-acting mutant form of cocaine esterase

A long-acting mutant form of a naturally occurring bacterial cocaine esterase (T172R/G173Q CocE; double mutant CocE (DM CocE)) has previously been shown to antagonize the reinforcing, convulsant, and lethal effects of cocaine in rodents. However, the effectiveness and therapeutic characteristics of DM CocE in nonhuman primates, in a more clinically relevant context, are unknown. The current studies were aimed at (1) characterizing the cardiovascular effects of cocaine in freely moving rhesus monkeys, (2) evaluating the capacity of DM CocE to ameliorate these cocaine-induced cardiovascular effects when administered 10 min after cocaine, and (3) assessing the immunological responses of monkeys to DM CocE following repeated administration. Intravenous administration of cocaine produced dose-dependent increases in mean arterial pressure (MAP) and heart rate (HR) that persisted throughout the 2-h observation period following a dose of 3.2 mg/kg cocaine. Cocaine failed to produce reliable changes in electrocardiograph (ECG) parameters, body temperature, and locomotor activity. DM CocE produced a rapid and dose-dependent amelioration of the cardiovascular effects, with saline-like MAP measures restored within 5-10 min, and saline-like HR measures restored within 20-40 min of DM CocE administration. Although administration of DM CocE produced increases in anti-CocE antibodies, they did not appear to have a neutralizing effect on the capacity of DM CocE to reverse the cardiovascular effects of cocaine. In conclusion, these findings in monkeys provide strong evidence to suggest that highly efficient cocaine esterases, such as DM CocE, can provide a potential therapeutic option for treatment of acute cocaine intoxication in humans.

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.

Effects of cocaine hydrolase on cocaine self-administration under a PR schedule and during extended access (escalation) in rats

RATIONALE AND OBJECTIVES

Previously, Albu-CocH, a cocaine hydrolase derived from human butyrylcholinesterase, blocked cocaine-induced reinstatement of drug seeking in rats. In the present study, rats were treated with Albu-CocH while self-administering cocaine under a progressive ratio (PR) schedule during 2-h sessions and under a fixed-ratio 1 (FR 1) schedule during 6-h sessions.

METHODS

In experiment 1, rats were treated with saline or Albu-CocH (2 or 4 mg/kg) before a single 2-h cocaine (0.2 mg/kg) self-administration (PR) session. In experiment 2, rats were treated with Albu-CocH or saline for the first seven of the 21-day 6-h sessions prior to cocaine (0.2 or 0.4 mg/kg) self-administration sessions (FR 1).

RESULTS

In experiment 1, Albu-CocH (vs saline) reduced cocaine infusions immediately following treatment compared with sessions pretreatment and posttreatment. In experiment 2, the Albu-CocH-treated groups (vs saline) showed an initial twofold to threefold increase in 0.2 and 0.4 mg/kg cocaine infusions over the 7 days of treatment, but they decreased to the infusion levels of saline controls by day 7. Cocaine (0.4 mg/kg) intake in the saline-treated group was elevated during the last 3 days of 6-h access compared with the first 3 days, indicating an escalation effect. Responding for 0.4 mg/kg (but not 0.2 mg/kg) cocaine during 2-h sessions after the 21 days of 6-h access was elevated in the saline groups (compared with 2-h sessions before long access) but not in the Albu-CocH-treated groups.

CONCLUSIONS

Albu-CocH decreased cocaine infusions under the PR schedule, indicating a reduced reward value of cocaine (experiment 1). However, Albu-CocH, compared with saline, temporarily increased cocaine infusions during long access. The post-long access 2-h cocaine intake was not increased in the Albu-CocH-treated groups as it was in the saline-treated groups. Albu-CocH is an effective agent for reducing cocaine reward under conditions of low cocaine exposure and chronic treatment.