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

Redox Properties of Flavin in BLUF and LOV Photoreceptor Proteins from Hybrid QM/MM Molecular Dynamics Simulation

Flavins play an important role in many oxidation and reduction processes in biological systems. For example, flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) are common cofactors found in enzymatic proteins that use the special redox properties of these flavin molecules for their catalytic or photoactive functions. The redox potential of the flavin is strongly affected by its (protein) environment; however, the underlying molecular interactions of this effect are still unknown. Using hybrid quantum mechanics/molecular mechanics (QM/MM) simulation techniques, we have studied the redox properties of flavin in the gas phase, aqueous solution, and two different protein environments, in particular, a BLUF and a LOV photoreceptor domain. By mapping the changes in electrostatic potential and solvent structure, we gain insight into how specific polarization of the flavin by its environment tunes the reduction potential. We find also that accurate calculation of the reduction potentials of these systems by using the hybrid QM/MM approach is hampered by a too limited sampling of the counterion configurations and by artifacts at the QM/MM boundary. We make suggestions for how these issues can be overcome.

Mechanistic insights into iron-sulfur clusters and flavin oxidation of a novel xanthine oxidoreductase from Sulfobacillus acidophilus TPY

Xanthine oxidoreductase (XOR) catalyzes the oxidation of purines (hypoxanthine and xanthine) to uric acid. XOR is widely used in various therapeutic and biotechnological applications. In this study, we characterized the biophysical and mechanistic properties of a novel bacterial XOR from Sulfobacillus acidophilus TPY (SaXOR). Our results showed that SaXOR is a heterotrimer consisting of three subunits, namely XoA, XoB, and XoC, which denote the molybdenum cofactor (Moco), 2Fe-2S, and FAD-binding domains, respectively. XoC was found to be stable when co-expressed with XoB, forming an XoBC complex. Furthermore, we prepared a fusion of XoB and XoC via a flexible linker (fusXoBC) and evaluated its function in comparison to that of XoBC. Spectroscopic analysis revealed that XoB harbors two 2Fe-2S clusters, whereas XoC bears a single-bound FAD cofactor. Electron transfer from reduced forms of XoC, XoBC, and fusXoBC to molecular oxygen (O2 ) during oxidative half-reaction yielded no flavin semiquinones, implying ultrafast single-electron transfer from 2Fe-2Sred to FAD. In the presence of XoA, XoBC and fusXoBC exhibited comparable XoA affinity and exploited a shared overall mechanism. Nonetheless, the linkage may accelerate the two-step, single-electron transfer cascade from 2Fe-2Sred to FAD while augmenting protein stability. Collectively, our findings provide novel insights into SaXOR properties and oxidation mechanisms divergent from prior mammalian and bacterial XOR paradigms.

Difference in the Charge-Separation Energetics between Distinct Conformers in the PixD Photoreceptor

Blue light using flavin (BLUF) domain proteins are photoreceptors in various organisms. The PixD BLUF domain can adopt two conformations, W91out and W91in, with Trp91 either proximal or distal to flavin (FMN). Using a quantum mechanical/molecular mechanical/polarizable continuum model approach, the energetics of charge-separated and biradical states in the two conformations were investigated. In the W91out conformation, the charge-separated state (FMN•-) is more stable than the photoexcited state (FMN*), whereas it is less stable due to an electrostatic repulsive interaction with the Ser28 side chain in the W91in conformation. This leads to a lower activation energy for the charge separation in the W91out conformation, resulting in a faster charge separation compared to that in the W91in conformation. In the W91out conformation, the radical state (FMNH•) is more stable than FMN•- and forms from FMN•-, leading to reorientation of the Gln50 side chain adjacent to FMN and formation of a hydrogen bond between Gln50 and FMN. Subsequently, a signaling state forms through charge recombination. In contrast, in the W91in conformation, FMN•- cannot proceed further, returning to the dark-adapted state, as FMNH• is less stable. Thus, formation of the signaling state exclusively occurs in the W91out conformation.

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.

Roles of Ferredoxin-NADP+ Oxidoreductase and Flavodoxin in NAD(P)H-Dependent Electron Transfer Systems

