Search

Flavin mono-nucleotides (FMNs) are cofactors that hold responsibility for carrying and transferring electrons in the electron transport chain stage of cellular respiration. Without being facilitated by FMNs, energy production is stagnant due to the interruption in most of the cellular processes. Investigation on FMN's functions, therefore, can gain holistic understanding about human diseases and molecular information on drug targets. We proposed a deep learning model using a two-dimensional convolutional neural network and position specific scoring matrices that could identify FMN interacting residues with the sensitivity of 83.7 percent, specificity of 99.2 percent, accuracy of 98.2 percent, and Matthews correlation coefficients of 0.85 for an independent dataset containing 141 FMN binding sites and 1,920 non-FMN binding sites. The proposed method outperformed other previous studies using similar evaluation metrics. Our positive outcome can also promote the utilization of deep learning in dealing with various problems in bioinformatics and computational biology.

Proteins like NADH:ubiquinone oxidoreductase (NQR), an essential enzyme and ion pump in the physiology of several pathogenic bacteria, tightly regulate the redox properties of their cofactors. Although flavin mononucleotide (FMN) is fully reduced in aqueous solution, FMN in subunits B and C of NQR exclusively undergo one-electron transitions during its catalytic cycle. Here, we perform ab initio calculations and molecular dynamics simulations to elucidate the mechanisms that regulate the redox state of FMN in NQR. QM/MM calculations show that binding site electrostatics disfavor anionic forms of FMNH(2), but permit a neutral form of the fully reduced flavin. The potential energy surface is unaffected by covalent bonding between FMN and threonine. Molecular dynamics simulations show that the FMN binding sites are inaccessible by water, suggesting that further reductions of the cofactors are limited or prohibited by the availability of water and other proton donors. These findings provide a deeper understanding of the mechanisms used by NQR to regulate electron transfer through the cofactors and perform its physiologic role. They also provide the first, to our knowledge, evidence of the simple concept that proteins regulate flavin redox states via water occlusion.

Dynamics of protein structures of five FMN binding isomers, WT, E13K, E13R, E13T and E13Q, were studied through time-dependent residue electrostatic energies (RESTs) of ionic residues and their correlation coefficients. The Arg86 and Glu119 have high positive correlation coefficients which are higher than 0.7 in most of FBP isomers. The Glu23 and FMN, Asp24 and FMN, and Asp61 and Arg63 pairs displayed relatively high negative coefficients. In these residues, the distributions of RESTs showed a few peaks. REST pairs have high positive coefficients when distribution of REST of one residue overlapped with that of REST of another residue, while the REST pairs displayed high negative coefficients when a peak distribution of REST of one residue overlapped with a minimum distribution of another residue. Correlations of REST pairs were not simple and displayed a few groups. Relationships were examined among the correlation coefficients, inter-residue distance (R) and signs of the charge. In the case of R ≤ 1 nm, the total number of the pairs with high negative coefficients in all FBP isomers was much greater when the pairs possess opposite signs of the charge, compared to that when the pairs possess same signs of the charge.

The nitroreductase superfamily of enzymes encompasses many flavin mononucleotide (FMN)-dependent catalysts promoting a wide range of reactions. All share a common core consisting of an FMN-binding domain, and individual subgroups additionally contain one to three sequence extensions radiating from defined positions within this core to support their unique catalytic properties. To identify the minimum structure required for activity in the iodotyrosine deiodinase subgroup of this superfamily, attention was directed to a representative from the thermophilic organism Thermotoga neapolitana (TnIYD). This representative was selected based on its status as an outlier of the subgroup arising from its deficiency in certain standard motifs evident in all homologues from mesophiles. We found that TnIYD lacked a typical N-terminal sequence and one of its two characteristic sequence extensions, neither of which was found to be necessary for activity. We also show that TnIYD efficiently promotes dehalogenation of iodo-, bromo-, and chlorotyrosine, analogous to related deiodinases (IYDs) from humans and other mesophiles. In addition, 2-iodophenol is a weak substrate for TnIYD as it was for all other IYDs characterized to date. Consistent with enzymes from thermophilic organisms, we observed that TnIYD adopts a compact fold and low surface area compared with IYDs from mesophilic organisms. The insights gained from our investigations on TnIYD demonstrate the advantages of focusing on sequences that diverge from conventional standards to uncover the minimum essentials for activity. We conclude that TnIYD now represents a superior starting structure for future efforts to engineer a stable dehalogenase targeting halophenols of environmental concern.

The biosynthesis of thioviridamide-like compounds has not been elucidated. Herein, we report that TvaF from the thioviridamide biosynthetic gene cluster is an FMN-dependent cysteine decarboxylase that transforms the C-terminal cysteine of precursor peptides into a thioenol motif and exhibits high substrate flexibility. We resolved the crystal structure of TvaF bound with FMN at 2.24 Å resolution. Key residues for FMN binding and catalytic activity of TvaF have been identified and evaluated by mutagenesis studies.

NADPH-cytochrome P450 reductase (CPR) plays a central role in the metabolism of insecticides. Numerous studies have shown that CPR is associated with insecticide resistance in insect. In this study, two transcripts of Spodoptera litura CPR (SlCPR-X1 and SlCPR-X2) were identified and cloned, and the deduced protein of SlCPR-X1 contains all the conserved CPR structural features (N-terminal membrane anchor, FMN, FAD and NADP binding domains, FAD binding motif, and catalytic residues). However, no N-terminal member anchor and a shorter FMN binding region have been identified in the deduced protein of SlCPR-X2. The specific expression patterns showed that SlCPR-X1 and SlCPR-X2 were detected in all tested developmental stages and tissues, but highly expressed in third-, fourth-, and fifth-instar larvae, and in midgut and fat body. In addition, compared with the susceptible strain, SlCPR-X1 and SlCPR-X2 were up-regulated and more inducible when treated with indoxacarb in the indoxacarb-resistant strain. However, the relative expression, up-regulation and induction of SlCPR-X1 were all higher than those of SlCPR-X2 in the indoxacarb-resistant strain. Furthermore, RNA interference and baculovirus expression system combined with MTT cytotoxicity assay demonstrated that only SlCPR-X1 with the N-terminal membrane anchor as the major CPR potentially involved in S. litura indoxacarb resistance. The outcome of this study further expands our understanding of the important role of insect CPR in xenobiotics detoxification and resistance development, and CPR could be a potential target for insecticide resistance management mediated by RNAi or CRISPR/Cas.

Cu(II)-enhanced microbial Cr(VI) reduction is common in the environment, yet its mechanism is unknown. The specific activity of chromate reductase, NfoR, from Staphylococcus aureus sp. LZ-01 was augmented 1.5-fold by Cu(II). Isothermal titration calorimetry and spectral data show that Cu(II) binds to NfoR nonspecifically. Further, Cu(II) stimulates the nitrobenzene reduction of NfoR, indicating that Cu(II) promotes electron transfer. The crystal structure of NfoR in complex with CuSO4 (1.46 A) was determined. The overall structure of NfoR-Cu(II) complex is a dimer that covalently binds with FMN and Cu(II)-binding pocket is located at the interface of the NfoR dimer. Structural superposition revealed that NfoR resembles the structure of class II chromate reductase. Site-directed mutagenesis revealed that Leu46 and Phe123 were involved in NADH binding, whereas Trp70 and Ser45 were the key residues for nitrobenzene binding. Furthermore, His100 and Asp171 were preferential affinity sites for Cu(II) and that Cys163 is an active site for FMN binding. Attenuation reductase activity in C163S can be partially restored to 54% wild type by increasing Cu(II) concentration. Partial restoration indicates dual-channel electron transfer of NfoR via Cu(II) and FMN. We propose a catalytic mechanism for Cu(II)-enhanced NfoR activity in which Cu(I) is formed transiently. Together, the current results provide an insight on Cu (II)-induced enhancement and benefit of Cr(VI) bioremediation.

