The cumulative electrostatic effect of aromatic stacking interactions and the negative electrostatic environment of the flavin mononucleotide binding site is a major determinant of the reduction potential for the flavodoxin from Desulfovibrio vulgaris [Hildenborough]

Retuning the potential of the electrochemical leaf

The electrochemical leaf enables the electrification and control of multi-enzyme cascades by exploiting two discoveries: (i) the ability to electrify the photosynthetic enzyme ferredoxin NADP+ reductase (FNR), driving it to catalyse the interconversion of NADP+/NADPH whilst it is entrapped in a highly porous, metal oxide electrode, and (ii) the evidence that additional enzymes can be co-entrapped in the electrode pores where, through one NADP(H)-dependent enzyme, extended cascades can be driven by electrical connection to FNR, via NADP(H) recycling. By changing a critical active-site tyrosine to serine, FNR's exclusivity for NADP(H) is swapped for unphosphorylated NAD(H). Here we present an electrochemical study of this variant FNR, and show that in addition to the intended inversion of cofactor preference, this change to the active site has altered FNR's tuning of the flavin reduction potential, making it less reductive. Exploiting the ability to monitor the variant's activity with NADP(H) as a function of potential has revealed a trapped intermediate state, relieved only by applying a negative overpotential, which allows catalysis to proceed. Inhibition by NADP+ (very tightly bound) with respect to NAD(H) turnover was also revealed and interestingly, this inhibition changes depending on the applied potential. These findings are of critical importance for future exploitation of the electrochemical leaf.

Tuning Active Site Electron Density for Enhanced Molecular Recognition and Catalysis

Enzymes have an extraordinary ability to utilize aromatic interactions for molecular recognition and catalysis. We here report molecularly imprinted nanoparticle receptors. The aromatic "wall" material in the imprinted binding site is used to enhance the molecular recognition of aromatic guests that have similar charges, shapes, and sizes but differ in π-electron density. Additionally, aromatic interactions are employed to activate an electron-rich aryl leaving group on a glycoside, mimicking the nucleoside hydrolase of the parasite Trypanosoma vivax.

Identifying a key spot for electron mediator-interaction to tailor CO dehydrogenase's affinity

Fe‒S cluster-harboring enzymes, such as carbon monoxide dehydrogenases (CODH), employ sophisticated artificial electron mediators like viologens to serve as potent biocatalysts capable of cleaning-up industrial off-gases at stunning reaction rates. Unraveling the interplay between these enzymes and their associated mediators is essential for improving the efficiency of CODHs. Here we show the electron mediator-interaction site on ChCODHs (Ch, Carboxydothermus hydrogenoformans) using a systematic approach that leverages the viologen-reactive characteristics of superficial aromatic residues. By enhancing mediator-interaction (R57G/N59L) near the D-cluster, the strategically tailored variants exhibit a ten-fold increase in ethyl viologen affinity relative to the wild-type without sacrificing the turn-over rate (kcat). Viologen-complexed structures reveal the pivotal positions of surface phenylalanine residues, serving as external conduits for the D-cluster to/from viologen. One variant (R57G/N59L/A559W) can treat a broad spectrum of waste gases (from steel-process and plastic-gasification) containing O2. Decoding mediator interactions will facilitate the development of industrially high-efficient biocatalysts encompassing gas-utilizing enzymes.

The Role of the si-Face Tyrosine of a Homodimeric Ferredoxin-NADP+ Oxidoreductase from Bacillus subtilis during Complex Formation and Redox Equivalent Transfer with NADP+/H and Ferredoxin

In the crystal structure of ferredoxin-NADP+ oxidoreductase from Bacillus subtilis (BsFNR), Tyr50 stacks on the si-face of the isoalloxazine ring portion of the FAD prosthetic group. This configuration is highly conserved among the homodimeric ferredoxin-NAD(P)+ oxidoreductases (FNR) from Gram-positive bacteria and photosynthetic bacteria. In this report, pre-steady state reactions of Tyr50 variants with NADP+/NADPH and ferredoxin from B. subtilis (BsFd) were examined with stopped-flow spectrophotometry to assess the effects of the mutation on the formation of FNR-substrate complexes and following redox equivalent transfer. Mixing oxidized BsFNRs with NADPH resulted in a rapid complex formation followed by a rate-limiting hydride transfer. The substitution substantially modulated the intensity of the charge transfer absorption band and decreased the observed hydride transfer rates compared to the wild type. Reduction of the Y50W mutant by NADPH proceeded in a monophasic manner, while the Y50G and Y50S mutants did in biphasic phases. The reduced Tyr50 mutants hardly promoted the reduction of NADP+. Mixing oxidized BsFNRs with reduced BsFd resulted in the reduction of the FNRs. The observed FNR reduction rates of the three variants were comparable, but in the Y50G and Y50S mutants, the amount of the reduced FNR at the rapid phase was decreased, and a slow FNR reduction phase was observed. The obtained results suggest that the replacements of Tyr50 with Gly and Ser permitted the conformational change in the reduced form, which induced an asymmetric kinetic behavior between the protomers of the homodimeric BsFNR.

Assessing the Performance of Non-Equilibrium Thermodynamic Integration in Flavodoxin Redox Potential Estimation

Flavodoxins are enzymes that contain the redox-active flavin mononucleotide (FMN) cofactor and play a crucial role in numerous biological processes, including energy conversion and electron transfer. Since the redox characteristics of flavodoxins are significantly impacted by the molecular environment of the FMN cofactor, the evaluation of the interplay between the redox properties of the flavin cofactor and its molecular surroundings in flavoproteins is a critical area of investigation for both fundamental research and technological advancements, as the electrochemical tuning of flavoproteins is necessary for optimal interaction with redox acceptor or donor molecules. In order to facilitate the rational design of biomolecular devices, it is imperative to have access to computational tools that can accurately predict the redox potential of both natural and artificial flavoproteins. In this study, we have investigated the feasibility of using non-equilibrium thermodynamic integration protocols to reliably predict the redox potential of flavodoxins. Using as a test set the wild-type flavodoxin from Clostridium Beijerinckii and eight experimentally characterized single-point mutants, we have computed their redox potential. Our results show that 75% (6 out of 8) of the calculated reaction free energies are within 1 kcal/mol of the experimental values, and none exceed an error of 2 kcal/mol, confirming that non-equilibrium thermodynamic integration is a trustworthy tool for the quantitative estimation of the redox potential of this biologically and technologically significant class of enzymes.

Magnetic coupling modulation in meta-nitroxide-functionalized isoalloxazine magnets with redox-active units as efficient side-modulators

