Catalytic Mechanism of Human Aldehyde Oxidase

DszA Catalyzes C-S Bond Cleavage through N5-Hydroperoxyl Formation

Due to its detrimental impact on human health and the environment, regulations demand ultralow sulfur levels on fossil fuels, in particular in diesel. However, current desulfurization techniques are expensive and cannot efficiently remove heteroaromatic sulfur compounds, which are abundant in crude oil and concentrate in the diesel fraction after distillation. Biodesulfurization via the four enzymes of the metabolic 4S pathway of the bacterium Rhodococcus erythropolis (DszA-D) is a possible solution. However, the 4S pathway needs to operate at least 500 times faster for industrial applicability, a goal currently pursued through enzyme engineering. In this work, we unveil the catalytic mechanism of the flavin monooxygenase DszA. Surprisingly, we found that this enzyme follows a recently proposed atypical mechanism that passes through the formation of an N5OOH intermediate at the re side of the cofactor, aided by a well-defined, predominantly hydrophobic O2 pocket. Besides clarifying the unusual chemical mechanism of the complex DszA enzyme, with obvious implications for understanding the puzzling chemistry of flavin-mediated catalysis, the result is crucial for the rational engineering of DszA, contributing to making biodesulfurization attractive for the oil refining industry.

Mechanism of AGT-Mediated Repair of dG-dC Cross-Links in the Drug Resistance to Chloroethylnitrosoureas: Molecular Docking, MD Simulation, and ONIOM (QM/MM) Investigation

Chloroethylnitrosoureas (CENUs) are important chemotherapies applied in the treatment of cancer. They exert anticancer activity by inducing DNA interstrand cross-links (ICLs) via the formation of two O6-alkylguanine intermediates, O6-chloroethylguanine (O6-ClEtG) and N1,O6-ethanoguanine (N1,O6-EtG). However, O6-alkylguanine-DNA alkyltransferase (AGT), a DNA-repair enzyme, can restore the O6-alkylguanine damages and thereby obstruct the formation of ICLs (dG-dC cross-link). In this study, the inhibitory mechanism of ICL formation was investigated to elucidate the drug resistance of CENUs mediated by AGT in detail. Based on the structures of the substrate-enzyme complexes obtained from docking and MD simulations, two ONIOM (QM/MM) models with different sizes of the QM region were constructed. The model with a larger QM region, which included the substrate (O6-ClEtG or N1,O6-EtG), a water molecule, and five residues (Tyr114, Cys145, His146, Lys165, and Glu172) in the active pocket of AGT, accurately described the repairing reaction and generated the results coinciding with the experimental outcomes. The repair process consists of two sequential steps: hydrogen transfer to form a thiolate anion on Cys145 and alkyl transfer from the O6 site of guanine (the rate-limiting step). The repair of N1,O6-EtG was more favorable than that of O6-ClEtG from both kinetics and thermodynamics aspects. Moreover, the comparison of the repairing process with the formation of dG-dC cross-link and the inhibition of AGT by O6-benzylguanine (O6-BG) showed that the presence of AGT could effectively interrupt the formation of ICLs leading to drug resistance, and the inhibition of AGT by O6-BG that was energetically more favorable than the repair of O6-ClEtG could not prevent the repair of N1,O6-EtG. Therefore, it is necessary to completely eliminate AGT activity before CENUs medication to enhance the chemotherapeutic effectiveness. This work provides reasonable explanations for the supposed mechanism of AGT-mediated drug resistance of CENUs and will assist in the development of novel CENU chemotherapies and their medication strategies.

Systems engineering of Escherichia coli for high-level glutarate production from glucose

