Role of the Somersault Rearrangement in the Oxidation Step for Flavin Monooxygenases (FMO). A Comparison between FMO and Conventional Xenobiotic Oxidation with Hydroperoxides

A multifunctional flavoprotein monooxygenase HspB for hydroxylation and C-C cleavage of 6-hydroxy-3-succinoyl-pyridine

Flavoprotein monooxygenases catalyze reactions, including hydroxylation and epoxidation, involved in the catabolism, detoxification, and biosynthesis of natural substrates and industrial contaminants. Among them, the 6-hydroxy-3-succinoyl-pyridine (HSP) monooxygenase (HspB) from Pseudomonas putida S16 facilitates the hydroxylation and C-C bond cleavage of the pyridine ring in nicotine. However, the mechanism for biodegradation remains elusive. Here, we refined the crystal structure of HspB and elucidated the detailed mechanism behind the oxidative hydroxylation and C-C cleavage processes. Leveraging structural information about domains for binding the cofactor flavin adenine dinucleotide (FAD) and HSP substrate, we used molecular dynamics simulations and quantum/molecular mechanics calculations to demonstrate that the transfer of an oxygen atom from the reactive FAD peroxide species (C4a-hydroperoxyflavin) to the C3 atom in the HSP substrate constitutes a rate-limiting step, with a calculated reaction barrier of about 20 kcal/mol. Subsequently, the hydrogen atom was rebounded to the FAD cofactor, forming C4a-hydroxyflavin. The residue Cys218 then catalyzed the subsequent hydrolytic process of C-C cleavage. Our findings contribute to a deeper understanding of the versatile functions of flavoproteins in the natural transformation of pyridine and HspB in nicotine degradation.IMPORTANCEPseudomonas putida S16 plays a pivotal role in degrading nicotine, a toxic pyridine derivative that poses significant environmental challenges. This study highlights a key enzyme, HspB (6-hydroxy-3-succinoyl-pyridine monooxygenase), in breaking down nicotine through the pyrrolidine pathway. Utilizing dioxygen and a flavin adenine dinucleotide cofactor, HspB hydroxylates and cleaves the substrate's side chain. Structural analysis of the refined HspB crystal structure, combined with state-of-the-art computations, reveals its distinctive mechanism. The crucial function of Cys218 was never discovered in its homologous enzymes. Our findings not only deepen our understanding of bacterial nicotine degradation but also open avenues for applications in both environmental cleanup and pharmaceutical development.

Environmental biotransformation mechanisms by flavin-dependent monooxygenase: A computational study

The enzyme-catalyzed metabolic biotransformation of xenobiotics plays a significant role in toxicology evolution and subsequently environmental health risk assessment. Recent studies noted that the phase I human flavin-dependent monooxygenase (e.g., FMO3) can catalyze xenobiotics into more toxic metabolites. However, details of the metabolic mechanisms are insufficient. To fill the mechanism in the gaps, the systemic density functional theory calculations were performed to elucidate diverse FMO-catalyzed oxidation reactions toward environmental pollutants, including denitrification (e.g., nitrophenol), N-oxidation (e.g., nicotine), desulfurization (e.g., fonofos), and dehalogenation (e.g., pentachlorophenol). Similar to the active center compound 0 of cytochrome P450, FMO mainly catalyzed reactions with the structure of the tricyclic isoalloxazine C-4a-hydroperoxide (FADHOOH). As will be shown, FMO-catalyzed pathways are more favorable with a concerted than stepwise mechanism; Deprotonation is necessary to initiate the oxidation reactions for phenolic substrates; The regioselectivity of nicotine by FMO prefers the N-oxidation other than N-demethylation pathway; Formation of the P-S-O triangle ring is the key step for desulfurization of fonofos by FMO. We envision that these fundamental mechanisms catalyzed by FMO with a computational method can be extended to other xenobiotics of similar structures, which may aid the high-throughput screening and provide theoretical predictions in the future.

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

