Denaturant dependence of equilibrium unfolding intermediates and denatured state structure of horse ferricytochrome c

Nuclear magnetic resonance spectroscopy is employed to characterize unfolding intermediates and the denatured state of horse ferricytochrome c in guanidine hydrochloride. Unfolded and partially unfolded species with non-native heme ligation are detected by analysis of hyperfine-shifted (1)H resonances. Two equilibrium unfolding intermediates with His-Lys heme axial ligation are detected, as are two unfolded species with bis-His heme ligation. These results are contrasted with previous results on horse ferricytochrome c denaturation by urea, for which only one unfolding intermediate and one unfolded species were detected by NMR spectroscopy. Urea and guanidine hydrochloride are often used interchangeably in protein denaturation studies, but these results and those of others indicate that unfolded and intermediate states in these two denaturants may have substantially different properties. Implications of these results for folding studies and the biological function of mitochondrial cytochromes c are discussed.     

Characterization of equilibrium intermediates in denaturant-induced unfolding of ferrous and ferric cytochromes c using magnetic circular dichroism, circular dichroism, and optical absorption spectroscopies

Protein unfolding during guanidine HCl denaturant titration of the reduced and oxidized forms of cytochrome c is monitored with magnetic circular dichroism (MCD), natural CD, and absorption of the heme bands and far-UV CD of the amide bands. Direct MCD spectral evidence is presented for bis-histidinyl heme ligation in the unfolded states of both the reduced and oxidized protein. For both redox states, the unfolding midpoints measured with MCD, which is an indicator of tertiary structure, are significantly lower than those measured with far-UV CD, an indicator of secondary structure. The disparate titration curves are interpreted in terms of a compound mechanism for denaturant-induced folding and unfolding involving a molten globulelike intermediate state (MG) with near-native secondary structure and nonnative tertiary structure and heme ligation. A comparison of the dependence of the free energy of formation of the MG intermediate on the redox state with the known contributions from heme ligation and solvation suggests that the heme is significantly more accessible to solvent in the MG intermediate than it is in the native state.     

