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

Effect of Dimer Structure and Inhomogeneous Broadening of Energy Levels on the Action of Flavomononucleotide in Rigid Polyvinyl Alcohol Films

The results of time-resolved fluorescence measurements of flavin mononucleotide (FMN) in rigid polyvinyl alcohol films (PVA) demonstrate that fluorescence intensity decays are strongly accelerated in the presence of fluorescent dimers and nonradiative energy transfer processes. The fluorescence decay originating both from H and J dimer states of FMN was experimentally observed for the first time. The mean fluorescence lifetimes for FMN dimers were obtained: τfl = 2.66 ns (at λ_exc = 445 nm) and τfl = 2.02 (at λ_exc = 487 nm) at λ_obs = 600 nm and T = 253 K from H and J state of dimers, respectively. We show that inhomogeneous orientational broadening of energy levels (IOBEL) affects the shape of the fluorescence decay and leads to the dependence of the average monomer fluorescence lifetime on excitation wavelength. IOBEL affected the nonradiative energy transfer and indicated that different flavin positioning in the protein pocket could (1) change the spectroscopic properties of flavins due to the existence of “blue” and “red” fluorescence centers, and (2) diminish the effectiveness of energy transfer between FMN molecules.

Folding of proteins with a flavodoxin-like architecture

The flavodoxin-like fold is a protein architecture that can be traced back to the universal ancestor of the three kingdoms of life. Many proteins share this α-β parallel topology and hence it is highly relevant to illuminate how they fold. Here, we review experiments and simulations concerning the folding of flavodoxins and CheY-like proteins, which share the flavodoxin-like fold. These polypeptides tend to temporarily misfold during unassisted folding to their functionally active forms. This susceptibility to frustration is caused by the more rapid formation of an α-helix compared to a β-sheet, particularly when a parallel β-sheet is involved. As a result, flavodoxin-like proteins form intermediates that are off-pathway to native protein and several of these species are molten globules (MGs). Experiments suggest that the off-pathway species are of helical nature and that flavodoxin-like proteins have a nonconserved transition state that determines the rate of productive folding. Folding of flavodoxin from Azotobacter vinelandii has been investigated extensively, enabling a schematic construction of its folding energy landscape. It is the only flavodoxin-like protein of which cotranslational folding has been probed. New insights that emphasize differences between in vivo and in vitro folding energy landscapes are emerging: the ribosome modulates MG formation in nascent apoflavodoxin and forces this polypeptide toward the native state.

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded proteins.

Comparison of chemical and thermal protein denaturation by combination of computational and experimental approaches. II

Chemical and thermal denaturation methods have been widely used to investigate folding processes of proteins in vitro. However, a molecular understanding of the relationship between these two perturbation methods is lacking. Here, we combined computational and experimental approaches to investigate denaturing effects on three structurally different proteins. We derived a linear relationship between thermal denaturation at temperature T(b) and chemical denaturation at another temperature T(u) using the stability change of a protein (ΔG). For this, we related the dependence of ΔG on temperature, in the Gibbs-Helmholtz equation, to that of ΔG on urea concentration in the linear extrapolation method, assuming that there is a temperature pair from the urea (T(u)) and the aqueous (T(b)) ensembles that produces the same protein structures. We tested this relationship on apoazurin, cytochrome c, and apoflavodoxin using coarse-grained molecular simulations. We found a linear correlation between the temperature for a particular structural ensemble in the absence of urea, T(b), and the temperature of the same structural ensemble at a specific urea concentration, T(u). The in silico results agreed with in vitro far-UV circular dichroism data on apoazurin and cytochrome c. We conclude that chemical and thermal unfolding processes correlate in terms of thermodynamics and structural ensembles at most conditions; however, deviations were found at high concentrations of denaturant.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.

Conformational selection or induced fit: a flux description of reaction mechanism

The mechanism of ligand binding coupled to conformational changes in macromolecules has recently attracted considerable interest. The 2 limiting cases are the "induced fit" mechanism (binding first) or "conformational selection" (conformational change first). Described here are the criteria by which the sequence of events can be determined quantitatively. The relative importance of the 2 pathways is determined not by comparing rate constants (a common misconception) but instead by comparing the flux through each pathway. The simple rules for calculating flux in multistep mechanisms are described and then applied to 2 examples from the literature, neither of which has previously been analyzed using the concept of flux. The first example is the mechanism of conformational change in the binding of NADPH to dihydrofolate reductase. The second example is the mechanism of flavodoxin folding coupled to binding of its cofactor, flavin mononucleotide. In both cases, the mechanism switches from being dominated by the conformational selection pathway at low ligand concentration to induced fit at high ligand concentration. Over a wide range of conditions, a significant fraction of the flux occurs through both pathways. Such a mixed mechanism likely will be discovered for many cases of coupled conformational change and ligand binding when kinetic data are analyzed by using a flux-based approach.

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.

Molecular crowding enhances native structure and stability of alpha/beta protein flavodoxin

To investigate the consequences of macromolecular crowding on the behavior of a globular protein, we performed a combined experimental and computational study on the 148-residue single-domain alpha/beta protein, Desulfovibrio desulfuricans apoflavodoxin. In vitro thermal unfolding experiments, as well as assessment of native and denatured structures, were probed by using far-UV CD in the presence of various amounts of Ficoll 70, an inert spherical crowding agent. Ficoll 70 has a concentration-dependent effect on the thermal stability of apoflavodoxin (DeltaT(m) of 20 degrees C at 400 mg/ml; pH 7). As judged by CD, addition of Ficoll 70 causes an increase in the amount of secondary structure in the native-state ensemble (pH 7, 20 degrees C) but only minor effects on the denatured state. Theoretical calculations, based on an off-lattice model and hard-sphere particles, are in good agreement with the in vitro data. The simulations demonstrate that, in the presence of 25% volume occupancy of spheres, native flavodoxin is thermally stabilized, and the free energy landscape shifts to favor more compact structures in both native and denatured states. The difference contact map reveals that the native-state compaction originates in stronger interactions between the helices and the central beta-sheet, as well as by less fraying in the terminal helices. This study demonstrates that macromolecular crowding has structural effects on the folded ensemble of polypeptides.

Survey of the year 2005: literature on applications of isothermal titration calorimetry

Isothermal titration calorimetry (ITC) can provide a full thermodynamic characterization of an interaction. Its usage does not suffer from constraints of molecular size, shape or chemical constitution. Neither is there any need for chemical modification or attachment to solid support. This ease of use has made it an invaluable instrumental resource and led to its appearance in many laboratories. Despite this, the value of the thermodynamic parameterization has, only quite recently, become widely appreciated. Although our understanding of the correlation between thermodynamic data and structural details continues to be somewhat naïve, a large number of publications have begun to improve the situation. In this overview of the literature for 2005, we have attempted to highlight works of interest and novelty. Furthermore, we draw attention to those works which we feel have provided a route to better analysis and increased our ability to understand the meaning of thermodynamic change on binding.

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.

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.

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.