Computational Survey of Recent Experimental Developments in the Hydroxylation Mechanism of Kynurenine 3-Monooxygenase

Recently, two new mechanistic proposals for the kynurenine 3-monooxygenase (KMO) catalyzed hydroxylation reaction of l-Kynurenine (l-Kyn) have been proposed. According to the first proposal, instead of the distal oxygen, the proximal oxygen of the hydroperoxide intermediate of flavin adenine dinucleotide (FAD) is transferred to the substrate ring. The second study proposes that l-Kyn participates in its base form in the reaction. To address these proposals, the reaction was reconsidered with a 386 atom quantum cluster model that is based on a recent X-ray structure (PDB id: 6FOX). The computations were carried out at the UB3LYP/6-311+G(2d,2p)//UB3LYP/6-31G(d,p) level with solvation (polarizable continuum model) and dispersion (DFT-D3(BJ)) corrections. To supplement the results of the density functional theory (DFT) calculations, molecular dynamics (MD) simulations of the protein-substrate complex were employed. The comparison of a proximal oxygen transfer mechanism to the distal oxygen transfer mechanism revealed that the former requires too high of a barrier energy while the latter validated our previous results. According to the MD simulations, the hydroperoxy moiety does not favor an alignment that might promote the proximal oxygen transfer mechanism. In the second part of the study, hydroxylation reaction with the base form of l-Kyn was sought. Although DFT calculations confirmed a much more facile reaction with the base form of l-Kyn, a mechanism which would allow the deprotonation of the l-Kyn before the oxygen transfer could not be determined with the X-ray-based positions. A concerted mechanism with both the oxygen transfer and the deprotonation required a high barrier energy. A stepwise mechanism involving the deprotonation of l-Kyn was found, starting from an MD frame. The overall barrier of the oxygen transfer step of this model was found to be in the range of that of with neutral l-Kyn. MD simulations supported the idea of ineffectiveness of the nearby shell surrounding the utilized active site core on the deprotonation of l-Kyn.

New frontiers in flavin-dependent monooxygenases

Flavin-dependent monooxygenases catalyze a wide variety of redox reactions in important biological processes and are responsible for the synthesis of highly complex natural products. Although much has been learned about FMO chemistry in the last ~80 years of research, several aspects of the reactions catalyzed by these enzymes remain unknown. In this review, we summarize recent advancements in the flavin-dependent monooxygenase field including aspects of flavin dynamics, formation and stabilization of reactive species, and the hydroxylation mechanism. Novel catalysis of flavin-dependent N-oxidases involving consecutive oxidations of amines to generate oximes or nitrones is presented and the biological relevance of the products is discussed. In addition, the activity of some FMOs have been shown to be essential for the virulence of several human pathogens. We also discuss the biomedical relevance of FMOs in antibiotic resistance and the efforts to identify inhibitors against some members of this important and growing family enzymes.

Predicting reactivity to drug metabolism: beyond P450s-modelling FMOs and UGTs

We present a study based on density functional theory calculations to explore the rate limiting steps of product formation for oxidation by Flavin-containing Monooxygenase (FMO) and glucuronidation by the UDP-glucuronosyltransferase (UGT) family of enzymes. FMOs are responsible for the modification phase of metabolism of a wide diversity of drugs, working in conjunction with Cytochrome P450 (CYP) family of enzymes, and UGTs are the most important class of drug conjugation enzymes. Reactivity calculations are important for prediction of metabolism by CYPs and reactivity alone explains around 70-85% of the experimentally observed sites of metabolism within CYP substrates. In the current work we extend this approach to propose model systems which can be used to calculate the activation energies, i.e. reactivity, for the rate-limiting steps for both FMO oxidation and glucuronidation of potential sites of metabolism. These results are validated by comparison with the experimentally observed reaction rates and sites of metabolism, indicating that the presented models are suitable to provide the basis of a reactivity component within generalizable models to predict either FMO or UGT metabolism.

