Common conformational changes in flavodoxins induced by FMN and anion binding: the structure of Helicobacter pylori apoflavodoxin

Flavodoxins, noncovalent complexes between apoflavodoxins and flavin mononucleotide (FMN), are useful models to investigate the mechanism of protein/flavin recognition. In this respect, the only available crystal structure of an apoflavodoxin (that from Anabaena) showed a closed isoalloxazine pocket and the presence of a bound phosphate ion, which posed many questions on the recognition mechanism and on the potential physiological role exerted by phosphate ions. To address these issues we report here the X-ray structure of the apoflavodoxin from the pathogen Helicobacter pylori. The protein naturally lacks one of the conserved aromatic residues that close the isoalloxazine pocket in Anabaena, and the structure has been determined in a medium lacking phosphate. In spite of these significant differences, the isoallozaxine pocket in H. pylori apoflavodoxin appears also closed and a chloride ion is bound at a native-like FMN phosphate site. It seems thus that it is a general characteristic of apoflavodoxins to display closed, non-native, isoalloxazine binding sites together with native-like, rather promiscuous, phosphate binding sites that can bear other available small anions present in solution. In this respect, both binding energy hot spots of the apoflavodoxin/FMN complex are initially unavailable to FMN binding and the specific spot for FMN recognition may depend on the dynamics of the two candidate regions. Molecular dynamics simulations show that the isoalloxazine binding loops are intrinsically flexible at physiological temperatures, thus facilitating the intercalation of the cofactor, and that their mobility is modulated by the anion bound at the phosphate site.     

Dissecting the energetics of the apoflavodoxin-FMN complex

Many flavoproteins are non-covalent complexes between FMN and an apoprotein. To understand better the stability of flavoproteins, we have studied the energetics of the complex between FMN and the apoflavodoxin from Anabaena PCC 7119 by a combination of site-directed mutagenesis, titration calorimetry, equilibrium binding constant determinations, and x-ray crystallography. Comparison of the strength of the wild type and mutant apoflavodoxin-FMN complexes and that of the complexes between wild type apoflavodoxin and shortened FMN analogues (riboflavin and lumiflavin) allows the dissection of the binding energy into contributions associated with the different parts of the FMN molecule. The estimated contribution of the phosphate is greatest, at 7 kcal mol(-1); that of the isoalloxazine is of around 5-6 kcal mol(-1) (mainly due to interaction with Trp-57 and Tyr-94 in the apoprotein) and the ribityl contributes least: around 1 kcal mol(-1). The stabilization of the complex is both enthalpic and entropic although the enthalpy contribution is dominant. Both the phosphate and the isoalloxazine significantly contribute to the enthalpy of binding. The ionic strength does not have a large effect on the stability of the FMN complex because, although it weakens the phosphate interactions, it strengthens the isoalloxazine-protein hydrophobic interactions. Phosphate up to 100 mM does not affect the strength of the riboflavin complex, which suggests the isoalloxazine and phosphate binding sites may be independent in terms of binding energy. Interestingly, we find crystallographic evidence of flexibility in one of the loops (57-62) involved in isoalloxazine binding.     

Mechanism of FMN Binding to the Apoflavodoxin from Helicobacter pylori

Flavodoxins are bacterial electron transport proteins whose redox competence is due to the presence of a tightly but noncovalently bound FMN molecule. While the thermodynamics of the complex are understood, the mechanism of association between the apoflavodoxin and the redox cofactor is not so clear. We investigate here the mechanism of FMN binding to the apoflavodoxin from Helicobacter pylori, an essential protein that is being used as a target to develop antimicrobials. This flavodoxin is structurally peculiar as it lacks the typical bulky residue interacting with the FMN re face but bears instead a small alanine. FMN binding is biphasic, regardless of the presence of phosphate molecules in solution, while riboflavin binding takes place in a single step, the rate constant of which coincides with the fast phase of FMN binding. A mutational study at the isoalloxazine and phosphate subsites for FMN binding clearly indicates that FMN association is always limited by interaction with the isoalloxazine subsite because mutating residues that interact with the phosphate moiety of FMN in the native complex hardly changes the observed rate constants and amplitudes. In contrast, replacing tyr92, which interacts with the isoalloxazine, greatly lowers the rate constants. Our analysis indicates that the two FMN binding phases observed are related neither with alternative or sequential interaction with the two binding subsites nor with the presence of bound phosphate. It is possible that they reflect the intrinsic conformational heterogeneity of the apoflavodoxin ensemble.     

