Mechanistic insights into the conversion of flavin adenine dinucleotide (FAD) to 8-formyl FAD in formate oxidase: a combined experimental and in-silico study

Abstract

Formate oxidase (FOx), which contains 8-formyl flavin adenine dinucleotide (FAD), exhibits a distinct advantage in utilizing ambient oxygen molecules for the oxidation of formic acid compared to other glucose-methanol-choline (GMC) oxidoreductase enzymes that contain only the standard FAD cofactor. The FOx-mediated conversion of FAD to 8-formyl FAD results in an approximate 10-fold increase in formate oxidase activity. However, the mechanistic details underlying the autocatalytic formation of 8-formyl FAD are still not well understood, which impedes further utilization of FOx. In this study, we employ molecular dynamics simulation, QM/MM umbrella sampling simulation, enzyme activity assay, site-directed mutagenesis, and spectroscopic analysis to elucidate the oxidation mechanism of FAD to 8-formyl FAD. Our results reveal that a catalytic water molecule, rather than any catalytic amino acids, serves as a general base to deprotonate the C8 methyl group on FAD, thus facilitating the formation of a quinone-methide tautomer intermediate. An oxygen molecule subsequently oxidizes this intermediate, resulting in a C8 methyl hydroperoxide anion that is protonated and dissociated to form OHC-RP and OH⁻. During the oxidation of FAD to 8-formyl FAD, the energy barrier for the rate-limiting step is calculated to be 22.8 kcal/mol, which corresponds to the required 14-hour transformation time observed experimentally. Further, the elucidated oxidation mechanism reveals that the autocatalytic formation of 8-formyl FAD depends on the proximal arginine and serine residues, R87 and S94, respectively. Enzymatic activity assay validates that the mutation of R87 to lysine reduces the k_cat value to 75% of the wild-type, while the mutation to histidine results in a complete loss of activity. Similarly, the mutant S94I also leads to the deactivation of enzyme. This dependency arises because the nucleophilic OH⁻ group and the quinone-methide tautomer intermediate are stabilized through the noncovalent interaction provided by R87 and S94. These findings not only explain the mechanistic details of each reaction step but also clarify the functional role of R87 and S94 during the oxidative maturation of 8-formyl FAD, thereby providing crucial theoretical support for the development of novel flavoenzymes with enhanced redox properties.

Keywords

Formate oxidase / 8-Formyl flavin adenine dinucleotide / Oxidative maturation / Molecular dynamics simulation / QM/MM umbrella sampling simulation

References

[1] Augustin P, Toplak M, Fuchs K, . Oxidation of the FAD cofactor to the 8-formyl-derivative in human electron-transferring flavoprotein. J Biol Chem, 2018, 293(8): 2829-2840.

[2] Berendsen HJC, Postma JPM, van Gunsteren WF, . Molecular dynamics with coupling to an external bath. J Chem Phys, 1984, 81(8): 3684-3690.

[3] Blay V, Pei D. Serine proteases: how did chemists tease out their catalytic mechanism?. Chem Texts, 2019, 5(4): 19.

[4] Darden T, York D, Pedersen L. Particle mesh Ewald: an N⋅log(N) method for Ewald sums in large systems. J Chem Phys, 1993, 98(12): 10089-10092.

[5] De Vivo M, Dal Peraro M, Klein ML. Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism. J Am Chem Soc, 2008, 130(33): 10955-10962.

[6] Doubayashi D, Ootake T, Maeda Y, . Formate oxidase, an enzyme of the glucose-methanol-choline oxidoreductase family, has a his-arg pair and 8-formyl-FAD at the catalytic site. Biosci Biotech Bioch, 2011, 75(9): 1662-1667.

[7] Doubayashi D, Oki M, Mikami B, . The microenvironment surrounding FAD mediates its conversion to 8-formyl-FAD in aspergillus oryzae RIB40 formate oxidase. J Biochem, 2019, 166(1): 67-75.

[8] Edmondson DE. Intramolecular hemiacetal formation in 8-formylriboflavine. Biochemistry, 1974, 13(14): 2817-2821.

[9] Elstner M, Porezag D, Jungnickel G, . Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys Rev B, 1998, 58(11): 7260-7268.

[10] Götz AW, Williamson MJ, Xu D, . Routine microsecond molecular dynamics simulations with AMBER on GPUs. 1. Generalized born. J Chem Theory Comput, 2012, 8(5): 1542-1555.

[11] Heath RS, Turner NJ. Recent advances in oxidase biocatalysts: enzyme discovery, cascade reactions and scale up. Curr Opin Green Sustain Chem, 2022, 38: 100693.

[12] Heuts DPHM, Scrutton NS, McIntire WS, . What’s in a covalent bond?. FEBS J, 2009, 276(13): 3405-3427.

[13] Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph, 1996, 14(1): 33-38.

[14] Kästner J, Senn HM, Thiel S, . QM/MM free-energy perturbation compared to thermodynamic integration and umbrella sampling: application to an enzymatic reaction. J Chem Theory Comput, 2006, 2(2): 452-461.

[15] Konjik V, Brünle S, Demmer U, . The crystal structure of RosB: insights into the reaction mechanism of the first member of a family of flavodoxin-like enzymes. Angew Chem Int Ed, 2017, 56(4): 1146-1151.

[16] Kräutler V, van Gunsteren WF, Hünenberger PH. A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations. J Comput Chem, 2001, 22(5): 501-508.

[17] Kumar S, Rosenberg JM, Bouzida D, . THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. J Comput Chem, 1992, 13(8): 1011-1021.

