1
|
Cruz-Hernandez A, Roney A, Goswami DG, Tewari-Singh N, Brown JM. A review of chemical warfare agents linked to respiratory and neurological effects experienced in Gulf War Illness. Inhal Toxicol 2022; 34:412-432. [PMID: 36394251 PMCID: PMC9832991 DOI: 10.1080/08958378.2022.2147257] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2022] [Accepted: 11/07/2022] [Indexed: 11/18/2022]
Abstract
Over 40% of veterans from the Persian Gulf War (GW) (1990-1991) suffer from Gulf War Illness (GWI). Thirty years since the GW, the exposure and mechanism contributing to GWI remain unclear. One possible exposure that has been attributed to GWI are chemical warfare agents (CWAs). While there are treatments for isolated symptoms of GWI, the number of respiratory and cognitive/neurological issues continues to rise with minimum treatment options. This issue does not only affect veterans of the GW, importantly these chronic multisymptom illnesses (CMIs) are also growing amongst veterans who have served in the Afghanistan-Iraq war. What both wars have in common are their regions and inhaled exposures. In this review, we will describe the CWA exposures, such as sarin, cyclosarin, and mustard gas in both wars and discuss the various respiratory and neurocognitive issues experienced by veterans. We will bridge the respiratory and neurological symptoms experienced to the various potential mechanisms described for each CWA provided with the most up-to-date models and hypotheses.
Collapse
Affiliation(s)
- Angela Cruz-Hernandez
- Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, The University of Colorado, Anschutz Medical Campus, Aurora, CO, USA
| | - Andrew Roney
- Department of Pharmacology and Toxicology, College of Osteopathic Medicine, Michigan State University, East Lansing, MI, USA
| | - Dinesh G Goswami
- Department of Pharmacology and Toxicology, College of Osteopathic Medicine, Michigan State University, East Lansing, MI, USA
| | - Neera Tewari-Singh
- Department of Pharmacology and Toxicology, College of Osteopathic Medicine, Michigan State University, East Lansing, MI, USA
| | - Jared M Brown
- Department of Pharmaceutical Sciences, Skaggs School of Pharmacy and Pharmaceutical Sciences, The University of Colorado, Anschutz Medical Campus, Aurora, CO, USA
| |
Collapse
|
2
|
De Boer D, Nguyen N, Mao J, Moore J, Sorin EJ. A Comprehensive Review of Cholinesterase Modeling and Simulation. Biomolecules 2021; 11:580. [PMID: 33920972 PMCID: PMC8071298 DOI: 10.3390/biom11040580] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2021] [Revised: 04/08/2021] [Accepted: 04/11/2021] [Indexed: 01/18/2023] Open
Abstract
The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.
Collapse
Affiliation(s)
- Danna De Boer
- Department of Chemistry & Biochemistry, California State University, Long Beach, CA 90840, USA;
| | - Nguyet Nguyen
- Department of Chemical Engineering, California State University, Long Beach, CA 90840, USA; (N.N.); (J.M.)
| | - Jia Mao
- Department of Chemical Engineering, California State University, Long Beach, CA 90840, USA; (N.N.); (J.M.)
| | - Jessica Moore
- Department of Biomedical Engineering, California State University, Long Beach, CA 90840, USA;
| | - Eric J. Sorin
- Department of Chemistry & Biochemistry, California State University, Long Beach, CA 90840, USA;
| |
Collapse
|
3
|
Enantioseparation, in vitro testing, and structural characterization of triple-binding reactivators of organophosphate-inhibited cholinesterases. Biochem J 2020; 477:2771-2790. [PMID: 32639532 DOI: 10.1042/bcj20200192] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Revised: 07/06/2020] [Accepted: 07/08/2020] [Indexed: 11/17/2022]
Abstract
The enantiomers of racemic 2-hydroxyimino-N-(azidophenylpropyl)acetamide-derived triple-binding oxime reactivators were separated, and tested for inhibition and reactivation of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibited with tabun (GA), cyclosarin (GF), sarin (GB), and VX. Both enzymes showed the greatest affinity toward the methylimidazole derivative (III) of 2-hydroxyimino-N-(azidophenylpropyl)acetamide (I). The crystal structure was determined for the complex of oxime III within human BChE, confirming that all three binding groups interacted with active site residues. In the case of BChE inhibited by GF, oximes I (kr = 207 M-1 min-1) and III (kr = 213 M-1 min-1) showed better reactivation efficiency than the reference oxime 2-PAM. Finally, the key mechanistic steps in the reactivation of GF-inhibited BChE with oxime III were modeled using the PM7R6 method, stressing the importance of proton transfer from Nε of His438 to Oγ of Ser203 for achieving successful reactivation.
Collapse
|
4
|
Hester KP, Bhattarai K, Jiang H, Agarwal PK, Pope C. Engineering Dynamic Surface Peptide Networks on Butyrylcholinesterase G117H for Enhanced Organophosphosphorus Anticholinesterase Catalysis. Chem Res Toxicol 2019; 32:1801-1810. [PMID: 31411024 DOI: 10.1021/acs.chemrestox.9b00146] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The single residue mutation of butyrylcholinesterase (BChEG117H) hydrolyzes a number of organophosphosphorus (OP) anticholinesterases. Whereas other BChE active site/proximal mutations have been investigated, none are sufficiently active to be prophylactically useful. In a fundamentally different computer simulations driven strategy, we identified a surface peptide loop (residues 278-285) exhibiting dynamic motions during catalysis and modified it via residue insertions. We evaluated these loop mutants using computer simulations, substrate kinetics, resistance to inhibition, and enzyme reactivation assays using both the choline ester and OP substrates. A slight but significant increase in reactivation was noted with paraoxon with one of the mutants, and changes in KM and catalytic efficiency were noted in others. Simulations suggested weaker interactions between OP versus choline substrates and the active site of all engineered versions of the enzyme. The results indicate that an improvement of OP anticholinesterase hydrolysis through surface loop engineering may be a more effective strategy in an enzyme with higher intrinsic OP compound hydrolase activity.
Collapse
Affiliation(s)
- Kirstin P Hester
- Department of Physiological Sciences , Oklahoma State University , Stillwater , Oklahoma 74078 , United States
| | - Krishna Bhattarai
- Department of Entomology and Plant Pathology , Oklahoma State University , Stillwater , Oklahoma 74078 , United States
| | - Haobo Jiang
- Department of Entomology and Plant Pathology , Oklahoma State University , Stillwater , Oklahoma 74078 , United States
| | - Pratul K Agarwal
- Department of Biochemistry & Cellular and Molecular Biology , University of Tennessee , Knoxville , Tennessee 37996 , United States.,Arium BioLabs , 2519 Caspian Drive , Knoxville , Tennessee 37932 , United States
| | - Carey Pope
- Department of Physiological Sciences , Oklahoma State University , Stillwater , Oklahoma 74078 , United States
| |
Collapse
|
5
|
[Organophosphorus poisoning: Towards enzymatic treatments]. ANNALES PHARMACEUTIQUES FRANÇAISES 2019; 77:349-362. [PMID: 31253354 DOI: 10.1016/j.pharma.2019.06.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Revised: 06/05/2019] [Accepted: 06/05/2019] [Indexed: 11/22/2022]
Abstract
Organophosphorus compounds (OP) are toxic molecules developed as insecticides and chemical warfare nerve agents (CWNAs). Most OP are neurotoxic and act as nervous system disruptors by blocking cholinergic transmission. They are therefore responsible for many poisonings worldwide. OP toxicity may result either from acute or chronic exposure, and their poisoning effect were evaluated using several animal models. These latter were also used for evaluating the efficacy of antidotes. Strategies based on enzymes that can trap (stoichiometric bioscavengers) or degrade (catalytic bioscavengers) OP, were particularly studied since they allow effective decontamination, without toxicity or environmental impact. This review summarizes the results obtained in vivo with enzymes through three levels: prophylaxis, treatment and external decontamination. The efficiency of enzymatic treatments in different animal models is presented and the relevance of these models is also discussed for a better extrapolation to humans.
Collapse
|
6
|
Lushchekina S, Masson P. Catalytic bioscavengers against organophosphorus agents: mechanistic issues of self-reactivating cholinesterases. Toxicology 2018; 409:91-102. [DOI: 10.1016/j.tox.2018.07.020] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 07/23/2018] [Accepted: 07/25/2018] [Indexed: 12/21/2022]
|
7
|
Wan X, Yao Y, Fang L, Liu J. Unexpected protonation state of Glu197 discovered from simulations of tacrine in butyrylcholinesterase. Phys Chem Chem Phys 2018; 20:14938-14946. [PMID: 29786716 DOI: 10.1039/c8cp01566j] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Butyrylcholinesterase (BChE) has been actively involved in drug discoveries from many fields for decades. In the crystal structure of the BChE-tacrine complex, there is an unanticipated formyl-proline molecule resolved very close to tacrine, raising an essential question on how reliable it is to apply the binding pose in a crystal structure to analyze related experimental observations, in which no formyl-proline is actually involved. In this study, by performing a series of 100 ns molecular dynamics simulations, we demonstrate that it is safe to employ the structural information from this crystal structure to analyze related experimental observations. Surprisingly, Glu197 needs to be protonated to have the structures simulated appropriately. It should be noted that Glu197 has been commonly considered as deprotonated in diverse analyses due to its low pKa in aqueous solution, for which some interpretations are inconsistent or unclear. Our further investigation shows that the protonated Glu197 plays a very important role in preserving His438 within the catalytic triad through stabilizing a highly conserved water molecule. Interestingly, the catalytic triad and Glu197 have been long recognized for possibly deviating largely from the crystal structure, which might be catalytically deficient and is generally considered to result from the difference between the crystal and aqueous environment. Herein, our results suggest that the large deviations of the catalytic triad and Glu197 from the crystal structure are caused by the inappropriate protonation state of Glu197. This finding shall provide an important clue that has been long missing for a better understanding of BChE-related puzzles or even reconsideration of some BChE-catalyzed reaction mechanisms.
Collapse
Affiliation(s)
- Xiao Wan
- School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei 430030, 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
|
Lushchekina SV, Schopfer LM, Grigorenko BL, Nemukhin AV, Varfolomeev SD, Lockridge O, Masson P. Optimization of Cholinesterase-Based Catalytic Bioscavengers Against Organophosphorus Agents. Front Pharmacol 2018; 9:211. [PMID: 29593539 PMCID: PMC5859046 DOI: 10.3389/fphar.2018.00211] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2017] [Accepted: 02/26/2018] [Indexed: 11/13/2022] Open
Abstract
Organophosphorus agents (OPs) are irreversible inhibitors of acetylcholinesterase (AChE). OP poisoning causes major cholinergic syndrome. Current medical counter-measures mitigate the acute effects but have limited action against OP-induced brain damage. Bioscavengers are appealing alternative therapeutic approach because they neutralize OPs in bloodstream before they reach physiological targets. First generation bioscavengers are stoichiometric bioscavengers. However, stoichiometric neutralization requires administration of huge doses of enzyme. Second generation bioscavengers are catalytic bioscavengers capable of detoxifying OPs with a turnover. High bimolecular rate constants (kcat/Km > 106 M−1min−1) are required, so that low enzyme doses can be administered. Cholinesterases (ChE) are attractive candidates because OPs are hemi-substrates. Moderate OP hydrolase (OPase) activity has been observed for certain natural ChEs and for G117H-based human BChE mutants made by site-directed mutagenesis. However, before mutated ChEs can become operational catalytic bioscavengers their dephosphylation rate constant must be increased by several orders of magnitude. New strategies for converting ChEs into fast OPase are based either on combinational approaches or on computer redesign of enzyme. The keystone for rational conversion of ChEs into OPases is to understand the reaction mechanisms with OPs. In the present work we propose that efficient OP hydrolysis can be achieved by re-designing the configuration of enzyme active center residues and by creating specific routes for attack of water molecules and proton transfer. Four directions for nucleophilic attack of water on phosphorus atom were defined. Changes must lead to a novel enzyme, wherein OP hydrolysis wins over competing aging reactions. Kinetic, crystallographic, and computational data have been accumulated that describe mechanisms of reactions involving ChEs. From these studies, it appears that introducing new groups that create a stable H-bonded network susceptible to activate and orient water molecule, stabilize transition states (TS), and intermediates may determine whether dephosphylation is favored over aging. Mutations on key residues (L286, F329, F398) were considered. QM/MM calculations suggest that mutation L286H combined to other mutations favors water attack from apical position. However, the aging reaction is competing. Axial direction of water attack is not favorable to aging. QM/MM calculation shows that F329H+F398H-based multiple mutants display favorable energy barrier for fast reactivation without aging.
Collapse
Affiliation(s)
- Sofya V Lushchekina
- Laboratory of Computer Modeling of Bimolecular Systems and Nanomaterials, N. M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, Moscow, Russia
| | - Lawrence M Schopfer
- Department of Biochemistry and Molecular Biology, Eppley Institute, University of Nebraska Medical Center, Omaha, NE, United States
| | - Bella L Grigorenko
- Laboratory of Computer Modeling of Bimolecular Systems and Nanomaterials, N. M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, Moscow, Russia.,Chemistry Department, Lomonosov State University, Moscow, Russia
| | - Alexander V Nemukhin
- Laboratory of Computer Modeling of Bimolecular Systems and Nanomaterials, N. M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, Moscow, Russia.,Chemistry Department, Lomonosov State University, Moscow, Russia
| | - Sergei D Varfolomeev
- Laboratory of Computer Modeling of Bimolecular Systems and Nanomaterials, N. M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, Moscow, Russia.,Chemistry Department, Lomonosov State University, Moscow, Russia
| | - Oksana Lockridge
- Department of Biochemistry and Molecular Biology, Eppley Institute, University of Nebraska Medical Center, Omaha, NE, United States
| | - Patrick Masson
- Neuropharmacology Laboratory, Kazan Federal University, Kazan, Russia
| |
Collapse
|
10
|
Vistoli G, Mazzolari A, Testa B, Pedretti A. Binding Space Concept: A New Approach To Enhance the Reliability of Docking Scores and Its Application to Predicting Butyrylcholinesterase Hydrolytic Activity. J Chem Inf Model 2017. [PMID: 28633528 DOI: 10.1021/acs.jcim.7b00121] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Docking simulations are very popular approaches able to assess the capacity of a given ligand to interact with a target. Docking simulations are usually focused on a single best complex even though many studies showed that ligands retain a significant mobility within a binding pocket by assuming different binding modes all of which may contribute to the monitored ligand affinity. The present study describes an innovative concept, the binding space, which allows an exploration of the ligand mobility within the binding pocket by simultaneously considering several ligand poses as generated by docking simulations. The multiple poses and the relative docking scores can then be analyzed by taking advantage of the same concepts already used in the property space analysis; hence the binding space can be parametrized by (a) mean scores, (b) score ranges, and (c) score sensitivity values. The first parameter represents a very simple procedure to account for the contribution of the often neglected alternative binding modes, while the last two descriptors encode the degree of mobility which a given ligand retains within the binding cavity (score range) as well as the ease with which a ligand explores such a mobility (score sensitivity). Here, the binding space concept is applied to the prediction of the hydrolytic activity of BChE by synergistically considering multiple poses and multiple protein structures. The obtained results shed light on the remarkable potential of the binding space concept, whose parameters allow a significant increase of the predictive power of the docking results as revealed by extended correlative analyses. Mean scores are the parameters affording the largest statistical improvement, and all the here proposed docking-based descriptors show enhancing effects in developing predictive models. Finally, the study describes a new score function (Contacts score) simply based on the number of surrounding residues which appears to be particularly productive in the framework of the binding space.
Collapse
Affiliation(s)
- Giulio Vistoli
- Dipartimento di Scienze Farmaceutiche, Università degli Studi di Milano , Via Mangiagalli, 25, I-20133 Milano, Italy
| | - Angelica Mazzolari
- Dipartimento di Scienze Farmaceutiche, Università degli Studi di Milano , Via Mangiagalli, 25, I-20133 Milano, Italy
| | - Bernard Testa
- Department of Pharmacy, University Hospital Centre (CHUV) , Rue du Bugnon, CH-1011 Lausanne, Switzerland
| | - Alessandro Pedretti
- Dipartimento di Scienze Farmaceutiche, Università degli Studi di Milano , Via Mangiagalli, 25, I-20133 Milano, Italy
| |
Collapse
|
11
|
Wei D, Huang X, Qiao Y, Rao J, Wang L, Liao F, Zhan CG. Catalytic Mechanisms for Cofactor-Free Oxidase-Catalyzed Reactions: Reaction Pathways of Uricase-Catalyzed Oxidation and Hydration of Uric Acid. ACS Catal 2017; 7:4623-4636. [PMID: 28890842 DOI: 10.1021/acscatal.7b00901] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
First-principles quantum mechanical/molecular mechanical (QM/MM)-free energy calculations have been performed to uncover how uricase catalyzes metabolic reactions of uric acid (UA), demonstrating that the entire reaction process of UA in uricase consists of two stages-oxidation followed by hydration. The oxidation consists of four steps: (1) chemical transformation from 8-hydroxyxythine to an anionic radical via a proton transfer along with an electron transfer, which is different from the previously proposed electron-transfer mechanism that involves a dianion intermediate (UA2-) during the catalytic reaction process; (2) proton transfer to the O2- anion (radical); (3) diradical recombination to form a peroxo intermediate; (4) dissociation of H2O2 to generate the dehydrourate. Hydration, for the most favorable pathway, is initiated by the nucleophilic attack of a water molecule on dehydrourate, along with a concerted proton transfer through residue Thr69 in the catalytic site. According to the calculated free energy profile, the hydration is the rate-determining step, and the corresponding free energy barrier of 16.2 kcal/mol is consistent with that derived from experimental kinetic data, suggesting that the computational insights into the catalytic mechanisms are reasonable. The mechanistic insights not only provide a mechanistic base for future rational design of uricase mutants with improved catalytic activity against uric acid as an improved enzyme therapy, but also are valuable for understanding a variety of other cofactor-free oxidase-catalyzed reactions involving an oxygen molecule.
Collapse
Affiliation(s)
- Donghui Wei
- College
of Chemistry and Molecular Engineering, Zhengzhou University, 100 Science Avenue, Zhengzhou, Henan 450001, P. R. China
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States,
| | - Xiaoqin Huang
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States,
- Center
for Theoretical Biological Physics, and Center for Research Computing, Rice University, Houston, Texas 77030, United States,
| | - Yan Qiao
- College
of Chemistry and Molecular Engineering, Zhengzhou University, 100 Science Avenue, Zhengzhou, Henan 450001, P. R. China
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States,
| | - Jingjing Rao
- Key
Laboratory of Medical Laboratory Diagnostics of the Education Ministry,
College of Laboratory Medicine, Chongqing Medical University, No.1, Yixueyuan Road, Chongqing 400016, China
| | - Lu Wang
- Key
Laboratory of Medical Laboratory Diagnostics of the Education Ministry,
College of Laboratory Medicine, Chongqing Medical University, No.1, Yixueyuan Road, Chongqing 400016, China
| | - Fei Liao
- Key
Laboratory of Medical Laboratory Diagnostics of the Education Ministry,
College of Laboratory Medicine, Chongqing Medical University, No.1, Yixueyuan Road, Chongqing 400016, China
| | - Chang-Guo Zhan
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States,
- Molecular
Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| |
Collapse
|
12
|
Yun H, Yu J, Kim S, Lee N, Lee J, Lee S, Kim ND, Yu C, Rho J. Expression and purification of biologically active recombinant human paraoxonase 1 from a Drosophila S2 stable cell line. Protein Expr Purif 2017; 131:34-41. [DOI: 10.1016/j.pep.2016.11.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Revised: 11/01/2016] [Accepted: 11/08/2016] [Indexed: 01/03/2023]
|
13
|
Molecular polymorphism of human enzymes as the basis of individual sensitivity to drugs. Supercomputer-assisted modeling as a tool for analysis of structural changes and enzymatic activity of proteins. Russ Chem Bull 2017. [DOI: 10.1007/s11172-016-1487-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
14
|
Abou-Donia MB, Siracuse B, Gupta N, Sobel Sokol A. Sarin (GB, O-isopropyl methylphosphonofluoridate) neurotoxicity: critical review. Crit Rev Toxicol 2016; 46:845-875. [PMID: 27705071 PMCID: PMC5764759 DOI: 10.1080/10408444.2016.1220916] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Sarin (GB, O-isopropyl methylphosphonofluoridate) is a potent organophosphorus (OP) nerve agent that inhibits acetylcholinesterase (AChE) irreversibly. The subsequent build-up of acetylcholine (ACh) in the central nervous system (CNS) provokes seizures and, at sufficient doses, centrally-mediated respiratory arrest. Accumulation of ACh at peripheral autonomic synapses leads to peripheral signs of intoxication and overstimulation of the muscarinic and nicotinic receptors, which is described as "cholinergic crisis" (i.e. diarrhea, sweating, salivation, miosis, bronchoconstriction). Exposure to high doses of sarin can result in tremors, seizures, and hypothermia. More seriously, build-up of ACh at neuromuscular junctions also can cause paralysis and ultimately peripherally-mediated respiratory arrest which can lead to death via respiratory failure. In addition to its primary action on the cholinergic system, sarin possesses other indirect effects. These involve the activation of several neurotransmitters including gamma-amino-butyric acid (GABA) and the alteration of other signaling systems such as ion channels, cell adhesion molecules, and inflammatory regulators. Sarin exposure is associated with symptoms of organophosphate-induced delayed neurotoxicity (OPIDN) and organophosphate-induced chronic neurotoxicity (OPICN). Moreover, sarin has been involved in toxic and immunotoxic effects as well as organophosphate-induced endocrine disruption (OPIED). The standard treatment for sarin-like nerve agent exposure is post-exposure injection of atropine, a muscarinic receptor antagonist, accompanied by an oxime, an AChE reactivator, and diazepam.
Collapse
Affiliation(s)
- Mohamed B Abou-Donia
- a Department of Pharmacology and Cancer Biology , Duke University , Durham , NC , USA
| | - Briana Siracuse
- a Department of Pharmacology and Cancer Biology , Duke University , Durham , NC , USA
| | - Natasha Gupta
- a Department of Pharmacology and Cancer Biology , Duke University , Durham , NC , USA
| | - Ashly Sobel Sokol
- a Department of Pharmacology and Cancer Biology , Duke University , Durham , NC , USA
| |
Collapse
|
15
|
Masson P, Lushchekina SV. Emergence of catalytic bioscavengers against organophosphorus agents. Chem Biol Interact 2016; 259:319-326. [DOI: 10.1016/j.cbi.2016.02.010] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2015] [Revised: 12/16/2015] [Accepted: 02/10/2016] [Indexed: 02/05/2023]
|
16
|
An Y, Zhu Y, Yao Y, Liu J. Is it possible to reverse aged acetylcholinesterase inhibited by organophosphorus compounds? Insight from the theoretical study. Phys Chem Chem Phys 2016; 18:9838-46. [PMID: 27000635 DOI: 10.1039/c5cp07991h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The main treatment for organophosphorus (OP) compound poisoning in clinics is to restore the activity of acetylcholinesterase (AChE) through oxime-induced reactivation of the phosphorylated OP-AChE adduct. It suffers from a competitive and irreversible aging reaction of the phosphorylated OP-AChE adduct, resulting in permanent inactivity of AChE. However, it was recently reported that N-methyl-2-methoxypyridinium species can act as methylating agents to methylate the methyl methane-phosphonate monoanion, in which the reaction mimics the reverse of the aging reaction of the phosphorylated OP-AChE adduct. If the aging reaction could be really reversed, the efficiency for the OP detoxification should be significantly improved, bringing up the possibility to develop an agent to reverse the aging process of the phosphorylated OP-AChE adduct. However, such a reaction with the N-methyl-2-methoxypyridinium species in the enzyme is still not reported so far. It is of great interest to know whether or not this reaction is observable in the enzyme, and more importantly, if it turns out to be not observable in the enzyme, why such a reaction proceeds quickly in aqueous solution but not in the enzyme. In the present study, we performed DFT calculations and quantum mechanical/molecular mechanical (QM/MM) calculations to reveal the fundamental mechanism for the methylation of both the methyl methane-phosphonate monoanion and the aged sarin-AChE adduct by N-methyl-2-methoxypyridinium species, respectively. The obtained results support the SN2 reaction mechanism, not the stepwise mechanism, for the methylation of the methyl methane-phosphonate monoanion by 9 reported N-methyl-2-methoxypyridinium compounds. The calculated free energy barriers are in good agreement with the experimental data. The methylation of the aged sarin-AChE adduct by one N-methyl-2-methoxypyridinium compound (labeled as compound 2) also employs the SN2 reaction mechanism with an extremely high free energy barrier of 30.4 ± 3.5 (or 26.6) kcal mol(-1), implying that this reaction in the enzyme hardly occurs. Our results clearly show that compound 2 forms a strong π-π stacking interaction with the aromatic ring of the W86 residue of AChE, making itself unable to approach sarin for the reverse of the aging process. On the basis of the structure and mechanism, several possible strategies have been suggested for designing methylating agents with higher activity against the aged sarin-AChE adduct.
Collapse
Affiliation(s)
- Yun An
- Institute of Theoretical and Simulational Chemistry, Academy of Fundamental and Interdisciplinary Science, Harbin Institute of Technology, 2 Yikuang Street, Harbin, 150080, P. R. China.
| | | | | | | |
Collapse
|
17
|
Wei D, Tang M, Zhan CG. Fundamental reaction pathway and free energy profile of proteasome inhibition by syringolin A (SylA). Org Biomol Chem 2016; 13:6857-65. [PMID: 26018983 DOI: 10.1039/c5ob00737b] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, molecular dynamics (MD) simulations and first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations have been performed to uncover the fundamental reaction pathway of proteasome with a representative inhibitor syringolin A (SylA). The calculated results reveal that the reaction process consists of three steps. The first step is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is a nucleophilic attack on the olefin carbon of SylA by the negatively charged Thr1-O(γ) atom. The last step is a proton transfer from Thr1-N(z) to another olefin carbon of SylA to complete the inhibition reaction process. The calculated free energy profile demonstrates that the second step should be the rate-determining step and has the highest free energy barrier of 24.6 kcal mol(-1), which is reasonably close to the activation free energy (∼22.4-23.0 kcal mol(-1)) derived from the available experimental kinetic data. In addition, our computational results indicate that no water molecule can assist the rate-determining step, since the second step is not involved in a proton transfer process. The obtained mechanistic insights should be valuable for understanding the inhibition process of proteasome by SylA and structurally related inhibitors at a molecular level, and thus provide a solid mechanistic base and valuable clues for future rational design of novel, more potent inhibitors of proteasome.
Collapse
Affiliation(s)
- Donghui Wei
- Department of Chemistry, Zhengzhou University, Daxue Road, Zhengzhou, Henan 450052, China
| | | | | |
Collapse
|
18
|
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
|
19
|
Nemukhin AV, Kulakova AM, Lushchekina SV, Ermilov AY, Varfolomeev SD. Modeling chemical transformations at the active sites of cholinesterases by quantum-based simulations. ACTA ACUST UNITED AC 2016. [DOI: 10.3103/s0027131415060061] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
20
|
Kulakova A, Lushchekina S, Grigorenko B, Nemukhin A. Modeling reactivation of the phosphorylated human butyrylcholinesterase by QM(DFTB)/MM calculations. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s0219633615500510] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Human butyrylcholinesterase (BChE) is a bioscavenger that protects the enzyme which is critical for the central nerve system, acetylcholinesterase, from poisoning by organophosphorus agents. Elucidating the details of the hydrolysis reaction mechanism is important to understand how the phosphorylated BChE can be reactivated. Application of the QM(DFTB)/MM(AMBER) method to construct the minimum energy pathways for the hydrolysis reaction of the diethylphosphorylated BChE allowed us to suggest a mechanism of reactivation of the wild-type and the G117H mutated enzyme. Unlike previous approaches assuming that either His438 or His117 serves as a general base in the catalysis, in our proposal the Glu197 residue is responsible for activation of the nucleophilic water molecule (Wat) leading to the chemical transformations that restore the catalytic Ser198 residue in BChE. In agreement with the experimental data, it is shown that the G117H mutation facilitates the reactivation of the inhibited enzyme.
Collapse
Affiliation(s)
- Anna Kulakova
- Chemistry Department, M. V. Lomonosov Moscow State University, 1-3 Leninskiye Gory, Moscow 119991, Russia
| | - Sofya Lushchekina
- N. M. Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119334, Russia
| | - Bella Grigorenko
- Chemistry Department, M. V. Lomonosov Moscow State University, 1-3 Leninskiye Gory, Moscow 119991, Russia
- N. M. Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119334, Russia
| | - Alexander Nemukhin
- Chemistry Department, M. V. Lomonosov Moscow State University, 1-3 Leninskiye Gory, Moscow 119991, Russia
- N. M. Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, 4 Kosygin Street, Moscow 119334, Russia
| |
Collapse
|
21
|
Reilly PJ, Rovira C. Computational Studies of Glycoside, Carboxylic Ester, and Thioester Hydrolase Mechanisms: A Review. Ind Eng Chem Res 2015. [DOI: 10.1021/acs.iecr.5b01312] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Peter J. Reilly
- Department
of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011-2230, United States
| | - Carme Rovira
- Departament de Química Orgànica
and Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain
- Institució
Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| |
Collapse
|
22
|
Qiao Y, Han K, Zhan CG. Reaction pathways and free energy profiles for cholinesterase-catalyzed hydrolysis of 6-monoacetylmorphine. Org Biomol Chem 2014; 12:2214-27. [PMID: 24595354 DOI: 10.1039/c3ob42464b] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
As the most active metabolite of heroin, 6-monoacetylmorphine (6-MAM) can penetrate into the brain for the rapid onset of heroin effects. The primary enzymes responsible for the metabolism of 6-MAM to the less potent morphine in humans are acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The detailed reaction pathways for AChE- and BChE-catalyzed hydrolysis of 6-MAM to morphine have been explored, for the first time, in the present study by performing first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the two enzymatic reaction processes follow similar catalytic reaction mechanisms, and the whole catalytic reaction pathway for each enzyme consists of four reaction steps. According to the calculated results, the second reaction step associated with the transition state TS2(a)/TS2(b) should be rate-determining for the AChE/BChE-catalyzed hydrolysis, and the free energy barrier calculated for the AChE-catalyzed hydrolysis (18.3 kcal mol(-1)) is 2.5 kcal mol(-1) lower than that for the BChE-catalyzed hydrolysis (20.8 kcal mol(-1)). The free energy barriers calculated for the AChE- and BChE-catalyzed reactions are in good agreement with the experimentally derived activation free energies (17.5 and 20.7 kcal mol(-1) for the AChE- and BChE-catalyzed reactions, respectively). Further structural analysis reveals that the aromatic residues Phe295 and Phe297 in the acyl pocket of AChE (corresponding to Leu286 and Val288 in BChE) contribute to the lower energy of TS2(a) relative to TS2(b). The obtained structural and mechanistic insights could be valuable for use in future rational design of a novel therapeutic treatment of heroin abuse.
Collapse
Affiliation(s)
- Yan Qiao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Zhongshan Road 457, Dalian 116023, P. R. China
| | | | | |
Collapse
|
23
|
Legler PM, Boisvert SM, Compton JR, Millard CB. Development of organophosphate hydrolase activity in a bacterial homolog of human cholinesterase. Front Chem 2014; 2:46. [PMID: 25077141 PMCID: PMC4100338 DOI: 10.3389/fchem.2014.00046] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2014] [Accepted: 06/19/2014] [Indexed: 01/10/2023] Open
Abstract
We applied a combination of rational design and directed evolution (DE) to Bacillus subtilis p-nitrobenzyl esterase (pNBE) with the goal of enhancing organophosphorus acid anhydride hydrolase (OPAAH) activity. DE started with a designed variant, pNBE A107H, carrying a histidine homologous with human butyrylcholinesterase G117H to find complementary mutations that further enhance its OPAAH activity. Five sites were selected (G105, G106, A107, A190, and A400) within a 6.7 Å radius of the nucleophilic serine Oγ. All 95 variants were screened for esterase activity with a set of five substrates: pNP-acetate, pNP-butyrate, acetylthiocholine, butyrylthiocholine, or benzoylthiocholine. A microscale assay for OPAAH activity was developed for screening DE libraries. Reductions in esterase activity were generally concomitant with enhancements in OPAAH activity. One variant, A107K, showed an unexpected 7-fold increase in its k cat/K m for benzoylthiocholine, demonstrating that it is also possible to enhance the cholinesterase activity of pNBE. Moreover, DE resulted in at least three variants with modestly enhanced OPAAH activity compared to wild type pNBE. A107H/A190C showed a 50-fold increase in paraoxonase activity and underwent a slow time- and temperature-dependent change affecting the hydrolysis of OPAA and ester substrates. Structural analysis suggests that pNBE may represent a precursor leading to human cholinesterase and carboxylesterase 1 through extension of two vestigial specificity loops; a preliminary attempt to transfer the Ω-loop of BChE into pNBE is described. Unlike butyrylcholinesterase and pNBE, introducing a G143H mutation (equivalent to G117H) did not confer detectable OP hydrolase activity on human carboxylesterase 1 (hCE1). We discuss the use of pNBE as a surrogate scaffold for the mammalian esterases, and the importance of the oxyanion-hole residues for enhancing the OPAAH activity of selected serine hydrolases.
Collapse
Affiliation(s)
- Patricia M. Legler
- Naval Research Laboratory, Center for Bio/Molecular Science and EngineeringWashington, DC, USA
| | | | | | - Charles B. Millard
- United States Army Medical Research and Materiel CommandFort Detrick, MD, USA
| |
Collapse
|
24
|
Qiao Y, Han K, Zhan CG. Fundamental reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of heroin. Biochemistry 2013; 52:6467-79. [PMID: 23992153 DOI: 10.1021/bi400709v] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The pharmacological function of heroin requires an activation process that transforms heroin into 6-monoacetylmorphine (6-MAM), which is the most active form. The primary enzyme responsible for this activation process in human plasma is butyrylcholinesterase (BChE). The detailed reaction pathway of the activation process via BChE-catalyzed hydrolysis has been explored computationally, for the first time, in this study via molecular dynamics simulation and first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the whole reaction process includes acylation and deacylation stages. The acylation consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the 3-acetyl group of heroin by the hydroxyl oxygen of the Ser198 side chain and the dissociation of 6-MAM. The deacylation also consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the acyl-enzyme intermediate by a water molecule and the dissociation of the acetic acid from Ser198. The calculated free energy profile reveals that the second transition state (TS2) should be rate-determining. The structural analysis reveals that the oxyanion hole of BChE plays an important role in the stabilization of rate-determining TS2. The free energy barrier (15.9 ± 0.2 or 16.1 ± 0.2 kcal/mol) calculated for the rate-determining step is in good agreement with the experimentally derived activation free energy (~16.2 kcal/mol), suggesting that the mechanistic insights obtained from this computational study are reliable. The obtained structural and mechanistic insights could be valuable for use in the future rational design of a novel therapeutic treatment of heroin abuse.
Collapse
Affiliation(s)
- Yan Qiao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science , Zhongshan Road 457, Dalian 116023, P. R. China
| | | | | |
Collapse
|
25
|
Li D, Huang X, Lin J, Zhan CG. Catalytic mechanism of cytochrome P450 for N-methylhydroxylation of nicotine: reaction pathways and regioselectivity of the enzymatic nicotine oxidation. Dalton Trans 2013; 42:3812-20. [PMID: 23303461 DOI: 10.1039/c2dt32106h] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The fundamental reaction mechanism of cytochrome P450 2A6 (CYP2A6)-catalyzed N-methylhydroxylation of (S)-(-)-nicotine and the free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical (QM/MM) reaction-coordinate calculations. In the CYP2A6-(S)-(-)-nicotine binding structures that allow for 5'-hydroxylation, the N-methyl group is also sufficiently close to the oxygen of Cpd I for the N-methylhydroxylation reaction to occur. It has been demonstrated that the CYP2A6-catalyzed N-methylhydroxylation reaction is a concerted process involving a hydrogen-transfer transition state on both the quartet and the doublet states. The N-methylhydroxylation reaction proceeds mainly in the doublet state, since the free energy barriers on the doublet state are lower than the corresponding ones on the quartet state. The calculated free energy barriers indicate that (S)-(-)-nicotine oxidation catalyzed by CYP2A6 proceeds with a high regioselective abstraction of the hydrogen at the 5'-position, rather than the hydrogen at the N-methyl group. The predicted regioselectivity of 93% is in agreement with the most recent experimentally reported regioselectivity of 95%. The binding mode of (S)-(-)-nicotine in the active site of CYP2A6 is an important determinant for the stereoselectivity of nicotine (S)-(-)-oxidation, whereas the regioselectivity of (S)-(-)-nicotine oxidation is determined mainly by the free energy barrier difference between the 5'-hydroxylation and N-methylhydroxylation reactions.
Collapse
Affiliation(s)
- Dongmei Li
- College of Pharmacy, State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071, China.
| | | | | | | |
Collapse
|