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

A quantitative LC-MS/MS method for investigation of polysubstance use involving heroin and cocaine

Concurrent use of heroin and cocaine (known as the "speedball") prevails among substance use disorder populations, especially in opioid-dependent individuals, with severe consequences and a high fatality rate. Little is known about the patterns and correlations of the concurrent use of heroin and cocaine. It is vital to investigate such a polydrug use in both humans and animals to uncover concomitant toxicity and the cause of fatal overdose (death). In this study, we aimed to shed some light on the role of cocaine in the etiology of heroin-related deaths in the context of molecular pharmacokinetics (PK). For the purpose, a high-performance liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) method for simultaneous determination of heroin, cocaine, and their metabolites in whole blood was developed and fully validated in accordance with the US Food and Drug Administration (FDA) and European Medicines Agency (EMA) guidelines. Then, this method was used to analyze heroin, cocaine, and their metabolites in blood from the rats intraperitoneally administered non-lethal 10 mg/kg heroin or 20 mg/kg cocaine alone, or their combination that is lethal with a proximal mortality of 33 %. The obtained results from the rats that experienced the lethal toxicity revealed that the concurrent use of heroin and cocaine significantly increased the risk of fatality from overdose. Heroin significantly slowed down the elimination of cocaine and its main metabolites in blood, while cocaine significantly enhanced heroin metabolism from 6-monoacetylmorphine (6-MAM) to morphine. Similar elimination half-lives for other heroin metabolites were observed. These findings are reported for the first time in this study, facilitating our understanding of the polysubstance metabolism and severe consequences produced by the polydrug use.

Discovery of a novel homocysteine thiolactone hydrolase and the catalytic activity of its natural variants

Homocysteine thiolactone (HTL), a toxic metabolite of homocysteine (Hcy) in hyperhomocysteinemia (HHcy), is known to modify protein structure and function, leading to protein damage through formation of N-Hcy-protein. HTL has been highly linked to HHcy-associated cardiovascular and neurodegenerative diseases. The protective role of HTL hydrolases against HTL-associated vascular toxicity and neurotoxicity have been reported. Although several endogeneous enzymes capable of hydrolyzing HTL have been identified, the primary enzyme responsible for its metabolism remains unclear. In this study, three human carboxylesterases were screened to explore new HTL hydrolase and human carboxylesterase 1 (hCES1) demonstrates the highest catalytic activity against HTL. Given the abundance of hCES1 in the liver and the clinical significance of its single-nucleotide polymorphisms (SNPs), six common hCES1 nonsynonymous coding SNP (nsSNPs) variants were examined and characterized for their kinetic parameters. Variants E220G and G143E displayed 7.3-fold and 13.2-fold lower catalytic activities than its wild-type counterpart. In addition, the detailed catalytic mechanism of hCES1 for HTL hydrolysis was computational investigated and elucidated by Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) method. The function of residues E220 and G143 in sustaining its hydrolytic activity of hCES1 was analyzed, and the calculated energy difference aligns well with experimental-derived results, supporting the validity of our computational insights. These findings provide insights into the potential protective role of hCES1 against HTL-associated toxicity, and warrant future studies on the possible association between specific genetic variants of hCES1 with impaired catalytic function and clinical susceptibility of HTL-associated cardiovascular and neurodegenerative diseases.

Role of carboxylesterase and arylacetamide deacetylase in drug metabolism, physiology, and pathology

Carboxylesterases (CES1 and CES2) and arylacetamide deacetylase (AADAC), which are expressed primarily in the liver and/or gastrointestinal tract, hydrolyze drugs containing ester and amide bonds in their chemical structure. These enzymes often catalyze the conversion of prodrugs, including the COVID-19 drugs remdesivir and molnupiravir, to their pharmacologically active forms. Information on the substrate specificity and inhibitory properties of these enzymes, which would be useful for drug development and toxicity avoidance, has accumulated. Recently,in vitroandin vivostudies have shown that these enzymes are involved not only in drug hydrolysis but also in lipid metabolism. CES1 and CES2 are capable of hydrolyzing triacylglycerol, and the deletion of their orthologous genes in mice has been associated with impaired lipid metabolism and hepatic steatosis. Adeno-associated virus-mediated human CES overexpression decreases hepatic triacylglycerol levels and increases fatty acid oxidation in mice. It has also been shown that overexpression of CES enzymes or AADAC in cultured cells suppresses the intracellular accumulation of triacylglycerol. Recent reports indicate that AADAC can be up- or downregulated in tumors of various organs, and its varied expression is associated with poor prognosis in patients with cancer. Thus, CES and AADAC not only determine drug efficacy and toxicity but are also involved in pathophysiology. This review summarizes recent findings on the roles of CES and AADAC in drug metabolism, physiology, and pathology.