MitoNEET is a mitochondrial outer membrane protein that hosts a redox active [2Fe-2S] cluster in the C-terminal cytosolic domain. Increasing evidence has shown that mitoNEET has an essential role in regulating energy metabolism in human cells. Previously, we reported that the [2Fe-2S] clusters in mitoNEET can be reduced by the reduced flavin mononucleotide (FMNH2) and oxidized by oxygen or ubiquinone-2, suggesting that mitoNEET may act as a novel redox enzyme catalyzing electron transfer from FMNH2 to oxygen or ubiquinone. Here, we explore the FMN binding site in mitoNEET by using FMN analogs and find that lumiflavin, like FMN, at nanomolar concentrations can mediate the redox transition of the mitoNEET [2Fe-2S] clusters in the presence of flavin reductase and NADH (100 μM) under aerobic conditions. The electron paramagnetic resonance (EPR) measurements show that both FMN and lumiflavin can dramatically change the EPR spectrum of the reduced mitoNEET [2Fe-2S] clusters and form a covalently bound complex with mitoNEET under blue light exposure, suggesting that FMN/lumiflavin has specific interactions with the [2Fe-2S] clusters in mitoNEET. In contrast, lumichrome, another FMN analog, fails to mediate the redox transition of the mitoNEET [2Fe-2S] clusters and has no effect on the EPR spectrum of the reduced mitoNEET [2Fe-2S] clusters under blue light exposure. Instead, lumichrome can effectively inhibit the FMNH2-mediated reduction of the mitoNEET [2Fe-2S] clusters, indicating that lumichrome may act as a potential inhibitor to block the electron transfer activity of mitoNEET.

NADPH cytochrome P450 oxidoreductase (CPR) is the obligatory electron supplier that sustains the activity of microsomal cytochrome P450 (CYP) enzymes. The variant nature of the isoform-specific proximal interface of microsomal CYPs indicates that CPR is capable of multiple degenerated interactions with CYPs for electron transfer, through different binding mechanisms, and which are still not well-understood. Recently, we showed that CPR dynamics allows formation of open conformations that can be sampled by its structurally diverse redox partners in a CYP-isoform dependent manner. To further investigate the role of the CPR FMN-domain in effective binding of CPR to its diverse acceptors and to clarify the underlying molecular mechanisms, five different CPR-FMN-domain random mutant libraries were created. These libraries were screened for mutants with increased activity when combined with specific CYP-isoforms. Seven CPR-FMN-domain mutants were identified, supporting a gain in activity for CYP1A2 (P117H, G144C, A229T), 2A6 (P117L/L125V, G175D, H183Y), or 3A4 (N151D). Effects were evaluated using extended enzyme kinetic analysis, cytochrome b5 competition, ionic strength effect on CYP activity, and structural analysis. Mutated residues were located either in or adjacent to several acidic amino acid stretches – formerly indicated to be involved in CPR:CYP interactions – or close to two tyrosine residues suggested to be involved in FMN binding. Several of the identified positions co-localize with mutations found in naturally occurring CPR variants that were previously shown to cause CYP-isoform-dependent effects. The mutations do not seem to significantly alter the geometry of the FMN-domain but are likely to cause very subtle alterations leading to improved interaction with a specific CYP. Overall, these data suggest that CYPs interact with CPR using an isoform specific combination of several binding motifs of the FMN-domain. publishersversion published

8 pags, 4 figs, 3 tabs. -- Nine figures are available athttp://www.biophysj.org/biophysj/supplemental/S0006-3495(18)30963-9., Flavin mononucleotide (FMN) and flavin-adenine dinucleotide (FAD) are essential flavoprotein cofactors. A riboflavin kinase (RFK) activity catalyzes riboflavin phosphorylation to FMN, which can then be transformed into FAD by an FMN:adenylyltransferase (FMNAT) activity. Two enzymes are responsible for each one of these activities in eukaryotes, whereas prokaryotes have a single bifunctional enzyme, FAD synthase (FADS). FADS folds in two independent modules: the C-terminal with RFK activity and the N-terminal with FMNAT activity. Differences in structure and chemistry for the FMNAT catalysis among prokaryotic and eukaryotic enzymes pointed to the FMNAT activity of prokaryotic FADS as a potential antimicrobial target, making the structural model of the bacterial FMNAT module in complex with substrates relevant to understand the FADS catalytic mechanism and to the discovery of antimicrobial drugs. However, such a crystallographic complex remains elusive. Here, we have used molecular docking and molecular dynamics simulations to generate energetically stable interactions of the FMNAT module of FADS from Corynebacterium ammoniagenes with ATP/Mg and FMN in both the monomeric and dimer-of-trimers assemblies reported for this protein. For the monomer, we have identified the residues that accommodate the reactive phosphates in a conformation compatible with catalysis. Interestingly, for the dimer-of-trimers conformation, we have found that the RFK module negatively influences FMN binding at the interacting FMNAT module. These results agree with calorimetric data of purified samples containing nearly 100% monomer or nearly 100% dimer-of-trimers, indicating that FMN binds to the monomer but not to the dimer-of-trimers. Such observations support regulation of flavin homeostasis by quaternary C. ammoniagenes FADS assemblies., This work has been supported by Departamento Administrativo de Ciencia,Tecnología e Innovación (Colciencias, programa ‘‘Es tiempo de volver’’ and ‘‘Max Planck agreement’’), Vicerrectoría de Investigación y Extensión Universidad Industrial de Santander (project 1818 to M.C.D), University of Antioquia (Colombia), the Max Planck Society (Germany), the Spanish Ministry of Economy, Industry and Competitiveness (BIO2016-75183-PAEI/FEDER, UE to M.M.) and the Government of Aragón-Fondo Europeo de Desarrollo Regional (E35_17R).

Cytochrome P450 reductase (CYPOR) provides electrons to all human microsomal cytochrome P450s (cyt P450s). The length and sequence of the “140s” FMN binding loop of CYPOR has been shown to be a key determinant of its redox potential and activity with cyt P450s. Shortening the “140s loop” by deleting glycine-141(ΔGly141) and by engineering a second mutant that mimics flavo-cytochrome P450 BM3 (ΔGly141/Glu142Asn) resulted in mutants that formed an unstable anionic semiquinone. In an attempt to understand the molecular basis of the inability of these mutants to support activity with cyt P450, we expressed, purified, and determined their ability to reduce ferric P450. Our results showed that the ΔGly141 mutant with a very mobile loop only reduced ~7% of cyt P450 with a rate similar to that of the wild type. On the other hand, the more stable loop in the ΔGly141/Glu142Asn mutant allowed for ~55% of the cyt P450 to be reduced ~60% faster than the wild type. Our results reveal that the poor activity of the ΔGly141 mutant is primarily accounted for by its markedly diminished ability to reduce ferric cyt P450. In contrast, the poor activity of the ΔGly141/Glu142Asn mutant is presumably a consequence of the altered structure and mobility of the “140s loop”.

Sulfur compounds in fossil fuels are a major source of environmental pollution, and microbial desulfurization has emerged as a promising technology for removing sulfur under mild conditions. The enzyme TdsC from the thermophile Paenibacillus sp. A11-2 is a two-component flavin-dependent monooxygenase that catalyzes the oxygenation of dibenzothiophene (DBT) to its sulfoxide (DBTO) and sulfone (DBTO2) during microbial desulfurization. The crystal structures of the apo and flavin mononucleotide (FMN)-bound forms of DszC, an ortholog of TdsC, were previously determined, although the structure of the ternary substrate–FMN–enzyme complex remains unknown. Herein, we report the crystal structures of the DBT–FMN–TdsC and DBTO–FMN–TdsC complexes. These ternary structures revealed many hydrophobic and hydrogen-bonding interactions with the substrate, and the position of the substrate could reasonably explain the two-step oxygenation of DBT by TdsC. We also determined the crystal structure of the indole-bound enzyme because TdsC, but not DszC, can also oxidize indole, and we observed that indole binding did not induce global conformational changes in TdsC with or without bound FMN. We also found that the two loop regions close to the FMN-binding site are disordered in apo-TdsC and become structured upon FMN binding. Alanine substitutions of Tyr-93 and His-388, which are located close to the substrate and FMN bound to TdsC, significantly decreased benzothiophene oxygenation activity, suggesting their involvement in supplying protons to the active site. Interestingly, these substitutions increased DBT oxygenation activity by TdsC, indicating that expanding the substrate-binding site can increase the oxygenation activity of TdsC on larger sulfur-containing substrates, a property that should prove useful for future microbial desulfurization applications.

