Free energy perturbation simulation on transition states and high-activity mutants of human butyrylcholinesterase for (-)-cocaine hydrolysis

In vitro and in vivo stability of a highly efficient long-acting cocaine hydrolase

It is recognized as a promising therapeutic strategy for cocaine use disorder to develop an efficient enzyme which can rapidly convert cocaine to physiologically inactive metabolites. We have designed and discovered a series of highly efficient cocaine hydrolases, including CocH5-Fc(M6) which is the currently known as the most efficient cocaine hydrolase with both the highest catalytic activity against (-)-cocaine and the longest biological half-life in rats. In the present study, we characterized the time courses of protein appearance, pH, structural integrity, and catalytic activity against cocaine in vitro and in vivo of a CocH5-Fc(M6) bulk drug substance produced in a bioreactor for its in vitro and in vivo stability after long-time storage under various temperatures (- 80, - 20, 4, 25, or 37 °C). Specifically, all the tested properties of the CocH5-Fc(M6) protein did not significantly change after the protein was stored at any of four temperatures including - 80, - 20, 4, and 25 °C for ~ 18 months. In comparison, at 37 °C, the protein was less stable, with a half-life of ~ 82 days for cocaine hydrolysis activity. Additionally, the in vivo studies further confirmed the linear elimination PK profile of CocH5-Fc(M6) with an elimination half-life of ~ 9 days. All the in vitro and in vivo data on the efficacy and stability of CocH5-Fc(M6) have consistently demonstrated that CocH5-Fc(M6) has the desired in vitro and in vivo stability as a promising therapeutic candidate for treatment of cocaine use disorder.

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.

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

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

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.

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.

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

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

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.

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

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

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

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

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.