Can cofactor-binding sites in proteins be flexible? Desulfovibrio desulfuricans flavodoxin binds FMN dimer

Flavodoxins catalyze redox reactions using the isoalloxazine moiety of the flavin mononucleotide (FMN) cofactor stacked between two aromatic residues located in two peptide loops. At high FMN concentrations that favor stacked FMN dimers in solution, isothermal titration calorimetric studies show that these dimers bind strongly to apo-flavodoxin from Desulfovibrio desulfuricans (30 degrees C, 20 mM Hepes, pH 7, K(D) = 5.8 microM). Upon increasing the temperature so the FMN dimers dissociate (as shown by (1)H NMR), only one-to-one (FMN-to-protein) binding is observed. Calorimetric titrations result in one-to-one binding also in the presence of phosphate or sulfate (30 degrees C, 13 mM anion, pH 7, K(D) = 0.4 microM). FMN remains dimeric in the presence of phosphate and sulfate, suggesting that specific binding of a divalent anion to the phosphate-binding site triggers ordering of the peptide loops so only one isoalloxazine can fit. Although the physiological relevance of FMN and other nucleotides as dimers has not been explored, our study shows that high-affinity binding to proteins of such dimers can occur in vitro. This emphasizes that the cofactor-binding site in flavodoxin is more flexible than previously expected.     

Phosphate group binding "cup" of PLP-dependent and non-PLP-dependent enzymes: leitmotif and variations

Pyridoxal-5'-phosphate (PLP) is widely used by many enzymes in reactions where amino acids are interconverted. Whereas the role of the pyridoxal ring in catalysis is well understood, the functional role of the single phosphate group in PLP has been less studied. Here we construct unambiguous connection diagrams that describe the interactions among the three non-ester phosphate oxygen atoms of PLP and surrounding atoms from the protein binding site and from water molecules, the so-called phosphate group binding "cup". These diagrams provide a simple means to identify common recognition motifs for the phosphate group in both similar and different protein folds. Diagrams were constructed and compared in the cases of five newly determined structures of PLP-dependent transferases (fold type I enzymes) and, additionally, two non-PLP protein complexes (indole-3-glycerol phosphate synthase (IGPS) with bound indole-3-glycerol phosphate (IGP) and old yellow enzyme (OYE) with bound flavin mononucleotide (FMN)). A detailed comparison of the diagrams shows that three positions out of ten in the structure of the phosphate group binding "cup" contain invariant atoms, while seven others are occupied by conserved atom types. This level of similarity was also observed in the fold type III (TIM beta/alpha-barrel) enzymes that bind three different ligands: PLP, IGP and FMN.     

Ligand-induced conformational changes in the capping subdomain of a bacterial old yellow enzyme homologue and conserved sequence fingerprints provide new insights into substrate binding