Design of prodrug for stereoisomers of omapatrilat to cross the blood-brain barrier using docking, homology modeling, MD, and QM/MM methods

In this study, we designed a suitable ester prodrug for omapatrilat to penetrate the blood-brain barrier and treat CNS diseases. Based on the ADMET properties, the methyl carboxylate ester of omapatrilat was chosen from among several prodrug structures. Sixteen methyl carboxylate esters were constructed for omapatrilat. The structure of brain carboxylesterase was derived via homology modeling, and molecular docking was used to determine the most potent stereoisomers against brain carboxylesterase. The top three stereoisomer complexes, and the apo form of the protein, were then considered using molecular dynamics simulation and MM/GBSA analysis. Following the simulation, structural analysis was performed using RMSD, RMSF, Rg, and hydrogen bond analysis tools. Our data demonstrated that the prodrug of RSSR is a suitable structure for crossing the blood-brain barrier and binding to brain carboxylesterase. In addition, we found via QM/MM calculation that the catalytic reaction of the prodrug of RSSR against brain carboxylesterase occurs via two steps, including acylation and diacylation steps. Based on our findings, we propose a clinical trial of a methyl carboxylate ester prodrug of omapatrilat's RSSR for the treatment of brain diseases.Communicated by Ramaswamy H. Sarma.

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.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

Going deeper into the toxicokinetics of synthetic cannabinoids: in vitro contribution of human carboxylesterases

Synthetic cannabinoids (SC) are new psychoactive substances known to cause intoxications and fatalities. One reason may be the limited data available concerning the toxicokinetics of SC, but toxicity mechanisms are insufficiently understood so far. Human carboxylesterases (hCES) are widely known to play a crucial role in the catalytic hydrolysis of drugs (of abuse). The aim of this study was to investigate the in vitro contribution of hCES to the metabolism of the 13 SC 3,5-AB-5F-FUPPYCA, AB-5F-P7AICA, A-CHMINACA, DMBA-CHMINACA, MBA-CHMINACA, MDMB-4F-BINACA, MDMB-4en-PINACA, MDMB-FUBICA, MDMB-5F-PICA, MMB-CHMICA, MMB-4en-PICA, MMB-FUBINACA, and MPhP-5F-PICA. The SC were incubated with recombinant hCES1b, hCES1c, or hCES2 and analyzed by liquid chromatography-ion trap mass spectrometry to assess amide or ester hydrolysis in an initial activity screening. Enzyme kinetic studies were performed if sufficient hydrolysis was observed. No hydrolysis of the amide linker was observed using those experimental conditions. Except for MDMB-5F-PICA, ester hydrolysis was always detected if an ester group was present in the head group. In general, SC with a terminal ester bearing a small alcohol part and a larger acyl part showed higher affinity to hCES1 isozymes. Due to the low hydrolysis rates, enzyme kinetics could not be modeled for the SC with a tert-leucine-derived moiety, but hydrolysis reactions of MPhP-5F-PICA and of those containing a valine-derived moiety followed classic Michaelis-Menten kinetics. In conclusion, drug-drug/drug-food interactions or hCES polymorphisms may prolong the half-life of SC and the current results help to estimate the risk of toxicity in the future after combining them with activity and clinical data.

Inhibition of Radix Scutellariae flavones on carboxylesterase mediated activations of prodrugs

AIMS

Carboxylesterase (CES) plays an essential role in the hydrolysis of ester prodrugs. Our study explored the inhibitions of Radix Scutellariae flavones, including baicalein (B), baicalin (BG), wogonin (W), wogonoside (WG), oroxylin A (OXA) and oroxylin A-7-O-glucuronide (OAG), on CES-mediated hydrolysis of seven prodrugs (capecitabine, clopidogrel, mycophenolate mofetil, dabigatran etexilate, acetylsalicylic acid, prasugrel and irinotecan).