Glutarate is a key monomer in polyester and polyamide production. The low efficiency of the current biosynthetic pathways hampers its production by microbial cell factories. Herein, through metabolic simulation, a lysine-overproducing E. coli strain Lys5 is engineered, achieving titer, yield, and productivity of 195.9 g/L, 0.67 g/g glucose, and 5.4 g/L·h, respectively. Subsequently, the pathway involving aromatic aldehyde synthase, monoamine oxidase, and aldehyde dehydrogenase (AMA pathway) is introduced into E. coli Lys5 to produce glutarate from glucose. To enhance the pathway's efficiency, rational mutagenesis on the aldehyde dehydrogenase is performed, resulting in the development of variant Mu5 with a 50-fold increase in catalytic efficiency. Finally, a glutarate tolerance gene cbpA is identified and genomically overexpressed to enhance glutarate productivity. With enzyme expression optimization, the glutarate titer, yield, and productivity of E. coli AMA06 reach 88.4 g/L, 0.42 g/g glucose, and 1.8 g/L·h, respectively. These findings hold implications for improving glutarate biosynthesis efficiency in microbial cell factories.

Transfer learning enhanced graph neural network for aldehyde oxidase metabolism prediction and its experimental application

Aldehyde oxidase (AOX) is a molybdoenzyme that is primarily expressed in the liver and is involved in the metabolism of drugs and other xenobiotics. AOX-mediated metabolism can result in unexpected outcomes, such as the production of toxic metabolites and high metabolic clearance, which can lead to the clinical failure of novel therapeutic agents. Computational models can assist medicinal chemists in rapidly evaluating the AOX metabolic risk of compounds during the early phases of drug discovery and provide valuable clues for manipulating AOX-mediated metabolism liability. In this study, we developed a novel graph neural network called AOMP for predicting AOX-mediated metabolism. AOMP integrated the tasks of metabolic substrate/non-substrate classification and metabolic site prediction, while utilizing transfer learning from 13C nuclear magnetic resonance data to enhance its performance on both tasks. AOMP significantly outperformed the benchmark methods in both cross-validation and external testing. Using AOMP, we systematically assessed the AOX-mediated metabolism of common fragments in kinase inhibitors and successfully identified four new scaffolds with AOX metabolism liability, which were validated through in vitro experiments. Furthermore, for the convenience of the community, we established the first online service for AOX metabolism prediction based on AOMP, which is freely available at https://aomp.alphama.com.cn.

Large Computational Survey of Intrinsic Reactivity of Aromatic Carbon Atoms with Respect to a Model Aldehyde Oxidase

Aldehyde oxidase (AOX) and other related molybdenum-containing enzymes are known to oxidize the C-H bonds of aromatic rings. This process contributes to the metabolism of pharmaceutical compounds and, therefore, is of vital importance to drug pharmacokinetics. The present work describes an automated computational workflow and its use for the prediction of intrinsic reactivity of small aromatic molecules toward a minimal model of the active site of AOX. The workflow is based on quantum chemical transition state searches for the underlying single-step oxidation reaction, where the automated protocol includes identification of unique aromatic C-H bonds, creation of three-dimensional reactant and product complex geometries via a templating approach, search for a transition state, and validation of reaction end points. Conformational search on the reactants, products, and the transition states is performed. The automated procedure has been validated on previously reported transition state barriers and was used to evaluate the intrinsic reactivity of nearly three hundred heterocycles commonly found in approved drug molecules. The intrinsic reactivity of more than 1000 individual aromatic carbon sites is reported. Stereochemical and conformational aspects of the oxidation reaction, which have not been discussed in previous studies, are shown to play important roles in accurate modeling of the oxidation reaction. Observations on structural trends that determine the reactivity are provided and rationalized.

Catalytic Reaction Mechanism of Bacterial GH92 α-1,2-Mannosidase: A QM/MM Metadynamics Study

The catalytic mechanism of aC a + 2 ${C{a}^{+2}}$ -dependent family 92 α ${{\rm \alpha }}$ -mannosidase, which is abundantly present in human gut flora and malfunctions leading to the lysosomal storage disease α-mannosidosis, has been investigated using quantum mechanics/molecular mechanics and metadynamics methods. Computational efforts show that the enzyme follows a conformational itinerary of and theC a + 2 ${C{a}^{+2}}$ ion serves a dual purpose, as it not only distorts the sugar ring but also plays a crucial role in orchestrating the arrangement of catalytic residues. This orchestration, in turn, contributes to the facilitation of O S 2 ${{{\rm \ }}^{{\rm O}}{{\rm S}}_{2}}$ conformers for the ensuing reaction. This mechanistic insight is well-aligned with the experimental predictions of the catalytic pathway, and the computed energies are of the same order of magnitude as the experimental estimations. Hence, our results extend the mechanistic understanding of glycosidases.