We have recently reported that Shewanella oneidensis, a Gram-negative gamma-proteobacterium with a rich arsenal of redox proteins, possesses four old yellow enzyme (OYE) homologues. Here, we report a series of high resolution crystal structures for one of these OYEs, Shewanella yellow enzyme 1 (SYE1), in its oxidized form at 1.4A resolution, which binds a molecule of PEG 400 in the active site, and in its NADH-reduced and p-hydroxybenzaldehyde- and p-hydroxyacetophenone-bound forms at 1.7A resolution. Although the overall structure of SYE1 reveals a monomeric enzyme based on the alpha(8)beta(8) barrel scaffold observed for other OYEs, the active site exhibits a unique combination of features: a strongly butterfly-bent FMN cofactor both in the oxidized and NADH-reduced forms, a collapsed and narrow active site tunnel, and a novel combination of conserved residues involved in the binding of phenolic ligands. Furthermore, we identify a second p-hydroxybenzaldehyde-binding site in a hydrophobic cleft next to the entry of the active site tunnel in the capping subdomain, formed by a restructuring of Loop 3 to an "open" conformation. This constitutes the first evidence to date for the entire family of OYEs that Loop 3 may indeed play a dynamic role in ligand binding and thus provides insights into the elusive NADH complex and into substrate binding in general. Structure-based sequence alignments indicate that the novelties we observe in SYE1 are supported by conserved residues in a number of structurally uncharacterized OYEs from the beta- and gamma-proteobacteria, suggesting that SYE1 represents a new subfamily of bacterial OYEs.     

Control of phosphorylase b conformation by a modified cofactor: crystallographic studies on R-state glycogen phosphorylase reconstituted with pyridoxal 5''-diphosphate.

Previous crystallographic studies on glycogen phosphorylase have described the different conformational states of the protein (T and R) that represent the allosteric transition and have shown how the properties of the 5'-phosphate group of the cofactor pyridoxal phosphate are influenced by these conformational states. The present work reports a study on glycogen phosphorylase b (GPb) complexed with a modified cofactor, pyridoxal 5'-diphosphate (PLPP), in place of the natural cofactor. Solution studies (Withers, S.G., Madsen, N.B., & Sykes, B.D., 1982, Biochemistry 21, 6716-6722) have shown that PLPP promotes R-state properties of the enzyme indicating that the cofactor can influence the conformational state of the protein. GPb complexed with pyridoxal 5'-diphosphate (PLPP) has been crystallized in the presence of IMP and ammonium sulfate in the monoclinic R-state crystal form and the structure refined from X-ray data to 2.8 A resolution to a crystallographic R value of 0.21. The global tertiary and quaternary structure in the vicinity of the Ser 14 and the IMP sites are nearly identical to those observed for the R-state GPb-AMP complex. At the catalytic site the second phosphate of PLPP is accommodated with essentially no change in structure from the R-state structure and is involved in interactions with the side chains of two lysine residues (Lys 568 and Lys 574) and the main chain nitrogen of Arg 569. Superposition of the T-state structure shows that were the PLPP to be incorporated into the T-state structure there would be a close contact with the 280s loop (residues 282-285) that would encourage the T to R allosteric transition. The second phosphate of the PLPP occupies a site that is distinct from other dianionic binding sites that have been observed for glucose-1-phosphate and sulfate (in the R state) and for heptulose-2-phosphate (in the T state). The results indicate mobility in the dianion recognition site, and the precise position is dependent on other linkages to the dianion. In the modified cofactor the second phosphate site is constrained by the covalent link to the first phosphate of PLPP. The observed position in the crystal suggests that it is too far from the substrate site to represent a site for catalysis.     

Insights into the mode of flavin mononucleotide binding and catalytic mechanism of bacterial chromate reductases: A molecular dynamics simulation study

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

Structural basis for recognition of polyglutamyl folates by thymidylate synthase.