Distinct isoforms of FAD-containing ferredoxin-NADP⁺ oxidoreductase (FNR) and ferredoxin (Fd) are involved in photosynthetic and non-photosynthetic electron transfer systems. The FNR (FAD)-Fd [2Fe-2S] redox pair complex switches between one- and two-electron transfer reactions in steps involving FAD semiquinone intermediates. In cyanobacteria and some algae, one-electron carrier Fd serves as a substitute for low-potential FMN-containing flavodoxin (Fld) during growth under low-iron conditions. This complex evolves into the covalent FNR (FAD)-Fld (FMN) pair, which participates in a wide variety of NAD(P)H-dependent metabolic pathways as an electron donor, including bacterial sulfite reductase, cytochrome P450 BM3, plant or mammalian cytochrome P450 reductase and nitric oxide synthase isoforms. These electron transfer systems share the conserved Ser-Glu/Asp pair in the active site of the FAD module. In addition to physiological electron acceptors, the NAD(P)H-dependent diflavin reductase family catalyzes a one-electron reduction of artificial electron acceptors such as quinone-containing anticancer drugs. Conversely, NAD(P)H: quinone oxidoreductase (NQO1), which shares a Fld-like active site, functions as a typical two-electron transfer antioxidant enzyme, and the NQO1 and UDP-glucuronosyltransfease/sulfotransferase pairs function as an antioxidant detoxification system. In this review, the roles of the plant FNR-Fd and FNR-Fld complex pairs were compared to those of the diflavin reductase (FAD-FMN) family. In the final section, evolutionary aspects of NAD(P)H-dependent multi-domain electron transfer systems are discussed.

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 (Em) 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 Em value, and when possible, we rationalized them on the basis of previous studies.

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.

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.

Molecular mechanism of metabolic NAD(P)H-dependent electron-transfer systems: The role of redox cofactors

NAD(P)H-dependent electron-transfer (ET) systems require three functional components: a flavin-containing NAD(P)H-dehydrogenase, one-electron carrier and metal-containing redox center. In principle, these ET systems consist of one-, two- and three-components, and the electron flux from pyridine nucleotide cofactors, NADPH or NADH to final electron acceptor follows a linear pathway: NAD(P)H → flavin → one-electron carrier → metal containing redox center. In each step ET is primarily controlled by one- and two-electron midpoint reduction potentials of protein-bound redox cofactors in which the redox-linked conformational changes during the catalytic cycle are required for the domain-domain interactions. These interactions play an effective ET reactions in the multi-component ET systems. The microsomal and mitochondrial cytochrome P450 (cyt P450) ET systems, nitric oxide synthase (NOS) isozymes, cytochrome b₅ (cyt b₅) ET systems and methionine synthase (MS) ET system include a combination of multi-domain, and their organizations display similarities as well as differences in their components. However, these ET systems are sharing of a similar mechanism. More recent structural information obtained by X-ray and cryo-electron microscopy (cryo-EM) analysis provides more detail for the mechanisms associated with multi-domain ET systems. Therefore, this review summarizes the roles of redox cofactors in the metabolic ET systems on the basis of one-electron redox potentials. In final Section, evolutionary aspects of NAD(P)H-dependent multi-domain ET systems will be discussed.

Electron transfer pathways in a multiheme cytochrome MtrF

In MtrF, an outer-membrane multiheme cytochrome, the 10 heme groups are arranged in heme binding domains II and IV along the pseudo-C₂ axis, forming the electron transfer (ET) pathways. Previous reports based on molecular dynamics simulations showed that the redox potential (E_m) values for the heme pairs located in symmetrical positions in domains II and IV were similar, forming bidirectional ET pathways [Breuer M, Zarzycki P, Blumberger J, Rosso KM (2012) J Am Chem Soc 134(24):9868-9871]. Here, we present the E_m values of the 10 hemes in MtrF, solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable residues and heme propionic groups. In contrast to previous studies, the E_m values indicated that the ET is more likely to be downhill from domain IV to II because of localization of acidic residues in domain IV. Reduction of hemes in MtrF lowered the E_m values, resulting in switching to alternative downhill ET pathways that extended to the flavin binding sites. These findings present an explanation of how MtrF serves as an electron donor to extracellular substrates.