Magnetic conversion can be accomplished in a variety of ways, as organic molecules with switchable magnetic characteristics offer numerous technological applications. It is crucial to find magnetism-switchable systems because, in the field of organic magnetic materials, the redox-induced magnetic reversal is very simple to achieve and shows significant applications. Herein, we computationally design isoalloxazine-based diradicals through oxidizing N10 and adding a nitroxide to C8 as the spin source (i.e. 8-nitroxide-isoalloxazine 10-oxide, an m-phenylene-like nitroxide diradical expanded with a redox unit as a side-modulator) and its N1/N5-hydrogenated/protonated diradical derivatives and introducing substituents (-OH, -NH₂, and -NO₂) to C6. We demonstrate that the basically modified structure exhibits ferromagnetic (FM) characteristics with a magnetic coupling constant (J) of 561.3 cm^-1 calculated at the B3LYP/6-311+G(d,p) level, obeying the meta-phenylene-mediated diradical character, and dihydrogenation can lead to an AFM diradical with considerably large J (-976.1 cm^-1). Surprisingly, protonation at N1 or N5 can lead to distinctly different magnetic variations (561.3 → -1602.9 cm^-1 at N1 versus 561.3 → 379.1 cm^-1 at N5). Analyses indicate that small singlet-triplet energy gaps and small energy gaps between the highest occupied and lowest unoccupied molecular orbitals (HOMO, LUMO) of the closed shell singlet state are the key features of these isoalloxazine diradicals, and aromaticity variations, significant spin delocalization from the π-conjugated structure and spin polarization from the non-Kekule structure induced by modification are responsible for the magnetic conversion. Furthermore, the spin alternation rule, the singly occupied molecular orbital (SOMO) effect, and the SOMO-SOMO energy splitting of the triplet state are used to analyze these distinct variations. This work provides a novel understanding of the structures and characteristics of modified isoalloxazine diradicals, as well as essential details for the intricate design and characterization of new isoalloxazine-based potential organic magnetic switches.

Machine Learning for Efficient Prediction of Protein Redox Potential: The Flavoproteins Case

Determining the redox potentials of protein cofactors and how they are influenced by their molecular neighborhoods is essential for basic research and many biotechnological applications, from biosensors and biocatalysis to bioremediation and bioelectronics. The laborious determination of redox potential with current experimental technologies pushes forward the need for computational approaches that can reliably predict it. Although current computational approaches based on quantum and molecular mechanics are accurate, their large computational costs hinder their usage. In this work, we explored the possibility of using more efficient QSPR models based on machine learning (ML) for the prediction of protein redox potential, as an alternative to classical approaches. As a proof of concept, we focused on flavoproteins, one of the most important families of enzymes directly involved in redox processes. To train and test different ML models, we retrieved a dataset of flavoproteins with a known midpoint redox potential (E_m) and 3D structure. The features of interest, accounting for both short- and long-range effects of the protein matrix on the flavin cofactor, have been automatically extracted from each protein PDB file. Our best ML model (XGB) has a performance error below 1 kcal/mol (∼36 mV), comparing favorably to more sophisticated computational approaches. We also provided indications on the features that mostly affect the E_m value, and when possible, we rationalized them on the basis of previous studies.

The encapsulin from Thermotoga maritima is a flavoprotein with a symmetry matched ferritin-like cargo protein

Bacterial nanocompartments, also known as encapsulins, are an emerging class of protein-based 'organelles' found in bacteria and archaea. Encapsulins are virus-like icosahedral particles comprising a ~ 25-50 nm shell surrounding a specific cargo enzyme. Compartmentalization is thought to create a unique chemical environment to facilitate catalysis and isolate toxic intermediates. Many questions regarding nanocompartment structure-function remain unanswered, including how shell symmetry dictates cargo loading and to what extent the shell facilitates enzymatic activity. Here, we explore these questions using the model Thermotoga maritima nanocompartment known to encapsulate a redox-active ferritin-like protein. Biochemical analysis revealed the encapsulin shell to possess a flavin binding site located at the interface between capsomere subunits, suggesting the shell may play a direct and active role in the function of the encapsulated cargo. Furthermore, we used cryo-EM to show that cargo proteins use a form of symmetry-matching to facilitate encapsulation and define stoichiometry. In the case of the Thermotoga maritima encapsulin, the decameric cargo protein with fivefold symmetry preferentially binds to the pentameric-axis of the icosahedral shell. Taken together, these observations suggest the shell is not simply a passive barrier-it also plays a significant role in the structure and function of the cargo enzyme.

Assignments of 19F NMR resonances and exploration of dynamics in a long-chain flavodoxin

Flavodoxin is a small protein that employs a non-covalently bound flavin to mediate single-electron transfer at low potentials. The long-chain flavodoxins possess a long surface loop that is proposed to interact with partner proteins. We have incorporated ¹⁹F-labeled tyrosine in long-chain flavodoxin from Rhodopseudomonas palustris to gain a probe of possible loop dynamics, exploiting the presence of a Tyr in the long loop in addition to Tyr residues near the flavin. We report ¹⁹F resonance assignments for all four Tyrs, and demonstration of a pair of resonances in slow exchange, both corresponding to a Tyr adjacent to the flavin. We also provide evidence for dynamics affecting the Tyr in the long loop. Thus, we show that ¹⁹F NMR of ¹⁹F-Tyr labeled flavodoxin holds promise for monitoring possible changes in conformation upon binding to partner proteins.

Flavins in the electron bifurcation process

The discovery of a new energy-coupling mechanism termed flavin-based electron bifurcation (FBEB) in 2008 revealed a novel field of application for flavins in biology. The key component is the bifurcating flavin endowed with strongly inverted one-electron reduction potentials (FAD/FAD•^- ≪ FAD•^-/FADH^-) that cooperatively transfers in its reduced state one low and one high-energy electron into different directions and thereby drives an endergonic with an exergonic reduction reaction. As energy splitting at the bifurcating flavin apparently implicates one-electron chemistry, the FBEB machinery has to incorporate prior to and behind the central bifurcating flavin 2e-to-1e and 1e-to-2e switches, frequently also flavins, for oxidizing variable medium-potential two-electron donating substrates and for reducing high-potential two-electron accepting substrates. The one-electron carriers ferredoxin or flavodoxin serve as low-potential (high-energy) electron acceptors, which power endergonic processes almost exclusively in obligate anaerobic microorganisms to increase the efficiency of their energy metabolism. In this review, we outline the global organization of FBEB enzymes, the functions of the flavins therein and the surrounding of the isoalloxazine rings by which their reduction potentials are specifically adjusted in a finely tuned energy landscape.

Structure and function of an unusual flavodoxin from the domain Archaea