The activity of the reversible decarboxylase enzyme Fdc1 is dependent on prenylated FMN (prFMN), a recently discovered cofactor. The oxidized prFMN supports a 1,3-dipolar cycloaddition mechanism that underpins reversible decarboxylation. Fdc1 is a distinct member of the UbiD family of enzymes, with the canonical UbiD catalyzing the (de)carboxylation of para-hydroxybenzoic acid-type substrates. Here we show that the Escherichia coli UbiD enzyme, which is implicated in ubiquinone biosynthesis, cannot be isolated in an active holoenzyme form despite the fact active holoFdc1 is readily obtained. Formation of holoUbiD requires reconstitution in vitro of the apoUbiD with reduced prFMN. Furthermore, although the Fdc1 apoenzyme can be readily reconstituted and activated, in vitro oxidation to the mature prFMN cofactor stalls at formation of a radical prFMN species in holoUbiD. Further oxidative maturation in vitro occurs only at alkaline pH, suggesting a proton-coupled electron transfer precedes formation of the fully oxidized prFMN. Crystal structures of holoUbiD reveal a relatively open active site potentially occluded from solvent through domain motion. The presence of a prFMN sulfite-adduct in one of the UbiD crystal structures confirms oxidative maturation does occur at ambient pH on a slow time scale. Activity could not be detected for a range of putative para-hydroxybenzoic acid substrates tested. However, the lack of an obvious hydrophobic binding pocket for the octaprenyl tail of the proposed ubiquinone precursor substrate does suggest UbiD might act on a non-prenylated precursor. Our data reveals an unexpected variation occurs in domain mobility, prFMN binding, and maturation by the UbiD enzyme family.

Hormesis is a biphasic concentration–response relationship. During the luminescence inhibition test of Vibrio qinghaiensis sp.-Q67 (Q67), some organic solvents display the hormesis phenomenon. However, the mechanism of hormesis with respect to organic solvents remains unclear. This study focuses on luciferase, which is the key factor in luminescent reactions, and explores its role in the mechanism of hormesis. Q67 luciferase has two subunits, α and β. Molecular docking and molecular dynamics simulations were carried out by taking organic solvents as ligand and the two subunits as receptor. In addition, the binding free energies of the complexes formed by the ligand and Q67 luciferase were calculated. The results showed that the organic solvent ligands exhibiting hormesis bind to the β subunit first, while those that do not exhibit hormesis bind more easily to the α subunit. The hormetic organic solvents bind to β subunit first at low concentration, and change the flexibility of residues Ser145–Arg165 located on the α subunit; this enables flavin mononucleotide (FMN) to bind to the α subunit, exhibiting the hormesis phenomenon. With the increasing concentration, redundant molecules start to bind to the α subunit and compete with/block FMN binding to the α subunit, resulting in inhibition.

17 pags, 9 figs, 1 tab. -- Supplementary materials can be found at http://www.mdpi.com/1422-0067/20/24/ 6203/s1, FAD synthase (FADS, or FMN:ATP adenylyl transferase) coded by the FLAD1 gene is the last enzyme in the pathway of FAD synthesis. The mitochondrial isoform 1 and the cytosolic isoform 2 are characterized by the following two domains: the C-terminal PAPS domain (FADSy) performing FAD synthesis and pyrophosphorolysis; the N-terminal molybdopterin-binding domain (FADHy) performing a Co++ /K+-dependent FAD hydrolysis. Mutations in FLAD1 gene are responsible for riboflavin responsive and non-responsive multiple acyl-CoA dehydrogenases and combined respiratory chain deficiency. In patients harboring frameshift mutations, a shorter isoform (hFADS6) containing the sole FADSy domain is produced representing an emergency protein. With the aim to ameliorate its function we planned to obtain an engineered more efficient hFADS6. Thus, the D238A mutant, resembling the D181A FMNAT “supermutant” of C. glabrata, was overproduced and purified. Kinetic analysis of this enzyme highlighted a general increase of Km, while the kcat was two-fold higher than that of WT. The data suggest that the FAD synthesis rate can be increased. Additional modifications could be performed to further improve the synthesis of FAD. These results correlate with previous data produced in our laboratory, and point towards the following proposals (i) FAD release is the rate limiting step of the catalytic cycle and (ii) ATP and FMN binding sites are synergistically connected., This work was supported by "Fondi di Ateneo" the Universita degli Studi di Bari (to M.B.); "Effetto di mutazioni di FLAD1 e di alterazioni dell'omeostasi delle flavine sullo stato redox e sulla biogenesi mitocondriale: uno studio integrato su fibroblasti umani" the Universita degli Studi di Bari (to M.B.); "Fondi di Ateneo" the Universita della Calabria (to C.I. and M.G.); the Spanish Agencia Estatal de Investigacion and Fondo Europeo de Desarrollo Regional BIO2016-75183-P AEI/FEDER, UE (to M.M.), and the Government of Aragon-FEDER Grupo de Referencia Biologia Estructural E35_17R (to M.M.).

Mr4511 from Methylobacterium radiotolerans is a 164 amino acid protein built of a flavin mononucleotide (FMN) binding, blue-light responsive LOV (Light, Oxygen, Voltage) core domain plus flanking regions. In contrast to the majority of LOV domains, Mr4511 lacks a tryptophan residue that was previously identified as a major quencher for the FMN triplet state in photosensitizers for singlet oxygen (SO) engineered from these photoreceptors. Here we show that for Mr4511 it is sufficient to only mutate the reactive cysteine responsible for the photocycle (Cys71) in the native protein to generate an efficient SO photosensitizer: both C71S and C71G variants exhibit SO quantum yields of formation, ΦΔ, around 0.2 in air-saturated solutions. Under oxygen saturated conditions, ΦΔ reaches ∼0.5 in deuterated buffer. The introduction of Trp112 in the canonical position for LOV domains dramatically lowers ΦΔ to values comparable to miniSOG, one of the early FMN binding proteins touted as a SO sensitizer. Besides its SO properties, Mr4511 is also exceedingly robust against denaturation with urea and is more photostable than free FMN.

The sulfur contents of fossil fuels have negative impacts on the environment and human health. The bio-catalytic desulfurization strategies and the biological refinement of fossil fuels are a cost-effective process compared to classical chemistry desulfurization. Rhodococcus erythropolis IGTS8 is able to metabolize the organic sulfur compound by the unique genes cluster (i.e. DszA, B, C and D genes) in the 4S metabolic pathway. The dszD gene codes a key enzyme for sulfur reduction in the gene cluster. In this study, the structure of the DszD enzyme was predicted and then the key residues toward FMN binding were identified which were Thr62, Ser63, Asn77, and Ala79. To investigate the effect of manipulation in key residues on the enzymatic activity of the DszD, different mutations were performed on key residues. The molecular docking simulation showed that A79I and A79N mutants have the lowest binding free energies compared to the wild-type enzyme in binding with FMN substrate. A 50 ns molecular dynamics (MD) simulation performed using GROMACS software. The RMSD and RMSF analysis showed that two mutants are more stable than the wild-type enzyme during MD simulation. The binding free energies between FMN substrate and complexes were calculated and analyzed by the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM-PBSA) method. The experimental results showed that the enzyme activity for the oxidoreductase process toward biodesulfurization increased 1.9 and 2.3 fold for A79I and A79N mutants, respectively.

Display Omitted Structural and dynamical properties of wild-type and mutated E13źK, E13R, E13T, and E13Q of FMN binding protein (FBP) dimers are investigated by means of molecular dynamics simulation (MDS).In crystal structures, subunit A (Sub A) and subunit B (Sub B) were almost completely equivalent in all of the five FBP dimers. However, in solution, structures of the two subunits were not equivalent as predicted by MDS.Water molecules were rarely accessible to Iso in all FBP dimers which are in contrast with other flavoenzymes. The structural and dynamical properties of five FMN binding protein (FBP) dimers, WT (wild type), E13K (Glu13 replaced by Lys), E13R (Glu13 replaced by Arg), E13T (Glu13 replaced by Thr) and E13Q (Glu13 replaced by Gln), were investigated using a method of molecular dynamics simulation (MDS). In crystal structures, subunit A (Sub A) and subunit B (Sub B) were almost completely equivalent in all of the five FBP dimers. However, the predicted MDS structures of the two subunits were not equivalent in solution, revealed by the distances and inter-planar angles between isoalloxazine (Iso) and aromatic amino acids (Trp32, Tyr35 and Trp106) as well as the hydrogen bonding pairs between Iso and nearby amino acids. Residue root of mean square fluctuations (RMSF) also displayed considerable differences between Sub A and Sub B and in the five FBP dimers. The dynamics of the whole protein structures were examined with the distance (RNN) between the peptide N atom of the N terminal (Met1) and the peptide N atom of the C terminal (Leu122). Water molecules were rarely accessible to Iso in all FBP dimers which are in contrast with other flavoenzymes.