Proton transfer reactions and hydrogen-bond networks in protein environments

In protein environments, proton transfer reactions occur along polar or charged residues and isolated water molecules. These species consist of H-bond networks that serve as proton transfer pathways; therefore, thorough understanding of H-bond energetics is essential when investigating proton transfer reactions in protein environments. When the pKa values (or proton affinity) of the H-bond donor and acceptor moieties are equal, significantly short, symmetric H-bonds can be formed between the two, and proton transfer reactions can occur in an efficient manner. However, such short, symmetric H-bonds are not necessarily stable when they are situated near the protein bulk surface, because the condition of matching pKa values is opposite to that required for the formation of strong salt bridges, which play a key role in protein-protein interactions. To satisfy the pKa matching condition and allow for proton transfer reactions, proteins often adjust the pKa via electron transfer reactions or H-bond pattern changes. In particular, when a symmetric H-bond is formed near the protein bulk surface as a result of one of these phenomena, its instability often results in breakage, leading to large changes in protein conformation.

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.

Redox modulation of flavin and tyrosine determines photoinduced proton-coupled electron transfer and photoactivation of BLUF photoreceptors

Photoinduced electron transfer in biological systems, especially in proteins, is a highly intriguing matter. Its mechanistic details cannot be addressed by structural data obtained by crystallography alone because this provides only static information on a given redox system. In combination with transient spectroscopy and site-directed manipulation of the protein, however, a dynamic molecular picture of the ET process may be obtained. In BLUF (blue light sensors using FAD) photoreceptors, proton-coupled electron transfer between a tyrosine and the flavin cofactor is the key reaction to switch from a dark-adapted to a light-adapted state, which corresponds to the biological signaling state. Particularly puzzling is the fact that, although the various naturally occurring BLUF domains show little difference in the amino acid composition of the flavin binding pocket, the reaction rates of the forward reaction differ quite largely from a few ps up to several hundred ps. In this study, we modified the redox potential of the flavin/tyrosine redox pair by site-directed mutagenesis close to the flavin C2 carbonyl and fluorination of the tyrosine, respectively. We provide information on how changes in the redox potential of either reaction partner significantly influence photoinduced proton-coupled electron transfer. The altered redox potentials allowed us furthermore to experimentally describe an excited state charge transfer intermediately prior to electron transfer in the BLUF photocycle. Additionally, we show that the electron transfer rate directly correlates with the quantum yield of signaling state formation.

Structure-activity characterization of sulfide:quinone oxidoreductase variants

Sulfide:quinone oxidoreductase (SQR) is a peripheral membrane protein that catalyzes the oxidation of sulfide species to elemental sulfur. The enzymatic reaction proceeds in two steps. The electrons from sulfides are transferred first to the enzyme cofactor, FAD, which, in turn, passes them onto the quinone pool in the membrane. Several wild-type SQR structures have been reported recently. However, the enzymatic mechanism of SQR has not been fully delineated. In order to understand the role of the catalytically essential residues in the enzymatic mechanism of SQR we produced a number of variants of the conserved residues in the catalytic site including the cysteine triad of SQR from the acidophilic, chemolithotrophic bacterium Acidithiobacillus ferrooxidans. These were structurally characterized and their activities for each reaction step were determined. In addition, the crystal structures of the wild-type SQR with sodium selenide and gold(I) cyanide have been determined. Previously we proposed a mechanism for the reduction of sulfides to elemental sulfur involving nucleophilic attack of Cys356 on C(4A) atom of FAD. Here we also consider an alternative anionic radical mechanism by direct electron transfer from Cys356 to the isoalloxazine ring of FAD.

Protein conformational gating of enzymatic activity in xanthine oxidoreductase

