1
|
Zhang Y, Mokkawes T, de Visser SP. Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides. Angew Chem Int Ed Engl 2023; 62:e202310785. [PMID: 37641517 DOI: 10.1002/anie.202310785] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Revised: 08/25/2023] [Accepted: 08/28/2023] [Indexed: 08/31/2023]
Abstract
Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.
Collapse
Affiliation(s)
- Yi Zhang
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M17DN, UK
- Department of Chemical Engineering, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Thirakorn Mokkawes
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M17DN, UK
- Department of Chemical Engineering, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Sam P de Visser
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M17DN, UK
- Department of Chemical Engineering, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| |
Collapse
|
2
|
Hermano Sampaio Dias A, Yadav R, Mokkawes T, Kumar A, Skaf MS, Sastri CV, Kumar D, de Visser SP. Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study. Inorg Chem 2023; 62:2244-2256. [PMID: 36651185 PMCID: PMC9923688 DOI: 10.1021/acs.inorgchem.2c03984] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C4-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.
Collapse
Affiliation(s)
- Artur Hermano Sampaio Dias
- Manchester
Institute of Biotechnology and Department of Chemical Engineering, The University of Manchester, 131 Princess Street, ManchesterM1 7DN, United Kingdom,Center
for Computing in Engineering & Sciences, University of Campinas, Rua Josué de Castro, s/n, Campinas13083-861, Brazil
| | - Rolly Yadav
- Department
of Chemistry, Indian Institute of Technology
Guwahati, Guwahati, Assam781039, India
| | - Thirakorn Mokkawes
- Manchester
Institute of Biotechnology and Department of Chemical Engineering, The University of Manchester, 131 Princess Street, ManchesterM1 7DN, United Kingdom
| | - Asheesh Kumar
- Department
of Physics, Babasaheb Bhimrao Ambedkar University, Lucknow, Uttar Pradesh (U.P.)226025, India
| | - Munir S. Skaf
- Center
for Computing in Engineering & Sciences, University of Campinas, Rua Josué de Castro, s/n, Campinas13083-861, Brazil
| | - Chivukula V. Sastri
- Department
of Chemistry, Indian Institute of Technology
Guwahati, Guwahati, Assam781039, India,
| | - Devesh Kumar
- Department
of Physics, Siddharth University, Kapilvastu, Siddharthnagar272202, India,
| | - Sam P. de Visser
- Manchester
Institute of Biotechnology and Department of Chemical Engineering, The University of Manchester, 131 Princess Street, ManchesterM1 7DN, United Kingdom,
| |
Collapse
|
3
|
Mokkawes T, Lim ZQ, de Visser SP. Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation? J Phys Chem B 2022; 126:9591-9606. [PMID: 36380557 PMCID: PMC9706573 DOI: 10.1021/acs.jpcb.2c07200] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C6-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C6-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.
Collapse
Affiliation(s)
- Thirakorn Mokkawes
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess
Street, Manchester M1 7DN, U.K.,Department
of Chemical Engineering, The University
of Manchester, Oxford
Road, Manchester M13 9PL, U.K.
| | - Ze Qing Lim
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess
Street, Manchester M1 7DN, U.K.,Department
of Chemical Engineering, The University
of Manchester, Oxford
Road, Manchester M13 9PL, U.K.
| | - Sam P. de Visser
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess
Street, Manchester M1 7DN, U.K.,Department
of Chemical Engineering, The University
of Manchester, Oxford
Road, Manchester M13 9PL, U.K.,
| |
Collapse
|
4
|
Peng W, Yan S, Zhang X, Liao L, Zhang J, Shaik S, Wang B. How Do Preorganized Electric Fields Function in Catalytic Cycles? The Case of the Enzyme Tyrosine Hydroxylase. J Am Chem Soc 2022; 144:20484-20494. [DOI: 10.1021/jacs.2c09263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Wei Peng
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| | - Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| | - Xuan Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| | - Langxing Liao
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| | - Jinyan Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| | - Sason Shaik
- Institute of Chemistry, The Hebrew University of Jerusalem, 9190407 Jerusalem, Israel
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005, People Republic of China
| |
Collapse
|