Thymidylate synthase (TS) catalyzes the final step in the de novo synthesis of thymidine. In vivo TS binds a polyglutamyl cofactor, polyglutamyl methylenetetrahydrofolate (CH2-H4folate), which serves as a carbon donor. Glutamate residues on the cofactor contribute as much as 3.7 kcal to the interaction between the cofactor, substrate, and enzyme. Because many ligand/receptor interactions appear to be driven largely by hydrophobic forces, it is surprising that the addition of hydrophilic, soluble groups such as glutamates increases the affinity of the cofactor for TS. The structure of a polyglutamyl cofactor analog bound in ternary complex with deoxyuridine monophosphate (dUMP) and Escherichia coli TS reveals how the polyglutamyl moiety is positioned in TS and accounts in a qualitative way for the binding contributions of the different individual glutamate residues. The polyglutamyl moiety is not rigidly fixed by its interaction with the protein except for the first glutamate residue nearest the p-aminobenzoic acid ring of folate. Each additional glutamate is progressively more disordered than the previous one in the chain. The position of the second and third glutamate residues on the protein surface suggests that the polyglutamyl binding site could be utilized by a new family of inhibitors that might fill the binding area more effectively than polyglutamate.     

Folding of Desulfovibrio desulfuricans flavodoxin is accelerated by cofactor fly-casting

Folding of cofactor-binding proteins involves ligand binding in addition to polypeptide folding. We here assess the kinetic folding/binding landscape for Desulfovibrio desulfuricans flavodoxin that coordinates an FMN cofactor. The apo-form folds in a two-step process involving a burst-phase intermediate. Studies on Tyr98Ala and Trp60Ala variants reveal that these aromatics-that stack with the FMN in the holo-form-are not participating in the apo-protein folding pathway. However, these residues are essential for FMN interactions with the unfolded protein during refolding of holo-flavodoxin. Unfolding of wild-type holo-flavodoxin is coupled to FMN dissociation whereas for Tyr98Ala and Trp60Ala holo-variants, FMN dissociates before polypeptide unfolding. Both variants refold as apo-proteins before FMN rebinds. In sharp contrast, refolding of unfolded wild-type holo-flavodoxin is over an order of magnitude faster than that of the apo-form, the pathway does not include a burst-phase intermediate, and the speed is independent of FMN excess ratio. These observations demonstrate that FMN binds rapidly to the unfolded polypeptide and guides folding straight to the native state. As this path to functional D. desulfuricans holo-flavodoxin is faster than if the cofactor binds to pre-folded apo-protein, this is one of few examples where molecular recognition via a "fly-casting" mechanism is kinetically favored.     

Crystal structure of enoyl-acyl carrier protein reductase (FabK) from Streptococcus pneumoniae reveals the binding mode of an inhibitor

Enoyl-acyl carrier protein (ACP) reductases are critical for bacterial type II fatty acid biosynthesis and thus are attractive targets for developing novel antibiotics. We determined the crystal structure of enoyl-ACP reductase (FabK) from Streptococcus pneumoniae at 1.7 A resolution. There was one dimer per asymmetric unit. Each subunit formed a triose phosphate isomerase (TIM) barrel structure, and flavin mononucleotide (FMN) was bound as a cofactor in the active site. The overall structure was similar to the enoyl-ACP reductase (ER) of fungal fatty acid synthase and to 2-nitropropane dioxygenase (2-ND) from Pseudomonas aeruginosa, although there were some differences among these structures. We determined the crystal structure of FabK in complex with a phenylimidazole derivative inhibitor to envision the binding site interactions. The crystal structure reveals that the inhibitor binds to a hydrophobic pocket in the active site of FabK, and this is accompanied by induced-fit movements of two loop regions. The thiazole ring and part of the ureido moiety of the inhibitor are involved in a face-to-face pi-pi stacking interaction with the isoalloxazine ring of FMN. The side-chain conformation of the proposed catalytic residue, His144, changes upon complex formation. Lineweaver-Burk plots indicate that the inhibitor binds competitively with respect to NADH, and uncompetitively with respect to crotonoyl coenzyme A. We propose that the primary basis of the inhibitory activity is competition with NADH for binding to FabK, which is the first step of the two-step ping-pong catalytic mechanism.     

Probing the nature of interactions in SH2 binding interfaces--evidence from electrospray ionization mass spectrometry.