MAIN METHODS

In vitro screenings were developed by incubating the flavones with prodrugs in rat plasma, intestine S9 and liver S9. Docking simulations were conducted using AMDock v1.5.2. In vivo evaluations were performed in rats co-administered with the selected flavone and prodrug via oral gavage/intravenous administration for five consecutive days.

KEY FINDINGS

The in vitro investigation showed that B and OXA demonstrated strongest inhibitions on the hydrolysis of irinotecan followed by dabigatran in rat plasma, intestine S9 and liver S9. Consistent results showed in the molecular docking analyses. Additionally, in rats receiving irinotecan, B/OXA intravenous and oral pre-treatments both led to reduction trends on the active metabolite SN-38 formation in plasma. Besides, significant decreases of SN-38/irinotecan plasma concentration ratios were found in the B/OXA oral pre-treatment group with quicker and stronger inhibition potential in OXA pre-treatment than that from B pre-treatment. OXA oral pre-treatment was also found to be able to significantly inhibit intestinal CES2 activities at 0.5 h and 5 h after irinotecan administration.

SIGNIFICANCE

Our current findings for the first time alert on potential CES-mediated HDIs between RS flavones and prodrugs, which provide a constructive information referring to rational drug combinations in clinical practice.

Rational Redesign of Enzyme via the Combination of Quantum Mechanics/Molecular Mechanics, Molecular Dynamics, and Structural Biology Study

Increasing demands for efficient and versatile chemical reactions have prompted innovations in enzyme engineering. A major challenge in engineering α-ketoglutarate-dependent oxygenases is to develop a rational strategy which can be widely used for directly evolving the desired mutant to generate new products. Herein, we report a strategy for rational redesign of a model enzyme, 4-hydroxyphenylpyruvate dioxygenase (HPPD), based on quantum mechanics/molecular mechanics (QM/MM) calculation and molecular dynamic simulations. This strategy enriched our understanding of the HPPD catalytic reaction pathway and led to the discovery of a series of HPPD mutants producing hydroxyphenylacetate (HPA) as the alternative product other than the native product homogentisate. The predicted HPPD-Fe(IV)═O-HPA intermediate was further confirmed by the crystal structure of Arabidopsis thaliana HPPD/S267W complexed with HPA. These findings not only provide a good understanding of the structure-function relationship of HPPD but also demonstrate a generally applicable platform for the development of biocatalysts.

Computational Design and Crystal Structure of a Highly Efficient Benzoylecgonine Hydrolase

Benzoylecgonine (BZE) is the major toxic metabolite of cocaine and is responsible for the long-term cocaine-induced toxicity owing to its long residence time in humans. BZE is also the main contaminant following cocaine consumption. Here, we identified the bacterial cocaine esterase (CocE) as a BZE-metabolizing enzyme (BZEase), which can degrade BZE into biological inactive metabolites (ecgonine and benzoic acid). CocE was redesigned by a reactant-state-based enzyme design theory. An encouraging mutant denoted as BZEase2, presented a >400-fold improved catalytic efficiency against BZE compared with wild-type (WT) CocE. In vivo, a single dose of BZEase2 (1 mg kg^-1 , IV) could eliminate nearly all BZE within only two minutes, suggesting the enzyme has the potential for cocaine overdose treatment and BZE elimination in the environment by accelerating BZE clearance. The crystal structure of a designed BZEase was also determined.

One- and Two-Proton Transfer Mechanisms Coexist in One Active Site

Acibenzolar-S-methyl (ASM) is one of the most successfully commercialized plant activators of the systemic acquired resistance (SAR). However, its activation (hydrolysis) mechanism catalyzed by the salicylic acid binding protein 2 (SABP2) remains elusive. The fundamental catalytic mechanism of the SABP2-catalyzed hydrolysis of the ASM had been investigated by extensive computational and experimental studies, including QM/MM simulations, charge transfer analysis, small-molecule synthesis, and biochemical assays. Here we report that the promiscuous SABP2 shows different catalytic mechanisms toward different substrates. To catalyze the ASM hydrolysis, the SABP2 uses a two-proton transfer mechanism, and the key intermediate is stabilized by the charge transfer effect; to catalyze the ethyl 1,2,3-benzothiadiazole-7-carboxylate (BTM, an ASM analogue) hydrolysis, the SABP2 applies the one-proton transfer mechanism, and the classic tetrahedral intermediate is stabilized by the electrostatic effect. The HPLC analyses of the SABP2 esterase activities toward the ASM and the BTM show comparable results with our computaional results, suggesting that the obtained computational mechanism insights are reasonable. The obtained mechanism is not only an important supplement to the theory of enzymes' catalytic promiscuity, but it also contributes a possible strategy for the design of next generation plant SAR activators.