[18] Lao YE, Heyerdahl F, Jacobsen D, . An enzymatic assay with formate oxidase for point-of-care diagnosis of methanol poisoning. Basic Clin Pharmacol, 2022, 131(6): 547-554.

[19] Leys D, Scrutton NS. Sweating the assets of flavin cofactors: new insight of chemical versatility from knowledge of structure and mechanism. Curr Opin Struc Biol, 2016, 41: 19-26.

[20] Li R, Zhang S, Li Q, . An atom-economic enzymatic cascade catalysis for high-throughput RAFT synthesis of ultrahigh molecular weight polymers. Angew Chem Int Ed, 2022, 61(46): e202213396.

[21] Linke JA, Rayat A, Ward JM. Production of indigo by recombinant bacteria. Bioresour Bioprocess, 2023, 10(1): 20.

[22] Maeda Y, Doubayashi D, Oki M, . The presence of a modified FAD in formate oxidase from Debaryomyces Vanjiriae MH201 and aspergillus oryzae RIB40. J Biosci Bioeng, 2009, 108: S106.

[23] Maeda Y, Doubayashi D, Oki M, . Expression in Escherichia coli of an unnamed protein gene from aspergillus oryzae RIB40 and cofactor analyses of the gene product as formate oxidase. Biosci Biotech Bioch, 2009, 73(12): 2645-2649.

[24] Maeda Y, Doubayashi D, Ootake T, . Crystallization and preliminary X-ray analysis of formate oxidase, an enzyme of the glucose-methanol-choline oxidoreductase family. Acta Cryst F, 2010, 66(9): 1064-1066.

[25] Marcos-Alcalde I, Setoain J, Mendieta-Moreno JI, . MEPSA: minimum energy pathway analysis for energy landscapes. Bioinformatics, 2015, 31(23): 3853-3855.

[26] Mark P, Nilsson L. Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J Phys Chem A, 2001, 105(43): 9954-9960.

[27] Niehaus TA, Elstner M, Frauenheim T, . Application of an approximate density-functional method to sulfur containing compounds. J Mol Struc-Theochem, 2001, 541(1): 185-194.

[28] Özpınar GA, Peukert W, Clark T. An improved generalized AMBER force field (GAFF) for urea. J Mol Model, 2010, 16(9): 1427-1440.

[29] Robbins JM, Souffrant MG, Hamelberg D, . Enzyme-mediated conversion of flavin adenine dinucleotide (FAD) to 8-Formyl FAD in formate oxidase results in a modified cofactor with enhanced catalytic properties. Biochemistry, 2017, 56(29): 3800-3807.

[30] Robbins JM, Bommarius AS, Gadda G. Mechanistic studies of formate oxidase from Aspergillus oryzae: a novel member of the glucose-methanol-choline oxidoreductase enzyme superfamily that oxidizes carbon acids. Arch Biochem Biophys, 2018, 643: 24-31.

[31] Robbins JM, Geng J, Barry BA, . Photoirradiation generates an ultrastable 8-formyl FAD semiquinone radical with unusual properties in formate oxidase. Biochemistry, 2018, 57(40): 5818-5826.

[32] Roe DR, Cheatham TE III. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. J Chem Theory Comput, 2013, 9(7): 3084-3095.

[33] Rosta E, Nowotny M, Yang W, . Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations. J Am Chem Soc, 2011, 133(23): 8934-8941.

[34] Scott AP, Radom L. Harmonic vibrational frequencies: an evaluation of Hartree-Fock, Møller-Plesset, quadratic configuration interaction, density functional theory, and semiempirical scale factors. J Phys Chem, 1996, 100(41): 16502-16513.

[35] Seabra GM, Walker RC, Elstner M, . Implementation of the SCC-DFTB method for hybrid QM/MM simulations within the Amber molecular dynamics package. J Phys Chem A, 2007, 111(26): 5655-5664.

[36] Srinivasan B. A guide to the Michaelis-Menten equation: steady state and beyond. FEBS J, 2022, 289(20): 6086-6098.

[37] Sun BY, Sui HL, Liu ZW, . Structure-guided engineering of a flavin-containing monooxygenase for the efficient production of indirubin. Bioresour Bioprocess, 2022, 9(1): 70.

[38] Tao Y, Zhao Q, Liu F, . Enzymes encapsulated in organic-inorganic hybrid nanoflower with spatial localization for sensitive and colorimetric detection of formate. J Colloid Interface Sci, 2024, 672: 97-106.

[39] Tian C, Kasavajhala K, Belfon KAA, . ff19SB: amino-acid-specific protein backbone parameters trained against quantum mechanics energy surfaces in solution. J Chem Theory Comput, 2020, 16(1): 528-552.

[40] Willot SJP, Hoang MD, Paul CE, . FOx news: towards methanol-driven biocatalytic oxyfunctionalisation reactions. ChemCatChem, 2020, 12(10): 2713-2716.

[41] Wongnate T, Chaiyen P. The substrate oxidation mechanism of pyranose 2-oxidase and other related enzymes in the glucose-methanol-choline superfamily. FEBS J, 2013, 280(13): 3009-3027.

[42] Yorita K, Matsuoka T, Misaki H, . Interaction of two arginine residues in lactate oxidase with the enzyme flavin: Conversion of FMN to 8-formyl-FMN. Proc Natl Acad Sci, 2000, 97(24): 13039-13044.

Funding

National Natural Science Foundation of China(32201030)

Science and Technology Department of Jilin Province(20230402041GH)

Education Department of Jilin Province(JJKH20220970KJ)

Development and Reform Commission of Jilin Province(2023C015)

Fundamental Research Funds of the Central Universities in China (2024-JCXK-11)