Does the Micronutrient Molybdenum Have a Role in Gestational Complications and Placental Health?

Molybdenum is an essential trace element for human health and survival, with molybdenum-containing enzymes catalysing multiple reactions in the metabolism of purines, aldehydes, and sulfur-containing amino acids. Recommended daily intakes vary globally, with molybdenum primarily sourced through the diet, and supplementation is not common. Although the benefits of molybdenum as an anti-diabetic and antioxidant inducer have been reported in the literature, there are conflicting data on the benefits of molybdenum for chronic diseases. Overexposure and deficiency can result in adverse health outcomes and mortality, although physiological doses remain largely unexplored in relation to human health. The lack of knowledge surrounding molybdenum intake and the role it plays in physiology is compounded during pregnancy. As pregnancy progresses, micronutrient demand increases, and diet is an established factor in programming gestational outcomes and maternal health. This review summarises the current literature concerning varied recommendations on molybdenum intake, the role of molybdenum and molybdoenzymes in physiology, and the contribution these play in gestational outcomes.

Functional Characterization and Protein Engineering of a Triterpene 3-/6-/2'-O-Glycosyltransferase Reveal a Conserved Residue Critical for the Regiospecificity

Engineering the function of triterpene glucosyltransferases (GTs) is challenging due to the large size of the sugar acceptors. In this work, we identified a multifunctional glycosyltransferase AmGT8 catalyzing triterpene 3-/6-/2'-O-glycosylation from the medicinal plant Astragalus membranaceus. To engineer its regiospecificity, a small mutant library was built based on semi-rational design. Variants A394F, A394D, and T131V were found to catalyze specific 6-O, 3-O, and 2'-O glycosylation, respectively. The origin of regioselectivity of AmGT8 and its A394F variant was studied by molecular dynamics and hydrogen deuterium exchange mass spectrometry. Residue 394 is highly conserved as A/G and is critical for the regiospecificity of the C- and O-GTs TcCGT1 and GuGT10/14. Finally, astragalosides III and IV were synthesized by mutants A394F, T131V and P192E. This work reports biocatalysts for saponin synthesis and gives new insights into protein engineering of regioselectivity in plant GTs.

Interrogating the Inhibition Mechanisms of Human Aldehyde Oxidase by X-ray Crystallography and NMR Spectroscopy: The Raloxifene Case

Human aldehyde oxidase (hAOX1) is mainly present in the liver and has an emerging role in drug metabolism, since it accepts a wide range of molecules as substrates and inhibitors. Herein, we employed an integrative approach by combining NMR, X-ray crystallography, and enzyme inhibition kinetics to understand the inhibition modes of three hAOX1 inhibitors-thioridazine, benzamidine, and raloxifene. These integrative data indicate that thioridazine is a noncompetitive inhibitor, while benzamidine presents a mixed type of inhibition. Additionally, we describe the first crystal structure of hAOX1 in complex with raloxifene. Raloxifene binds tightly at the entrance of the substrate tunnel, stabilizing the flexible entrance gates and elucidating an unusual substrate-dependent mechanism of inhibition with potential impact on drug-drug interactions. This study can be considered as a proof-of-concept for an efficient experimental screening of prospective substrates and inhibitors of hAOX1 relevant in drug discovery.

THE CATALYTIC MECHANISM OF THE RETAINING GLYCOSYLTRANSFERASE MANNOSYLGLYCERATE SYNTHASE