Mechanism of Kynurenine 3-Monooxygenase-Catalyzed Hydroxylation Reaction: A Quantum Cluster Approach

The mechanism of the hydroxylation reaction between l-Kyn and model flavin adenine dinucleotide (FAD)-hydroperoxide was investigated via density functional theory (DFT) calculations in the absence and in the presence of the kynurenine 3-monooxygenase (KMO) enzyme by considering possible pathways that can lead to the product 3-hydroxykynurenine (3-HK). Crystal structure (pdb code: 5NAK )-based calculations involved a quantum cluster model in which the active site of the enzyme with the substrate l-Kyn was represented with 348 atoms. According to the deduced mechanism, KMO-catalyzed hydroxylation reaction takes place with four transformations. In the initial transition state, FAD delivers its peroxy hydroxyl to the l-Kyn ring, creating an sp³-hybridized carbon center. Then, the hydrogen on the hydroxyl moiety is immediately transferred back to the proximal oxygen that remained on FAD. These consequent transformations are in line with the somersault rearrangement previously described for similar enzymatic systems. The second step corresponds to a hydride shift from the sp³-hybridized carbon of the substrate ring to its adjacent carbon, producing the keto form of 3-HK. Then, keto-3-HK is transformed into its enol form (3-HK) with a water-assisted tautomerization. Lastly, FAD is oxidized with a water-assisted dehydration, which also involves 3-HK as a catalyst. In the proposed pathway, Asn54, Pro318, and a crystal water molecule were seen to play significant roles in the proton relays. The energies obtained via the cluster approach were calculated at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d,p) level with solvation (polarizable continuum model) and dispersion (DFT-D3(BJ)) corrections.

Pathways for Arene Oxidation in Non-Heme Diiron Enzymes: Lessons from Computational Studies on Benzoyl Coenzyme A Epoxidase

Oxygenation of aromatic rings using O₂ is catalyzed by several non-heme carboxylate-bridged diiron enzymes. In order to provide a general mechanistic description for these reactions, computational studies were carried out at the ONIOM(B3LYP/BP86/Amber) level on the non-heme diiron enzyme benzoyl coenzyme A epoxidase, BoxB. The calculations revealed four possible pathways for attacking the aromatic ring: (a) electrophilic (2e^-) attack by a bis(μ-oxo)-diiron(IV) species (Q pathway); (b) electrophilic (2e^-) attack via the σ* orbital of a μ-η²:η²-peroxo-diiron(III) intermediate (Pσ* pathway); (c) radical (1e^-) attack via the π*-orbital of a superoxo-diiron(II,III) species (Pπ* pathway); (d) radical (1e^-) attack of a partially quenched bis(μ-oxo)-diiron(IV) intermediate (Q' pathway). The results allowed earlier work of de Visser on olefin epoxidation by diiron complexes and QM-cluster studies of Liao and Siegbahn on BoxB to be put into a broader perspective. Parallels with epoxidation using organic peracids were also examined. Specifically for the BoxB enzyme, the Q pathway was found to be the most preferred, but the corresponding bis(μ-oxo)-diiron(IV) species is significantly destabilized and not expected to be directly observable. Epoxidation via the Pσ* pathway represents an energetically somewhat higher lying alternative; possible strategies for experimental discrimination are discussed. The selectivity toward epoxidation is shown to stem from a combination of inherent electronic properties of the thioacyl substituent and enzymatic constraints. Possible implications of the results for toluene monooxygenases are considered as well.

Insights into the catalytic mechanism of chlorophenol 4-monooxygenase: a quantum mechanics/molecular mechanics study

The degradation mechanism of chlorophenol 4-monooxygenase toward pollutants 2,4,5-trichlorophenol, 2,4,6-trichlorophenol, and 2,5-dichloro-p-hydroquinone was studied by QM/MM investigations.

Use of density functional theory in drug metabolism studies

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.