Evaluation the synergistic antitumor effect of methotrexate-camptothecin codelivery prodrug from self-assembly process to acid-catalyzed both drugs release: A comprehensive theoretical study

Therapeutic efficiency of amphiphilic methotrexate-camptothecin (MTX-CPT) prodrug compared to free drug mixture (MTX/CPT) has been investigated using all-atom molecular dynamics simulation and first principles density functional theory calculations. This comparison revealed that MTX-CPT prodrug tends to form spherical self-assembled nanoparticle (NP), while free MTX/CPT mixture forms rod-shape NP. These observations are attributed to a structural defect in the MTX-CPT prodrug and solvation free energies of MTX, CPT and MTX-CPT molecules. The results provided evidence that noncovalent interactions (NCIs) among the pharmaceutical drugs play a very important role in anticancer agents aggregation process, leading to enhanced stability of the self-assembled NPs. It is found that the stability of MTX-CPT self-assembled NP is greater than the MTX/CPT NP due to the synergistic effect of hydrogen bonding between monomers and solvent (water). Moreover, the noncatalyzed as well as catalyzed hydrolysis reactions of MTX-CPT prodrug are theoretically studied at the PCM(water)//M06-2X/6-31G(d,p) computational level to shed additional light on the role of acidic condition in tumor tissues. We found that the ester hydrolysis in mild acidic solutions is a concerted reaction. In an agreement between theory and experiment, we also confirmed that the activation energies of the catalyzed-hydrolysis steps are much lower than the activation energies of the corresponding steps in the noncatalyzed reaction. Thus, the MTX-CPT prodrug reveals very promising properties as a pH-controlled drug delivery system.

Effects of Cebranopadol on Cocaine-induced Hyperactivity and Cocaine Pharmacokinetics in Rats

Cebranopadol is known as a highly potent analgesic. Recent studies also demonstrated that administration of cebranopadol significantly decreased cocaine self-administration and significantly reduced cue-induced cocaine-seeking behaviors in rats. However, it was unclear whether these interesting behavioral observations are related to any potential effects of cebranopadol on cocaine pharmacokinetics or cocaine-induced hyperactivity. In principle, a promising therapeutic candidate for cocaine dependence treatment may alter the cocaine pharmacokinetics and/or attenuate cocaine-induced reward and hyperactivity and, thus, decrease cocaine self-administration and reduce cue-induced cocaine-seeking behaviors. In this study, we examined possible effects of cebranopadol on cocaine pharmacokinetics and cocaine-induced hyperactivity for the first time. According to our animal data in rats, cebranopadol did not significantly alter the pharmacokinetics of cocaine. According to our more extensive locomotor activity testing data, cebranopadol itself also dose-dependently induced hyperactivity in rats at doses higher than 50 µg/kg. Cebranopadol at a low dose of 25 µg/kg (p.o.) did not induce significant hyperactivity itself, but significantly potentiated cocaine-induced hyperactivity on Days 4 to 7 after the repeated daily dosing of the drug.

Cryptophycin-55/52 based antibody-drug conjugates: Synthesis, efficacy, and mode of action studies