Bacterial riboswitches are non-coding RNA structural elements that direct gene expression in numerous metabolic pathways. The key regulatory roles of riboswitches, and the urgent need for new classes of antibiotics to treat multi-drug resistant bacteria, has led to efforts to develop small-molecules that mimic natural riboswitch ligands to inhibit metabolic pathways and bacterial growth. Recently, we reported the results of a phenotypic screen targeting the riboflavin biosynthesis pathway in the Gram-negative bacteria Escherichia coli that led to the identification of ribocil, a small molecule inhibitor of the flavin mononucleotide (FMN) riboswitch controlling expression of this biosynthetic pathway. Although ribocil is structurally distinct from FMN, ribocil functions as a potent and highly selective synthetic mimic of the natural ligand to repress riboswitch-mediated ribB gene expression and inhibit bacterial growth both in vitro and in vivo. Herein, we expand our analysis of ribocil; including mode of binding in the FMN binding pocket of the riboswitch, mechanisms of resistance and structure-activity relationship guided efforts to generate more potent analogs.

In this report, the main contributions of FMN were employed in the reductive cleavage reaction of AzrC protein (as a member of azoreductase family). Molecular dynamics simulations of three models in the presence and absence of FMN and ligand were performed to gather information about the dynamic nature of active site residues of AzrC. Combination of pairwise decomposition and alanine scanning calculations provides critical information about the FMN binding sites. The MD results analyzed by alanine scanning method revealed the high negative scores for N 10 (A) A, N 12 (A) A, S 17 (A) A and Y 151 (A) A mutations, which were in agreement with pairwise decomposition analyses. Hydrogen bond analyses indicated that these residues play critical roles in establishing appropriate hydrogen bonds between AzrC and FMN. Negative energy results for nonpolar residues such as W 103 (A), M 102 (A) and F 105 (A) and binding free energy analyses of three complexes indicate that the VDW interactions could be regarded as some favorable contribution in FMN and AzrC protein and confirmed the critical role of FMN in ligand binding (35.84 %), in addition to its catalytic function. This information could be used for future experimental investigations.

Enzymes from natural sources protect the environment via complex biological mechanisms, which aid in reductive immobilization of toxic metals including chromium. Nevertheless, progress was being made in elucidating high-resolution crystal structures of reductases and their binding with flavin mononucleotide (FMN) to understand the underlying mechanism of chromate reduction. Therefore, herein, we employed molecular dynamics (MD) simulations, principal component analysis (PCA), and binding free energy calculations to understand the dynamics behavior of these enzymes with FMN. Six representative chromate reductases in monomeric and dimeric forms were selected to study the mode, dynamics, and energetic component that drive the FMN binding process. As evidenced by MD simulation, FMN prefers to bind the cervix formed between the catalytic domain surrounded by strong conserved hydrogen bonding, electrostatic, and hydrophobic contacts. The slight movement and reorientation of FMN resulted in breakage of some crucial H-bonds and other nonbonded contacts, which were well compensated with newly formed H-bonds, electrostatic, and hydrophobic interactions. The critical residues aiding in tight anchoring of FMN within dimer were found to be strongly conserved in the bacterial system. The molecular mechanics combined with the Poisson-Boltzmann surface area binding free energy of the monomer portrayed that the van der Waals and electrostatic energy contribute significantly to the total free energy, where, the polar solvation energy opposes the binding of FMN. The proposed proximity relationships between enzyme and FMN binding site presented in this study will open up better avenues to engineer enzymes with optimized chromate reductase activity for sustainable bioremediation of heavy metals.

The human gut microbiota encodes β-glucuronidases (GUSs) that play key roles in health and disease via the metabolism of glucuronate-containing carbohydrates and drugs. Hundreds of putative bacterial GUS enzymes have been identified by metagenomic analysis of the human gut microbiome, but less than 10% have characterized structures and functions. Here we describe a set of unique gut microbial GUS enzymes that bind flavin mononucleotide (FMN). First, we show using mass spectrometry, isothermal titration calorimetry, and x-ray crystallography that a purified GUS from the gut commensal microbe Faecalibacterium prausnitzii binds to FMN on a surface groove located 30 A away from the active site. Second, utilizing structural and functional data from this FMN-binding GUS, we analyzed the 279 unique GUS sequences from the Human Microbiome Project database and identified 14 putative FMN-binding GUSs. We characterized four of these hits and solved the structure of two, the GUSs from Ruminococcus gnavus and Roseburia hominis, which confirmed that these are FMN binders. Third, binding and kinetic analysis of the FMN-binding site mutants of these five GUSs show that they utilize a conserved site to bind FMN that is not essential for GUS activity, but can affect KM. Lastly, a comprehensive structural review of the PDB reveals that the FMN-binding site employed by these enzymes is unlike any structurally characterized FMN binders to date. These findings reveal the first instance of an FMN-binding glycoside hydrolase and suggest a potential link between FMN and carbohydrate metabolism in the human gut microbiota.

Endothelial and neuronal nitric oxide synthases (eNOS, nNOS) are important signal generators in a number of processes including angiogenesis and neurotransmission. The homologous inducible isoform (iNOS) occupies a multitude of conformational states in a catalytic cycle, including subnanosecond input and output states and a distribution of ‘open’ conformations with average lifetimes of ~4.3 ns. In this study, fluorescence lifetime spectroscopy was used to probe conformational states of purified eNOS and nNOS in the presence of chaotropes, calmodulin, NADP+ and NADPH. Two-domain FMN/oxygenase constructs of nNOS were also examined with respect to calmodulin effects. Optical biosensing was used to analyze calmodulin binding in the presence of NADP+ and NADPH. Calmodulin binding induced a shift of the population away from the input and to the open and output states of NOS. NADP+ shifted the population towards the input state. The oxygenase domain, lacking the input state, provided a measure of calmodulin-induced open-output transitions. A mechanism for regulation by calmodulin and an elucidation of the catalytic mechanism are suggested by a ‘conformational lockdown’ model. Calmodulin speeds transitions between input and open and between open and output states, effectively reducing the conformational manifold, speeding catalysis. Conformational control of catalysis involves reorientation of the FMN binding domain, of which fluorescence lifetime is an indicator. The approach described herein is a new tool for biophysical and structural analysis of NOS enzymes, regulatory events and other homologous reductase-containing enzymes.A note to the readerThis manuscript has over the past several years been submitted to and rejected by several journals, usually on the basis of reviewer opinion that it was not an important enough result to merit inclusion in the journal. Owing to the passing of the first author and the loss of his expertise in fluorescence lifetime spectroscopy, it has become too onerous a task to continually revise the manuscript to suit the whims of reviewers who nevertheless still reject the work. We are thus simply releasing the final form of the manuscript to BioRxiv in the hopes that it finds a readership who will find, as we do, that the results are of value to the field.

Dibenzothiophene (DBT) and their derivatives, accounting for the major part of the sulfur components in crude oil, make one of the most significant pollution sources. The DBT sulfone monooxygenase BdsA, one of the key enzymes in the “4S” desulfurization pathway, catalyzes the oxidation of DBT sulfone to 2′-hydroxybiphenyl 2-sulfonic acid (HBPSi). Here, we determined the crystal structure of BdsA from Bacillus subtilis WU-S2B, at the resolution of 2.2 Å, and the structure of the BdsA-FMN complex at 2.4 Å. BdsA and the BdsA-FMN complex exist as tetramers. DBT sulfone was placed into the active site by molecular docking. Seven residues (Phe12, His20, Phe56, Phe246, Val248, His316, and Val372) are found to be involved in the binding of DBT sulfone. The importance of these residues is supported by the study of the catalytic activity of the active site variants. Structural analysis and enzyme activity assay confirmed the importance of the right position and orientation of FMN and DBT sulfone, as well as the involvement of Ser139 as a nucleophile in catalysis. This work combined with our previous structure of DszC provides a systematic structural basis for the development of engineered desulfurization enzymes with higher efficiency and stability.

