1
|
Wei H, LeSaint JE, Jin Z, Zhan CG, Zheng F. Long-lasting blocking of interoceptive effects of cocaine by a highly efficient cocaine hydrolase in rats. Sci Rep 2024; 14:927. [PMID: 38195724 PMCID: PMC10776848 DOI: 10.1038/s41598-023-50678-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Accepted: 12/22/2023] [Indexed: 01/11/2024] Open
Abstract
Cocaine dependence is a serious world-wide public health problem without an FDA-approved pharmacotherapy. We recently designed and discovered a highly efficient long-acting cocaine hydrolase CocH5-Fc(M6). The present study examined the effectiveness and duration of CocH5-Fc(M6) in blocking interoceptive effects of cocaine by performing cocaine discrimination tests in rats, demonstrating that the duration of CocH5-Fc(M6) in blocking cocaine discrimination was dependent on cocaine dose and CocH5-Fc(M6) plasma concentration. Particularly, a dose of 3 mg/kg CocH5-Fc(M6) effectively attenuated discriminative stimulus effects of 10 mg/kg cocaine, cumulative doses of 10 and 32 mg/kg cocaine, and cumulative doses of 10, 32 and 56 mg/kg cocaine by ≥ 20% for 41, 19, and 10 days, and completely blocked the discriminative stimulus effects for 30, 13, and 5 days with corresponding threshold plasma CocH5-Fc(M6) concentrations of 15.9, 72.2, and 221 nM, respectively, under which blood cocaine concentration was negligible. Additionally, based on the data obtained, cocaine discrimination model is more sensitive than the locomotor activity to reveal cocaine effects and that CocH5-Fc(M6) itself has no long-term toxicity regarding behavioral activities such as lever pressing and food consumption in rats, further demonstrating that CocH5-Fc(M6) has the desired properties as a promising therapeutic candidate for prevenance of cocaine dependence.
Collapse
Affiliation(s)
- Huimei Wei
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
| | - Johnathan E LeSaint
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
| | - Zhenyu Jin
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
| |
Collapse
|
2
|
Zheng F, Hou S, Xue L, Yang W, Zhan CG. Human Butyrylcholinesterase Mutants for (-)-Cocaine Hydrolysis: A Correlation Relationship between Catalytic Efficiency and Total Hydrogen Bonding Energy with an Oxyanion Hole. J Phys Chem B 2023; 127:10723-10729. [PMID: 38063500 DOI: 10.1021/acs.jpcb.3c06392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
A combined computational and experimental study has been carried out to explore and test a quantitative correlation relationship between the relative catalytic efficiency (RCE) of human butyrylcholinesrase (BChE) mutant-catalyzed hydrolysis of substrate (-)-cocaine and the total hydrogen bonding energy (tHBE) of the carbonyl oxygen of the substrate with the oxyanion hole of the enzyme in the modeled transition-state structure (TS1), demonstrating a satisfactory linear correlation relationship between ln(RCE) and tHBE. The satisfactory correlation relationship has led us to computationally predict and experimentally confirm new human BChE mutants that have a further improved catalytic activity against (-)-cocaine, including the most active one (the A199S/F227S/S287G/A328W/Y332G mutant) with a 2790-fold improved catalytic efficiency (kcat/KM = 2.5 × 109 min-1 M-1) compared to the wild-type human BChE. Compared to the reference mutant (the A199S/S287G/A328W/Y332G mutant) tested in the reported clinical development of an enzyme therapy for cocaine dependence treatment, this new mutant (with a newly predicted additional F227S mutation) has an improved catalytic efficiency against (-)-cocaine by ∼2.6-fold. The good agreement between the computational and experimental ln(RCE) values suggests that the obtained correlation relationship is robust for computational prediction. A similar correlation relationship could also be explored in studying BChE or other serine hydrolases/esterases with an oxyanion hole stabilizing the carbonyl oxygen in the rate-determining reaction step of the enzymatic hydrolysis of other substrates.
Collapse
Affiliation(s)
- Fang Zheng
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Shurong Hou
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Liu Xue
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Wenchao Yang
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| |
Collapse
|
3
|
Catalytic activities of a highly efficient cocaine hydrolase for hydrolysis of biologically active cocaine metabolites norcocaine and benzoylecgonine. Sci Rep 2023; 13:640. [PMID: 36635293 PMCID: PMC9837138 DOI: 10.1038/s41598-022-27280-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 12/29/2022] [Indexed: 01/13/2023] Open
Abstract
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 (kcat = 9,210 min-1, KM = 20.9 μM, and kcat/KM = 1.87 × 105 min-1 M-1) and benzoylecgonine (kcat = 158 min-1, KM = 286 μM, and kcat/KM = 5.5 × 105 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.
Collapse
|
4
|
Galvis CEP, Kouznetsov VV. Cocaine: from a controlled drug by the DEA to an approved drug by the FDA. CURR ORG CHEM 2022. [DOI: 10.2174/1385272826666220221102840] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Abstract:
NumbrinoTM, a nasal anesthetic solution based on cocaine hydrochloride, was recently approved by the FDA for diagnostic procedures and surgery – now, what will be the role of the controlling agencies and the legislation of producing countries to face the imminent medicinal use of cocaine? Would chemistry be the protagonist in proposing alternative approaches to help and overcome this issue? The debate is open with this commentary
Collapse
Affiliation(s)
- Carlos E. Puerto Galvis
- Laboratorio de Química Orgánica y Biomolecular, CMN, Universidad Industrial de Santander. Parque Tecnológico Guati-guará, Km 2 Vía Refugio Piedecuesta 681011, Colombia
| | - Vladimir V. Kouznetsov
- Laboratorio de Química Orgánica y Biomolecular, CMN, Universidad Industrial de Santander. Parque Tecnológico Guati-guará, Km 2 Vía Refugio Piedecuesta 681011, Colombia
| |
Collapse
|
5
|
Arabnejad H, Bombino E, Colpa DI, Jekel PA, Trajkovic M, Wijma HJ, Janssen DB. Computational Design of Enantiocomplementary Epoxide Hydrolases for Asymmetric Synthesis of Aliphatic and Aromatic Diols. Chembiochem 2020; 21:1893-1904. [PMID: 31961471 PMCID: PMC7383614 DOI: 10.1002/cbic.201900726] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2019] [Revised: 01/16/2020] [Indexed: 12/13/2022]
Abstract
The use of enzymes in preparative biocatalysis often requires tailoring enzyme selectivity by protein engineering. Herein we explore the use of computational library design and molecular dynamics simulations to create variants of limonene epoxide hydrolase that produce enantiomeric diols from meso‐epoxides. Three substrates of different sizes were targeted: cis‐2,3‐butene oxide, cyclopentene oxide, and cis‐stilbene oxide. Most of the 28 designs tested were active and showed the predicted enantioselectivity. Excellent enantioselectivities were obtained for the bulky substrate cis‐stilbene oxide, and enantiocomplementary mutants produced (S,S)‐ and (R,R)‐stilbene diol with >97 % enantiomeric excess. An (R,R)‐selective mutant was used to prepare (R,R)‐stilbene diol with high enantiopurity (98 % conversion into diol, >99 % ee). Some variants displayed higher catalytic rates (kcat) than the original enzyme, but in most cases KM values increased as well. The results demonstrate the feasibility of computational design and screening to engineer enantioselective epoxide hydrolase variants with very limited laboratory screening.
Collapse
Affiliation(s)
- Hesam Arabnejad
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Elvira Bombino
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Dana I. Colpa
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Peter A. Jekel
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Milos Trajkovic
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Hein J. Wijma
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| | - Dick B. Janssen
- Biotransformation and Biocatalysis, Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 49747 AGGroningenThe Netherlands
| |
Collapse
|
6
|
Cai Y, Zhou S, Jin Z, Wei H, Shang L, Deng J, Zhan CG, Zheng F. Reengineering of Albumin-Fused Cocaine Hydrolase CocH1 (TV-1380) to Prolong Its Biological Half-Life. AAPS JOURNAL 2019; 22:5. [PMID: 31754920 DOI: 10.1208/s12248-019-0377-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2019] [Accepted: 08/29/2019] [Indexed: 11/30/2022]
Abstract
Therapeutic treatment of cocaine toxicity or addiction is a grand medical challenge. As a promising therapeutic strategy for treatment of cocaine toxicity and addiction to develop a highly efficient cocaine hydrolase (CocH) capable of accelerating cocaine metabolism to produce physiologically/biologically inactive metabolites, our previously designed A199S/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), known as cocaine hydrolase-1 (CocH1), possesses the desirably high catalytic activity against cocaine. The C-terminus of CocH1, truncated after amino acid #529, was fused to human serum albumin (HSA) to extend the biological half-life. The C-terminal HSA-fused CocH1 (CocH1-HSA), known as Albu-CocH1, Albu-CocH, AlbuBChE, Albu-BChE, or TV-1380 in literature, has shown favorable preclinical and clinical profiles. However, the actual therapeutic value of TV-1380 for cocaine addiction treatment is still limited by the short half-life. In this study, we designed and tested a new type of HSA-fused CocH1 proteins, i.e., N-terminal HSA-fused CocH1, with or without a linker between the HSA and CocH1 domains. It has been demonstrated that the catalytic activity of these new fusion proteins against cocaine is similar to that of TV-1380. However, HSA-CocH1 (without a linker) has a significantly longer biological half-life (t1/2 = 14 ± 2 h) compared to the corresponding C-terminal HSA-fused CocH1, i.e., CocH1-HSA (TV-1380 with t1/2 = 5-8 h), in rats. Further, the N-terminal HSA-fused CocH1 proteins with a linker have further prolonged biological half-lives: t1/2 = 17 ± 2 h for both HSA-EAAAK-CocH1 and HSA-PAPAP-CocH1, and t1/2 = 18 ± 3 h for HSA-(PAPAP)2-CocH1. These N-terminal HSA-fused CocH1 proteins may serve as more promising protein drug candidates for cocaine addiction treatment.
Collapse
Affiliation(s)
- Yingting Cai
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Shuo Zhou
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Zhenyu Jin
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Huimei Wei
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Linyue Shang
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Jing Deng
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA.
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky, 40536, USA.
| |
Collapse
|
7
|
Liu SY, Xiong H, Yang JQ, Yang SH, Li Y, Yang WC, Yang GF. Discovery of Butyrylcholinesterase-Activated Near-Infrared Fluorogenic Probe for Live-Cell and In Vivo Imaging. ACS Sens 2018; 3:2118-2128. [PMID: 30203965 DOI: 10.1021/acssensors.8b00697] [Citation(s) in RCA: 56] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Butyrylcholinesterase (BChE) is widely distributed in various tissues and highly implicated in several important human diseases, especially Alzheimer's disease (AD). However, the role of BChE in AD is still controversial, which may be partially attributed to the lack of a direct tool for real-time and noninvasive monitoring of BChE in in vivo. Here, we report three rationally designed near-infrared fluorogenic probes that possess excellent discrimination for butyrylcholinesterase (BChE) over the related enzyme acetylcholinesterase (AChE). The refined probe, BChE-NIRFP, not only functions as an exquisite substrate for BChE in in vitro assays but also represents a superb "signal-on" imaging tool to real-time track BChE levels in human cells, zebrafish, and a mouse model of AD. A further application of BChE-NIRFP to identify the cellular mechanism reveals that Aβ fibrils and insulin resistance may be important contributors to the abnormally elevated BChE levels observed during AD progression. Based on the results from the present study, this new probe is a valuable tool for basic and clinical research designed to obtain a complete understanding of the physiological roles of BChE in diverse human diseases, particularly AD.
Collapse
Affiliation(s)
- Shi-Yu Liu
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
| | - Hao Xiong
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
| | - Jia-Qian Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
| | - Shu-Hou Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
| | - Yingfu Li
- Department of Biochemistry and Biomedical Sciences, McMaster University, 1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada
| | - Wen-Chao Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
| | - Guang-Fu Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, International Joint Research Center for Intelligent Biosensor Technology and Health, and Chemical Biology Center, College of Chemistry, Central China Normal University, Wuhan 430079, P.R. China
- Collaborative Innovation Center of Chemical Science and Engineering, Tianjin 30071, P.R. China
| |
Collapse
|
8
|
Yao J, Chen X, Zheng F, Zhan CG. Catalytic Reaction Mechanism for Drug Metabolism in Human Carboxylesterase-1: Cocaine Hydrolysis Pathway. Mol Pharm 2018; 15:3871-3880. [PMID: 30095924 DOI: 10.1021/acs.molpharmaceut.8b00354] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Carboxylesterase-1 (CE-1) is a crucial enzyme responsible for metabolism/activation/inactivation of xenobiotics (therapeutic agents, prodrugs, abused drugs, and organophosphorus nerve agents etc.) and also involved in many other biological processes. In this study, we performed extensive computational modeling and simulations to understand the fundamental reaction mechanism of cocaine hydrolysis catalyzed by CE-1, revealing that CE-1-catalyzed cocaine hydrolysis follows a novel reaction pathway with only two reaction steps: a single-step acylation process and a single-step deacylation process. In the transition states of both single-step processes, the cocaine NH group joins the oxyanion hole to form an additional hydrogen bond with the negatively charged carbonyl oxygen atom of the cocaine. Thus, the transition states are stabilized by both intermolecular and intramolecular hydrogen bonds with the methyl ester of cocaine, specifically the carbonyl oxygen atom. The rate-limiting transition state is associated with the acylation process, and the activation free energy barrier was predicted to be 20.1 kcal/mol. Further, in vitro experimental kinetic analysis was performed for human CE-1-catalyzed cocaine hydrolysis. For CE-1-catalyzed cocaine hydrolysis, the computationally predicted free energy barrier (20.1 kcal/mol) is reasonably close to the experimentally derived turnover number ( kcat = 0.058 min-1), indicating the reasonability of the computational results. The obtained novel mechanistic insights are expected to benefit not only CE-1 related rational drug discovery but also future research on the catalytic mechanism of other esterases.
Collapse
Affiliation(s)
- Jianzhuang Yao
- School of Biological Science and Techonology , University of Jinan , Jinan 250022 , China
| | - Xiabin Chen
- School of Medicine , Hangzhou Normal University , Hangzhou 311121 , China
| | | | | |
Collapse
|
9
|
Garcia-Borràs M, Houk KN, Jiménez-Osés G. Computational Design of Protein Function. COMPUTATIONAL TOOLS FOR CHEMICAL BIOLOGY 2017. [DOI: 10.1039/9781788010139-00087] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.
Collapse
Affiliation(s)
- Marc Garcia-Borràs
- Department of Chemistry and Biochemistry, University of California Los Angeles California CA 90095-1569 USA
| | - Kendall N. Houk
- Department of Chemistry and Biochemistry, University of California Los Angeles California CA 90095-1569 USA
| | - Gonzalo Jiménez-Osés
- Departamento de Química, Centro de Investigación en Síntesis Química Universidad de La Rioja 26006 Logroño La Rioja Spain
| |
Collapse
|
10
|
Plant-expressed cocaine hydrolase variants of butyrylcholinesterase exhibit altered allosteric effects of cholinesterase activity and increased inhibitor sensitivity. Sci Rep 2017; 7:10419. [PMID: 28874829 PMCID: PMC5585256 DOI: 10.1038/s41598-017-10571-z] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2017] [Accepted: 08/10/2017] [Indexed: 12/04/2022] Open
Abstract
Butyrylcholinesterase (BChE) is an enzyme with broad substrate and ligand specificities and may function as a generalized bioscavenger by binding and/or hydrolyzing various xenobiotic agents and toxicants, many of which target the central and peripheral nervous systems. Variants of BChE were rationally designed to increase the enzyme’s ability to hydrolyze the psychoactive enantiomer of cocaine. These variants were cloned, and then expressed using the magnICON transient expression system in plants and their enzymatic properties were investigated. In particular, we explored the effects that these site-directed mutations have over the enzyme kinetics with various substrates of BChE. We further compared the affinity of various anticholinesterases including organophosphorous nerve agents and pesticides toward these BChE variants relative to the wild type enzyme. In addition to serving as a therapy for cocaine addiction-related diseases, enhanced bioscavenging against other harmful agents could add to the practicality and versatility of the plant-derived recombinant enzyme as a multivalent therapeutic.
Collapse
|
11
|
Yang SH, Sun Q, Xiong H, Liu SY, Moosavi B, Yang WC, Yang GF. Discovery of a butyrylcholinesterase-specific probe via a structure-based design strategy. Chem Commun (Camb) 2017; 53:3952-3955. [DOI: 10.1039/c7cc00577f] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We report herein the structure-based design and application of a fluorogenic butyrylcholinesterase probe that could discriminate butyrylcholinesterase from acetylcholinesterase.
Collapse
Affiliation(s)
- Shu-Hou Yang
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Qi Sun
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Hao Xiong
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Shi-Yu Liu
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Behrooz Moosavi
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Wen-Chao Yang
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| | - Guang-Fu Yang
- Key Laboratory of Pesticide and Chemical Biology
- Ministry of Education
- College of Chemistry
- Central China Normal University
- Wuhan 430079
| |
Collapse
|
12
|
Zhang Y, Huang X, Han K, Zheng F, Zhan CG. Free energy profiles of cocaine esterase-cocaine binding process by molecular dynamics and potential of mean force simulations. Chem Biol Interact 2016; 259:142-147. [PMID: 27163853 DOI: 10.1016/j.cbi.2016.05.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2016] [Revised: 04/29/2016] [Accepted: 05/05/2016] [Indexed: 11/29/2022]
Abstract
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 KM 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.
Collapse
Affiliation(s)
- Yuxin Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Zhongshan Road 457, Dalian 116023, PR China; Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA; Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Xiaoqin Huang
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Keli Han
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Zhongshan Road 457, Dalian 116023, PR China
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA; Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA; Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.
| |
Collapse
|
13
|
Yao J, Yuan Y, Zheng F, Zhan CG. Unexpected Reaction Pathway for butyrylcholinesterase-catalyzed inactivation of "hunger hormone" ghrelin. Sci Rep 2016; 6:22322. [PMID: 26922910 PMCID: PMC4770301 DOI: 10.1038/srep22322] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2015] [Accepted: 02/12/2016] [Indexed: 11/16/2022] Open
Abstract
Extensive computational modeling and simulations have been carried out, in the present study, to uncover the fundamental reaction pathway for butyrylcholinesterase (BChE)-catalyzed hydrolysis of ghrelin, demonstrating that the acylation process of BChE-catalyzed hydrolysis of ghrelin follows an unprecedented single-step reaction pathway and the single-step acylation process is rate-determining. The free energy barrier (18.8 kcal/mol) calculated for the rate-determining step is reasonably close to the experimentally-derived free energy barrier (~19.4 kcal/mol), suggesting that the obtained mechanistic insights are reasonable. The single-step reaction pathway for the acylation is remarkably different from the well-known two-step acylation reaction pathway for numerous ester hydrolysis reactions catalyzed by a serine esterase. This is the first time demonstrating that a single-step reaction pathway is possible for an ester hydrolysis reaction catalyzed by a serine esterase and, therefore, one no longer can simply assume that the acylation process must follow the well-known two-step reaction pathway.
Collapse
Affiliation(s)
- Jianzhuang Yao
- Molecular Modeling and Biopharmaceutical Center and College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.,Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Yaxia Yuan
- Molecular Modeling and Biopharmaceutical Center and College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.,Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center and College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.,Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center and College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.,Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| |
Collapse
|
14
|
Yao Y, Liu J, Zheng F, Zhan CG. Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine. Theor Chem Acc 2016; 135. [PMID: 28250715 DOI: 10.1007/s00214-015-1788-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yuan Yao
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536; The Academy of Fundamental and Interdisciplinary Sciences, Harbin Institute of Technology, Harbin 150080, P.R. China
| | - Junjun Liu
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536; Tongji School of Pharmacy, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei 430030, P.R. China
| | - Fang Zheng
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center and Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| |
Collapse
|
15
|
Otzen M, Crismaru CG, Postema CP, Wijma HJ, Heberling MM, Szymanski W, de Wildeman S, Janssen DB. Metabolism of β-valine via a CoA-dependent ammonia lyase pathway. Appl Microbiol Biotechnol 2015; 99:8987-98. [PMID: 26004802 PMCID: PMC4619459 DOI: 10.1007/s00253-015-6551-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2014] [Revised: 02/24/2015] [Accepted: 03/17/2015] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Marleen Otzen
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Ciprian G Crismaru
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Christiaan P Postema
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Hein J Wijma
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Matthew M Heberling
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Wiktor Szymanski
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands
| | - Stefaan de Wildeman
- DSM Pharmaceutical Products, Geleen, The Netherlands.,BioBased Materials, Faculty of Humanities and Sciences, Maastricht University, Chemelot, The Netherlands
| | - Dick B Janssen
- Department of Biochemistry, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG, Groningen, The Netherlands.
| |
Collapse
|
16
|
Yuan Y, Huang X, Midde NM, Quizon PM, Sun WL, Zhu J, Zhan CG. Molecular mechanism of HIV-1 Tat interacting with human dopamine transporter. ACS Chem Neurosci 2015; 6:658-665. [PMID: 25695767 DOI: 10.1021/acschemneuro.5b00001] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Yaxia Yuan
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Xiaoqin Huang
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Narasimha M. Midde
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Pamela M. Quizon
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Wei-Lun Sun
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Jun Zhu
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center and ‡Department of
Pharmaceutical Sciences,
College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department of Drug Discovery and Biomedical Sciences, South Carolina
College of Pharmacy, University of South Carolina, Columbia, South Carolina 29208, United States
| |
Collapse
|
17
|
Sharir-Ivry A, Varatharaj R, Shurki A. Valence bond and enzyme catalysis: a time to break down and a time to build up. Chemistry 2015; 21:7159-69. [PMID: 25808731 DOI: 10.1002/chem.201406236] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2014] [Indexed: 11/12/2022]
Abstract
Understanding enzyme catalysis and developing ability to control of it are two great challenges in biochemistry. A few successful examples of computational-based enzyme design have proved the fantastic potential of computational approaches in this field, however, relatively modest rate enhancements have been reported and the further development of complementary methods is still required. Herein we propose a conceptually simple scheme to identify the specific role that each residue plays in catalysis. The scheme is based on a breakdown of the total catalytic effect into contributions of individual protein residues, which are further decomposed into chemically interpretable components by using valence bond theory. The scheme is shown to shed light on the origin of catalysis in wild-type haloalkane dehalogenase (wt-DhlA) and its mutants. Furthermore, the understanding gained through our scheme is shown to have great potential in facilitating the selection of non-optimal sites for catalysis and suggesting effective mutations to enhance the enzymatic rate.
Collapse
Affiliation(s)
- Avital Sharir-Ivry
- Institute for Drug Design, School of Pharmacy, The Lise Meitner-Minerva Center for Computational Quantum Chemistry, The Hebrew University of Jerusalem, Jerusalem 91120 (Israel), Fax: (+972) 2-675-7076
| | | | | |
Collapse
|
18
|
Zheng F, Xue L, Hou S, Liu J, Zhan M, Yang W, Zhan CG. A highly efficient cocaine-detoxifying enzyme obtained by computational design. Nat Commun 2014; 5:3457. [PMID: 24643289 PMCID: PMC3996704 DOI: 10.1038/ncomms4457] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2013] [Accepted: 02/17/2014] [Indexed: 11/09/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Liu Xue
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Shurong Hou
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Junjun Liu
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Max Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Wenchao Yang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| |
Collapse
|
19
|
Zhan CG. Novel pharmacological approaches to treatment of drug overdose and addiction. Expert Rev Clin Pharmacol 2014; 2:1-4. [PMID: 21072135 DOI: 10.1586/17512433.2.1.1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Chang-Guo Zhan
- University of Kentucky, Department of Pharmaceutical Sciences, College of Pharmacy, 789 South Limestone Street, Lexington, KY 40536, USA TEL: 859-323-3943 FAX: 859-323-3575
| |
Collapse
|
20
|
Hou S, Xue L, Yang W, Fang L, Zheng F, Zhan CG. Substrate selectivity of high-activity mutants of human butyrylcholinesterase. Org Biomol Chem 2013; 11:7477-85. [PMID: 24077614 PMCID: PMC3836059 DOI: 10.1039/c3ob41713a] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shurong Hou
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA.
| | | | | | | | | | | |
Collapse
|
21
|
Hao GF, Tan Y, Yang SG, Wang ZF, Zhan CG, Xi Z, Yang GF. Computational and experimental insights into the mechanism of substrate recognition and feedback inhibition of protoporphyrinogen oxidase. PLoS One 2013; 8:e69198. [PMID: 23935953 PMCID: PMC3720618 DOI: 10.1371/journal.pone.0069198] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2013] [Accepted: 06/05/2013] [Indexed: 11/20/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ge-Fei Hao
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P. R. China
| | - Ying Tan
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P. R. China
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin, P. R. China
| | - Sheng-Gang Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P. R. China
| | - Zhi-Fang Wang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P. R. China
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin, P. R. China
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, Lexington, Kentucky, United States of America
| | - Zhen Xi
- State Key Laboratory of Elemento-Organic Chemistry, Nankai University, Tianjin, P. R. China
| | - Guang-Fu Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P. R. China
| |
Collapse
|
22
|
van der Kamp MW, Mulholland AJ. Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. Biochemistry 2013; 52:2708-28. [PMID: 23557014 DOI: 10.1021/bi400215w] [Citation(s) in RCA: 399] [Impact Index Per Article: 36.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
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.
Collapse
Affiliation(s)
- Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK.
| | | |
Collapse
|
23
|
Yang H, Wong MW. Oxyanion Hole Stabilization by C–H···O Interaction in a Transition State—A Three-Point Interaction Model for Cinchona Alkaloid-Catalyzed Asymmetric Methanolysis of meso-Cyclic Anhydrides. J Am Chem Soc 2013; 135:5808-18. [DOI: 10.1021/ja4005893] [Citation(s) in RCA: 130] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hui Yang
- Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543
| | - Ming Wah Wong
- Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543
| |
Collapse
|
24
|
Xue L, Hou S, Yang W, Fang L, Zheng F, Zhan CG. Catalytic activities of a cocaine hydrolase engineered from human butyrylcholinesterase against (+)- and (-)-cocaine. Chem Biol Interact 2013; 203:57-62. [PMID: 22917637 PMCID: PMC3527670 DOI: 10.1016/j.cbi.2012.08.003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2012] [Revised: 07/31/2012] [Accepted: 08/02/2012] [Indexed: 11/27/2022]
Abstract
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.
Collapse
Affiliation(s)
| | | | - Wenchao Yang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Lei Fang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| |
Collapse
|
25
|
Liu J, Zhan CG. Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine. J Chem Theory Comput 2012; 8:1426-1435. [PMID: 23066354 PMCID: PMC3469271 DOI: 10.1021/ct200810d] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Junjun Liu
- Tongji School of Pharmacy, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei 430030, P.R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
26
|
Huang X, Zhao X, Zheng F, Zhan CG. Cocaine esterase-cocaine binding process and the free energy profiles by molecular dynamics and potential of mean force simulations. J Phys Chem B 2012; 116:3361-8. [PMID: 22385120 DOI: 10.1021/jp2111605] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Xiaoqin Huang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, USA
| | | | | | | |
Collapse
|
27
|
Chen X, Fang L, Liu J, Zhan CG. Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine. Biochemistry 2012; 51:1297-305. [PMID: 22304234 DOI: 10.1021/bi201786s] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
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).
Collapse
Affiliation(s)
- Xi Chen
- Key Laboratory of Pesticide & Chemical Biology of the Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, PR China
| | | | | | | |
Collapse
|
28
|
Yang X, Yuan Y, Zhan CG, Liao F. Uricases as therapeutic agents to treat refractory gout: Current states and future directions. Drug Dev Res 2011; 73:66-72. [PMID: 22665944 DOI: 10.1002/ddr.20493] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Xiaolan Yang
- Unit for Analytical Probes and Protein Biotechnology, Key Laboratory of Medical Laboratory Diagnosis of the Education Ministry of China, College of Laboratory Medicine, Chongqing Medical University, Chongqing 400016, China
| | | | | | | |
Collapse
|
29
|
Huang X, Zheng F, Zhan CG. Human butyrylcholinesterase-cocaine binding pathway and free energy profiles by molecular dynamics and potential of mean force simulations. J Phys Chem B 2011; 115:11254-60. [PMID: 21902185 DOI: 10.1021/jp2047807] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present study, we have performed combined molecular dynamics and potential of mean force (PMF) simulations to determine the enzyme-substrate (ES) binding pathway and the corresponding free energy profiles for wild-type butyrylcholinesterase (BChE) binding with (-)/(+)-cocaine and for the A328W/Y332G mutant binding with (-)-cocaine. According to the PMF simulations, for each ES binding system, the substrate first binds with the enzyme at a peripheral anionic site around the entrance of the active-site gorge to form the first ES complex (ES1-like) during the binding process. Further evolution from the ES1-like complex to the nonprereactive ES complex is nearly barrierless, with a free energy barrier lower than 1.0 kcal/mol. So, the nonprereactive ES binding process should be very fast. The rate-determining step of the entire ES binding process is the subsequent evolution from the nonprereactive ES complex to the prereactive ES complex. Further accounting for the entire ES binding process, the PMF-based simulations qualitatively reproduced the relative order of the experimentally derived binding free energies (ΔG(bind)), although the simulations systematically overestimated the magnitude of the binding affinity and systematically underestimated the differences between the ΔG(bind) values. The obtained structural and energetic insights into the entire ES binding process provide a valuable base for future rational design of high-activity mutants of BChE as candidates for an enzyme therapy for cocaine overdose and abuse.
Collapse
Affiliation(s)
- Xiaoqin Huang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, USA
| | | | | |
Collapse
|
30
|
Fang L, Pan Y, Muzyka JL, Zhan CG. Active site gating and substrate specificity of butyrylcholinesterase and acetylcholinesterase: insights from molecular dynamics simulations. J Phys Chem B 2011; 115:8797-805. [PMID: 21682268 PMCID: PMC3135420 DOI: 10.1021/jp112030p] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Lei Fang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536
| | - Yongmei Pan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536
| | | | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536
| |
Collapse
|
31
|
Recent progress in protein drug design and discovery with a focus on novel approaches to the development of anti-cocaine medications. Future Med Chem 2011; 1:515-28. [PMID: 20161378 DOI: 10.4155/fmc.09.20] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
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.
Collapse
|
32
|
Liu J, Zhao X, Yang W, Zhan CG. Reaction mechanism for cocaine esterase-catalyzed hydrolyses of (+)- and (-)-cocaine: unexpected common rate-determining step. J Phys Chem B 2011; 115:5017-25. [PMID: 21486046 PMCID: PMC3087188 DOI: 10.1021/jp200975v] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
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.
Collapse
Affiliation(s)
| | | | - Wenchao Yang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
33
|
Huang X, Pan Y, Zheng F, Zhan CG. 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. J Phys Chem B 2011; 114:13545-54. [PMID: 20883001 DOI: 10.1021/jp106539w] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Xiaoqin Huang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 S. Limestone, Lexington, Kentucky 40536, USA
| | | | | | | |
Collapse
|
34
|
Xue L, Ko MC, Tong M, Yang W, Hou S, Fang L, Liu J, Zheng F, Woods JH, Tai HH, Zhan CG. Design, preparation, and characterization of high-activity mutants of human butyrylcholinesterase specific for detoxification of cocaine. Mol Pharmacol 2011; 79:290-7. [PMID: 20971807 PMCID: PMC3033707 DOI: 10.1124/mol.110.068494] [Citation(s) in RCA: 70] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2010] [Accepted: 10/22/2010] [Indexed: 11/22/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Liu Xue
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street Lexington, KY 40536, USA
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
35
|
Abstract
“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.”
Collapse
Affiliation(s)
- Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536, USA
| |
Collapse
|
36
|
Chen X, Fang L, Liu J, Zhan CG. Reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of acetylcholine. J Phys Chem B 2010; 115:1315-22. [PMID: 21175195 DOI: 10.1021/jp110709a] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Xi Chen
- Key Laboratory of Pesticide & Chemical Biology of the Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, P R China
| | | | | | | |
Collapse
|
37
|
Zheng F, Yang W, Xue L, Hou S, Liu J, Zhan CG. Design of high-activity mutants of human butyrylcholinesterase against (-)-cocaine: structural and energetic factors affecting the catalytic efficiency. Biochemistry 2010; 49:9113-9. [PMID: 20886866 DOI: 10.1021/bi1011628] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The present study was aimed to explore the correlation between the protein structure and catalytic efficiency of butyrylcholinesterase (BChE) mutants against (-)-cocaine by modeling the rate-determining transition state (TS1), i.e., the transition state for the first step of chemical reaction process, of (-)-cocaine hydrolysis catalyzed by various mutants of human BChE in comparison with the wild type. Molecular modeling of the TS1 structures revealed that mutations on certain nonactive site residues can indirectly affect the catalytic efficiency of the enzyme against (-)-cocaine through enhancing or weakening the overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. Computational insights and predictions were supported by the catalytic activity data obtained from wet experimental tests on the mutants of human BChE, including five new mutants reported for the first time. The BChE mutants with at least ∼1000-fold improved catalytic efficiency against (-)-cocaine compared to the wild-type BChE are all associated with the TS1 structures having stronger overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. The combined computational and experimental data demonstrate a reasonable correlation relationship between the hydrogen-bonding distances in the TS1 structure and the catalytic efficiency of the enzyme against (-)-cocaine.
Collapse
Affiliation(s)
- Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, Kentucky 40536, United States
| | | | | | | | | | | |
Collapse
|
38
|
Yang W, Xue L, Fang L, Chen X, Zhan CG. Characterization of a high-activity mutant of human butyrylcholinesterase against (-)-cocaine. Chem Biol Interact 2010; 187:148-52. [PMID: 20060817 PMCID: PMC2891156 DOI: 10.1016/j.cbi.2010.01.004] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2009] [Revised: 12/30/2009] [Accepted: 01/04/2010] [Indexed: 11/30/2022]
Abstract
Cocaine addiction and overdose are a well-known public health problem. There is no approved medication available for cocaine abuse treatment. Our recently designed and discovered high-activity mutant (A199S/S287G/A328W/Y332G) of human butyrylcholinesterase (BChE) has been recognized to be worth exploring for clinical application in humans as a potential anti-cocaine medication. The catalytic rate constant (k(cat)) and Michaelis-Menten constant (K(M)) for (-)-cocaine hydrolysis catalyzed by A199S/S287G/A328W/Y332G BChE (without fusion with any other peptide) have been determined to be 3,060 min(-1) and 3.1 microM, respectively, in the present study. The determined kinetic parameters reveal that the un-fused A199S/S287G/A328W/Y332G mutant has a approximately 1,080-fold improved catalytic efficiency (k(cat)/K(M)) against (-)-cocaine compared to the wild-type BChE. The approximately 1,080-fold improvement in the catalytic efficiency of the un-fused A199S/S287G/A328W/Y332G mutant is very close to the previously reported the approximately 1,000-fold improvement in the catalytic efficiency of the A199S/S287G/A328W/Y332G mutant fused with human serum albumin. These results suggest that the albumin fusion did not significantly change the catalytic efficiency of the BChE mutant while extending the plasma half-life. In addition, we have also examined the catalytic activities of the A199S/S287G/A328W/Y332G mutant against two other substrates, acetylthiocholine (ATC) and butyrylthiocholine (BTC). It has been shown that the A199S/S287G/A328W/Y332G mutations actually decreased the catalytic efficiencies of BChE against ATC and BTC, while considerably improving the catalytic efficiency of BChE against (-)-cocaine.
Collapse
Affiliation(s)
- Wenchao Yang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, USA
| | - Liu Xue
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, USA
| | - Lei Fang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, USA
| | - Xi Chen
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, USA
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536, USA
| |
Collapse
|
39
|
Yang W, Pan Y, Fang L, Gao D, Zheng F, Zhan CG. Free energy perturbation simulation on transition states and high-activity mutants of human butyrylcholinesterase for (-)-cocaine hydrolysis. J Phys Chem B 2010; 114:10889-96. [PMID: 20677742 PMCID: PMC2930763 DOI: 10.1021/jp104989b] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A unified computational approach based on free energy perturbation (FEP) simulations of transition states has been employed to calculate the mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state for (-)-cocaine hydrolysis catalyzed by the currently most promising series of mutants of human butyrylcholinesterase (BChE) that contain the A199S/A328W/Y332G mutations. The FEP simulations were followed by Michaelis-Menten kinetics analysis determining the individual k(cat) and K(M) values missing for the A199S/F227A/A328W/Y332G mutant in this series. The calculated mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state are in good agreement with the experimental kinetic data, demonstrating that the unified computational approach based on the FEP simulations of the transition states may be valuable for future computational design of new BChE mutants with a further improved catalytic efficiency against (-)-cocaine.
Collapse
Affiliation(s)
- Wenchao Yang
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P. R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Yongmei Pan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Lei Fang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Daquan Gao
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone, Lexington, KY 40536
| |
Collapse
|
40
|
Computational simulations of structural role of the active-site W374C mutation of acetyl-coenzyme-A carboxylase: Multi-drug resistance mechanism. J Mol Model 2010; 17:495-503. [DOI: 10.1007/s00894-010-0742-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2010] [Accepted: 05/04/2010] [Indexed: 10/19/2022]
|
41
|
Zhao PL, Wang L, Zhu XL, Huang X, Zhan CG, Wu JW, Yang GF. Subnanomolar inhibitor of cytochrome bc1 complex designed by optimizing interaction with conformationally flexible residues. J Am Chem Soc 2010; 132:185-94. [PMID: 19928849 DOI: 10.1021/ja905756c] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cytochrome bc(1) complex (EC 1.10.2.2, bc(1)), an essential component of the cellular respiratory chain and the photosynthetic apparatus in photosynthetic bacteria, has been identified as a promising target for new drugs and agricultural fungicides. X-ray diffraction structures of the free bc(1) complex and its complexes with various inhibitors revealed that the phenyl group of Phe274 in the binding pocket exhibited significant conformational flexibility upon different inhibitors binding to optimize respective pi-pi interactions, whereas the side chains of other hydrophobic residues showed conformational stability. Therefore, in the present study, a strategy of optimizing the pi-pi interaction with conformationally flexible residues was proposed to design and discover new bc(1) inhibitors with a higher potency. Eight new compounds were designed and synthesized, among which compound 5c, with a K(i) value of 570 pM, was identified as the most promising drug or fungicide candidate, significantly more potent than the commercially available bc(1) inhibitors, including azoxystrobin (AZ), kresoxim-methyl (KM), and pyraclostrobin (PY). To our knowledge, this is the first bc(1) inhibitor discovered from structure-based design with a potency of subnanomolar K(i) value. For all of the compounds synthesized and assayed, the calculated binding free energies correlated reasonably well with the binding free energies derived from the experimental K(i) values, with a correlation coefficient of r(2) = 0.89. The further inhibitory kinetics studies revealed that 5c is a noncompetitive inhibitor with respect to substrate cytochrome c, but it is a competitive inhibitor with respect to substrate ubiquinol. Due to its subnanomolar K(i) potency and slow dissociation rate constant (k(-0) = 0.00358 s(-1)), 5c could be used as a specific probe for further elucidation of the mechanism of bc(1) function and as a new lead compound for future drug discovery.
Collapse
Affiliation(s)
- Pei-Liang Zhao
- Key Laboratory of Pesticide & Chemical Biology, College of Chemistry, Ministry of Education, Central China Normal University, Wuhan 430079, PR China
| | | | | | | | | | | | | |
Collapse
|
42
|
Huang X, Gu HH, Zhan CG. Mechanism for cocaine blocking the transport of dopamine: insights from molecular modeling and dynamics simulations. J Phys Chem B 2010; 113:15057-66. [PMID: 19831380 DOI: 10.1021/jp900963n] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Molecular modeling and dynamics simulations have been performed to study how cocaine inhibits dopamine transporter (DAT) for the transport of dopamine. The computationally determined DAT-ligand binding mode is totally different from the previously proposed overlap binding mode in which cocaine- and dopamine-binding sites are the same (Beuming, T.; et al. Nat. Neurosci. 2008, 11, 780-789). The new cocaine-binding site does not overlap with, but is close to, the dopamine-binding site. Analysis of all results reveals that when cocaine binds to DAT, the initial binding site is likely the one modeled in this study because this binding site can naturally accommodate cocaine. Then cocaine may move to the dopamine-binding site after DAT makes some necessary conformational change and expands the binding site cavity. It has been demonstrated that cocaine may inhibit the transport of dopamine through both blocking the initial DAT-dopamine binding and reducing the kinetic turnover of the transporter following the DAT-dopamine binding. The relative contributions to the phenomenological inhibition of the transport of dopamine from blocking the initial binding and reducing the kinetic turnover can be different in different types of assays. The obtained general structural and mechanistic insights are consistent with available experimental data and could be valuable for guiding future studies toward understanding cocaine's inhibiting of other transporters.
Collapse
Affiliation(s)
- Xiaoqin Huang
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536, USA
| | | | | |
Collapse
|
43
|
Zhu XL, Ge-Fei H, Zhan CG, Yang GF. Computational simulations of the interactions between acetyl-coenzyme-A carboxylase and clodinafop: resistance mechanism due to active and nonactive site mutations. J Chem Inf Model 2009; 49:1936-43. [PMID: 19594140 DOI: 10.1021/ci900174d] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Grass weed populations resistant to acetyl-CoA carboxylase-inhibiting (ACCase; EC 6.4.1.2) herbicides represent a major problem for the sustainable development of modern agriculture. In the present study, extensive computational simulations, including homology modeling, molecular dynamics (MD) simulations, and molecular mechanics-Poisson-Boltzmann surface area (MM/PBSA) calculations, have been carried out to uncover the detailed molecular mechanism of Alopecurus myosuroides resistance to clodinafop, a commercial herbicide targeting ACCase. All the computational model and energetic results indicated that W374C, I388N, D425G, and G443A mutations have great effects on the conformational change of the binding pocket and the hydrogen-bonding interactions. The pi-pi interaction between ligand and the residue of Phe377 and Tyr161', playing an important contribution to the binding affinity, were decreased after mutations. In addition, the hydrogen-bonding interactions between clodinafop and the residues (Ile158' and Ala54') disappeared or decreased significantly upon mutation. As a result, the mutant-type ACCase has a lower affinity for the inhibitor binding than the wild-type enzyme, which accounts for the molecular basis of herbicidal resistance. The structural and mechanistic insights obtained from the present study will provide a valuable clue for future designing of a promising inhibitor to reduce drug resistance associated with both active and nonactive site mutations.
Collapse
Affiliation(s)
- Xiao-Lei Zhu
- Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, PR China
| | | | | | | |
Collapse
|
44
|
Liu J, Hamza A, Zhan CG. Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase. J Am Chem Soc 2009; 131:11964-75. [PMID: 19642701 PMCID: PMC2738781 DOI: 10.1021/ja903990p] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The fundamental reaction mechanism of cocaine esterase (CocE)-catalyzed hydrolysis of (-)-cocaine and the corresponding free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, the entire hydrolysis reaction consists of four reaction steps, including the nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by the hydroxyl group of Ser117, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally dissociation between the (-)-cocaine benzoyl group and Ser117 of CocE. The third reaction step involving the nucleophilic attack of a water molecule was found to be rate-determining, which is remarkably different from (-)-cocaine hydrolysis catalyzed by wild-type butyrylcholinesterase (BChE; where the formation of the prereactive BChE-(-)-cocaine complex is rate-determining) or its mutants containing Tyr332Gly or Tyr332Ala mutation (where the first chemical reaction step is rate-determining). Besides, the role of Asp259 in the catalytic triad of CocE does not follow the general concept of the "charge-relay system" for all serine esterases. The free energy barrier calculated for the rate-determining step of CocE-catalyzed hydrolysis of (-)-cocaine is 17.9 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 16.2 kcal/mol. In the present study, where many sodium ions are present, the effects of counterions are found to be significant in determining the free energy barrier. The finding of the significant effects of counterions on the free energy barrier may also be valuable in guiding future mechanistic studies on other charged enzymes.
Collapse
Affiliation(s)
- Junjun Liu
- Key Laboratory of Pesticide and Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P. R. China
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| | - Adel Hamza
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, KY 40536
| |
Collapse
|
45
|
Bioactive conformation analysis of cyclic imides as protoporphyrinogen oxidase inhibitor by combining DFT calculations, QSAR and molecular dynamic simulations. Bioorg Med Chem 2009; 17:4935-42. [DOI: 10.1016/j.bmc.2009.06.003] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2009] [Revised: 05/31/2009] [Accepted: 06/02/2009] [Indexed: 11/20/2022]
|
46
|
Pan Y, Muzyka JL, Zhan CG. Model of human butyrylcholinesterase tetramer by homology modeling and dynamics simulation. J Phys Chem B 2009; 113:6543-52. [PMID: 19402731 DOI: 10.1021/jp8114995] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A mutant of human butyrylcholinesterase (BChE) with high activity against cocaine would be highly promising as a drug for therapeutic treatment of cocaine abuse and overdose. It is desirable to design a recombinant BChE mutant with a long half-life in human circulation. Studies showed that BChE subunits can be assembled by a peptide containing the proline-rich attachment domain (PRAD) to form a stable tetramer. The models of BChE tetramer complexed with PRAD with various sequences have been constructed, in the present study, on the basis of homology modeling and molecular dynamics simulation of explicit water-solvated systems. The 3D models enable us to understand how the BChE subunits are arranged in the tetramer and how the tetramerization domain of BChE is associated with PRAD to form a stable tetramer of human BChE. It has been shown that the six conserved hydrophobic residues located on the C-terminal of BChE are responsible for the key electrostatic and hydrophobic interactions between the tetramerization domain of BChE and PRAD. The simulated tetramer structures suggest that mutation of three residues, i.e., Phe547, Met554, and Phe561, to other hydrophobic residues may be beneficial for increasing the binding between the tetramerization domain of BChE and PRAD. Thus, the detailed structural insights obtained from this study may be valuable for rational design of a recombinant BChE tetramer with a longer residence time in circulation.
Collapse
Affiliation(s)
- Yongmei Pan
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536, USA
| | | | | |
Collapse
|
47
|
Hao GF, Zhu XL, Ji FQ, Zhang L, Yang GF, Zhan CG. Understanding the mechanism of drug resistance due to a codon deletion in protoporphyrinogen oxidase through computational modeling. J Phys Chem B 2009; 113:4865-75. [PMID: 19284797 DOI: 10.1021/jp807442n] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Protoporphyrinogen oxidase (PPO; EC 1.3.3.4) is the last common enzyme for the enzymatic transformation of protoporphyrinogen-IX to protoporphyrin-IX, which is the key common intermediate leading to heme and chlorophyll. Hence, PPO has been identified as one of the most importance action targets for the treatment of some important diseases including cancer and variegated porphyria (VP). In the agricultural field, PPO inhibitors have been used as herbicides for many years. Recently, a unique drug resistance was found to be associated with a nonactive site residue (Gly210) deletion rather than substitution in A. tuberculatus PPO. In the present study, extensive computational simulations, including homology modeling, molecular dynamics (MD) simulations, and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations, have been carried out to uncover the detailed molecular mechanism of drug resistance associated with Gly210 deletion. Although Gly210 in the wild-type A. tuberculatus PPO has no direct interaction with the inhibitors, all the computational models and energetic results indicated that Gly210 deletion has great effects on the hydrogen-bonding network and the conformational change of the binding pocket. An interchain hydrogen bond between Gly210 with Ser424, playing an important role in stabilizing the local conformation of the wild-type enzyme, disappeared after Gly210 deletion. As a result, the mutant-type PPO has a lower affinity than the wild-type enzyme, which accounts for the molecular mechanism of drug resistance. The structural and mechanistic insights obtained from the present study provide a new starting point for future rational design of novel PPO inhibitors to overcome drug resistance associated with Gly210 deletion.
Collapse
Affiliation(s)
- Ge-Fei Hao
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, PR China
| | | | | | | | | | | |
Collapse
|
48
|
Yang W, Pan Y, Zheng F, Cho H, Tai HH, Zhan CG. Free-energy perturbation simulation on transition states and redesign of butyrylcholinesterase. Biophys J 2009; 96:1931-8. [PMID: 19254552 DOI: 10.1016/j.bpj.2008.11.051] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2008] [Accepted: 11/24/2008] [Indexed: 11/28/2022] Open
Abstract
It is recognized that an ideal anti-cocaine treatment is to accelerate cocaine metabolism by producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., butyrylcholinesterase (BChE)-catalyzed hydrolysis of cocaine. BChE mutants with a higher catalytic activity against (-)-cocaine are highly desired for use as an exogenous enzyme in humans. To develop a rational design for high-activity mutants, we carried out free-energy perturbation (FEP) simulations on various mutations of the transition-state structures in addition to the corresponding free-enzyme structures by using an extended FEP procedure. The FEP simulations on the mutations of both the free-enzyme and transition-state structures allowed us to calculate the mutation-caused shift of the free-energy change from the free enzyme (BChE) to the transition state, and thus to theoretically predict the mutation-caused shift of the catalytic efficiency (k(cat)/K(M)). The computational predictions are supported by the kinetic data obtained from the wet experiments, demonstrating that the FEP-based computational design approach is promising for rational design of high-activity mutants of an enzyme. One of the BChE mutants designed and discovered in this study has an approximately 1800-fold improved catalytic efficiency against (-)-cocaine compared to wild-type BChE. The high-activity mutant may be therapeutically valuable.
Collapse
Affiliation(s)
- Wenchao Yang
- Key Laboratory of Pesticide and Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan, PR China
| | | | | | | | | | | |
Collapse
|
49
|
Design and synthesis of N-2,6-difluorophenyl-5-methoxyl-1,2,4-triazolo[1,5-a]-pyrimidine-2-sulfonamide as acetohydroxyacid synthase inhibitor. Bioorg Med Chem 2009; 17:3011-7. [DOI: 10.1016/j.bmc.2009.03.018] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2009] [Revised: 03/09/2009] [Accepted: 03/10/2009] [Indexed: 11/21/2022]
|
50
|
Zheng F, Yang W, Ko MC, Liu J, Cho H, Gao D, Tong M, Tai HH, Woods JH, Zhan CG. Most efficient cocaine hydrolase designed by virtual screening of transition states. J Am Chem Soc 2008; 130:12148-55. [PMID: 18710224 DOI: 10.1021/ja803646t] [Citation(s) in RCA: 143] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cocaine is recognized as the most reinforcing of all drugs of abuse. There is no anticocaine medication available. The disastrous medical and social consequences of cocaine addiction have made the development of an anticocaine medication a high priority. It has been recognized that an ideal anticocaine medication is one that accelerates cocaine metabolism producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., cocaine hydrolysis catalyzed by plasma enzyme butyrylcholinesterase (BChE). However, wild-type BChE has a low catalytic efficiency against the abused cocaine. Design of a high-activity enzyme mutant is extremely challenging, particularly when the chemical reaction process is rate-determining for the enzymatic reaction. Here we report the design and discovery of a high-activity mutant of human BChE by using a novel, systematic computational design approach based on transition-state simulations and activation energy calculations. The novel computational design approach has led to discovery of the most efficient cocaine hydrolase, i.e., a human BChE mutant with an approximately 2000-fold improved catalytic efficiency, promising for therapeutic treatment of cocaine overdose and addiction as an exogenous enzyme in human. The encouraging discovery resulted from the computational design not only provides a promising anticocaine medication but also demonstrates that the novel, generally applicable computational design approach is promising for rational enzyme redesign and drug discovery.
Collapse
Affiliation(s)
- Fang Zheng
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 725 Rose Street, Lexington, Kentucky 40536, USA
| | | | | | | | | | | | | | | | | | | |
Collapse
|