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

Quantum chemical modeling of enantioselective sulfoxidation and epoxidation reactions by indole monooxygenase VpIndA1

Indole monooxygenases (IMOs) are enzymes from the family of Group E monooxygenases, requiring flavin adenine dinucleotide (FAD) for their activities. IMOs play important roles in both sulfoxidation and epoxidation reactions. The broad substrate range and high selectivity of IMOs make them promising biocatalytic tools for synthesizing chiral compounds. In the present study, quantum chemical calculations using the cluster approach were performed to investigate the reaction mechanism and the enantioselectivity of the IMO from Variovorax paradoxus EPS (VpIndA1). The sulfoxidation of methyl phenyl sulfide (MPS) and the epoxidation of indene were chosen as the representative reactions. The calculations confirmed that the FADOOH intermediate is the catalytic species in the VpIndA1 reactions. The oxidation of MPS adopts a one-step mechanism involving the direct oxygen-transfer from FADOOH to the substrate and the proton transfer from the -OH group back to FAD, while the oxidation of indene follows a stepwise mechanism involving a carbocation intermediate. It was computationally predicted that VpIndA1 prefers the formation of (S)-product for the MPS sulfoxidation and (1S,2R)-product for the indene epoxidation, consistent with the experimental observations. Importantly, the factors controlling the stereo-preference of the two reactions are identified. The findings in the present study provide valuable insights into the VpIndA1-catalyzed reactions, which are essential for the rational design of this enzyme and other IMOs for industrial applications. It is also worth emphasizing that the quantum chemical cluster approach is again demonstrated to be powerful in studying the enantioselectivity of enzymatic reactions.

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.

Mechanism of Orbital Interactions in the Sharpless Epoxidation with Ti(IV) Peroxides: A DFT Study

The M06-2X DFT functional has been employed to examine monomeric titanium(IV) hydroperoxo catalysts that model the individual steps in the dimeric titanium(IV)-catalyzed Sharpless reaction. This is the first example of a transition structure for titanium(IV) tert-butyl hydroperoxide-catalyzed epoxidation that describes the molecular motion required for oxygen atom transfer. These epoxidation catalysts have been examined for both bimolecular reactions with E-2-butene and the intramolecular epoxidation of allyl alcohol. The transition structure for the bimolecular peroxyacetic acid epoxidation of E-2-butene has been shown to be spiro in nature, and likewise, the intramolecular epoxidation of allyl alcohol is also nearly spiro. The significance of the O-C-C═C dihedral angle of allyl alcohol is examined for the Ti(IV) tert-butyl hydroperoxide epoxidation mechanism. Evidence is presented that supports a hexacoordinate titanium peroxo environment that exists in the dimeric form of the Sharpless catalyst. The mechanism for a 1,3-rearrangement of the alkoxide ligand in a titanium hydroperoxide to the Ti center in concert with oxygen atom transfer of the proximal oxygen to the C═C bond of the substrate is presented. The dimerization of Ti(IV)-(R,R)-diethyl tartrate-diisopropoxide and its hydrolysis have been calculated. The mechanism for rapid ligand exchange with alkyl hydroperoxides involving the Ti(O-i-Pr)₄ precursor is examined to show how the active epoxidation catalyst is produced.

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.

Pharmacophore-Based Virtual Screening of Novel Competitive Inhibitors of the Neurodegenerative Disease Target Kynurenine-3-Monooxygenase

The pathogenesis of several neurodegenerative diseases such as Alzheimer’s or Huntington’s disease has been associated with metabolic dysfunctions caused by imbalances in the brain and cerebral spinal fluid levels of neuroactive metabolites. Kynurenine monooxygenase (KMO) is considered an ideal therapeutic target for the regulation of neuroactive tryptophan metabolites. Despite significant efforts, the known KMO inhibitors lack blood–brain barrier (BBB) permeability and upon the mimicking of the substrate binding mode, are subject to produce reactive oxygen species as a side reaction. The computational drug design is further complicated by the absence of complete crystal structure information for human KMO (hKMO). In the current work, we performed virtual screening of readily available compounds using several protein–ligand complex pharmacophores. Each of the pharmacophores accounts for one of three distinct reported KMO protein-inhibitor binding conformations. As a result, six novel KMO inhibitors were discovered based on an in vitro fluorescence assay. Compounds VS1 and VS6 were predicted to be BBB permeable and avoid the hydrogen peroxide production dilemma, making them valuable, novel hit compounds for further drug property optimization and advancement in the drug design pipeline.

Importance of substituents in ring opening: a DFT study on a model reaction of thiazole to thioamide