Organohalide respiration (OHR) is the energy metabolism of anaerobic bacteria able to use halogenated organic compounds as terminal electron acceptors. While the terminal enzymes in OHR, so-called reductive dehalogenases, are well-characterized, the identity of proteins potentially involved in electron transfer to the terminal enzymes remains elusive. Among the accessory genes identified in OHR gene clusters, the C subunit (rdhC) could well code for the missing redox protein between the quinol pool and the reductive dehalogenase, although it was initially proposed to act as transcriptional regulator. RdhC sequences are characterised by the presence of multiple transmembrane segments, a flavin mononucleotide (FMN) binding motif and two conserved CX3CP motifs. Based on these features, we propose a curated selection of RdhC proteins identified in general sequence databases. Beside the Firmicutes from which RdhC sequences were initially identified, sequences belong to three additional phyla, the Chloroflexi, the Proteobacteria and the Bacteriodetes. The diversity of RdhC sequences mostly respects the phylogenetic distribution, suggesting that rdhC genes emerged relatively early in the evolution of the OHR metabolism. PceC, the C subunit of the tetrachloroethene (PCE) reductive dehalogenase is encoded by the conserved pceABCT gene cluster identified in Dehalobacter restrictus PER-K23 and in several strains of Desulfitobacterium hafniense. Surfaceome analysis of D. restrictus cells confirmed the predicted topology of the FMN-binding domain (FBD) of PceC that is the exocytoplasmic face of the membrane. Starting from inclusion bodies of a recombinant FBD protein, strategies for successful assembly of the FMN cofactor and refolding were achieved with the use of the flavin-trafficking protein from D. hafniense TCE1. Mass spectrometry analysis and site-directed mutagenesis of rFBD revealed that threonine-168 of PceC is binding FMN covalently. Our results suggest that PceC, and more generally RdhC proteins, may play a role in electron transfer in the metabolism of organohalide respiration.

The bacterial 2-nitroreductase NbaA is the primary enzyme initiating the degradation of 2-nitrobenzoate (2-NBA), and its activity is controlled by posttranslational modifications. To date, the structure of NbaA remains to be elucidated. In this study, the crystal structure of a Cys194Ala NbaA mutant was determined to a 1.7-Å resolution. The substrate analog 2-NBA methyl ester was used to decipher the substrate binding site by inhibition of the wild-type NbaA protein. Tandem mass spectrometry showed that 2-NBA methyl ester produced a 2-NBA ester bond at the Tyr193 residue in the wild-type NbaA but not residues in the Tyr193Phe mutant. Moreover, covalent binding of the 2-NBA methyl ester to Tyr193 reduced the reactivity of the Cys194 residue on the peptide link. The Tyr193 hydroxyl group was shown to be essential for enzyme catalysis, as a Tyr193Phe mutant resulted in fast dissociation of flavin mononucleotide (FMN) from the protein with the reduced reactivity of Cys194. FMN binding to NbaA varied with solution NaCl concentration, which was related to the catalytic activity but not to cysteine reactivity. These observations suggest that the Cys194 reactivity is negatively affected by a posttranslational modification of the adjacent Tyr193 residue, which interacts with FMN and the substrate in the NbaA catalytic site.

Flavodoxins in combination with the flavin mononucleotide (FMN) cofactor play important roles for electron transport in prokaryotes. Here, novel insights into the FMN-binding mechanism to flavodoxins-4 were obtained from the NMR structures of the apo-protein from Lactobacillus acidophilus (YP_193882.1) and comparison of its complex with FMN. Extensive reversible conformational changes were observed upon FMN binding and release. The NMR structure of the FMN complex is in agreement with the crystal structure (PDB ID: 3EDO) and exhibits the characteristic flavodoxin fold, with a central five-stranded parallel β-sheet and five α-helices forming an α/β-sandwich architecture. The structure differs from other flavoproteins in that helix α2 is oriented perpendicular to the β-sheet and covers the FMN-binding site. This helix reversibly unfolds upon removal of the FMN ligand, which represents a unique structural rearrangement among flavodoxins.

Microbes can reduce hexavalent chromium Cr (VI) to the less toxic and soluble trivalent Cr (III). Copper stimulates microbial reduction of Cr (VI) by the Bacillus, Ochrobactrum, and Gluconobacter species; however, the mechanism remains unclear. In our study, the rate of Cr (VI) reduction by Staphylococcus aureus LZ-01 was increased by 210 % when supplemented with 60 μM Cu (II). A putative NAD(P)H-flavin oxidoreductase gene (nfoR) was upregulated under Cr (VI) stress. NfoR-knockout mutant displayed impaired reduction of Cr (VI) and Cu (II)-enhanced Cr (VI) reduction by nfoR isogenic mutant was attenuated in the presence of Cu (II). In vitro tests showed an increased V max value of 25.22 μM min−1 mg−1 NfoR in the presence of Cu (II). Together, these results indicate that NfoR is responsible for Cu (II) enhancement. Isothermal titration calorimetry (ITC) assays confirmed the interaction of NfoR with Cu (II) at the dissociation constant of 85.5 μM. Site-directed mutagenesis indicates that His100, His128, and Met165 residues may be important for Cu (II) binding, while Cys163 is necessary for the FMN binding of NfoR. These findings show that Cu (II)-enhanced NfoR belongs to a new branch of Cr (VI) reductases and profoundly influences Cr (VI) reduction.

All cytochromes P450s in the endoplasmic reticulum rely on P450 oxidoreductase (POR) for their catalytic activities. Mutations in POR cause metabolic disorders of steroid hormone biosynthesis and affect certain drug metabolizing P450 activities. We studied mutations A115V, T142A, Q153R identified in the flavin mononucleotide (FMN) binding domain of POR that interacts with partner proteins and P284L located in the hinge region that is required for flexibility and domain movements in POR. Human wild type and mutant POR as well as CYP3A4 and CYP19A1 proteins in recombinant form were expressed in bacteria, and purified proteins were reconstituted in liposomes for enzyme kinetic assays. Quality of POR protein was checked by cytochrome c reduction assay as well as flavin content measurements. We found that proteins carrying mutations A115V, T142A located close to the FMN binding site had reduced flavin content compared to wild type POR and lost almost all activity to metabolize androstenedione via CYP19A1 and showed reduced CYP3A4 activity. The variant P284L identified from apparently normal subjects also had severe loss of both CYP19A1 and CYP3A4 activities, indicating this to be a potentially disease causing mutation. The mutation Q153R initially identified in a patient with disordered steroidogenesis showed remarkably increased activities of both CYP19A1 and CYP3A4 without any significant change in flavin content, indicating improved protein-protein interactions between POR Q153R and some P450 proteins. These results indicate that effects of mutations on activities of individual cytochromes P450 can be variable and a detailed analysis of each variant with different partner proteins is necessary to accurately determine the genotype-phenotype correlations of POR variants.

Effect of the charge (negative, positive or neutral) of amino acid residue-13 on the photoinduced electron transfer (ET) from Trp32, Tyr35 and Trp106 to the excited isoalloxazine was evaluated in the flavin mononucleotide-binding protein from Desulfovibrio vulgaris isolate Miyazaki F (DvFBP). The protein structures of the wild type and the four isoforms where glutamic acid-13 is replaced with lysine (E13K), arginine (E13R), threonine (E13T) and glutamine (E13Q) in aqueous solution were obtained by molecular dynamics simulation. The distances between the amino acid residue-13 and isoalloxazine (Iso), and between the amino acid residue-13 and the ET donors were longer than 1 nm. The ET rates were evaluated with the Kakitani and Mataga model (KM theory) from their ultrafast fluorescence dynamics by means of a non-linear least squares method. Electrostatic (ES) energies between the photo-products and other ionic groups in the proteins markedly varied among ET donors and among the DvFBP isoforms, while the oth...

In recent years, RNA has been shown to fulfil a number of cellular functions. This has led to much interest in elucidation of the structure of functional RNA molecules, and thus, in the preparation of suitably functionalized RNAs. The chemical synthesis of RNAs allows for the site-specific modification; however, is limited to sequences of about 60–70 nucleotides in length. At the example of the flavine mononucleotide (FMN) responsive aptamer of the ypaA riboswitch from B. subtilis, we demonstrate the highly efficient preparation of site-specifically modified long-mer RNAs. Our strategy consists of the chemical synthesis of fragments followed by enzymatic or chemical ligation. Splint ligation with T4 RNA ligase turned out to be most successful among the enyzymatic protocols. Highly efficient chemical ligation was performed by azide–alkyne cycloaddition of suitably modified RNA fragments. Wild-type and 2-aminopurine (2-AP)-modified variants of the ypaA aptamer were prepared. FMN binding to all synthesized ypaA aptamer variants is demonstrated. However, dissociation of FMN from its binding site by reduction of the isoalloxazin unit as demonstrated before for a small-hairpin-derived aptazyme could not be shown. This implies that either FMN is less accessible to reduction when it is bound to its natural aptamer; that reduced FMN remains bound to the aptamer; or that FMN upon reduction indeed is released from its binding site, without the aptamer folding back in the natural ligand-free state. The results of this study are of general interest to the preparation of site-specifically modified RNAs for investigation into structure and function.

