1
|
Jaladanki CK, Gahlawat A, Rathod G, Sandhu H, Jahan K, Bharatam PV. Mechanistic studies on the drug metabolism and toxicity originating from cytochromes P450. Drug Metab Rev 2020; 52:366-394. [DOI: 10.1080/03602532.2020.1765792] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Chaitanya K. Jaladanki
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| | - Anuj Gahlawat
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| | - Gajanan Rathod
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| | - Hardeep Sandhu
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| | - Kousar Jahan
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| | - Prasad V. Bharatam
- Department of Medicinal Chemistry and Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research (NIPER), SAS Nagar, Punjab, India
| |
Collapse
|
2
|
Neves Cruz J, Santana de Oliveira M, Gomes Silva S, Pedro da Silva Souza Filho A, Santiago Pereira D, Lima E Lima AH, de Aguiar Andrade EH. Insight into the Interaction Mechanism of Nicotine, NNK, and NNN with Cytochrome P450 2A13 Based on Molecular Dynamics Simulation. J Chem Inf Model 2019; 60:766-776. [PMID: 31622091 DOI: 10.1021/acs.jcim.9b00741] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Tobacco smoke contains various cancer-causing toxic substances, including nicotine and nitrosamines 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) and N'-nitrosonornicotine (NNN). The cytochrome 2A13 is involved in nicotine metabolism and in the activation of the pro-carcinogenic agents NNK and NNN, by means of α-hydroxylation reactions. Despite the significance of cytochrome 2A13 in the biotransformation of these molecules, its conformational mechanism and the molecular basis involved in the process are not fully understood. In this study, we used molecular dynamics and principal component analysis simulations for an in-depth analysis of the essential protein motions involved in the interaction of cytochrome 2A13 with its substrates. We also evaluated the interaction of these substrates with the amino acid residues in the binding pocket of cytochrome 2A13. Furthermore, we quantified the nature of these chemical interactions from free energy calculations using the Molecular Mechanics/Generalized Born Surface Area method. The ligands remained favorably oriented toward compound I (cytochrome P450 O═FeIV state), to undergo α-hydroxylation. The hydrogen bond with asparagine 297 was essential to maintaining the substrates in a favorable catalytic orientation. The plot of first principal motion vs second principal motion revealed that the enzyme's interaction with nicotine and NNK involved different conformational subgroups, whereas the conformational subgroups in the interaction with NNN are more similar. These results provide new mechanistic insights into the mode of interaction of the substrates with the active site of cytochrome 2A13, in the presence of compound I, which is essential for α-hydroxylation.
Collapse
Affiliation(s)
- Jorddy Neves Cruz
- Adolpho Ducke Laboratory , Emílio Goeldi Paraense Museum , Belém , Pará 66040-170, Brazil.,Laboratory of Agro-Industry , Embrapa Eastern Amazon , Belém , Pará 66040-170, Brazil
| | | | - Sebastião Gomes Silva
- Adolpho Ducke Laboratory , Emílio Goeldi Paraense Museum , Belém , Pará 66040-170, Brazil
| | | | | | - Anderson Henrique Lima E Lima
- Laboratory of Planning and Development of Pharmaceuticals , Federal University of Pará , Belém , Pará 70770-901, Brazil
| | | |
Collapse
|
3
|
Wei WJ, Qian HX, Wang WJ, Liao RZ. Computational Understanding of the Selectivities in Metalloenzymes. Front Chem 2018; 6:638. [PMID: 30622942 PMCID: PMC6308299 DOI: 10.3389/fchem.2018.00638] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Accepted: 12/07/2018] [Indexed: 01/26/2023] Open
Abstract
Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.
Collapse
Affiliation(s)
| | | | | | - Rong-Zhen Liao
- Key Laboratory of Material Chemistry for Energy Conversion and Storage, Ministry of Education, Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica, Hubei Key Laboratory of Materials Chemistry and Service Failure, School of Chemistry and Chemical Engineering, Huazhong University of Science and Technology, Wuhan, China
| |
Collapse
|
4
|
Theoretical Insights into Imidazolidine Oxidation of Imidacloprid by Cytochrome P450 3A4. J Mol Graph Model 2018; 80:173-181. [DOI: 10.1016/j.jmgm.2018.01.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Revised: 01/09/2018] [Accepted: 01/10/2018] [Indexed: 12/18/2022]
|
5
|
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
|
6
|
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
|
7
|
Olsen L, Oostenbrink C, Jørgensen FS. Prediction of cytochrome P450 mediated metabolism. Adv Drug Deliv Rev 2015; 86:61-71. [PMID: 25958010 DOI: 10.1016/j.addr.2015.04.020] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2014] [Revised: 03/30/2015] [Accepted: 04/27/2015] [Indexed: 10/23/2022]
Abstract
Cytochrome P450 enzymes (CYPs) form one of the most important enzyme families involved in the metabolism of xenobiotics. CYPs comprise many isoforms, which catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be formed. However, it is often hard to rationalize what metabolites these enzymes generate. In recent years, many different in silico approaches have been developed to predict binding or regioselective product formation for the different CYP isoforms. These comprise ligand-based methods that are trained on experimental CYP data and structure-based methods that consider how the substrate is oriented in the active site or/and how reactive the part of the substrate that is accessible to the heme group is. We will review key aspects for various approaches that are available to predict binding and site of metabolism (SOM), what outcome can be expected from the predictions, and how they could potentially be improved.
Collapse
|
8
|
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
|
9
|
Hirao H, Thellamurege N, Zhang X. Applications of density functional theory to iron-containing molecules of bioinorganic interest. Front Chem 2014; 2:14. [PMID: 24809043 PMCID: PMC4010748 DOI: 10.3389/fchem.2014.00014] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2013] [Accepted: 03/10/2014] [Indexed: 12/29/2022] Open
Abstract
The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.
Collapse
Affiliation(s)
- Hajime Hirao
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological UniversitySingapore, Singapore
| | | | | |
Collapse
|
10
|
Cheng CL, Zhang MZ, Zhao GJ. Mechanical stability and thermal conductivity of β-barrel in green fluorescent protein by steered molecular dynamics. RSC Adv 2014. [DOI: 10.1039/c3ra42679c] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
11
|
Rydberg P, Jørgensen FS, Olsen L. Use of density functional theory in drug metabolism studies. Expert Opin Drug Metab Toxicol 2013; 10:215-27. [PMID: 24295134 DOI: 10.1517/17425255.2014.864278] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
INTRODUCTION The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes. AREAS COVERED We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom. EXPERT OPINION DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.
Collapse
Affiliation(s)
- Patrik Rydberg
- University of Copenhagen, Department of Drug Design and Pharmacology , Denmark
| | | | | |
Collapse
|
12
|
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
|
13
|
Hirao H, Chuanprasit P, Cheong YY, Wang X. How is a metabolic intermediate formed in the mechanism-based inactivation of cytochrome P450 by using 1,1-dimethylhydrazine: hydrogen abstraction or nitrogen oxidation? Chemistry 2013; 19:7361-9. [PMID: 23592585 DOI: 10.1002/chem.201300689] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2013] [Indexed: 11/10/2022]
Abstract
A precise understanding of the mechanism-based inactivation of cytochrome P450 enzymes (P450s) at the quantum mechanical level should allow more reliable predictions of drug-drug interactions than those currently available. Hydrazines are among the molecules that act as mechanism-based inactivators to terminate the function of P450s, which are essential heme enzymes responsible for drug metabolism in the human body. Despite its importance, the mechanism explaining how a metabolic intermediate (MI) is formed from hydrazine is not fully understood. We used density functional theory (DFT) calculations to compare four possible mechanisms underlying the reaction between 1,1-dimethylhydrazine (or unsymmetrical dimethylhydrazine, UDMH) and the reactive compound I (Cpd I) intermediate of P450. Our DFT calculations provided a clear view on how an aminonitrene-type MI is formed from UDMH. In the most favorable pathway, hydrogen is spontaneously abstracted from the N2 atom of UDMH by Cpd I, followed by a second hydrogen abstraction from the N2 atom by Cpd II. Nitrogen oxidation of nitrogen atoms and hydrogen abstraction from the C-H bond of the methyl group were found to be less favorable than the hydrogen abstraction from the N-H bond. We also found that the reaction of protonated UDMH with Cpd I is rather sluggish. The aminonitrene-type MI binds to the ferric heme more strongly than a water molecule. This is consistent with the notion that the catalytic cycle of P450 is impeded when such an MI is produced through the P450-catalyzed reaction.
Collapse
Affiliation(s)
- Hajime Hirao
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371.
| | | | | | | |
Collapse
|