Apoflavodoxin: structure, stability, and FMN binding

Flavodoxin in a binary surfactant system consisting of the nonionic 1-decanoyl-rac-glycerol and the zwitterionic lauryldimethylamine-N-oxide: molecular dynamics simulation approach

Due to the short time constant of the spin-spin relaxation process, there is a limitation in the preparation of NMR sample solution for large proteins. To overcome this problem, reverse micelle systems are used.  Here, molecular dynamics simulation was used to study the structure of flavodoxin in a quaternary mixture of 1-decanoyl-rac-glycerol, lauryldimethylamine-N-oxide, pentane and hexanol.  Hexanol was used as co-solvent. Simulations were performed at three different co-solvent concentrations.  The proportion of components in the mixture was selected according to experimental conditions.  For comparison, simulation of flavodoxin in water was also performed.  The simulation results show that the C$$\alpha$$-RMSD for the protein in water is less than for the surfactant mixture.  Also, the radius of gyration of flavodoxin increased in the presence of surfactants.  The distance between the two residues trp-57 and phe-94, as a measure of protein activity, was obtained from the simulations.  The results showed that in the surfactant mixtures this distance increases.  Analysis of the secondary structure of the protein shows that the N-terminal part of the flavodoxin is more affected by surfactants.  The flavodoxin diffusion coefficient in the surfactant mixture decreased in relation to its diffusion coefficient in water.

New Frontiers in Druggability

A powerful early approach to evaluating the druggability of proteins involved determining the hit rate in NMR-based screening of a library of small compounds. Here, we show that a computational analog of this method, based on mapping proteins using small molecules as probes, can reliably reproduce druggability results from NMR-based screening and can provide a more meaningful assessment in cases where the two approaches disagree. We apply the method to a large set of proteins. The results show that, because the method is based on the biophysics of binding rather than on empirical parametrization, meaningful information can be gained about classes of proteins and classes of compounds beyond those resembling validated targets and conventionally druglike ligands. In particular, the method identifies targets that, while not druggable by druglike compounds, may become druggable using compound classes such as macrocycles or other large molecules beyond the rule-of-five limit.

Prediction of flavin mono-nucleotide binding sites using modified PSSM profile and ensemble support vector machine

Flavin mono-nucleotide (FMN) closely evolves in many biological processes. In this study, a computational method was proposed to identify FMN binding sites based on amino acid sequences of proteins only. A modified Position Specific Score Matrix was used to characterize the local environmental sequence information, and a visible improvement of performance was obtained. Also, the ensemble SVM was applied to solve the imbalanced data problem. Additionally, an independent dataset was built to evaluate the practical performance of the method, and a satisfactory accuracy of 87.87% was achieved. It demonstrates that the method is effective in predicting FMN-binding sites.

Underexposed polar residues and protein stabilization

Increasing protein stability is interesting for practical reasons and because it tests our understanding of protein energetics. We explore here the feasibility of stabilizing proteins by replacing underexposed polar residues by apolar ones of similar size and shape. We have compared the stability of wild-type apoflavodoxin with that of a few carefully selected mutants carrying Y → F, Q → L, T → V or K → M replacements. Although a clear inverse correlation between native solvent exposures of replaced polar residues and stability of mutants is observed, most mutations fail to stabilize the protein. The promising exceptions are the two Q → L mutations tested, which characteristically combine the greatest reduction in polar burial with the greatest increase in apolar burial relative to wild type. Analysis of published stability data corresponding to a variety of mutant proteins confirms that, unlike Y → F or T → V replacements, Q → L mutations tend to be stabilizing, and it suggests that N → L mutations might be stabilizing as well. On the other hand, we show that the stability changes associated to the apoflavodoxin mutations can be rationalized in terms of differential polar and apolar burials upon folding plus a generic destabilizing penalty term. Simple equations combining these contributions predict stability changes in a large data set of 113 mutants (Y → F, Q → L or T → V) similarly well as more complex algorithms available on the Internet.

Pseudosymmetry, high copy number and twinning complicate the structure determination of Desulfovibrio desulfuricans (ATCC 29577) flavodoxin