Cryptophycin-52 (CR52), a tubulin inhibitor, exhibits promising antitumor activity in vitro (picomolar level) and in mouse xenograft models. However, the narrow therapeutic window in clinical trials limits its further development. Antibody-drug conjugate (ADC), formed by coupling cytotoxic compound (payload) to an antibody via a linker, can deliver drug to tumor locations in a targeted manner by antibody, enhancing the therapeutic effects and reducing toxic and side effects. In this study, we aim to explore the possibility of CR52-based ADC for tumor targeted therapy. Due to the lack of a coupling site in CR52, its prodrug cryptophycin-55 (CR55) containing a free hydroxyl was synthesized and conjugated to the model antibody trastuzumab (anti-HER2 antibody drug approved by FDA for breast cancer therapy) via the linkers based on Mc-NHS and Mc-Val-Cit-PAB-PNP. The average drug-to-antibody ratios (DARs) of trastuzumab-CR55 conjugates (named T-L1-CR55, T-L2-CR55, and T-L3-CR55) were 3.50, 3.29, and 3.35, respectively. These conjugates exhibited potent cytotoxicity in HER2-positive tumor cell lines with IC₅₀ values at low nanomolar levels (0.58-1.19 nM). Further, they displayed significant antitumor activities at the doses of 10 mg/kg in established ovarian cancer (SKOV3) and gastric cancer (NCI-N87) xenograft models without overt toxicities. Finally, the drug releases were analyzed and the results indicated that T-L3-CR55 was able to effectively release CR55 and further epoxidized to CR52, which may be responsible for its best performance in antitumor activities. In conclusion, our results demonstrated that these conjugates have the potential for tumor targeted therapy, which provides insights to further research the CR55/CR52-based ADC for tumor therapy.

Catalytic Hydrolysis Mechanism of Cocaine by Human Carboxylesterase 1: An Orthoester Intermediate Slows Down the Reaction

Human carboxylesterase 1 (hCES1) is a major carboxylesterase in the human body and plays important roles in the metabolism of a wide variety of substances, including lipids and drugs, and therefore is attracting more and more attention from areas including lipid metabolism, pharmacokinetics, drug–drug interactions, and prodrug activation. In this work, we studied the catalytic hydrolysis mechanism of hCES1 by the quantum mechanics computation method, using cocaine as a model substrate. Our results support the four-step theory of the esterase catalytic hydrolysis mechanism, in which both the acylation stage and the deacylation stage include two transition states and a tetrahedral intermediate. The roles and cooperation of the catalytic triad, S221, H468, and E354, were also analyzed in this study. Moreover, orthoester intermediates were found in hCES1-catalyzed cocaine hydrolysis reaction, which significantly elevate the free energy barrier and slow down the reaction. Based on this finding, we propose that hCES1 substrates with β-aminocarboxylester structure might form orthoester intermediates in hCES1-catalyzed hydrolysis, and therefore prolong their in vivo half-life. Thus, this study helps to clarify the catalytic mechanism of hCES1 and elucidates important details of its catalytic process, and furthermore, provides important insights into the metabolism of hCES1 substrates and drug designing.

Human Acetyl-CoA Carboxylase 1 Is an Isomerase: Carboxyl Transfer Is Activated by Catalytic Effect of Isomerization

Obesity and its related diseases such as cancer and diabetes are leading life-threatening issues in the modern world. Thus, new drugs toward obesity and obesity-caused diseases are highly desired. Human acetyl-CoA carboxylase 1 (hACC1) in charge of the rate-limiting step of the human fatty acid synthesis was recognized as an attractive target for rational drug design. The fundamental reaction mechanism and nature of the transition state of hACC1 remain unclear. In this study, combined quantum mechanics and molecular mechanics (QM/MM), molecular dynamics (MD), and free-energy simulations were performed to investigate the catalytic mechanism of the hACC1-catalyzed carboxyl-transfer reaction. Our computational results show a three-step mechanism for carboxyl transferase (CT)-catalyzed reaction, including isomerization of carboxybiotin, proton-transfer from acetyl-CoA to carboxybiotin, and carboxylation of acetyl-CoA enolate. Interestingly, isomerization of carboxybiotin is the rate-limiting step of the entire reaction pathway, indicating hACC1 has the catalytic effect of isomerization and thus might be an isomerase also. The activation free-energy barrier of carboxyl-transfer catalyzed by hACC1 was calculated to be 16.4 kcal/mol, in excellent agreement with the experimental result (16.7 kcal/mol). The obtained reaction mechanism together with the nature of the transition state provides helpful knowledge not only for future investigation of other ACCs but also for rational design of hACC1 inhibitors, such as TS analogue. The catalytic effect of hACC1 isomerization is discussed.

The catalytic mechanism of S-acyltransferases: acylation is triggered on by a loose transition state and deacylation is turned off by a tight transition state

The reaction of S-acyltransferase is characterized by a loose transition state.