In mammals, xanthine oxidoreductase can exist as xanthine dehydrogenase (XDH) and xanthine oxidase (XO). The two enzymes possess common redox active cofactors, which form an electron transfer (ET) pathway terminated by a flavin cofactor. In spite of identical protein primary structures, the redox potential difference between XDH and XO for the flavin semiquinone/hydroquinone pair (E(sq/hq)) is ~170 mV, a striking difference. The former greatly prefers NAD(+) as ultimate substrate for ET from the iron-sulfur cluster FeS-II via flavin while the latter only accepts dioxygen. In XDH (without NAD(+)), however, the redox potential of the electron donor FeS-II is 180 mV higher than that for the acceptor flavin, yielding an energetically uphill ET. On the basis of new 1.65, 2.3, 1.9, and 2.2 Å resolution crystal structures for XDH, XO, the NAD(+)- and NADH-complexed XDH, E(sq/hq) were calculated to better understand how the enzyme activates an ET from FeS-II to flavin. The majority of the E(sq/hq) difference between XDH and XO originates from a conformational change in the loop at positions 423-433 near the flavin binding site, causing the differences in stability of the semiquinone state. There was no large conformational change observed in response to NAD(+) binding at XDH. Instead, the positive charge of the NAD(+) ring, deprotonation of Asp429, and capping of the bulk surface of the flavin by the NAD(+) molecule all contribute to altering E(sq/hq) upon NAD(+) binding to XDH.

Proton-binding sites of acid-sensing ion channel 1

Acid-sensing ion channels (ASICs) are proton-gated cation channels that exist throughout the mammalian central and peripheral nervous systems. ASIC1 is the most abundant of all the ASICs and is likely to modulate synaptic transmission. Identifying the proton-binding sites of ASCI1 is required to elucidate its pH-sensing mechanism. By using the crystal structure of ASIC1, the protonation states of each titratable site of ASIC1 were calculated by solving the Poisson-Boltzmann equation under conditions wherein the protonation states of all these sites are simultaneously in equilibrium. Four acidic-acidic residue pairs—Asp238-Asp350, Glu220-Asp408, Glu239-Asp346, and Glu80-Glu417—were found to be highly protonated. In particular, the Glu80-Glu417 pair in the inner pore was completely protonated and possessed 2 H⁺, implying its possible importance as a proton-binding site. The pK_a of Glu239, which forms a pair with a possible pH-sensing site Asp346, differs among each homo-trimer subunit due to the different H-bond pattern of Thr237 in the different protein conformations of the subunits. His74 possessed a pK_a of ≈6–7. Conservation of His74 in the proton-sensitive ASIC3 that lacks a residue corresponding to Asp346 may suggest its possible pH-sensing role in proton-sensitive ASICs.

High-resolution crystal structures of the flavoprotein NrdI in oxidized and reduced states--an unusual flavodoxin. Structural biology

The small flavoprotein NrdI is an essential component of the class Ib ribonucleotide reductase system in many bacteria. NrdI interacts with the class Ib radical generating protein NrdF. It is suggested to be involved in the rescue of inactivated diferric centres or generation of active dimanganese centres in NrdF. Although NrdI bears a superficial resemblance to flavodoxin, its redox properties have been demonstrated to be strikingly different. In particular, NrdI is capable of two-electron reduction, whereas flavodoxins are exclusively one-electron reductants. This has been suggested to depend on a lesser destabilization of the negatively-charged hydroquinone state than in flavodoxins. We have determined the crystal structures of NrdI from Bacillus anthracis, the causative agent of anthrax, in the oxidized and semiquinone forms, at resolutions of 0.96 and 1.4 Å, respectively. These structures, coupled with analysis of all curated NrdI sequences, suggest that NrdI defines a new structural family within the flavodoxin superfamily. The conformational behaviour of NrdI in response to FMN reduction is very similar to that of flavodoxins, involving a peptide flip in a loop near the N5 atom of the flavin ring. However, NrdI is much less negatively charged than flavodoxins, which is expected to affect its redox properties significantly. Indeed, sequence analysis shows a remarkable spread in the predicted isoelectric points of NrdIs, from approximately pH 4-10. The implications of these observations for class Ib ribonucleotide reductase function are discussed.