Flavin mononucleotide (FMN)-binding proteins (FBP) from Desulfovibrio vulgaris , Miyazaki F, contain FMN as a cofactor. Upon photo-excitation of the FBP, photoinduced electron transfer (ET) takes place from tryptophan (Trp)32, tyrosine (Tyr)35 and Trp105 residues to the excited isoalloxazine ring (Iso ∗ ) of FMN. The ultrafast fluorescence dynamics of the wild type, and the E13K, E13R, E13T and E13Q substitution isoforms of FBP were simultaneously analyzed with molecular dynamics simulation structures and Kakitani-Mataga ET theory, to obtain the time-dependent ET rates. A bell-shaped behavior of the logarithmic ET rates (In k ET ) from Trp32 and Trp106 to Iso ∗ versus the centre-to-centre donor–acceptor distances ( Rc ) was obtained in the time domain of sub-ps to ps in all five FBP isoforms. The In k ET from Tyr35 to Iso ∗ linearly decreased with the Rc distance, which were much slower than those of Trp32 and Trp105. These findings suggest that the nuclear term in the ET theory is important in the ultrafast time domain, whereas the electronic coupling term may be dominant in slower ET rates.

The ubiX gene (UniProtKB code Q489U8) of Colwellia psychrerythraea strain 34H has been annotated as a putative flavin mononucleotide (FMN)-dependent aromatic acid decarboxylase. Based on previous studies of homologous proteins, CpsUbiX is thought to catalyze the decarboxylation of 3-octaprenyl-4-hydroxybenzoate to produce 2-polyprenylphenol in the ubiquinone-biosynthesis pathway using a noncovalently bound FMN molecule as a cofactor. However, the detailed mechanisms of this important enzyme are not yet clear and need to be further elucidated. In this study, it was found that the V47S single mutation resulted in a loss of FMN binding, resulting in the production of FMN-free CpsUbiX protein. This mutation is likely to destabilize FMN–protein interactions without affecting the overall structural folding. To fully characterize the conformational changes upon FMN binding and the enzymatic mechanism at the molecular level, the wild-type (FMN-bound) and V47S mutant (FMN-free) CpsUbiX proteins were purified and crystallized using the sitting-drop vapour-diffusion method. Furthermore, complete diffraction data sets for FMN-bound (space group C2221) and FMN-free (space group P23) forms were obtained to 2.0 and 1.76 Å resolution, respectively.

Flavin mononucleotide (FMN)-binding proteins (FBPs) play an important role in the electron transport process in bacteria. In this study, the structures of the FBP from Desulfovibrio vulgaris (DvFBP) (Miyazaki F) were compared between those obtained experimentally by nuclear magnetic resonance (NMR) spectroscopy and those derived from molecular dynamics simulations (MDSs). A high-residue root of mean square deviation (RMSD) was observed in residues located at both sides of the wings (Gly22, Glu23, Asp24, Ala59, Arg60, Asp61, Glu62, Gly75, Arg76, Asn77, Gly78 and Pro79), while a low-residue RMSD was found in residues located in a hollow of the structure (Asn12, Glu13, Gly14, Val15, Val16, Asn30, Thr31, Trp32, Asn33, Ser34, Gly69, Ser70, Arg71 and Lys72). Inter-planar angles between the Phe7 and Iso and between the Phe7 and Trp106 residues were remarkably different between the MDS- and NMR-derived DvFBP structures. Distribution of the torsion angles around the covalent bonds in the aliphatic chain of FMN was...

Flavodoxins are a family of FMN binding proteins widely distributed in prokaryotes. They involve in various electron transfer reactions using the non-covalently bound FMN cofactor as the redox center. The Escherichia coli yqcA gene was identified to encode a short-chain favodoxin based on sequence information. However, the structure of YqcA protein is unknown and its exact biological function in cell is yet to be investigated. Herein, we report the resonance assignments of 1H, 13C and 15N atoms of E. coli YqcA in both the apo and holo states.

Nitric oxide synthases (NOSs) synthesize nitric oxide (NO), a signaling molecule, from l-arginine, utilizing electrons from NADPH. NOSs are flavo-hemo proteins, with two flavin molecules (FAD and FMN) and one heme per monomer, which require the binding of calcium/calmodulin (Ca(2+)/CaM) to produce NO. It is therefore important to understand the molecular factors influencing CaM binding from a structure/function perspective. A crystal structure of the CaM-bound iNOS FMN-binding domain predicted a salt bridge between R536 of human iNOS and E47 of CaM. To characterize the interaction between the homologous Arg of rat nNOS (R753) and murine iNOS (R530) with CaM, the Arg was mutated to Ala and, in iNOS, to Glu. The mutation weakens the interaction between nNOS and CaM, decreasing affinity by ~3-fold. The rate of electron transfer from FMN is greatly attenuated; however, little effect on electron transfer from FAD is observed. The mutated proteins showed reduced FMN binding, from 20% to 60%, suggesting an influence of this residue on FMN incorporation. The weakened FMN binding may be due to conformational changes caused by the arginine mutation. Our data show that this Arg residue plays an important role in CaM binding and influences FMN binding.

Relationship was studied between rates of photoinduced electron transfer (PET) from individual aromatic amino acid and hydrogen bonding (H-bond) chain of Tyr-Gln-isoalloxazine (Iso) among seven flavin photoreceptors (FPR); four TePixDs [wild type (WT), Q50N (amino Gln50 was replaced by Asn), Q50A (Gln50 replaced by Ala), and Y8F (Tyr8 replaced by Phe)] and three AppAs [WT, W92F (Trp92 replaced by Phe), and Y9F (Tyr9 replaced by Phe)]. The PET rates were obtained from the reported fluorescence dynamics of the seven FPRs with atomic coordinates determined by molecular dynamic simulations (MDS) and Kakitani and Mataga (KM) electron transfer theory. In TePixD isoforms, the fastest PET rate was obtained to be of 4.0 × 10 −2 ps −1 from Tyr8 to Iso* in WT, while the rates were 8.7 × 10 −3 ps −1 from Trp91 in Q50N, 7.0 × 10 −2 ps −1 from Trp91 in Q50A and 1.1 × 10 −3 ps −1 from Trp91 in Y8F. Similarly, the fastest PET rates in AppAs were 4.4 × 10 −3 and 1.1 × 10 −3 ps −1 from Tyr9 to Iso* in both WT and W92F. Formations of H-bonds between Tyr and Iso were examined by MDS snapshots. Tyr8-Gln50-Iso H-bond chain in TePixD was found only in WT among the four isoforms, whereas Tyr9-Gln51-Iso H-bond chain in AppA was found in both WT and W92F. It was found that the PET rates from Tyr8 in TePixD and Tyr9 in AppA were fastest when Tyr, Gln and Iso form H-bond chains. The rates from the Tyrs were much faster than one from Tyr with similar donor-acceptor distance in non-FPR flavoprotein of FMN binding protein. The reasons for it were discussed with PET parameters contained in KM theory.

A putative chromate ion binding site was identified proximal to a rigidly bound FMN from electron densities in the crystal structure of the quinone reductase from Gluconacetobacter hansenii (Gh-ChrR) (3s2y.pdb). To clarify the location of the chromate binding site, and to understand the role of FMN in the NADPH-dependent reduction of chromate, we have expressed and purified four mutant enzymes involving the site-specific substitution of individual side chains within the FMN binding pocket that form non-covalent bonds with the ribityl phosphate (i.e., S15A and R17A in loop 1 between β1 sheet and α1 helix) or the isoalloxanzine ring (E83A or Y84A in loop 4 between the β3 sheet and α4 helix). Mutations that selectively disrupt hydrogen bonds between either the N3 nitrogen on the isoalloxanzine ring (i.e., E83) or the ribitylphos- phoate (i.e., S15) respectively result in 50% or 70% reductions in catalytic rates of chromate reduction. In comparison, mutations that disrupt π-π ring stacking interactions with the isoal-loxanzine ring (i.e., Y84) or a salt bridge with the ribityl phosphate result in 87% and 97% inhibittion. In all cases there are minimal alterations in chromate binding affinities. Collectively, these results support the hypothesis that chromate binds proximal to FMN, and implicate a structural role for FMN positioning for optimal chromate reduction rates. As side chains proximal to the β3/α4 FMN binding loop 4 contribute to both NADH and metal ion binding, we propose a model in which structural changes around the FMN binding pocket couples to both chromate and NADH binding sites.