We have adopted nanoflow electrospray ionization mass spectrometry (ESI-MS) and isothermal titration calorimetry (ITC) to probe the mechanism of peptide recognition by the SH2 domain from the Src family tyrosine kinase protein, Fyn. This domain is involved in the mediation of intracellular signal transduction pathways by interaction with proteins containing phosphorylated tyrosine (Y*) residues. The binding of tyrosyl phosphopeptides can mimic these interactions. Specificity in these interactions has been attributed to the interaction of the Y* and residues proximal and C-terminal to it. Previous studies have established that for specific binding with Fyn, the recognition sequence consists of pTyr-Glu-Glu-Ile. The specific interactions involve the binding of Y* with the ionic, and the Y* + 3 Ile residue with the hydrophobic binding pockets on the surface of the Fyn SH2 domain. In this work, a variation in the Y* + 3 residue of this high-affinity sequence was observed to result in changes in the relative binding affinities as determined in solution (ITC) and in the gas phase (nanoflow ESI-MS). X-ray analysis shows that a feature of the Src family SH2 domains is the involvement of water molecules in the peptide binding site. Under the nanoflow ESI conditions, water molecules appear to be maintained in the Fyn SH2-ligand complex. Compelling evidence for these molecules being incorporated in the SH2-peptide interface is provided by the prevalence of the peaks assigned to water-bound over the water-free complex at high-energy conditions. Thus, the stability of water protein-ligand complex appears to be intimately linked to the presence of water.     

Towards a new therapeutic target: Helicobacter pylori flavodoxin

Helicobacter pylori flavodoxin is the electronic acceptor of the pyruvate-oxidoreductase complex (POR) that catalyzes pyruvate oxidative decarboxilation. Inactivation of this metabolic route precludes bacterial survival. Because flavodoxin is not present in the human host, substances interfering electronic transport from POR might be well suited for eradication therapies against the bacterium. H. pylori flavodoxin presents a peculiar cofactor (FMN) binding site, compared to other known flavodoxins, where a conserved aromatic residue is replaced by alanine. A cavity thus appears under the cofactor that can be filled with small organic molecules. We have cloned H. pylori fldA gene, expressed the protein in Escherichia coli and characterized the purified flavodoxin. Thermal up-shift assays of flavodoxin with different concentrations of benzylamine, as well as fluorescence titration experiments indicate benzylamine binds in the pocket near the FMN binding site. It seems thus that low affinity inhibitors of H. pylori flavodoxin can be easily found that, after improvement, may give rise to leads.     

Crystal structures of human prostatic acid phosphatase in complex with a phosphate ion and alpha-benzylaminobenzylphosphonic acid update the mechanistic picture and offer new insights into inhibitor design

The X-ray crystal structure of human prostatic acid phosphatase (PAP) in complex with a phosphate ion has been determined at 2.4 A resolution. This structure offers a snapshot of the final intermediate in the catalytic mechanism and does not support the role of Asp 258 as a proton donor in catalysis. A total of eight hydrogen bonds serve to strongly bind the phosphate ion within the active site. Bound PEG molecules from the crystallization matrix have allowed the identification of a channel within the molecule that likely plays a role in molecular recognition and in macromolecular substrate selectivity. Additionally, the structure of PAP in complex with a phosphate derivative, alpha-benzylaminobenzylphosphonic acid, a potent inhibitor (IC(50) = 4 nM), has been determined to 2.9 A resolution. This structure gives new insight into the determinants of binding hydrophobic ligands within the active site and allows us to explain PAP's preference for aromatic substrates.     

Characterization of the acetate binding pocket in the Methanosarcina thermophila acetate kinase