The crystal structure of oxidized flavodoxin from Desulfovibrio desulfuricans (ATCC 29577) was determined by molecular replacement in two crystal forms, P3(1)21 and P4(3), at 2.5 and 2.0 A resolution, respectively. Structure determination in space group P3(1)21 was challenging owing to the presence of pseudo-translational symmetry and a high copy number in the asymmetric unit (8). Initial phasing attempts in space group P3(1)21 by molecular replacement using a poor search model (46% identity) and multi-wavelength anomalous dispersion were unsuccessful. It was necessary to solve the structure in a second crystal form, space group P4(3), which was characterized by almost perfect twinning, in order to obtain a suitable search model for molecular replacement. This search model with complementary approaches to molecular replacement utilizing the pseudo-translational symmetry operators determined by analysis of the native Patterson map facilitated the selection and manual placement of molecules to generate an initial solution in the P3(1)21 crystal form. During the early stages of refinement, application of the appropriate twin law, (-h, -k, l), was required to converge to reasonable R-factor values despite the fact that in the final analysis the data were untwinned and the twin law could subsequently be removed. The approaches used in structure determination and refinement may be applicable to other crystal structures characterized by these complicating factors. The refined model shows flexibility of the flavin mononucleotide coordinating loops indicated by the isolation of two loop conformations and provides a starting point for the elucidation of the mechanism used for protein-partner recognition.

Detection of structural changes in a cofactor binding protein by using a wheat germ cell-free protein synthesis system coupled with unnatural amino acid probing

A cell-free protein synthesis system is a powerful tool with which unnatural amino acids can be introduced into polypeptide chains. Here, the authors describe unnatural amino acid probing in a wheat germ cell-free translation system as a method for detecting the structural changes that occur in a cofactor binding protein on a conversion of the protein from an apo-form to a holo-form. The authors selected the FMN-binding protein from Desulfovibrio vulgaris as a model protein. The apo-form of the protein was synthesized efficiently in the absence of FMN. The purified apo-form could be correctly converted to the holo-form. Thus, the system could synthesize the active apo-form. Gel filtration chromatography, analytical ultracentrifugation, and circular dichroism-spectra studies suggested that the FMN-binding site of the apo-form is open as compared with the holo-form. To confirm this idea, the unnatural amino acid probing was performed by incorporating 3-azido-L-tyrosine at the Tyr35 residue in the FMN-binding site. The authors optimized three steps in their system. The introduced 3-azido-L-tyrosine residue was subjected to specific chemical modification by a fluorescein-triarylphosphine derivative. The initial velocity of the apo-form reaction was 20 fold faster than that of the holo-form, demonstrating that the Tyr35 residue in the apo-form is open to solvent.

Energetics of aliphatic deletions in protein cores

Although core residues can sometimes be replaced by shorter ones without introducing significant changes in protein structure, the energetic consequences are typically large and destabilizing. Many efforts have been devoted to understand and predict changes in stability from analysis of the environment of mutated residues, but the relationships proposed for individual proteins have often failed to describe additional data. We report here 17 apoflavodoxin large-to-small mutations that cause overall protein destabilizations of 0.6-3.9 kcal.mol(-1). By comparing two-state urea and three-state thermal unfolding data, the overall destabilizations observed are partitioned into effects on the N-to-I and on the I-to-U equilibria. In all cases, the equilibrium intermediate exerts a "buffering" effect that reduces the impact of the overall destabilization on the N-to-I equilibrium. The performance of several structure-energetics relationships, proposed to explain the energetics of hydrophobic shortening mutations, has been evaluated by using an apoflavodoxin data set consisting of 14 mutations involving branching-conservative aliphatic side-chain shortenings and a larger data set, including similar mutations implemented in seven model proteins. Our analysis shows that the stability changes observed for any of the different types of mutations (LA, IA, IV, and VA) in either data set are best explained by a combination of differential hydrophobicity and of the calculated volume of the modeled cavity (as previously observed for LA and IA mutations in lysozyme T4). In contrast, sequence conservation within the flavodoxin family, which is a good predictor for charge-reversal stabilizing mutations, does not perform so well for aliphatic shortening ones.

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.

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.

FMN binding and unfolding of Desulfovibrio desulfuricans flavodoxin: "hidden" intermediates at low denaturant concentrations

The flavin mononucleotide (FMN) cofactor in Desulfovibrio desulfuricans flavodoxin stays associated with the polypeptide upon guanidine hydrochloride (GuHCl) induced unfolding. Using isothermal titration calorimetry (ITC), we determined the affinity of FMN for the flavodoxin polypeptide as a function of both urea and GuHCl concentrations (pH 7, 25 degrees C). The FMN affinity for folded and GuHCl-unfolded flavodoxin differs 10-fold, which is in agreement with the difference in thermodynamic stability between the apo- and holo-forms. In contrast, the urea-unfolded protein does not interact with FMN and equilibrium unfolding of holo-flavodoxin in urea results in FMN dissociation prior to polypeptide unfolding. ANS-binding, near-UV circular dichroism (CD), acrylamide quenching and FMN-emission experiments reveal the presence of native-like intermediates, not detected by far-UV CD and aromatic fluorescence detection methods, in low concentrations of both denaturants. Time-resolved experiments show that FMN binding is fastest at GuHCl concentrations where the native-like intermediate species is populated.