Mechanism-based tuning of a LOV domain photoreceptor

Phototropin-like LOV domains form a cysteinyl-flavin adduct in response to blue light but show considerable variation in output signal and the lifetime of the photo-adduct signaling state. Mechanistic studies of the slow-cycling fungal LOV photoreceptor Vivid (VVD) reveal the importance of reactive cysteine conformation, flavin electronic environment and solvent accessibility for adduct scission and thermal reversion. Proton inventory, pH effects, base catalysis and structural studies implicate flavin N(5) deprotonation as rate-determining for recovery. Substitutions of active site residues Ile74, Ile85, Met135 and Met165 alter photoadduct lifetimes by over four orders of magnitude in VVD, and similar changes in other LOV proteins show analogous effects. Adduct state decay rates also correlate with changes in conformational and oligomeric properties of the protein necessary for signaling. These findings link natural sequence variation of LOV domains to function and provide a means to design broadly reactive light-sensitive probes.

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.

Light-induced hydrogen bonding pattern and driving force of electron transfer in AppA BLUF domain photoreceptor

The AppA BLUF (blue light sensing using FAD) domain from Rhodobacter sphaeroides serves as a blue light-sensing photoreceptor. The charge separation process between Tyr-21 and flavin plays an important role in the light signaling state by transforming the dark state conformation to the light state one. By solving the linearized Poisson-Boltzmann equation, I calculated E(m) for Tyr-21, flavin, and redox-active Trp-104 and revealed the electron transfer (ET) driving energy. Rotation of the Gln-63 side chain that converts protein conformation from the dark state to the light state is responsible for the decrease of 150 mV in E(m) for Tyr-21, leading to the significantly larger ET driving energy in the light state conformation. The pK(a) values of protonation for flavin anions are essentially the same in both dark and light state crystal structures. In contrast to the ET via Tyr-21, formation of the W state results in generation of only the dark state conformation (even if the initial conformation is in the light state); this could explain why Trp-104-mediated ET deactivates the light-sensing yield and why the activity of W104A mutant is similar to that of the light-adapted native BLUF.

Redox potential difference between Desulfovibrio vulgaris and Clostridium beijerinckii flavodoxins

The redox potential of the flavin mononucleotide (FMN) hydroquinones for one-electron reduction in the Desulfovibrio vulgaris ( D. vulgaris) flavodoxin ( E sq/hq for FMNH (*)/FMNH (-)) was calculated using the crystal structure of the relevant hydroquinone form and compared to the results of the Clostridium beijerinckii ( C. beijerinckii) flavodoxin. In D. vulgaris and C. beijerinckii flavodoxins, the protein side chain causes significant downshifts of 170 and 240 mV in E sq/hq, respectively. In the C. beijerinckii flavodoxin, the E sq/hq downshift because of the protein side chain is essentially compensated by the counter influence of the protein backbone ( E sq/hq upshift of 260 mV). However, in the D. vulgaris flavodoxin, the corresponding protein backbone influence on E sq/hq is significantly small, i.e., less than half of that in the C. beijerinckii flavodoxin. In particular, there is a significant difference in the influence of the protein backbone of the so-called 60s loop region between the two flavodoxins. The E sq/hq difference can be best explained by the lower compensation of the side chain influence by the backbone influence in the D. vulgaris flavodoxin than in the C. beijerinckii flavodoxin.

Modulation of the protein environment in the hydrophilic pore of the ammonia transporter protein AmtB upon GlnK protein binding

The conduction of ammonia/ammonium (NH3/NH4(+)) through the channel protein AmtB is inhibited by the binding of the signal transduction protein GlnK. In the AmtB-GlnK binding interface, there exists an NH3/NH4(+) binding site--Am6. The calculated pK(a) values at the Am6 sites in both the AmtB-GlnK complex and isolated AmtB implies the dominance of an uncharged NH3 state. The GlnK protein binding causes a significant downshift in the Am6 pK(a) value of the AmtB. However, this downshift is perfectly compensated by the reorientation of the protein backbone (carbonyl group of Cys312 from the AmtB part) upon AmtB-GlnK complex formation.