Folding of horse cytochrome c in the reduced state

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10 degrees C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded DeltaG(H2O)=18.8(+/-1.45) kcal mol(-1), and C(m)=5.1(+/-0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of -164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U<-->N) of 80-90 % of unfolded molecules with a time constant in the range 190-550 micros estimated by linear extrapolation and model calculations. The remaining 10-20 % of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I<-->U<-->N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10-20 % of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pK(a) value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is approximately 300(+/-125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I<-->U<-->N scheme have yielded the microscopic rate constants.     

Heme ligation and conformational plasticity in the isolated c domain of cytochrome cd1 nitrite reductase

The heme ligation in the isolated c domain of Paracoccus pantotrophus cytochrome cd(1) nitrite reductase has been characterized in both oxidation states in solution by NMR spectroscopy. In the reduced form, the heme ligands are His69-Met106, and the tertiary structure around the c heme is similar to that found in reduced crystals of intact cytochrome cd1 nitrite reductase. In the oxidized state, however, the structure of the isolated c domain is different from the structure seen in oxidized crystals of intact cytochrome cd1, where the c heme ligands are His69-His17. An equilibrium mixture of heme ligands is present in isolated oxidized c domain. Two-dimensional exchange NMR spectroscopy shows that the dominant species has His69-Met106 ligation, similar to reduced c domains. This form is in equilibrium with a high-spin form in which Met106 has left the heme iron. Melting studies show that the midpoint of unfolding of the isolated c domain is 320.9 +/- 1.2 K in the oxidized and 357.7 +/- 0.6 K in the reduced form. The thermally denatured forms are high-spin in both oxidation states. The results reveal how redox changes modulate conformational plasticity around the c heme and show the first key steps in the mechanism that lead to ligand switching in the holoenzyme. This process is not solely a function of the properties of the c domain. The role of the d1 heme in guiding His17 to the c heme in the oxidized holoenzyme is discussed.     

Optical spectroscopic differentiation of various equilibrium denatured states of horse cytochrome c

Thermal unfolding of cytochrome c (cyt c) from several states has been studied using equilibrium spectroscopic techniques. CD in the uv, vibrational circular dichroism, infrared, and uv-vis absorption spectra measured at various temperatures, pHs, salt concentrations, and GuHCl concentrations are used to show the conformational as well as heme structural differences between native and various denatured states. The difference in thermal denaturation behaviors of cyt c starting from acid denatured, molten globule (MG), and the A and native states are explored. Different final high temperature states were observed for cytochrome c unfolding from four different initial states (native, MG, A, and acid denatured state) by electronic CD, Fourier transform infrared (FTIR), and vibrational CD (VCD). Consistent with this, different thermal unfolding pathways for the MG and A states are suggested by the FTIR and VCD data for this process.     

Fast folding of cytochrome c.

Native iso-2 cytochrome c contains two residues (His 18, Met 80) coordinated to the covalently attached heme. On unfolding of iso-2, the His 18 ligand remains coordinated to the heme iron, whereas Met 80 is displaced by a non-native heme ligand, His 33 or His 39. To test whether non-native His-heme ligation slows folding, we have constructed a double mutant protein in which the non-native ligands are replaced by asparagine and lysine, respectively (H33N,H39K iso-2). The double mutant protein, which cannot form non-native histidine-heme coordinate bonds, folds significantly faster than normal iso-2 cytochrome c: gamma = 14-26 ms for H33N,H39K iso-2 versus gamma = 200-1,100 ms for iso-2. These results with iso-2 cytochrome c strongly support the hypothesis that non-native His-heme ligation results in a kinetic barrier to fast folding of cytochrome c. Assuming that the maximum rate of a conformational search is about 10(11) s-1, the results imply that the direct folding pathway of iso-2 involves passage through on the order of 10(9) or fewer partially folded conformers.     

A model for the misfolded bis-His intermediate of cytochrome c: the 1-56 N-fragment

We have characterized the ferric and ferrous forms of the heme-containing (1-56 residues) N-fragment of horse heart cytochrome c (cyt c) at different pH values and low ionic strength by UV-visible absorption and resonance Raman (RR) scattering. The results are compared with native cyt c in the same experimental conditions as this may provide a deeper insight into the cyt c unfolding-folding process. Folding of cyt c leads to a state having the heme iron coordinated to a histidine (His18) and a methionine (Met80) as axial ligands. At neutral pH the N-fragment (which lacks Met80) shows absorption and RR spectra that are consistent with the presence of a bis-His low spin heme, like several non-native forms of the parental protein. In particular, the optical spectra are identical to those of cyt c in the presence of a high concentration of denaturants; this renders the N-fragment a suitable model to study the heme pocket microenvironment of the misfolded (His-His) intermediate formed during folding of cyt c. Acid pH affects the ligation state in both cyt c and the N-fragment. Data obtained as a function of pH allow a correlation between the structural properties in the heme pocket of the N-fragment and those of non-native forms of cyt c. The results underline that the (57-104 residues) segment under native-like conditions imparts structural stability to the protein by impeding solvent access into the heme pocket.     

Y97V substitution in the horse cytochrome c causes accumulation of the equilibrium intermediate

Equilibrium unfolding experiments on several mutant forms of horse heart cytochrome c were performed. By means of absorbance spectroscopy, the accumulation of an equilibrium intermediate was revealed upon unfolding of Y97V mutant protein, and its structural properties were characterized. The data obtained allow one to conclude that the equilibrium intermediate corresponds to the earliest kinetic intermediate Ic in cytochrome c folding reaction. A comparative analysis of spectral properties of unfolded states of cytochrome c induced by urea or guanidine hydrochloride is presented.     

pH-induced conformation changes and equilibrium unfolding in yeast iso-2 cytochrome c

The relationship between pH-induced conformational changes in iso-2 cytochrome c from Saccharomyces cerevisiae and the guanidine hydrochloride induced unfolding transition has been investigated. Comparison of equilibrium unfolding transitions at acid, neutral, and alkaline pH shows that stability toward guanidine hydrochloride denaturation is decreased at low pH but increased at high pH. In the acid range the decrease in stability of the folded protein is correlated with changes in the visible spectrum, which indicate conversion to a high-spin heme state--probably involving the loss of heme ligands. The increase in stability at high pH is correlated with a pH-induced conformational change with an apparent pK near 8. As in the case of homologous cytochromes c, this transition involves the loss of the 695-nm absorbance band with only minor changes in other optical parameters. For the unfolded protein, optical spectroscopy and 1H NMR spectroscopy are consistent with a random coil unfolded state in which amino acid side chains serve as (low-spin) heme ligands at both neutral and alkaline pH. However, the paramagnetic region of the proton NMR spectrum of unfolded iso-2 cytochrome c indicates a change in the (low-spin) heme-ligand complex at high pH. Apparently, the folded and unfolded states of the (inactive) alkaline form differ from the corresponding states of the less stable native protein.     

Ligand exchange during unfolding of cytochrome c.

The productive folding pathway of cytochrome c passes through an obligatory HW intermediate in which the heme is coordinated by a solvent water molecule and a native ligand, His-18, prior to the formation of the folded HM state with both the native His-18 and Met-80 heme coordination. Two off pathway intermediates, a five-coordinated state (5C) and a bis-histidine state (HH), were also identified during the folding reaction. In the present work, the thermodynamics and the kinetics of the unfolding reaction of cytochrome c were investigated with resonance Raman scattering, tryptophan fluorescence spectroscopy, and circular dichroism. The objective of these experiments was to determine if the protein opens up and diverges into the differing heme ligation states through a many pathway mechanism or if it passes through intermediate states analogous to those observed during the folding reaction. Equilibrium unfolding results indicate that, in contrast to 5C, the stability of HH with respect to HW decreases as the concentration of GdnHCl increases. The difference in their response to the denaturant indicates that the polypeptide structure of 5C is relatively loose as compared with HH in which the polypeptide is misfolded. Time-resolved resonance Raman measurements show that strikingly similar ligand exchange reactions occur during unfolding as were observed during folding. Combined with fluorescence data, a kinetic model is proposed in which local structural rearrangements controlled by heme ligand exchange reactions appear prior to the global relaxation of the polypeptide chain.     

Surface-enhanced resonance Raman spectroscopy and spectroscopy study of redox-induced conformational equilibrium of cytochrome c adsorbed on DNA-modified metal electrode

The redox-induced conformational equilibrium of cytochrome c (cyt c) adsorbed on DNA-modified metal electrode and the interaction mechanism of DNA with cyt c have been studied by electrochemical, spectroscopic and spectroelectrochemical techniques. The results indicate that the external electric field induces potential-dependent coordination equilibrium of the adsorbed cyt c between its oxidized state (with native six-coordinate low-spin and non-native five-coordinate high-spin heme configuration) and its reduced state (with native six-coordinate low-spin heme configuration) on DNA-modified metal electrode. The strong interactions between DNA and cyt c induce the self-aggregation of cyt c adsorbed on DNA. The orientational distribution of cyt c adsorbed on DNA-modified metal electrode is potential-dependent, which results in the deviation from an ideal Nernstian behavior of the adsorbed cyt c at high electrode potentials. The electric-field-induced increase in the activation barrier of proton-transfer steps attributed to the rearrangement of the hydrogen bond network and the self-aggregation of cyt c upon adsorption on DNA-modified electrode strongly decrease the interfacial electron transfer rate. In addition, the strongly Coulombic interactions between DNA and cyt c only disturb the microenvironment of the heme, and do not affect the states of heme ligation and spin. The secondary structure of the adsorbed cyt c is retained, while the conformation of DNA is changed from the B form DNA to A form DNA.     

The conformational manifold of ferricytochrome c explored by visible and far-UV electronic circular dichroism spectroscopy

The oxidized state of cytochrome c is a subject of continuous interest, owing to the multitude of conformations which the protein can adopt in solution and on surfaces of artificial and cell membranes. The structural diversity corresponds to a variety of functions in electron transfer, peroxidase and apoptosis processes. In spite of numerous studies, a comprehensive analysis and comparison of native and non-native states of ferricytochrome c has thus far not been achieved. This results in part from the fact that the influence of solvent conditions (i.e., ionic strength, anion concentration, temperature dependence of pH values) on structure, function and equilibrium thermodynamics has not yet been thoroughly assessed. The current study is a first step in this direction, in that it provides the necessary experimental data to compare different non-native states adopted at high temperature and alkaline pH. To this end, we employed visible electronic circular dichroism (ECD) and absorption spectroscopy to probe structural changes of the heme environment in bovine and horse heart ferricytochrome c as a function of temperature between 278 and 363 K at different neutral and alkaline pH values. A careful selection of buffers enabled us to monitor the partial unfolding of the native state at room temperature while avoiding a change to an alkaline state at high temperatures. We found compelling evidence for the existence of a thermodynamic intermediate of the thermal unfolding/folding process, termed III h, which is structurally different from the alkaline states, IV 1 and IV 2, contrary to current belief. At neutral or slightly acidic pH, III h is populated in a temperature region between 320 and 345 K. The unfolded state of the protein becomes populated at higher temperatures. The ECD spectra of the B-bands of bovine and horse heart cytochrome c (pH 7.0) exhibit a pronounced couplet that is maintained below 343 K, before protein unfolding replaces it by a rather strong positive Cotton band. A preliminary vibronic analysis of the B-band profile reveals that the couplet reflects a B-band splitting of 350 cm (-1), which is mostly of electronic origin, due to the internal electric field in the heme cavity. Our results suggest that the conformational transition from the native state, III, into a thermally activated intermediate state, III h, does not substantially affect the internal electric field and causes only moderate rearrangements of the heme pocket, which involves changes, rather than a rupture, of the Fe (3+)-M80 linkage. In the unfolded state, as well as in the alkaline states IV and V, the band splitting is practically eliminated, but the positive Cotton effect observed for the B-band suggests that the proximal environment, encompassing H18 and the two cysteine residues 14 and 17, is most likely still intact and covalently bound to the heme chromophore. Both alkaline states IV and V were found to melt via intermediate states. Unfolded states probed at neutral and alkaline pH can be discriminated, owing to the different intensities of the Cotton bands of the respective B-band transitions. Differences between the ECD intensities of the B-bands of the different unfolded states and alkaline states most likely reflect different degrees of openness of the corresponding heme crevice.     

NMR study of the conformational transition of cytochrome c upon the displacement of Met80 by exogenous ligand: structural and magnetic characterization of azidoferricytochrome c

As the exogenous ligand-cytochrome c complexes were purported to represent models for the unfolding intermediate of cytochrome c, NMR spectroscopy has been utilized to study the azide adduct of horse heart cytochrome c. The structure of azidoferricytochrome c was modeled by restrained energy minimization using paramagenetic pseudocontact shifts as constraints. The bound azide moiety was found to be tilted approximately 15 degrees from the heme normal. The displacement of Met80 by the exogenous azide molecule causes large structural rearrangement in the distal cavity. Furthermore, the conformation transition associated with the swing out of the loop containing Met80 and the shift of the 50s-helix increases the solvent accessibility of the heme group. To elucidate the heme electronic structure of the complex, the paramagnetic 13C shifts were analyzed in terms of a model based on the pi molecular orbitals of the heme under perturbed D(4) symmetry. It turned out that the His-Fe bonding provides the protein constraint that orients the in-plane anisotropy in the complex. The electronic properties are in accordance with the calculated magnetic susceptibility anisotropy and the structural information.     

Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: a comparative study of the folding reactions of alpha-lactalbumin and lysozyme

The refolding kinetics of alpha-lactalbumin at different concentrations of guanidine hydrochloride have been investigated by means of kinetic circular dichroism and stopped-flow absorption measurements. The refolding reaction consists of at least two stages, the instantaneous accumulation of the transient intermediate that has peptide secondary structure and the subsequent slow process associated with formation of tertiary structure. The transient intermediate is compared with the well-characterized equilibrium intermediate observed during the denaturant-induced unfolding. Stabilities of the secondary structures against the denaturant, affinities for Ca2+, and tryptophan absorption properties of the transient and equilibrium intermediates were investigated. In all of these respects, the transient intermediate is identical with the equilibrium one, demonstrating the validity of the use of the equilibrium intermediate as a model of the folding intermediate. Essentially the same transient intermediate was also detected in the folding of lysozyme, the protein known to be homologous to alpha-lactalbumin but whose equilibrium unfolding is represented as a two-state reaction. The stability and cooperativity of the secondary structure of the intermediate of lysozyme are compared with those of alpha-lactalbumin. The results show that the protein folding occurring via the intermediate is not limited to the proteins that show equilibrium intermediates. Although the unfolding equilibria of most proteins are well approximated as a two-state reaction, the two-state hypothesis may not be applicable to the folding reaction under the native condition. Two models of protein folding, intermediate-controlled folding model and multiple-pathway folding model, which are different in view of the role of the intermediate in determining the pathway of folding, are also discussed.     

Many residues in cytochrome c populate alternative states under equilibrium conditions

A curved temperature dependence of an amide proton NMR chemical shift indicates that it explores discrete alternative conformations at least 1% of the time; that is, it accesses conformations that lie within 5 kcal/mol(-1) of the ground state. The simulations presented show how curvature varies with the nature of the alternative state, and are compared to experimental results. From studies in different denaturant concentrations, it is concluded that at least 25% of residues in reduced horse cytochrome c, covering most of the protein, with the exception of the center of the N- and C-terminal helices, visit alternative states under equilibrium conditions. The conformational ensemble of the protein therefore has high structural entropy. The density of alternative states is particularly high near the heme ligand Met80, which is of interest because both redox change and the first identified stage in unfolding are associated with change in Met80 ligation. By combining theoretical and experimental approaches, it is concluded that the alternative states each comprise approximately five residues, have in general less structure than the native state, and are accessed independently. They are therefore locally unfolded structures. The locations of the alternative states match what is known of the global unfolding pathway of cytochrome c, suggesting that they may determine the pathway.     

Characterization of Hydrogenobacter thermophilus cytochromes c 552 expressed in the cytoplasm and periplasm of Escherichia coli

Hydrogenobacter thermophilus cytochrome c(552) ( Ht cyt c(552)) is a small monoheme protein in the cytochrome c(551) family. Ht cyt c(552) is unique because it is hypothesized to undergo spontaneous cytoplasmic maturation (covalent heme attachment) when expressed in Escherichia coli. This is in contrast to the usual maturation route for bacterial cytochromes c that occurs in the cellular periplasm, where maturation factors direct heme attachment. Here, the expression of Ht cyts c(552) in the periplasm as well as the cytoplasm of E. coli is reported. The products are characterized by absorption, circular dichroism, and NMR spectroscopy as well as mass spectrometry, proteolysis, and denaturation studies. The periplasmic product's properties are found to be indistinguishable from those reported for protein isolated from Ht cells, while the major cytoplasmic product exhibits structural anomalies in the region of the N-terminal helix. These anomalies are shown to result from the retention of the N-terminal methionine in the cytoplasmic product, and not from heme attachment errors. The (1)H NMR chemical shifts of the heme methyls of the oxidized ( S=1/2) expression products display a unique pattern not previously reported for a cytochrome c with histidine-methionine axial ligation, although they are consistent with native-like heme ligation. These results support the hypothesis that proper heme attachment can occur spontaneously in the E. coli cytoplasm for Ht cyt c(552).     

The unfolding of the cytochromes c in methanol and acid.

The cytochromes c are a family of hemoproteins that share a number of structural features: a thioether linkage between the protein and the heme, histidine and methionine as the fifth and sixth iron ligands, and a tertiary structure known as the "cytochrome fold." These proteins follow a common mechanism of equilibrium unfolding in methanol and acid, differing only in their reactivity to the denaturing conditions. The reduced cytochromes c exhibit an increased conformational stability which is consistent with the presence of a strengthened iron-methionine linkage in the reduced state.     

Trichloroacetic acid and trifluoroacetic acid-induced unfolding of cytochrome c: stabilization of a native-like folded intermediate(1)

A systematic investigation of trichloroacetic acid (TCA) and trifluoroacetic acid (TFA)-induced equilibrium unfolding of native horse cytochrome c has been carried out using a combination of optical spectroscopy and electrospray ionization mass spectroscopy (ESI MS). In the presence of an increasing concentration of TCA the native cytochrome c does not undergo significant unfolding but stabilization of a partially folded intermediate is observed. This TCA-induced partially folding intermediate of cytochrome c had an enhanced secondary structure and slightly disrupted tertiary structure compared to native protein and undergoes extensive unfolding in the presence of TFA. However, in the presence of an increasing concentration of TFA, cytochrome c was found to undergo extensive unfolding characterized by a significant breakdown of the secondary and tertiary structure of protein. The TFA-unfolded cytochrome c was found to undergo folding in the presence of TCA and low guanidine hydrochloride (GdmCl) resulting in the stabilization of the partially folded intermediate. The effectiveness of TCA as compared to TFA in the stabilization of intermediates was further supported by the observation that low concentrations of TCA were found to induce refolding of HCl-denatured cytochrome c whereas, under similar concentrations of acid, no significant effect on the unfolded structure of protein was observed in the presence of TFA. ESI MS studies indicated that the trichloroacetate anion has a greater affinity for cytochrome c compared to trifluoroacetate anion, which might be the reason for the stabilization of the native-like folded intermediate during TCA-induced denaturation of cytochrome c as compared to extensive unfolding observed in the presence of TFA.     

Last in, first out: the role of cofactor binding in flavodoxin folding

Although many proteins require the binding of a ligand to be functional, the role of ligand binding during folding is scarcely investigated. Here, we have reported the influence of the flavin mononucleotide (FMN) cofactor on the global stability and folding kinetics of Azotobacter vinelandii holoflavodoxin. Earlier studies have revealed that A. vinelandii apoflavodoxin kinetically folds according to the four-state mechanism: I(1) <=> unfolded apoflavodoxin <=> I(2) <=> native apoflavodoxin. I(1)an off-pathway molten globule-like is intermediate that populates during denaturant-induced equilibrium unfolding; I(2) is a high energy on-pathway folding intermediate that never populates to a significant extent. Here, we have presented extensive denaturant-induced equilibrium unfolding data of holoflavodoxin, holoflavodoxin with excess FMN, and apoflavodoxin as well as kinetic folding and unfolding data of holoflavodoxin. All folding data are excellently described by a five-state mechanism: I(1) + FMN <=> unfolded apoflavodoxin + FMN <=> I(2) + FMN <=> native apoflavodoxin + FMN<=> holoflavodoxin. The last step in flavodoxin folding is thus the binding of FMN to native apoflavodoxin. I(1),I(2), and unfolded apoflavodoxin do not interact to a significantextent with FMN. The autonomous formation of native apoflavodoxin is essential during holoflavodoxin folding. Excess FMN does not accelerate holoflavodoxin folding, and FMN does not act as a nucleation site for folding. The stability of holoflavodoxin is so high that even under strongly denaturing conditions FMN needs to be released first before global unfolding of the protein can occur.