To protect their intracellular proteins, extremophile microorganisms synthesize molecules called compatible solutes. These molecules are the result of the attachment of a small negatively charged molecule to a sugar molecule. It has been found that these molecules, not only protect the microorganism against osmotic stress, as initially thought, but also against other extreme conditions. The observation that these molecules can confer protection against extreme conditions to isolated enzymes from different organisms made them an exciting prospect for potential biotechnological applications. One of the most widespread compatible solute in hyperthermophile organisms is the molecule 2-O-α-D-mannosyl-D-glycerate (MG). In addition to confer protection to proteins against extreme conditions, MG was found to prevent Alzheimer's β-amyloid aggregation and reduce α-synuclein fibril formation in Parkinson's disease.  In this work we studied, using computational methods, the catalytic mechanism of the synthesis of MG by the enzyme mannosylglycerate synthase (MGS) from the thermophilic bacteria  Rhodothermus marinus . MGS is a promiscuous enzyme, accepting a variety of sugar donors and acceptors. This feature can be used to synthesize other molecules with potential biotechnological applications beyond MG. The unravelling of the catalytic mechanism with atomistic resolution and the associated free energies and electrostatic profiles of the stationary states obtained in the present work will help future investigations to full explore the potential of MGS.

Structural, Mechanistic, and Functional Insights into an Arthrobacter nicotinovorans Molybdenum Hydroxylase Involved in Nicotine Degradation

Arthrobacter nicotinovorans decomposes nicotine through the pyridine pathway. 6-hydroxypseudooxynicotine 2-oxidoreductase (also named ketone dehydrogenase, Kdh) is an important enzyme in nicotine degradation pathway of A. nicotinovorans, and is responsible for the second hydroxylation of nicotine. Kdh belongs to the molybdenum hydroxylase family, and catalyzes the oxidation of 6-hydroxy-pseudooxynicotine (6-HPON) to 2,6-dihydroxy-pseudooxynicotine (2,6-DHPON). We determined the crystal structure of the Kdh holoenzyme from A. nicotinovorans, with its three subunits KdhL, KdhM, and KdhS, and their associated cofactors molybdopterin cytosine dinucleotide (MCD), two iron-sulfur clusters (Fe₂S₂), and flavin adenine dinucleotide (FAD), respectively. In addition, we obtained a structural model of the substrate 6-HPON-bound Kdh through molecular docking, and performed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations to unveil the catalytic mechanism of Kdh. The residues Glu345, Try551, and Glu748 of KdhL were found to participate in substrate binding, and Phe269 and Arg383 of KdhL were found to contribute to stabilize the MCD conformation. Furthermore, site-directed mutagenesis and enzymatic activity assays were performed to support our structural and computational results, which also revealed a trend of increasing catalytic efficiency with the increase in the buffer pH. Lastly, our electrochemical results demonstrated electron transfer among the various cofactors of Kdh. Therefore, our work provides a comprehensive structural, mechanistic, and functional study on the molybdenum hydroxylase Kdh in the nicotine degradation pathway of A. nicotinovorans.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

A comprehensive review of in silico approaches for the prediction and modulation of aldehyde oxidase-mediated drug metabolism: The current features, challenges and future perspectives

The importance of aldehyde oxidase (AOX) in drug metabolism necessitates the development and application of the in silico rational drug design methods as an integral part of drug discovery projects for the early prediction and modulation of AOX-mediated metabolism. The current study represents an up-to-date and thorough review of in silico studies of AOX-mediated metabolism and modulation methods. In addition, the challenges and the knowledge gap that should be covered have been discussed. The importance of aldehyde oxidase (AOX) in drug metabolism is a hot topic in drug discovery. Different strategies are available for the modulation of the AOX-mediated metabolism of drugs. Application of the rational drug design methods as an integral part of drug discovery projects is necessary for the early prediction of AOX-mediated metabolism. The current study represents a comprehensive review of AOX molecular structure, AOX-mediated reactions, AOX substrates, AOX inhibition, approaches to modify AOX-mediated metabolism, prediction of AOX metabolism/substrates/inhibitors, and the AOX related structure-activity relationship (SAR) studies. Furthermore, an up-to-date and thorough review of in silico studies of AOX metabolism has been carried out. In addition, the challenges and the knowledge gap that should be covered in the scientific literature have been discussed in the current review.

The complete catalytic mechanism of xanthine oxidase: a computational study

In this article, quantum mechanical/molecular mechanical (QM/MM) methods were used to study the full catalytic mechanism of xanthine oxidase (XO).

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.