Acetate kinase catalyzes the reversible magnesium-dependent synthesis of acetyl phosphate by transfer of the ATP gamma-phosphoryl group to acetate. Inspection of the crystal structure of the Methanosarcina thermophila enzyme containing only ADP revealed a solvent-accessible hydrophobic pocket formed by residues Val(93), Leu(122), Phe(179), and Pro(232) in the active site cleft, which identified a potential acetate binding site. The hypothesis that this was a binding site was further supported by alignment of all acetate kinase sequences available from databases, which showed strict conservation of all four residues, and the recent crystal structure of the M. thermophila enzyme with acetate bound in this pocket. Replacement of each residue in the pocket produced variants with K(m) values for acetate that were 7- to 26-fold greater than that of the wild type, and perturbations of this binding pocket also altered the specificity for longer-chain carboxylic acids and acetyl phosphate. The kinetic analyses of variants combined with structural modeling indicated that the pocket has roles in binding the methyl group of acetate, influencing substrate specificity, and orienting the carboxyl group. The kinetic analyses also indicated that binding of acetyl phosphate is more dependent on interactions of the phosphate group with an unidentified residue than on interactions between the methyl group and the hydrophobic pocket. The analyses also indicated that Phe(179) is essential for catalysis, possibly for domain closure. Alignments of acetate kinase, propionate kinase, and butyrate kinase sequences obtained from databases suggested that these enzymes have similar catalytic mechanisms and carboxylic acid substrate binding sites.     

Crystal structures of human mitochondrial branched chain aminotransferase reaction intermediates: ketimine and pyridoxamine phosphate forms

The three-dimensional structures of the isoleucine ketimine and the pyridoxamine phosphate forms of human mitochondrial branched chain aminotransferase (hBCATm) have been determined crystallographically at 1.9 A resolution. The hBCATm-catalyzed transamination can be described in molecular terms together with the earlier solved pyridoxal phosphate forms of the enzyme. The active site lysine, Lys202, undergoes large conformational changes, and the pyridine ring of the cofactor tilts by about 18 degrees during catalysis. A major determinant of the enzyme's substrate and stereospecificity for L-branched chain amino acids is a group of hydrophobic residues that form three hydrophobic surfaces and lock the side chain in place. Short-chain aliphatic amino acid side chains are unable to interact through van der Waals contacts with any of the surfaces whereas bulky aromatic side chains would result in significant steric hindrance. As shown by modeling, and in agreement with previous biochemical data, glutamate but not aspartate can form hydrogen bond interactions. The carboxylate group of the bound isoleucine is on the same side as the phosphate group of the cofactor. These active site interactions are largely retained in a model of the human cytosolic branched chain aminotransferase (hBCATc), suggesting that residues in the second tier of interactions are likely to determine the specificity of hBCATc for the drug gabapentin. Finally, the structures reveal a unique role for cysteine residues in the mammalian BCAT. Cys315 and Cys318, which immediately follow a beta-turn (residues 311-314) and are located just outside the active site, form an unusual thiol-thiolate hydrogen bond. This beta-turn positions Thr313 for its interaction with the pyridoxal phosphate oxygens and substrate alpha-carboxylate group.     

Trapping conformational states along ligand-binding dynamics of peptide deformylase: the impact of induced fit on enzyme catalysis

For several decades, molecular recognition has been considered one of the most fundamental processes in biochemistry. For enzymes, substrate binding is often coupled to conformational changes that alter the local environment of the active site to align the reactive groups for efficient catalysis and to reach the transition state. Adaptive substrate recognition is a well-known concept; however, it has been poorly characterized at a structural level because of its dynamic nature. Here, we provide a detailed mechanism for an induced-fit process at atomic resolution. We take advantage of a slow, tight binding inhibitor-enzyme system, actinonin-peptide deformylase. Crystal structures of the initial open state and final closed state were solved, as well as those of several intermediate mimics captured during the process. Ligand-induced reshaping of a hydrophobic pocket drives closure of the active site, which is finally "zipped up" by additional binding interactions. Together with biochemical analyses, these data allow a coherent reconstruction of the sequence of events leading from the encounter complex to the key-lock binding state of the enzyme. A "movie" that reconstructs this entire process can be further extrapolated to catalysis.     