The Pden_2689 gene encoding FerA, an NADH:flavin oxidoreductase required for growth of Paracoccus denitrificans under iron limitation, was cloned and overexpressed as a C-terminally His6-tagged derivative. The binding of substrates and products was detected and quantified by isothermal titration calorimetry and fluorometric titration. FerA binds FMN and FAD with comparable affinity in an enthalpically driven, entropically opposed process. The reduced flavin is bound more loosely than the oxidized one, which was confirmed by a negative shift in the redox potential of FMN after addition of FerA. Initial velocity and substrate analogs inhibition studies showed that FerA follows a random-ordered sequence of substrate (NADH and FMN) binding. The primary kinetic isotope effects from stereospecifically deuterated nicotinamide nucleotides demonstrated that hydride transfer occurs from the pro-S position and contributes to rate limitation for the overall reaction. The crystal structure of FerA revealed a twisted seven-stranded antiparallel β-barrel similar to that of other short chain flavin reductases. Only minor structural changes around Arg106 took place upon FMN binding. The solution structure FerA derived from small angle X-ray scattering (SAXS) matched the dimer assembly predicted from the crystal structure. Site-directed mutagenesis pinpointed a role of Arg106 and His146 in binding of flavin and NADH, respectively. Pull down experiments performed with cytoplasmic extracts resulted in a negative outcome indicating that FerA might physiologically act without association with other proteins. Rapid kinetics experiments provided evidence for a stabilizing effect of another P. denitrificans protein, the NAD(P)H:acceptor oxidoreducase FerB, against spontaneous oxidation of the FerA-produced dihydroflavin.

Streptococcus pneumoniae (Sp) strain TIGR4 is a virulent, encapsulated serotype that causes bacteremia, otitis media, meningitis and pneumonia. Increased bacterial resistance and limited efficacy of the available vaccine to some serotypes complicate the treatment of diseases associated to this microorganism. Flavodoxins are bacterial proteins involved in several important metabolic pathways. The Sp flavodoxin (Spfld) gene was recently reported to be essential for the establishment of meningitis in a rat model, which makes SpFld a potential drug target. To facilitate future pharmacological studies, we have cloned and expressed SpFld in E. coli and we have performed an extensive structural and biochemical characterization of both the apo form and its active complex with the FMN cofactor. SpFld is a short-chain flavodoxin containing 146 residues. Unlike the well-characterized long-chain apoflavodoxins, the Sp apoprotein displays a simple two-state thermal unfolding equilibrium and binds FMN with moderate affinity. The X-ray structures of the apo and holo forms of SpFld differ at the FMN binding site, where substantial rearrangement of residues at the 91-100 loop occurs to permit cofactor binding. This work will set up the basis for future studies aiming at discovering new potential drugs to treat S. pneumoniae diseases through the inhibition of SpFld., We acknowledge financial support from BFU2010-16297 and BFU2010-19504 [Ministerio de Ciencia e Innovación Spain], BFU2013-47064-P, BIO2014-57314-REDT and CTQ2013-44367-C2-2-P [Ministerio de Economía y Competitividad, Spain], and DGA (Protein Targets B89). We also thank synchrotron radiation sources DLS (Oxford), and in particular beamline I04-1 (experiment number MX8035-3 and MX8035-11). The research leading to these results has also received funding from the FP7 (2007–2013) under BIOSTRUCTX-7687. A.R.C. was funded by a Banco Santander Central Hispano/Universidad de Zaragoza predoctoral fellowship. M. C-G was recipient of a predoctoral fellowship from the Government of Aragón.

The Rhodococcus rhodochrous strain 11Y XplA enzyme is an unusual cytochrome P450-flavodoxin fusion enzyme that catalyzes reductive denitration of the explosive hexahydro-1,3,5-trinitro-1,3,5-triazene (RDX). We show by light scattering that XplA is a monomeric enzyme. XplA has high affinity for imidazole (K(d) = 1.6 μM), explaining previous reports of a red-shifted XplA Soret band in pure enzyme. The true Soret maximum of XplA is at 417 nm. Similarly, unusually weak XplA flavodoxin FMN binding (K(d) = 1.09 μM) necessitates its purification in the presence of the cofactor to produce hallmark flavin contributions absent in previously reported spectra. Structural and ligand-binding data reveal a constricted active site able to accommodate RDX and small inhibitory ligands (e.g. 4-phenylimidazole and morpholine) while discriminating against larger azole drugs. The crystal structure also identifies a high affinity imidazole binding site, consistent with its low K(d), and shows active site penetration by PEG, perhaps indicative of an evolutionary lipid-metabolizing function for XplA. EPR studies indicate heterogeneity in binding mode for RDX and other ligands. The substrate analog trinitrobenzene does not induce a substrate-like type I optical shift but creates a unique low spin EPR spectrum due to influence on structure around the distal water heme ligand. The substrate-free heme iron potential (-268 mV versus NHE) is positive for a low spin P450, and the elevated potential of the FMN semiquinone/hydroquinone couple (-172 mV) is also an adaptation that may reflect (along with the absence of a key Thr/Ser residue conserved in oxygen-activating P450s) the evolution of XplA as a specialized RDX reductase catalyst.

The novel cytochrome P450/redox partner fusion enzyme CYP116B1 from Cupriavidus metallidurans was expressed in and purified from Escherichia coli. Isolated CYP116B1 exhibited a characteristic Fe(II)CO complex with Soret maximum at 449 nm. EPR and resonance Raman analyses indicated low-spin, cysteinate-coordinated ferric haem iron at both 10 K and ambient temperature, respectively, for oxidized CYP116B1. The EPR of reduced CYP116B1 demonstrated stoichiometric binding of a 2Fe-2S cluster in the reductase domain. FMN binding in the reductase domain was confirmed by flavin fluorescence studies. Steady-state reduction of cytochrome c and ferricyanide were supported by both NADPH/NADH, with NADPH used more efficiently (K(m[NADPH]) = 0.9 ± 0.5 μM and K(m[NADH]) = 399.1 ± 52.1 μM). Stopped-flow studies of NAD(P)H-dependent electron transfer to the reductase confirmed the preference for NADPH. The reduction potential of the P450 haem iron was -301 ± 7 mV, with retention of haem thiolate ligation in the ferrous enzyme. Redox potentials for the 2Fe-2S and FMN cofactors were more positive than that of the haem iron. Multi-angle laser light scattering demonstrated CYP116B1 to be monomeric. Type I (substrate-like) binding of selected unsaturated fatty acids (myristoleic, palmitoleic and arachidonic acids) was shown, but these substrates were not oxidized by CYP116B1. However, CYP116B1 catalysed hydroxylation (on propyl chains) of the herbicides S-ethyl dipropylthiocarbamate (EPTC) and S-propyl dipropylthiocarbamate (vernolate), and the subsequent N-dealkylation of vernolate. CYP116B1 thus has similar thiocarbamate-oxidizing catalytic properties to Rhodoccocus erythropolis CYP116A1, a P450 involved in the oxidative degradation of EPTC.

LOV domains function as blue light-sensing modules in various photoreceptors in plants, fungi, algae, and bacteria. A LOV/LOV protein (LLP) has been found from Arabidopsis thaliana (AtLLP) as a two LOV domain-containing protein. However, its function remains unknown. We isolated cDNA clones coding for an LLP homolog from tomato (Solanum lycopersicum) and two homologs from the moss Physcomitrella patens. The tomato LLP (SlLLP) contains two LOV domains (LOV1 and LOV2 domains), as in AtLLP. Most of the amino acids required for association with chromophore are conserved in both LOV domains, except that the amino acid at the position equivalent to the cysteine essential for cysteinyl adduct formation is glycine in the LOV1 domain as in AtLLP. When expressed in Escherichia coli, SlLLP binds FMN and undergoes a self-contained photocycle upon irradiation of blue light. Analyses using mutant SlLLPs revealed that SlLLP binds FMN in both LOV domains, although the LOV1 domain does not show spectral changes on irradiation. However, when Gly(66) in the LOV1 domain, which is located at the position equivalent to the essential cysteine of LOV domains, is replaced by cysteine, the mutated LOV1 domain shows light-induced spectral changes. In addition, all four LOV domains of P. patens LLPs (PpLLP1 and PpLLP2) show the typical features of LOV domains, including the reactive cysteine in each. This study shows that plants have a new LOV domain-containing protein family with the typical biochemical and photochemical properties of other LOV domain-containing proteins such as the phototropins.