A double-deletion method to quantifying incremental binding energies in proteins from experiment: example of a destabilizing hydrogen bonding pair

The contribution of a specific hydrogen bond in apoflavodoxin to protein stability is investigated by combining theory, experiment and simulation. Although hydrogen bonds are major determinants of protein structure and function, their contribution to protein stability is still unclear and widely debated. The best method so far devised to estimate the contribution of side-chain interactions to protein stability is double mutant cycle analysis, but the interaction energies so derived are not identical to incremental binding energies (the energies quantifying net contributions of two interacting groups to protein stability). Here we introduce double-deletion analysis of 'isolated' residue pairs as a means to precisely quantify incremental binding. The method is exemplified by studying a surface-exposed hydrogen bond in a model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these residues by alanines slightly destabilizes the protein due to a decrease in hydrophobic surface burial. Subtraction of this effect, however, clearly indicates that the hydrogen-bonded groups in fact destabilize the native conformation. In addition, molecular dynamics simulations and classic double mutant cycle analysis explain quantitatively that, due to frustration, the hydrogen bond must form in the native structure because when the two groups get approximated upon folding their binding becomes favorable. We would like to remark that 1), this is the first time the contribution of a specific hydrogen bond to protein stability has been measured by experiment; and 2), more hydrogen bonds need to be analyzed to draw general conclusions on protein hydrogen bond energetics. To that end, the double-deletion method should be of help.

The Long and Short Flavodoxins

Flavodoxins are classified in two groups according to the presence or absence of a approximately 20-residue loop of unknown function. In the accompanying paper (36), we have shown that the differentiating loop from the long-chain Anabaena PCC 7119 flavodoxin is a peripheral structural element that can be removed without preventing the proper folding of the apoprotein. Here we investigate the role played by the loop in the stability and folding mechanism of flavodoxin by comparing the equilibrium and kinetic behavior of the full-length protein with that of loop-lacking, shortened variants. We show that, when the loop is removed, the three-state equilibrium thermal unfolding of apoflavodoxin becomes two-state. Thus, the loop is responsible for the complexity shown by long-chain apoflavodoxins toward thermal denaturation. As for the folding reaction, both shortened and wild type apoflavodoxins display three-state behavior but their folding mechanisms clearly differ. Whereas the full-length protein populates an essentially off-pathway transient intermediate, the additional state observed in the folding of the shortened variant analyzed seems to be simply an alternative native conformation. This finding suggests that the long loop may also be responsible for the accumulation of the kinetic intermediate observed in the full-length protein. Most revealing, however, is that the influence of the loop on the overall conformational stability of apoflavodoxin is quite low and the natively folded shortened variant Delta(120-139) is almost as stable as the wild type protein. The fact that the loop, which is not required for a proper folding of the polypeptide, does not even play a significant role in increasing the conformational stability of the protein supports our proposal (36) that the differentiating loop of long-chain flavodoxins may be related to a recognition function, rather than serving a structural purpose.

The Long and Short Flavodoxins

Flavodoxins are well known one-domain alpha/beta electron-transfer proteins that, according to the presence or absence of a approximately 20-residue loop splitting the fifth beta-strand of the central beta-sheet, have been classified in two groups: long and short-chain flavodoxins, respectively. Although the flavodoxins have been extensively used as models to study electron transfer, ligand binding, protein stability and folding issues, the role of the loop has not been investigated. We have constructed two shortened versions of the long-chain Anabaena flavodoxin in which the split beta-strand has been spliced to remove the original loop. The two variants have been carefully analyzed using various spectroscopic and hydrodynamic criteria, and one of them is clearly well folded, indicating that the long loop is a peripheral element of the structure of long flavodoxins. However, the removal of the loop (which is not in contact with the cofactor in the native structure) markedly decreases the affinity of the apoflavodoxin-FMN complex. This seems related to the fact that, in long flavodoxins, the adjacent tyrosine-bearing FMN binding loop (which is longer and thus more flexible than in short flavodoxins) is stabilized in its competent conformation by interactions with the excised loop. The modest role played by the long loop of long flavodoxins in the structure of these proteins (and in its conformational stability, see Lopez-Llano, J., Maldonado, S., Jain, S., Lostao, A., Godoy-Ruiz, R., Sanchez-Ruiz, Cortijo, M., Fernandez-Recio, J., and Sancho, J. (2004) J. Biol. Chem. 279, 47184-47191) opens the possibility that its conservation in so many species is related to a functional role yet to be discovered. In this respect, we discuss the possibility that the long loop is involved in the recognition of some flavodoxin partners. In addition, we report on a structural feature of flavodoxins that could indicate that the short flavodoxins derive from the long ones.