Mechanism of flavin mononucleotide cofactor binding to the Desulfovibrio vulgaris flavodoxin. 1. Kinetic evidence for cooperative effects associated with the binding of inorganic phosphate and the 5'-phosphate moiety of the cofactor

The pathway(s) by which the flavin cofactor binds to the apoflavoprotein is the subject of some debate. The crystal and NMR structures of several different flavodoxins have provided some insight, although there is disagreement about the location of the initial interaction between the flavin mononucleotide (FMN) and the apoflavodoxin and the degree of protein conformational change associated with cofactor binding [Genzor, C. G., Perales-Alcon, A., Sancho, J., and Romero, A. (1996) Nat. Struct. Biol. 3, 329-332; Steensma, E., and van Mierlo, C. P. M. (1998) J. Mol. Biol. 282, 653-666]. Binding kinetics using stopped-flow spectrofluorimetry and phosphate competition studies were used to develop a model for flavin binding to the flavodoxin from Desulfovibrio vulgaris. In the presence of phosphate, the time course of fluorescence quenching associated with FMN binding to apoflavodoxin was biphasic, whereas riboflavin, which lacks the 5'-phosphate group of FMN, displayed monophasic binding kinetics. When the concentration of phosphate in solution was increased, the FMN binding rates of the two phases behaved differently; the rate of one phase decreased, while the rate of the other increased. A similar increase in the single phase associated with riboflavin binding was also observed. This has led to the following model. The binding of the flavin isoalloxazine ring to its subsite is dependent on the presence of a phosphate group in the phosphate-binding subsite. When phosphate is in the buffer solution, FMN can bind in either of two ways: by the initial insertion of the 5'-phosphate group followed by ring binding or, when inorganic phosphate from solution is bound, the insertion of the isoalloxazine ring first. Riboflavin, which lacks the phosphate moiety of FMN, binds only in the presence of inorganic phosphate, presumably due to the binding of this group in the phosphate-binding subsite. These results suggest that cooperative interactions exist between the phosphate subsite and the ring-binding region in the D. vulgaris flavodoxin that are necessary for isoalloxazine ring binding.     

Effect of inorganic phosphate on FMN binding and loop flexibility in Desulfovibrio desulfuricans apo-flavodoxin

The complex between flavin mononucleotide (FMN) and apo-flavodoxin is dominated by isoalloxazine-stacking interactions and 5'-phosphate hydrogen bonds. We show here that FMN binding to Desulfovibrio desulfuricans apo-flavodoxin is faster and the affinity is higher in the presence of inorganic phosphate as compared to in its absence (I=110 mM, pH 7, 20 degrees C). The transition-state of complex formation was investigated by phi-value analysis using Trp60Ala and Tyr98Ala apo-flavodoxin variants. We find that Tyr98 is highly involved in the FMN/protein transition state independent of inorganic phosphate, whereas the participation of Trp60 is modulated by inorganic phosphate. The phi-value for Trp60 is higher without phosphate, implying that at this condition stronger binding of Trp60 is required in the transition state to assure successful complex formation. Consistent with the experimental data, all-atom molecular dynamic simulations reveal that the presence of an anion in the phosphate subsite restricts the mobility of the Trp60-containing loop in terms of both backbone and side-chain movements, but has no effect on the Tyr98-containing loop. The overall thermodynamic stability of apo-flavodoxin is higher in the presence of inorganic phosphate as compared to in its absence (I=110 mM, pH 7, 20 degrees C). Kinetic experiments reveal that the additional stability originates in slower unfolding. The combined experimental and computational observations demonstrate that phosphate has an ordering effect on the Trp60-containing loop, which positions Trp60 favorably for FMN binding and increases the barrier for protein unfolding.     

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.