Determination of the role of the Carboxyl-terminal leucine-122 in FMN-binding protein by mutational and structural analysis

Physical quantity of residue electrostatic energy in flavin mononucleotide binding protein dimer

The electrostatic (ES) energy of each residue was for the first time quantitatively evaluated in a flavin mononucleotide binding protein (FBP). A residue electrostatic energy (RES) was obtained as the sum of the ES energies between atoms in each residue and all other atoms in the FBP dimer using atomic coordinates obtained by a molecular dynamics (MD) simulation. ES is one of the most important energies among the interaction energies in a protein. It is determined from the RES, the residues which mainly contribute to stabilize the structure of each subunit, and the binding energy between two subunits can be estimated. The RES of all residues in subunit A (Sub A) and subunit B (Sub B) were attractive forces, even though the residues contain net negative or positive charges. This reveals that the ES energies of any of the residues can contribute to stabilize the protein structure. The total binding ES energy over all residues among the subunits was distributed between -0.2 to -1.2 eV (mean = -0.67 eV) from the MD simulation time.

Flavin Binding to the Deca-heme Cytochrome MtrC: Insights from Computational Molecular Simulation

Certain dissimilatory bacteria have the remarkable ability to use extracellular metal oxide minerals instead of oxygen as terminal electron sinks, using a process known as “extracellular respiration”. Specialized multiheme cytochromes located on the outer membrane of the microbe were shown to be crucial for electron transfer from the cell surface to the mineral. This process is facilitated by soluble, biogenic flavins secreted by the organism for the purpose of acting as an electron shuttle. However, their interactions with the outer-membrane cytochromes are not established on a molecular scale. Here, we study the interaction between the outer-membrane deca-heme cytochrome MtrC from Shewanella oneidensis and flavin mononucleotide (FMN in fully oxidized quinone form) using computational docking. We find that interaction of FMN with MtrC is significantly weaker than with known FMN-binding proteins, but identify a mildly preferred interaction site close to heme 2 with a dissociation constant (K_d) = 490 μM, in good agreement with recent experimental estimates, K_d = 255 μM. The weak interaction with MtrC can be qualitatively explained by the smaller number of hydrogen bonds that the planar headgroup of FMN can form with this protein compared to FMN-binding proteins. Molecular dynamics simulation gives indications for a possible conformational switch upon cleavage of the disulphide bond of MtrC, but without concomitant increase in binding affinities according to this docking study. Overall, our results suggest that binding of FMN to MtrC is reversible and not highly specific, which may be consistent with a role as redox shuttle that facilitates extracellular respiration.

Conformational difference between two subunits in flavin mononucleotide binding protein dimers from Desulfovibrio vulgaris (MF): molecular dynamics simulation

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

Experimental Evidence for a Revision in the Annotation of Putative Pyridoxamine 5'-Phosphate Oxidases P(N/M)P from Fungi

Pyridoxinamine 5'-phosphate oxidases (P(N/M)P oxidases) that bind flavin mononucleotide (FMN) and oxidize pyridoxine 5'-phosphate or pyridoxamine 5'-phosphate to form pyridoxal 5'-phosphate (PLP) are an important class of enzymes that play a central role in cell metabolism. Failure to generate an adequate supply of PLP is very detrimental to most organisms and is often clinically manifested as a neurological disorder in mammals. In this study, we analyzed the function of YLR456W and YPR172W, two homologous genes of unknown function from S. cerevisiae that have been annotated as putative P(N/M)P oxidases based on sequence homology. Different experimental approaches indicated that neither protein catalyzes PLP formation nor binds FMN. On the other hand, our analysis confirmed the enzymatic activity of Pdx3, the S. cerevisiae protein previously implicated in PLP biosynthesis by genetic and structural characterization. After a careful sequence analysis comparing the putative and confirmed P(N/M)P oxidases, we found that the protein domain (PF01243) that led to the YLR456W and YPR172W annotation is a poor indicator of P(N/M)P oxidase activity. We suggest that a combination of two Pfam domains (PF01243 and PF10590) present in Pdx3 and other confirmed P(N/M)P oxidases would be a stronger predictor of this molecular function. This work exemplifies the importance of experimental validation to rectify genome annotation and proposes a revision in the annotation of at least 400 sequences from a wide variety of fungal species that are homologous to YLR456W and are currently misrepresented as putative P(N/M)P oxidases.

Effects of the disappearance of one charge on ultrafast fluorescence dynamics of the FMN binding protein

Crystal structures of E13T (Glu13 was replaced by Thr13) and E13Q (Glu13 was replaced by Gln13) FMN binding proteins (FMN-bp) from Desulfovibrio vulgaris, strain Miyazaki F, were determined by the X-ray diffraction method. Geometrical factors related to photoinduced electron transfer from Trp32, Tyr35, and Trp106 to the excited isoalloxazine (Iso*) were compared among the three forms of FMN-bp. The rate of ET is considered to be fastest from Trp32 to Iso* in FMN-bp and then from Tyr35 and Trp106. The distances between Iso and Trp32 did not change appreciably (0.705-0.712 nm) among WT, E13T, and E13Q FMN-bps, though the distances between Iso and Tyr35 or Trp106 became a little shorter by ca. 0.01 nm in both mutated FMN-bps. The distances between the residue at 13 and the ET donors or acceptor in the mutated proteins, however, changed markedly, compared to WT. Hydrogen bonding pairs and distances between Iso and surrounding amino acids were not modified when Glu13 was replaced by Thr13 or Gln13. Effects of elimination of ionic charge at Glu13 on the ultrafast fluorescence dynamics in E13T and E13Q were investigated comparing to WT, by means of a fluorescence up-conversion method. Fluorescence lifetimes were tau(1) = 107 fs (alpha(1) = 0.86), tau(2) = 475 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.02) in E13T and tau(1) = 134 fs (alpha(1) = 0.85), alpha(2) = 746 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.03) in E13Q, which are compared to the reported lifetimes in WT, tau(1) = 168 fs (alpha(1) = 0.95) and alpha(2) = 1.4 ps (alpha(2) = 0.05). Average lifetimes (tau(AV) = Sigma(i=1)(2or3)alpha(i)tau(i)) were 0.75 ps in E13T, 1.10 ps in E13Q, and 0.23 ps in WT, which implies that tau(AV) was 3.3 times longer in E13T and 4.8 times longer in E13Q, compared to WT. The ultrafast fluorescence dynamics of WT did not change when solvent changed from H(2)O to D(2)O. Static ET rates (inverse of average lifetimes) were analyzed with static structures of the three systems of FMN-bp. Net electrostatic (ES) energies of Iso and Trp32, on which ET rates depend, were 0.0263 eV in WT, 0.322 eV in E13T, and 0.412 eV in E13Q. The calculated ET rates were in excellent agreement with the observed ones in all systems.