Inhibition of pea ferredoxin-NADP(H) reductase by Zn-ferrocyanide

Ferredoxin-NADP(H) reductases (FNRs) represent a prototype of enzymes involved in numerous metabolic pathways. We found that pea FNR ferricyanide diaphorase activity was inhibited by Zn2+ (Ki 1.57 microM). Dichlorophenolindophenol diaphorase activity was also inhibited by Zn2+ (Ki 1.80 microM), but the addition of ferrocyanide was required, indicating that the inhibitor is an arrangement of both ions. Escherichia coli FNR was also inhibited by Zn-ferrocyanide, suggesting that inhibition is a consequence of common structural features of these flavoenzymes. The inhibitor behaves in a noncompetitive manner for NADPH and for artificial electron acceptors. Analysis of the oxidation state of the flavin during catalysis in the presence of the inhibitor suggests that the electron-transfer process between NADPH and the flavin is not significantly altered, and that the transfer between the flavin and the second substrate is mainly affected. Zn-ferrocyanide interacts with the reductase, probably increasing the accessibility of the prosthetic group to the solvent. Ferredoxin reduction was also inhibited by Zn-ferrocyanide in a noncompetitive manner, but the observed Ki was about nine times higher than those for the diaphorase reactions. The electron transfer to Anabaena flavodoxin was not affected by Zn-ferrocyanide. Binding of the apoflavodoxin to the reductase was sufficient to overcome the inhibition by Zn-ferrocyanide, suggesting that the interaction of FNRs with their proteinaceous electron partners may induce a conformational change in the reductase that alters or completely prevents the inhibitory effect.

Native hydrogen bonds in a molten globule: the apoflavodoxin thermal intermediate

The structure and energetics of protein-folding intermediates are poorly understood. We have identified, in the thermal unfolding of the apoflavodoxin from Anabaena PCC 7119, an equilibrium intermediate with spectroscopic properties of a molten globule and substantial enthalpy and heat capacity of unfolding. The structure of the intermediate is probed by mutagenesis (and phi analysis) of polar residues involved in surface-exposed hydrogen bonds connecting secondary-structure elements in the native protein. All hydrogen bonds analysed are formed in the molten globule intermediate, either with native strength or debilitated. This suggests the overall intermediate's topology and surface tertiary interactions are close to native, and indicates that hydrogen bonding may contribute significantly to shape the conformation and energetics of folding intermediates.

Stabilization of apoflavodoxin by replacing hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or Gln

Knowledge of protein stability principles provides a means to increase protein stability in a rational way. Here we explore the feasibility of stabilizing proteins by replacing solvent-exposed hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or GLN: The rationale behind this is a previous observation that, in some cases, neutral hydrogen bonds may be more stable that charged ones. We identified, in the apoflavodoxin from Anabaena PCC 7119, three surface-exposed aspartate or glutamate residues involved in hydrogen bonding with a single partner and we mutated them to asparagine or glutamine, respectively. The effect of the mutations on apoflavodoxin stability was measured by both urea and temperature denaturation. We observed that the three mutant proteins are more stable than wild-type (on average 0.43 kcal/mol from urea denaturation and 2.8 degrees C from a two-state analysis of fluorescence thermal unfolding data). At high ionic strength, where potential electrostatic repulsions in the acidic apoflavodoxin should be masked, the three mutants are similarly more stable (on average 0.46 kcal/mol). To rule out further that the stabilization observed is due to removal of electrostatic repulsions in apoflavodoxin upon mutation, we analysed three control mutants and showed that, when the charged residue mutated to a neutral one is not hydrogen bonded, there is no general stabilizing effect. Replacing hydrogen-bonded charged Asp or Glu residues by Asn or Gln, respectively, could be a straightforward strategy to increase protein stability.