NADPH-cytochrome P450 reductase (CPR) transfers two reducing equivalents derived from NADPH via FAD and FMN to microsomal P450 monooxygenases in one-electron transfer steps. The crystal structure of yeast CPR (yCPR) contains a surface-exposed FMN binding site (FMN2 site) at the interface of the FMN binding and connecting domains, in addition to the single buried site that has been observed in rat CPR. This finding provides a testable hypothesis of how intramolecular (between FAD and FMN) and intermolecular (between FMN and P450) electron transfer may occur in CPR. To verify that occupancy of the FMN2 site is not an artifact of crystallization, a surface plasmon resonance (SPR) biosensor technique has been applied to probe the selectivity of this site under functional conditions. A series of kinetic and equilibrium binding experiments involving yCPR immobilized on different sensor chip surfaces was performed using FMN and FAD, as well as FMN-derived compounds, including riboflavin, dimethylalloxazine, and alloxazine, and other molecules which resemble the planar isoalloxazine ring structure. Only FMN and FAD showed stoichiometric binding responses. Binding affinity for FMN was in the submicromolar range, thirty times higher than that for FAD. Association kinetic rates for the yCPR/FMN complex were up to 60-fold higher than for the yCPR/FAD complex. Taken together, these data indicate that (i) the surface-exposed site in yCPR is highly selective toward binding flavins, (ii) binding of FMN in this site is notably favored, and finally, (iii) both the phosphate group and the isoalloxazine ring of FMN are essential for binding.

Old Yellow Enzymes (OYEs) are NAD(P)H dehydrogenases of not fully resolved physiological roles that are widespread among bacteria, plants, and fungi and have a great potential for biotechnological applications. We determined the apo form crystal structure of a member of the OYE class, glycerol trinitrate reductase XdpB, from Agrobacterium bohemicum R89-1 at 2.1 Å resolution. In agreement with the structures of the related bacterial OYEs, the structure revealed the TIM barrel fold with an N-terminal β-hairpin lid, but surprisingly, the structure did not contain its cofactor FMN. Its putative binding site was occupied by a pentapeptide TTSDN from the C-terminus of a symmetry related molecule. Biochemical experiments confirmed a specific concentration-dependent oligomerization and a low FMN content. The blocking of the FMN binding site can exist in vivo and regulates enzyme activity. Our bioinformatic analysis indicated that a similar self-inhibition could be expected in more OYEs which we designated as subgroup OYE C1. This subgroup is widespread among G-bacteria and can be recognized by the conserved sequence GxxDYP in proximity of the C termini. In proteobacteria, the C1 subgroup OYEs are typically coded in one operon with short-chain dehydrogenase. This operon is controlled by the tetR-like transcriptional regulator. OYEs coded in these operons are unlikely to be involved in the oxidative stress response as the other known members of the OYE family because no upregulation of XdpB was observed after exposing A. bohemicum R89-1 to oxidative stress.

Photoinduced electron transfer (PET) from Tyr21 to isoalloxazine (Iso) in the excited state (Iso*) is considered to be an initial step of the photosensing function of the blue-light sensing using flavin adenine dinucleotide (BLUF) component of the anti-repressor of the photosynthetic regulation (AppA). The PET mechanism was investigated via fluorescence dynamics of AppA and Kakitani and Mataga (KM) theories as well as by molecular dynamic (MD) simulation. The local structures of both the Y21F and W104F mutant AppAs around the Iso binding sites were quite different from those of the wild type (WT) AppA. The distances between Iso and Trp104 in Y21F, and between Iso and Tyr21 in W104F were shorter by 0.06 nm and 0.02 nm, respectively, compared to the WT. The frequency factor, ν 0 , in Tyr21 was 1.15-fold greater than that in Trp104. The critical distance between adiabatic and non-adiabatic PET processes, R 0 , was found to be very long in the AppA Tyr21. The large values of ν 0 and R 0 for Tyr21 of AppA compared to those in a non photosensing flavoprotein, FMN binding protein (FBP), were elucidated by hydrogen bond (H bond) chain between Tyr21 and Iso through Gln63. Interaction energies among Iso*, Trp104, Tyr21 and Gln63 in WT were calculated using the semi-empirical PM3 method. The amount of the transferred charge from Trp104 to Iso* in the WT exhibited a maximum at an interaction energy of around −20 kcal/mol, but decreased as the interaction energy (absolute value) increased.

Patients with congenital adrenal hyperplasia, exhibiting combined CYP17 and CYP21 deficiency, were shown by Arlt et al. (2004) to harbor a 541T→G mutation in exon 5 of POR (encoding NADPH-cytochrome P450 reductase, CYPOR), which resulted in a Y181D substitution that obliterated electron transfer capacity. Using bacterial expression models, we examined catalytic and physical properties of the human CYPOR Y181D variant. As purified, Y181D lacked flavin mononucleotide (FMN) and NADPH-cytochrome c reductase (NCR) activity but retained normal flavin adenine dinucleotide binding and NADPH utilization. Titration of the purified protein with FMN restored 64 of wild-type (WT) NCR activity in Y181D with an activation constant of ∼2 μM. As determined by FMN fluorescence quenching, Y181D had KdFMN = 7.3 μM. Biplasmid coexpression of CYPOR and CYP1A2, at the physiological ratio of ∼1:10 in the engineered MK_1A2_POR Escherichia coli strain, showed the compromised capacity of Y181D to support CYP1A2-catalyzed metabolism of the procarcinogens 2-aminoanthracene, 2-amino-3-methylimidazo(4,5-f)quinoline, and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone. Isolated MK1A2_POR membranes confirmed FMN stimulation of Y181D NCR activity with a 1.6 μM activation constant. CYP1A2 ethoxyresorufin-O-dealkylase activity of the MK1A2_PORY181D membranes, undetectable in the absence of added FMN, increased to 37% of MK1A2_PORWT membranes with a 1.2 μM FMN activation constant. Therefore, we conclude that compromised FMN binding is the specific molecular defect causing POR deficiency in patients with Y181D mutation and that this defect, in large part, can be overcome in vitro by FMN addition.

Mutations in pyridoxine 5'-phosphate oxidase are known to cause neonatal epileptic encephalopathy. This disorder has no cure or effective treatment and is often fatal. Pyridoxine 5'-phosphate oxidase catalyzes the oxidation of pyridoxine 5'-phosphate to pyridoxal 5'-phosphate, the active cofactor form of vitamin B(6) required by more than 140 different catalytic activities, including enzymes involved in amino acid metabolism and biosynthesis of neurotransmitters. Our aim is to elucidate the mechanism by which a homozygous missense mutation (R229W) in the oxidase, linked to neonatal epileptic encephalopathy, leads to reduced oxidase activity. The R229W variant is approximately 850-fold less efficient than the wild-type enzyme due to an approximately 192-fold decrease in pyridoxine 5'-phosphate affinity and an approximately 4.5-fold decrease in catalytic activity. There is also an approximately 50-fold reduction in the affinity of the R229W variant for the FMN cofactor. A 2.5 A crystal structure of the R229W variant shows that the substitution of Arg-229 at the FMN binding site has led to a loss of hydrogen-bond and/or salt-bridge interactions between FMN and Arg-229 and Ser-175. Additionally, the mutation has led to an alteration of the configuration of a beta-strand-loop-beta-strand structure at the active site, resulting in loss of two critical hydrogen-bond interactions involving residues His-227 and Arg-225, which are important for substrate binding and orientation for catalysis. These results provide a molecular basis for the phenotype associated with the R229W mutation, as well as providing a foundation for understanding the pathophysiological consequences of pyridoxine 5'-phosphate oxidase mutations.

Two previously reported holoprotein crystal forms of the flavodoxin-like E. coli protein WrbA, diffracting to 2.6 and 2.0 A resolution, and new crystals of WrbA apoprotein diffracting to 1.85 A, are refined and analysed comparatively through the lens of flavodoxin structures. The results indicate that differences between apo- and holoWrbA crystal structures are manifested on many levels of protein organization as well as in the FMN-binding sites. Evaluation of the influence of crystal contacts by comparison of lattice packing reveals the protein's global response to FMN binding. Structural changes upon cofactor binding are compared with the monomeric flavodoxins. Topologically non-equivalent residues undergo remarkably similar local structural changes upon FMN binding to WrbA or to flavodoxin, despite differences in multimeric organization and residue types at the binding sites. Analysis of the three crystal structures described here, together with flavodoxin structures, rationalizes functional similarities and differences of the WrbAs relative to flavodoxins, leading to a new understanding of the defining features of WrbAs. The results suggest that WrbAs are not a remote and unusual branch of the flavodoxin family as previously thought but rather a central member with unifying structural features.