Flavodoxins, electron transfer proteins essential for diverse metabolisms in microbes from the domain Bacteria, are extensively characterized. Remarkably, although genomic annotations of flavodoxins are widespread in microbes from the domain Archaea, none have been isolated and characterized. Herein is described the structural, biochemical, and physiological characterization of an unusual flavodoxin (FldA) from Methanosarcina acetivorans, an acetate-utilizing methane-producing microbe of the domain Archaea In contrast to all flavodoxins, FldA is homodimeric, markedly less acidic, and stabilizes an anionic semiquinone. The crystal structure reveals an flavin mononucleotide (FMN) binding site unique from all other flavodoxins that provides a rationale for stabilization of the anionic semiquinone and a remarkably low reduction potentials for both the oxidized/semiquinone (-301 mV) and semiquinone/hydroquinone couples (-464 mV). FldA is up-regulated in acetate-grown versus methanol-grown cells and shown here to substitute for ferredoxin in mediating the transfer of low potential electrons from the carbonyl of acetate to the membrane-bound electron transport chain that generates ion gradients driving ATP synthesis. FldA offers potential advantages over ferredoxin by (i) sparing iron for abundant iron-sulfur proteins essential for acetotrophic growth and (ii) resilience to oxidative damage.

Functional group surface modifications for enhancing the formation and performance of exoelectrogenic biofilms on the anode of a bioelectrochemical system

Various new energy technologies have been developed to reduce reliance on fossil fuels. The bioelectrochemical system (BES), an integrated microbial-electrochemical energy conversion process, is projected to be a sustainable and environmentally friendly energy technology. However, low power density is still one of the main limiting factors restricting the practical application of BESs. To enhance power output, functional group modification on anode surfaces has been primarily developed to improve the bioelectrochemical performances of BESs in terms of startup, power density, chemical oxygen demand (COD) removal and coulombic efficiency (CE). This modification could change the anode surface characteristics: roughness, hydrophobicity, biocompatibility, chemical bonding and electrochemically active surface area. This will facilitate bacterial adhesion, biofilm formation and extracellular electron transfer (EET). Additionally, some antibacterial functional groups are applied on air cathodes in order to suppress aerobic biofilms and enhance cathodic oxygen reduction reactions (ORRs). Various modification strategies such as: soaking, heat treatment and plasma modification have been reported to introduce functional groups typically as O-, N- and S-containing groups. In this review, the effects of anode functional groups on electroactive bacteria through the whole biofilm formation process are summarized. In addition, the application of those modification technologies to improve bioelectricity generation, resource recovery, bioelectrochemical analysis and the production of value-added chemicals and biofuels is also discussed. Accordingly, this review aims to help scientists select the most appropriate functional groups and up-to-date methods to improve biofilm formation.

Structural insight into the high reduction potentials observed for Fusobacterium nucleatum flavodoxin

Flavodoxins are small flavin mononucleotide (FMN)-containing proteins that mediate a variety of electron transfer processes. The primary sequence of flavodoxin from Fusobacterium nucleatum, a pathogenic oral bacterium, is marked with a number of distinct features including a glycine to lysine (K13) substitution in the highly conserved phosphate-binding loop (T/S-X-T-G-X-T), variation in the aromatic residues that sandwich the FMN cofactor, and a more even distribution of acidic and basic residues. The E_ox/sq (oxidized/semiquinone; -43 mV) and E_sq/hq (semiquinone/hydroquinone; -256 mV) are the highest recorded reduction potentials of known long-chain flavodoxins. These more electropositive values are a consequence of the apoprotein binding to the FMN hydroquinone anion with ~70-fold greater affinity compared to the oxidized form of the cofactor. Inspection of the FnFld crystal structure revealed the absence of a hydrogen bond between the protein and the oxidized FMN N5 atom, which likely accounts for the more electropositive E_ox/sq . The more electropositive E_sq/hq is likely attributed to only one negatively charged group positioned within 12 Å of the FMN N1. We show that natural substitutions of highly conserved residues partially account for these more electropositive reduction potentials.

Redox control of iodotyrosine deiodinase

The redox chemistry of flavoproteins is often gated by substrate and iodotyrosine deiodinase (IYD) has the additional ability to switch between reaction modes based on the substrate. Association of fluorotyrosine (F-Tyr), an inert substrate analog, stabilizes single electron transfer reactions of IYD that are not observed in the absence of this ligand. The co-crystal of F-Tyr and a T239A variant of human IYD have now been characterized to provide a structural basis for control of its flavin reactivity. Coordination of F-Tyr in the active site of this IYD closely mimics that of iodotyrosine and only minor perturbations are observed after replacement of an active site Thr with Ala. However, loss of the side chain hydroxyl group removes a key hydrogen bond from flavin and suppresses the formation of its semiquinone intermediate. Even substitution of Thr with Ser decreases the midpoint potential of human IYD between its oxidized and semiquinone forms of flavin by almost 80 mV. This decrease does not adversely affect the kinetics of reductive dehalogenation although an analogous Ala variant exhibits a 6.7-fold decrease in its kcat /Km . Active site ligands lacking the zwitterion of halotyrosine are not able to induce closure of the active site lid that is necessary for promoting single electron transfer and dehalogenation. Under these conditions, a basal two-electron process dominates catalysis as indicated by preferential reduction of nitrophenol rather than deiodination of iodophenol.

Electrostatic Modification for Promotion of Flavin-Mediated Oxidation of a Probe for Flavin Detection