Thiazole ring is an important active molecular skeleton of drugs. Thiazole in natural products and drugs are usually harmlessly eliminated. However, hepatotoxic reactions may occur due to the biological activation of thiazole to produce reactive thioamide. A typical example is hepatotoxic sudoxicam and safety meloxicam. The only structural difference between them is a methyl group on C5 position of thiazole in meloxicam. The molecular basis for the difference remains unknown and the bioactivation mechanism of the thiazole ring is still obscure. Quantum chemical calculations were performed to elucidate the activation mechanism of the thiazole ring under P450 catalysis, and the influence of the substituents on the activation pathways of thiazole ring was also studied. The calculated results show that the activation of thiazole is closely related to the substituents on the thiazole and spin state of Cpd I. The thiazole and substituted thiazole directly open the ring when catalyzed by doublet spin state Cpd I that catalyzed by the quartet spin state Cpd I can open the ring directly or indirectly, which is related to the substituents. Thiazoles modified with electron-donating substituents mainly undergo direct ring opening, while thiazoles modified with electron-withdrawing groups or weak electron-donating groups mainly undergo indirect ring-opening process accompanied by intermediate formation. The research results laid the foundation for the design of thiazole ring drugs, and also laid a theoretical foundation for the study of reducing the toxicity of thiazole ring drugs.

In silico methods predict new blood-brain barrier permeable structure for the inhibition of kynurenine 3-monooxygenase

Kynurenine 3-monooxygenase (KMO) regulates the levels of bioactive substances in the kynurenine pathway of tryptophan catabolism and its activity is tied to so many diseases that finding an appropriate inhibitor for KMO has become an urgent task. This especially proved to be difficult for the central nervous system related diseases due to the requirement that the supposed inhibitor should be both blood brain barrier permeable and should not cause hydrogen peroxide as a harmful side product. In this in silico study, we present our step-wise approach, whose starting point is based on the important experimental observations. To tackle the problem, a library of 7561938 structures was obtained from Zinc15 database utilizing the tranche browser. From this library, a subset of 501777 structures was determined with the considerations of their functional groups that constrain their applicability. Then, the binding affinity ranking of this set of structures was determined via virtual screening. Starting from the structures whose affinities are the highest among this subset, the ADMET properties were checked through in silico methods and the binding properties of the selected inhibitor candidates were further investigated via molecular dynamics simulations and MM/GBSA calculations. According to the computational results of this study, ZINC_71915355 has passed all the evaluations and is a potentially BBB permeable structure that can inhibit KMO. Additionally, ZINC_19827377 was identified as a new potential KMO inhibitor which may be more suitable for peripheral administration. From the in silico study presented herein, ZINC_71915355 and ZINC_19827377 structures, which showed high binding affinity without harmful H₂O₂ production, along with the tailored properties can now serve as powerful candidates for KMO inhibition and these hits are worth of further experimental validation.

Hydrogen movements in the oxidative half-reaction of kynurenine 3-monooxygenase from Pseudomonas fluorescens reveal the mechanism of hydroxylation

Kynurenine 3-monoxygenase (KMO) catalyzes the conversion of l-kynurenine (L-Kyn) to 3-hydroxykynurenine (3-OHKyn) in the pathway for tryptophan catabolism. We have investigated the effects of pH and deuterium substitution on the oxidative half-reaction of KMO from P. fluorescens (PfKMO). The three phases observed during the oxidative half reaction are formation of the hydroperoxyflavin, hydroxylation and product release. The measured rate constants for these phases proved largely unchanging with pH, suggesting that the KMO active site is insulated from exchange with solvent during catalysis. A solvent inventory study indicated that a solvent isotope effect of 2-3 is observed for the hydroxylation phase and that two or more protons are in flight during this step. An inverse isotope effect of 0.84 ± 0.01 on the rate constant for the hydroxylation step with ring perdeutero-L-Kyn as a substrate indicates a shift from sp² to sp³ hybridization in the transition state leading to the formation of a non-aromatic intermediate. The pH dependence of transient state data collected for the substrate analog meta-nitrobenzoylalanine indicate that groups proximal to the hydroperoxyflavin are titrated in the range pH 5-8.5 and can be described by a pKa of 8.8. That higher pH values do not slow the rate of hydroxylation precludes that the pKa measured pertains to the proton of the hydroperoxflavin. Together, these observations indicate that the C4a-hydroperoxyflavin has a pKa ≫ 8.5, that a non-aromatic species is the immediate product of hydroxylation and that at least two solvent derived protons are in-flight during oxygen insertion to the substrate aromatic ring. A unifying mechanistic proposal for these observations is proposed.