Electrostatic effects on the redox photochemistry of synthetic probes (1, 2, and 1-Zn) are examined by adjusting the thermodynamic driving force of their oxidation reactions. The redox photochemistry was simply controlled by introducing a zinc binding site (2,2'-dipicolylamine (DPA)) on the coumarin moiety of probe 2. Zinc complexation produced a positively charged environment on the coumarin (1-Zn), which lowered the electron density of a nearby 9 H-xanthene ring, attenuating the auto-oxidation of 1-Zn by 45 % compared with that of probe 1 at 298 K. The positive net charge of 1-Zn also provided an attractive Coulombic force toward the phosphate of flavin mononucleotide and flavin adenine dinucleotide, which lowered the reduction potential of the electron acceptor (isoalloxazine) and improved intermolecular electron transfer from the 9 H-xanthene ring to isoalloxazine. The flavin-mediated oxidation rate of 1-Zn was increased to 1.5 times that of probe 2. Probe 1-Zn showed highly selective sensing behaviour toward flavins, producing an intense brightness (ϵΦ_F =2.80×10³  m^-1  cm^-1 ) in the long-wavelength regions (λ_max =588 nm) upon flavin-mediated oxidation. Furthermore, probes 1-Zn and 2 were successfully applied to eosinophil imaging and the differential diagnosis of eosinophilia; this demonstrates their use as diagnostic tools.

Hydrophilicity boosted extracellular electron transfer in Shewanella loihica PV-4

A superhydrophilic electrode enables the drastically boosted bacterial EET activity ofShewanella loihicaPV-4. It is proposed that a hydrophilic electrode favors the reduced state of OMCs, and consequently both the EET activity and cell proliferation are highly facilitated.

Role of tryptophan residues of Erv1: Trp95 and Trp183 are important for its folding and oxidase function

Erv1 (essential for respiration and viability 1) is a FAD-dependent sulphydryl oxidase with a tryptophan-rich catalytic domain. We show that Trp95 and Trp183 are important for stabilizing the folding, FAD-binding, and function of Erv1, whilst other four tryptophan residues are not functionally important.

Erv1 is an FAD-dependent thiol oxidase of the ERV (essential for respiration and viability)/ALR (augmenter of liver regeneration) sub-family and an essential component of the mitochondrial import and assembly pathway. Erv1 contains six tryptophan residues, which are all located in the highly conserved C-terminal FAD-binding domain. Though important structural roles were predicted for the invariable Trp⁹⁵, no experimental study has been reported. In the present study, we investigated the structural and functional roles of individual tryptophan residues of Erv1. Six single tryptophan-to-phenylalanine yeast mutant strains were generated and their effects on cell viability were tested at various temperatures. Then, the mutants were purified from Escherichia coli. Their effects on folding, FAD-binding and Erv1 activity were characterized. Our results showed that Erv1^W95F has the strongest effect on the stability and function of Erv1 and followed by Erv1^W183F. Erv1^W95F results in a decrease in the T_m of Erv1 by 23°C, a significant loss of the oxidase activity and thus causing cell growth defects at both 30°C and 37°C. Erv1^W183F induces changes in the oligomerization state of Erv1, along with a pronounced effect on the stability of Erv1 and its function at 37°C, whereas the other mutants had no clear effect on the function of Erv1 including the highly conserved Trp¹⁵⁷ mutant. Finally, computational analysis indicates that Trp⁹⁵ plays a key role in stabilizing the isoalloxazine ring to interact with Cys¹³³. Taken together, the present study provided important insights into the molecular mechanism of how thiol oxidases use FAD in catalysing disulfide bond formation.

Flavin-Protein Complexes: Aromatic Stacking Assisted by a Hydrogen Bond

Enzyme-catalyzed reactions often rely on a noncovalently bound cofactor whose reactivity is tuned by its immediate environment. Flavin cofactors, the most versatile catalyst encountered in biology, are often maintained within the protein throughout numbers of complex ionic and aromatic interactions. Here, we have investigated the role of π-π stacking and hydrogen bond interactions between a tyrosine and the isoalloxazine moiety of the flavin adenine dinucleotide (FAD) in an FAD-dependent RNA methyltransferase. Combining several static and time-resolved spectroscopies as well as biochemical approaches, we showed that aromatic stacking is assisted by a hydrogen bond between the phenol group and the amide of an adjacent active site loop. A mechanism of recognition and binding of the redox cofactor is proposed.

A switch between one- and two-electron chemistry of the human flavoprotein iodotyrosine deiodinase is controlled by substrate

Reductive dehalogenation is not typical of aerobic organisms but plays a significant role in iodide homeostasis and thyroid activity. The flavoprotein iodotyrosine deiodinase (IYD) is responsible for iodide salvage by reductive deiodination of the iodotyrosine derivatives formed as byproducts of thyroid hormone biosynthesis. Heterologous expression of the human enzyme lacking its N-terminal membrane anchor has allowed for physical and biochemical studies to identify the role of substrate in controlling the active site geometry and flavin chemistry. Crystal structures of human IYD and its complex with 3-iodo-l-tyrosine illustrate the ability of the substrate to provide multiple interactions with the isoalloxazine system of FMN that are usually provided by protein side chains. Ligand binding acts to template the active site geometry and significantly stabilize the one-electron-reduced semiquinone form of FMN. The neutral form of this semiquinone is observed during reductive titration of IYD in the presence of the substrate analog 3-fluoro-l-tyrosine. In the absence of an active site ligand, only the oxidized and two-electron-reduced forms of FMN are detected. The pH dependence of IYD binding and turnover also supports the importance of direct coordination between substrate and FMN for productive catalysis.

Cell-secreted flavins bound to membrane cytochromes dictate electron transfer reactions to surfaces with diverse charge and pH

The variety of solid surfaces to and from which microbes can deliver electrons by extracellular electron transport (EET) processes via outer-membrane c-type cytochromes (OM c-Cyts) expands the importance of microbial respiration in natural environments and industrial applications. Here, we demonstrate that the bifurcated EET pathway of OM c-Cyts sustains the diversity of the EET surface in Shewanella oneidensis MR-1 via specific binding with cell-secreted flavin mononucleotide (FMN) and riboflavin (RF). Microbial current production and whole-cell differential pulse voltammetry revealed that RF and FMN enhance EET as bound cofactors in a similar manner. Conversely, FMN and RF were clearly differentiated in the EET enhancement by gene-deletion of OM c-Cyts and the dependency of the electrode potential and pH. These results indicate that RF and FMN have specific binding sites in OM c-Cyts and highlight the potential roles of these flavin-cytochrome complexes in controlling the rate of electron transfer to surfaces with diverse potential and pH.

Origin and evolution of the sodium -pumping NADH: ubiquinone oxidoreductase

The sodium -pumping NADH: ubiquinone oxidoreductase (Na⁺-NQR) is the main ion pump and the primary entry site for electrons into the respiratory chain of many different types of pathogenic bacteria. This enzymatic complex creates a transmembrane gradient of sodium that is used by the cell to sustain ionic homeostasis, nutrient transport, ATP synthesis, flagellum rotation and other essential processes. Comparative genomics data demonstrate that the nqr operon, which encodes all Na⁺-NQR subunits, is found in a large variety of bacterial lineages with different habitats and metabolic strategies. Here we studied the distribution, origin and evolution of this enzymatic complex. The molecular phylogenetic analyses and the organizations of the nqr operon indicate that Na⁺-NQR evolved within the Chlorobi/Bacteroidetes group, after the duplication and subsequent neofunctionalization of the operon that encodes the homolog RNF complex. Subsequently, the nqr operon dispersed through multiple horizontal transfer events to other bacterial lineages such as Chlamydiae, Planctomyces and α, β, γ and δ -proteobacteria. Considering the biochemical properties of the Na⁺-NQR complex and its physiological role in different bacteria, we propose a detailed scenario to explain the molecular mechanisms that gave rise to its novel redox- dependent sodium -pumping activity. Our model postulates that the evolution of the Na⁺-NQR complex involved a functional divergence from its RNF homolog, following the duplication of the rnf operon, the loss of the rnfB gene and the recruitment of the reductase subunit of an aromatic monooxygenase.

Electron transferases

The flavin isoalloxazine ring in electron transferases functions in a redox capacity, being able to take up electrons from a donor to subsequently deliver them to an acceptor. The main characteristics of these flavoproteins, including their unique ability to mediate obligatory processes of two-electron transfers with those involving single-electron transfer, are here described. To illustrate the versatility of these proteins, the acquired knowledge of the function of the two electron transferases involved in the cyanobacterial photosynthetic electron transfer from photosystem I to NADP(+) is presented. Many aspects of their biochemistry and biophysics have been extensively characterized using site-directed mutagenesis, steady-state and transient kinetics, spectroscopy, calorimetry, X-ray crystallography, electron paramagnetic resonance, and computational methods.

Crystal structure of cindoxin, the P450cin redox partner

The crystal structure of the flavin mononucleotide (FMN)-containing redox partner to P450cin, cindoxin (Cdx), has been determined to 1.3 Å resolution. The overall structure is similar to that of the FMN domain of human cytochrome P450 reductase. A Brownian dynamics–molecular dynamics docking method was used to produce a model of Cdx with its redox partner, P450cin. This Cdx–P450cin model highlights the potential importance of Cdx Tyr96 in bridging the FMN and heme cofactors as well P450cin Arg102 and Arg346. Each of the single-site Ala mutants exhibits ∼10% of the wild-type activity, thus demonstrating the importance of these residues for binding and/or electron transfer. In the well-studied P450cam system, redox partner binding stabilizes the open low-spin conformation of P450cam and greatly decreases the stability of the oxy complex. In sharp contrast, Cdx does not shift P450cin to a low-spin state, although the stability of oxy-P450cin is decreased 10-fold in the presence of Cdx. This indicates that Cdx may have a modest effect on the open–closed equilibrium in P450cin compared to that in P450cam. It has been postulated that part of the effector role of Pdx on P450cam is to promote a significant structural change that makes available a proton relay network involving Asp251 required for O₂ activation. The structure around the corresponding Asp in P450cin, Asp241, provides a possible structural reason for why P450cin is less dependent on its redox partner for functionally important structural changes.

Structure of dihydromethanopterin reductase, a cubic protein cage for redox transfer

Dihydromethanopterin reductase (Dmr) is a redox enzyme that plays a key role in generating tetrahydromethanopterin (H4MPT) for use in one-carbon metabolism by archaea and some bacteria. DmrB is a bacterial enzyme understood to reduce dihydromethanopterin (H2MPT) to H4MPT using flavins as the source of reducing equivalents, but the mechanistic details have not been elucidated previously. Here we report the crystal structure of DmrB from Burkholderia xenovorans at a resolution of 1.9 Å. Unexpectedly, the biological unit is a 24-mer composed of eight homotrimers located at the corners of a cubic cage-like structure. Within a homotrimer, each monomer-monomer interface exhibits an active site with two adjacently bound flavin mononucleotide (FMN) ligands, one deeply buried and tightly bound and one more peripheral, for a total of 48 ligands in the biological unit. Computational docking suggested that the peripheral site could bind either the observed FMN (the electron donor for the overall reaction) or the pterin, H2MPT (the electron acceptor for the overall reaction), in configurations ideal for electron transfer to and from the tightly bound FMN. On this basis, we propose that DmrB uses a ping-pong mechanism to transfer reducing equivalents from FMN to the pterin substrate. Sequence comparisons suggested that the catalytic mechanism is conserved among the bacterial homologs of DmrB and partially conserved in archaeal homologs, where an alternate electron donor is likely used. In addition to the mechanistic revelations, the structure of DmrB could help guide the development of anti-obesity drugs based on modification of the ecology of the human gut.

The conjugated oligoelectrolyte DSSN+ enables exceptional coulombic efficiency via direct electron transfer for anode-respiring Shewanella oneidensis MR-1—a mechanistic study

Biofilm electrochemistry reveals that DSSN+ increases coulombic efficiency by enhancing the native direct electron transfer pathway of S. oneidensis MR-1.

Rate enhancement of bacterial extracellular electron transport involves bound flavin semiquinones

Extracellular redox-active compounds, flavins and other quinones, have been hypothesized to play a major role in the delivery of electrons from cellular metabolic systems to extracellular insoluble substrates by a diffusion-based shuttling two-electron-transfer mechanism. Here we show that flavin molecules secreted by Shewanella oneidensis MR-1 enhance the ability of its outer-membrane c-type cytochromes (OM c-Cyts) to transport electrons as redox cofactors, but not free-form flavins. Whole-cell differential pulse voltammetry revealed that the redox potential of flavin was reversibly shifted more than 100 mV in a positive direction, in good agreement with increasing microbial current generation. Importantly, this flavin/OM c-Cyts interaction was found to facilitate a one-electron redox reaction via a semiquinone, resulting in a 10(3)- to 10(5)-fold faster reaction rate than that of free flavin. These results are not consistent with previously proposed redox-shuttling mechanisms but suggest that the flavin/OM c-Cyts interaction regulates the extent of extracellular electron transport coupled with intracellular metabolic activity.

Distant residues mediate picomolar binding affinity of a protein cofactor

Numerous proteins require cofactors to be active. Computer simulations suggest that cooperative interaction networks achieve optimal cofactor binding. There is a need for the experimental identification of the residues crucial for stabilizing these networks and thus for cofactor binding. Here we investigate the electron transporter flavodoxin, which contains flavin mononucleotide as non-covalently bound cofactor. We show that after binding flavin mononucleotide with nanomolar affinity, the protein relaxes extremely slowly (time constant ~5 days) to an energetically more favourable state with picomolar-binding affinity. Rare small-scale openings of this state are revealed through H/D exchange of N(3)H of flavin. We find that H/D exchange can pinpoint amino acids that cause tight cofactor binding. These hitherto unknown residues are dispersed throughout the structure, and many are located distantly from the flavin and seem irrelevant to flavodoxin's function. Quantification of the thermodynamics of ligand binding is important for understanding, engineering, designing and evolving ligand-binding proteins.

Flavodoxin requires tight binding of its FMN cofactor to be active, but the residues involved are unknown. In this biophysical study, FMN binding is shown to change from nanomolar to picomolar affinity on extremely slow protein relaxation and the residues responsible for cofactor binding are identified.

NADPH-cytochrome P450 oxidoreductase: prototypic member of the diflavin reductase family

NADPH-cytochrome P450 oxidoreductase (CYPOR) and nitric oxide synthase (NOS), two members of the diflavin oxidoreductase family, are multi-domain enzymes containing distinct FAD and FMN domains connected by a flexible hinge. FAD accepts a hydride ion from NADPH, and reduced FAD donates electrons to FMN, which in turn transfers electrons to the heme center of cytochrome P450 or NOS oxygenase domain. Structural analysis of CYPOR, the prototype of this enzyme family, has revealed the exact nature of the domain arrangement and the role of residues involved in cofactor binding. Recent structural and biophysical studies of CYPOR have shown that the two flavin domains undergo large domain movements during catalysis. NOS isoforms contain additional regulatory elements within the reductase domain that control electron transfer through Ca(2+)-dependent calmodulin (CaM) binding. The recent crystal structure of an iNOS Ca(2+)/CaM-FMN construct, containing the FMN domain in complex with Ca(2+)/CaM, provided structural information on the linkage between the reductase and oxgenase domains of NOS, making it possible to model the holo iNOS structure. This review summarizes recent advances in our understanding of the dynamics of domain movements during CYPOR catalysis and the role of the NOS diflavin reductase domain in the regulation of NOS isozyme activities.

Understanding the FMN cofactor chemistry within the Anabaena Flavodoxin environment

The chemical versatility of flavin cofactors within the flavoprotein environment allows them to play main roles in the bioenergetics of all type of organisms, particularly in energy transformation processes such as photosynthesis or oxidative phosphorylation. Despite the large diversity of properties shown by flavoproteins and of the biological processes in which they are involved, only two flavin cofactors, FMN and FAD (both derived from the 7,8-dimethyl-10-(1'-D-ribityl)-isoalloxazine), are usually found in these proteins. Using theoretical and experimental approaches we have carried out an evaluation of the effects introduced upon substituting the 7- and/or 8-methyls of the isoalloxazine ring in the chemical and oxido-reduction properties of the different atoms of the ring on free flavins and on the photosynthetic Anabaena Flavodoxin (a flavoprotein that replaces Ferredoxin as electron carrier from Photosystem I to Ferredoxin-NADP(+) reductase). In Anabaena Flavodoxin both the protein environment and the redox state contribute to modulate the chemical reactivity of the isoalloxazine ring. Anabaena apoflavodoxin is shown to be designed to stabilise/destabilise each one of the FMN redox states (but not of the analogues produced upon substitution of the 7- and/or 8-methyls groups) in the adequate proportions to provide Flavodoxin with the particular properties required for the functions in which it is involved in vivo. The 7- and/or 8-methyl groups of the ixoalloxazine can be discarded as the gate for electrons exchange in Anabaena Fld, but a key role in this process is envisaged for the C6 atom of the flavin and the backbone atoms of Asn58.

A partial convergence in action of methylene blue and artemisinins: antagonism with chloroquine, a reversal with verapamil, and an insight into the antimalarial activity of chloroquine

Artemisinins rapidly oxidize leucomethylene blue (LMB) to methylene blue (MB); they also oxidize dihydroflavins such as the reduced conjugates RFH₂ of riboflavin (RF), and FADH₂ of the cofactor flavin adenine dinucleotide (FAD), to the corresponding flavins. Like the artemisinins, MB oxidizes FADH₂, but unlike artemisinins, it also oxidizes NAD(P)H. Like MB, artemisinins are implicated in the perturbation of redox balance in the malaria parasite by interfering with parasite flavoenzyme disulfide reductases. The oxidation of LMB by artemisinin is inhibited by chloroquine (CQ), an inhibition that is abruptly reversed by verapamil (VP). CQ also inhibits artemisinin-mediated oxidation of RFH₂ generated from N-benzyl-1,4-dihydronicotinamide (BNAH)-RF, or FADH₂ generated from NADPH or NADPH-Fre, an effect that is also modulated by verapamil. The inhibition likely proceeds by the association of LMB or dihydroflavin with CQ, possibly involving donor-acceptor or π complexes that hinder oxidation by artemisinin. VP competitively associates with CQ, liberating LMB or dihydroflavin from their respective CQ complexes. The observations explain the antagonism between CQ-MB and CQ-artemisinins in vitro, and are reconcilable with CQ perturbing intraparasitic redox homeostasis. They further suggest that a VP-CQ complex is a means by which VP reverses CQ resistance, wherein such a complex is not accessible to the putative CQ-resistance transporter (PfCRT).

15N solid-state NMR as a probe of flavin H-bonding

Flavins mediate a wide variety of chemical reactions in biology. To learn how one cofactor can be made to execute different reactions in different enzymes, we are developing solid-state NMR (SSNMR) to probe the flavin electronic structure, via the (15)N chemical shift tensor principal values (δ(ii)). We find that SSNMR has superior responsiveness to H-bonds, compared to solution NMR. H-bonding to a model of the flavodoxin active site produced an increase of 10 ppm in the δ(11) of N5, although none of the H-bonds directly engage N5, and solution NMR detected only a 4 ppm increase in the isotropic chemical shift (δ(iso)). Moreover SSNMR responded differently to different H-bonding environments, as H-bonding with water caused δ(11) to decrease by 6 ppm, whereas δ(iso) increased by less than 1 ppm. Our density functional theoretical (DFT) calculations reproduce the observations, validating the use of computed electronic structures to understand how H-bonds modulate the flavin's reactivity.

On the midpoint potential of the FAD chromophore in a BLUF-domain containing photoreceptor protein

The redox-midpoint potential of the FAD chromophore in the BLUF domain of anti-transcriptional regulator AppA from Rhodobacter sphaeroides equals ∼-260mV relative to the calomel electrode. Altering the structure of its chromophore-binding pocket through site-directed mutagenesis brings this midpoint potential closer to that of free flavin in aqueous solution. The redox-midpoint potential of this BLUF domain is intermediate between those of LOV domains and Cryptochromes, which may rationalize the primary photochemistry observed in these three flavin-containing photoreceptor families. These results also imply that LOV domains, among the flavin-containing photosensory receptors, are least sensitive to intracellular chemical reduction in the dark.

Aromatic stacking between nucleobase and enzyme promotes phosphate ester hydrolysis in dUTPase

Aromatic interactions are well-known players in molecular recognition but their catalytic role in biological systems is less documented. Here, we report that a conserved aromatic stacking interaction between dUTPase and its nucleotide substrate largely contributes to the stabilization of the associative type transition state of the nucleotide hydrolysis reaction. The effect of the aromatic stacking on catalysis is peculiar in that uracil, the aromatic moiety influenced by the aromatic interaction is relatively distant from the site of hydrolysis at the alpha-phosphate group. Using crystallographic, kinetics, optical spectroscopy and thermodynamics calculation approaches we delineate a possible mechanism by which rate acceleration is achieved through the remote π–π interaction. The abundance of similarly positioned aromatic interactions in various nucleotide hydrolyzing enzymes (e.g. most families of ATPases) raises the possibility of the reported phenomenon being a general component of the enzymatic catalysis of phosphate ester hydrolysis.

Calculating chemically accurate redox potentials for engineered flavoproteins from classical molecular dynamics free energy simulations

The tricyclic isoalloxazine nucleus of the redox cofactors flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) acts as an electron sink in life-sustaining biological electron transfer (eT). The functional diversity of flavin-containing proteins (flavoproteins) transcends that of free flavins. A large body of experimental evidence attributes natural control of flavoprotein-mediated eT to tuning of the thermodynamic driving force by the protein environment. Understanding and engineering such modulation by the protein environment of the flavin redox potential (DeltaE(o)) is valuable in biotechnology and device design. In this study we employed classical molecular dynamics free energy simulations (MDFES), within a thermodynamic integration (TI) formalism, to calculate the change in FMN first reduction potential (DeltaDeltaE(o)(ox/sq)) imparted by 6 flavoprotein active site mutations. The combined performance of the AMBER ff03 (protein) and GAFF (cofactor) force fields was benchmarked against experimental data for mutations close to the isoalloxazine re- and si-faces that perturb the wild-type DeltaE(o)(ox/sq) value in Anabaena flavodoxin. The classical alchemical approach used in this study overestimates the magnitude of DeltaE(o) values, in common with other studies. Nevertheless, chemically accurate DeltaDeltaE(o) values--calculated to within 1 kcal mol(-1) of the experimental value--were obtained for five of the six mutations studied. We have shown that this approach is practical for quantitative in silico screening of the effect of mutations on the first reduction potential where experimental values and structural data are available for the wild-type flavoprotein. This approach promises to be useful as an integral part of future interdisciplinary strategies to engineer desired thermodynamic properties in flavoproteins of biotechnological interest.

Structural and mechanistic aspects of flavoproteins: photosynthetic electron transfer from photosystem I to NADP+

This minireview covers the research carried out in recent years into different aspects of the function of the flavoproteins involved in cyanobacterial photosynthetic electron transfer from photosystem I to NADP(+), flavodoxin and ferredoxin-NADP(+) reductase. Interactions that stabilize protein-flavin complexes and tailor the midpoint potentials in these proteins, as well as many details of the binding and electron transfer to protein and ligand partners, have been revealed. In addition to their role in photosynthesis, flavodoxin and ferredoxin-NADP(+) reductase are ubiquitous flavoenzymes that deliver NAD(P)H or low midpoint potential one-electron donors to redox-based metabolisms in plastids, mitochondria and bacteria. They are also the basic prototypes for a large family of diflavin electron transferases with common functional and structural properties. Understanding their mechanisms should enable greater comprehension of the many physiological roles played by flavodoxin and ferredoxin-NADP(+) reductase, either free or as modules in multidomain proteins. Many aspects of their biochemistry have been extensively characterized using a combination of site-directed mutagenesis, steady-state and transient kinetics, spectroscopy and X-ray crystallography. Despite these considerable advances, various key features of the structural-function relationship are yet to be explained in molecular terms. Better knowledge of these systems and their particular properties may allow us to envisage several interesting applications of these proteins beyond their physiological functions.

Effect of the Insertion of a Glycine Residue into the Loop Spanning Residues 536-541 on the Semiquinone State and Redox Properties of the Flavin Mononucleotide-Binding Domain of Flavocytochrome P450BM-3 from Bacillus megaterium

Despite sharing sequence and structural similarities with other diflavin reductases such as NADPH-cytochrome P450 reductase (CPR) and nitric oxide synthase, flavocytochrome P450BM-3 displays some unique redox and electron transferring properties, including the inability to thermodynamically stabilize the neutral semiquinone (SQ) state of the flavin mononucleotide (FMN) cofactor. Rather, the anionic SQ species is only transiently formed during rapid reduction. Why is this? The absence of a conserved glycine residue and, as a consequence, the shorter and less flexible cofactor-binding loop in P450BM-3 represents a notable difference from other diflavin reductases and the structurally related flavodoxin. This difference may facilitate the formation of a strong hydrogen bond between backbone amide NH group of Asn537 and N5 of the oxidized FMN, an interaction not found in the other proteins. In the flavodoxin, the conserved glycine residue plays a crucial role in a redox-linked conformational change that contributes to the thermodynamic stabilization of the neutral SQ species of the FMN through the formation of a hydrogen bond with the N5H group of the flavin. In this study, a glycine residue was inserted after Tyr536 in the loop within the isolated FMN-binding domain as well as the diflavin reductase domain of P450BM-3, a position equivalent to Gly141 in human CPR. As a result, the insertion variant was observed to accumulate the neutral form of the FMN SQ species much like CPR. The midpoint potential for the SQ/HQ couple decreased by 68 mV, while that for the OX/SQ couple remained unchanged. (15)N NMR data provide evidence of the disruption of the hydrogen bond between the backbone amide group of Asn537 and the N5 atom in the oxidized state of the FMN. Molecular models suggest that the neutral FMN SQ could be stabilized through hydrogen bonding with the backbone carbonyl group of the inserted glycine residue in a manner similar to that of CPR and the flavodoxin. The insertion of the glycine at the same location within the diflavin domain resulted in a purified protein that retained nearly stoichiometric levels of bound FAD but tended to lose the FMN cofactor. This preparation retained one-third of the ferricyanide reductase activity but <1% of the cytochrome c reductase activity of the wild type. However, the insertion variant reconstituted with FMN regained nearly half of the wild-type cytochrome c reductase activity. These results demonstrate the importance of the unique structural characteristics of the shorter loop in P450BM-3 in establishing the unique redox properties of the FMN in this protein but not its general cytochrome reductase activity.

Flavodoxin: a compromise between efficiency and versatility in the electron transfer from Photosystem I to Ferredoxin-NADP(+) reductase

Under iron-deficient conditions Flavodoxin (Fld) replaces Ferredoxin in Anabaena as electron carrier from Photosystem I (PSI) to Ferredoxin-NADP(+) reductase (FNR). Several residues modulate the Fld interaction with FNR and PSI, but no one appears as specifically critical for efficient electron transfer (ET). Fld shows a strong dipole moment, with its negative end directed towards the flavin ring. The role of this dipole moment in the processes of interaction and ET with positively charged surfaces exhibited by PSI and FNR has been analysed by introducing single and multiple charge reversal mutations on the Fld surface. Our data confirm that in this system interactions do not rely on a precise complementary surface of the reacting molecules. In fact, they indicate that the initial orientation driven by the alignment of dipole moment of the Fld molecule with that of the partner contributes to the formation of a bunch of alternative binding modes competent for the efficient ET reaction. Additionally, the fact that Fld uses different interaction surfaces to dock to PSI and to FNR is confirmed.

Redox properties and crystal structures of a Desulfovibrio vulgaris flavodoxin mutant in the monomeric and homodimeric forms

The mutant S64C of the short-chain flavodoxin from Desulfovibrio vulgaris has been designed to introduce an accessible and reactive group on the protein surface. Crystals have been obtained of both the monomeric and homodimeric forms of the protein, with the cofactor FMN in either the oxidized or the one electron-reduced (semiquinone) state, and the structures have been determined to high resolution. The redox properties of the different species have been investigated and the variations observed with respect to wild type have been related to the structural changes induced by the mutation and S-S bridge formation.

On the oxygen reactivity of flavoprotein oxidases: an oxygen access tunnel and gate in brevibacterium sterolicum cholesterol oxidase

The flavoprotein cholesterol oxidase from Brevibacterium sterolicum (BCO) possesses a narrow channel that links the active center containing the flavin to the outside solvent. This channel has been proposed to serve for the access of dioxygen; it contains at its "bottom" a Glu-Arg pair (Glu-475-Arg-477) that was found by crystallographic studies to exist in two forms named "open" and "closed," which in turn was suggested to constitute a gate functioning in the control of oxygen access. Most mutations of residues that flank the channel have minor effects on the oxygen reactivity. Mutations of Glu-311, however, cause a switch in the basic kinetic mechanism of the reaction of reduced BCO with dioxygen; wild-type BCO and most mutants show a saturation behavior with increasing oxygen concentration, whereas for Glu-311 mutants a linear dependence is found that is assumed to reflect a "simple" second order process. This is taken as support for the assumption that residue Glu-311 finely tunes the Glu-475-Arg-477 pair, forming a gate that functions in modulating the access/reactivity of dioxygen.

Tuning of the FMN binding and oxido-reduction properties by neighboring side chains in Anabaena flavodoxin

Contribution of three regions (phosphate-binding, 50's and 90's loops) of Anabaena apoflavodoxin to FMN binding and reduction potential was studied. Thr12 and Glu16 did not influence FMN redox properties, but Thr12 played a role in FMN binding. Replacement of Trp57 with Glu, Lys or Arg moderately shifted E(ox/sq) and E(sq/hq) and altered the energetic of the FMN redox states binding profile. Our data indicate that the side chain of position 57 does not modulate E(ox/sq) by aromatic stacking or solvent exclusion, but rather by influencing the relative strength of the H-bond between the N(5) of the flavin and the Asn58-Ile59 bond. A correlation was observed between the isoalloxazine increase in solvent accessibility and less negative E(sq/hq). Moreover, E(sq/hq) became less negative as positively charged residues were added near to the isoalloxazine. Ile59 and Ile92 were simultaneously mutated to Ala or Glu. These mutations impaired FMN binding, while shifting E(sq/hq) to less negative values and E(ox/sq) to more negative. These effects are discussed on the bases of the X-ray structures of some of the Fld mutants, suggesting that in Anabaena Fld the structural control of both electron transfer steps is much more subtle than in other Flds.

Influence of the protein environment on the redox potentials of flavodoxins from Clostridium beijerinckii

The flavin mononucleotide (FMN) quinones in flavodoxin have two characteristic redox potentials, namely, Em(FMNH./FMNH-) for the one-electron reduction of the protonated FMN (E1) and Em(FMN/FMNH.) for the proton-coupled one-electron reduction (E2). These redox potentials in native and mutant flavodoxins obtained from Clostridium beijerinckii were calculated by considering the protonation states of all titratable sites as well as the energy contributed at the pKa value of FMN during protonation at the N5 nitrogen (pKa(N5)). E1 is sensitive to the subtle differences in the protein environments in the proximity of FMN. The protein dielectric volume that prevents the solvation of charged FMN quinones is responsible for the downshift of 130-160 mV of the E1 values with respect to that in an aqueous solution. The influence of the negatively charged 5'-phosphate group of FMN quinone on E1 could result in a maximum shift of 90 mV. A dramatic difference of 130 mV in the calculated E2 values of FMN quinone of the native and G57T mutant flavodoxins is due to the difference in the pKa(N5) values. This is due to the difference in the influence exerted by the carbonyl group of the protein backbone at residue 57.

15N solid-state NMR provides a sensitive probe of oxidized flavin reactive sites

Flavins are central to the reactivity of a wide variety of enzymes and electron transport proteins. There is great interest in understanding the basis for the different reactivities displayed by flavins in different protein contexts. We propose solid-state nuclear magnetic resonance (SS-NMR) as a tool for directly observing reactive positions of the flavin ring and thereby obtaining information on their frontier orbitals. We now report the SS-NMR signals of the redox-active nitrogens N1 and N5, as well as that of N3. The chemical shift tensor of N5 is over 720 ppm wide, in accordance with the predictions of theory and our calculations. The signal of N3 can be distinguished on the basis of coupling to 1H absent for N1 and N5, as well as the shift tensor span of only 170 ppm, consistent with N3's lower aromaticity and lack of a nonbonding lone pair. The isotropic shifts and spans of N5 and N1 reflect two opposite extremes of the chemical shift range for "pyridine-type" N's, consistent with their electrophilic and nucleophilic chemical reactivities, respectively. Upon flavin reduction, N5's chemical shift tensor contracts dramatically to a span of less than 110 ppm, and the isotropic chemical shift changes by approximately 300 ppm. Both are consistent with loss of N5's nonbonding lone pair and decreased aromaticity, and illustrate the responsiveness of the 15N chemical shift principal values to electronic structure. Thus. 15N chemical shift principal values promise to be valuable tools for understanding electronic differences that underlie variations in flavin reactivity, as well as the reactivities of other heterocyclic cofactors.

Native-specific stabilization of flavodoxin by the FMN cofactor: structural and thermodynamical explanation

Flavodoxins are useful models to investigate protein/cofactor interactions. The binding energy of the apoflavodoxin-FMN complex is high and therefore the holoflavodoxin is expected to be more stable than the apoprotein. This expectation has been challenged by reports on the stability of Desulfovibrio desulfuricans flavodoxin indicating that FMN binds to the unfolded polypeptide with similar affinity as to the native state, thus causing no net effect on protein stability. In previous work, we have analyzed in detail the stability of the apoflavodoxin from Anabaena PCC 7119 and the energetics of its functional complex with FMN. Here, we use the Anabaena holoprotein to directly investigate the contribution of the bound cofactor to protein stability through a detailed analysis of the chemical and thermal denaturation equilibria. Our data clearly shows that FMN binding largely stabilizes the protein towards both chemical and thermal denaturation, and that the stabilization observed at 25 degrees C in low ionic strength conditions is precisely the one expected if full release of the cofactor takes place upon flavodoxin unfolding. On the other hand, the binding of FMN to the native polypeptide is shown to simplify the thermal unfolding so that, while apoflavodoxin follows a three-state mechanism, the holoprotein unfolds in a two-state fashion. Comparison of the X-ray structure of native apoflavodoxin with the phi-structure of the thermal intermediate indicates that the increase in cooperativity driven by the cofactor originates in its preferential binding to the native state, which is a consequence of the disorganization in the intermediate of the FMN binding loops and of an adjacent longer loop.