Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Do proteins always benefit from a stability increase? Relevant and residual stabilisation in a three-state protein by charge optimisation

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.     

Equilibrium phi-analysis of a molten globule: the 1-149 apoflavodoxin fragment

The apoflavodoxin fragment comprising residues 1-149 that can be obtained by chemical cleavage of the C-terminal alpha-helix of the full-length protein is known to populate a molten globule conformation that displays a cooperative behaviour and experiences two-state urea and thermal denaturation. Here, we have used a recombinant form of this fragment to investigate molten globule energetics and to derive structural information by equilibrium Phi-analysis. We have characterized 15 mutant fragments designed to probe the persistence of native interactions in the molten globule and compared their conformational stability to that of the equivalent full-length apoflavodoxin mutants. According to our data, most of the mutations analysed modify the stability of the molten globule fragment following the trend observed when the same mutations are implemented in the full-length protein. However, the changes in stability observed in the molten globule are much smaller and the Phi-values calculated are (with a single exception) below 0.4. This is consistent with an overall and significant debilitation of the native structure. Nevertheless, the fact that the molten globule fragment can be stabilised using as a guide the native structure of the full-length protein (by increasing helix propensity, optimising charge interactions and filling small cavities) suggests that the overall structure of the molten globule is still quite close to native, in spite of the lowered stability observed.     

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.     

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.     

Salt-induced stabilization of apoflavodoxin at neutral pH is mediated through cation-specific effects

Electrostatic contributions to the conformational stability of apoflavodoxin were studied by measurement of the proton and salt-linked stability of this highly acidic protein with urea and temperature denaturation. Structure-based calculations of electrostatic Gibbs free energy were performed in parallel over a range of pH values and salt concentrations with an empirical continuum method. The stability of apoflavodoxin was higher near the isoelectric point (pH 4) than at neutral pH. This behavior was captured quantitatively by the structure-based calculations. In addition, the calculations showed that increasing salt concentration in the range of 0 to 500 mM stabilized the protein, which was confirmed experimentally. The effects of salts on stability were strongly dependent on cationic species: K(+), Na(+), Ca(2+), and Mg(2+) exerted similar effects, much different from the effect measured in the presence of the bulky choline cation. Thus cations bind weakly to the negatively charged surface of apoflavodoxin. The similar magnitude of the effects exerted by different cations indicates that their hydration shells are not disrupted significantly by interactions with the protein. Site-directed mutagenesis of selected residues and the analysis of truncation variants indicate that cation binding is not site-specific and that the cation-binding regions are located in the central region of the protein sequence. Three-state analysis of the thermal denaturation indicates that the equilibrium intermediate populated during thermal unfolding is competent to bind cations. The unusual increase in the stability of apoflavodoxin at neutral pH affected by salts is likely to be a common property among highly acidic proteins.     

The equilibrium unfolding of Azotobacter vinelandii apoflavodoxin II occurs via a relatively stable folding intermediate.

A flavodoxin from Azotobacter vinelandii is chosen as a model system to study the folding of alpha/beta doubly wound proteins. The guanidinium hydrochloride induced unfolding of apoflavodoxin is demonstrated to be reversible. Apoflavodoxin thus can fold in the absence of the FMN cofactor. The unfolding curves obtained for wild-type, C69A and C69S apoflavodoxin as monitored by circular dichroism and fluorescence spectroscopy do not coincide. Apoflavodoxin unfolding occurs therefore not via a simple two-state mechanism. The experimental data can be described by a three-state mechanism of apoflavodoxin equilibrium unfolding in which a relatively stable intermediate is involved. The intermediate species lacks the characteristic tertiary structure of native apoflavodoxin as deduced from fluorescence spectroscopy, but has significant secondary structure as inferred from circular dichroism spectroscopy. Both spectroscopic techniques show that thermally-induced unfolding of apoflavodoxin also proceeds through formation of a similar molten globule-like species. Thermal unfolding of apoflavodoxin is accompanied by anomalous circular dichroism characteristics: the negative ellipticity at 222 nM increases in the transition zone of unfolding. This effect is most likely attributable to changes in tertiary interactions of aromatic side chains upon protein unfolding. From the presented results and hydrogen/deuterium exchange data, a model for the equilibrium unfolding of apoflavodoxin is presented.     

Anabaena apoflavodoxin hydrogen exchange: on the stable exchange core of the alpha/beta(21345) flavodoxin-like family

An important issue in modern protein biophysics is whether structurally homologous proteins share common stability and/or folding features. Flavodoxin is an archetypal alpha/beta protein organized in three layers: a central beta-sheet (strand order 21345) flanked by helices 1 and 5 on one side and helices 2, 3, and 4 on the opposite side. The backbone internal dynamics of the apoflavodoxin from Anabaena is analyzed here by the hydrogen exchange method. The hydrogen exchange rates indicate that 46 amide protons, distributed throughout the structure of apoflavodoxin, exchange relatively slowly at pH 7.0 (k(ex) < 10(-1) min(-1)). According to their distribution in the structure, protein stability is highest on the beta-sheet, helix 4, and on the layer formed by helices 1 and 5. The exchange kinetics of Anabaena apoflavodoxin was compared with those of the apoflavodoxin from Azotobacter, with which it shares a 48% sequence identity, and with Che Y and cutinase, two other alpha/beta (21345) proteins with no significant sequence homology with flavodoxins. Both similarities and differences are observed in the cores of these proteins. It is of interest that a cluster of a few structurally equivalent residues in the central beta-strands and in helix 5 is common to the cores.     

A general approach for detecting folding intermediates from steady-state and time-resolved fluorescence of single-tryptophan-containing proteins

During denaturant-induced equilibrium (un)folding of wild-type apoflavodoxin from Azotobacter vinelandii, a molten globule-like folding intermediate is formed. This wild-type protein contains three tryptophans. In this study, we use a general approach to analyze time-resolved fluorescence and steady-state fluorescence data that are obtained upon denaturant-induced unfolding of a single-tryptophan-containing variant of apoflavodoxin [i.e., W74/F128/F167 (WFF) apoflavodoxin]. The experimental data are assembled in matrices, and subsequent singular-value decomposition of these matrices (i.e., based on either steady-state or time-resolved fluorescence data) shows the presence of three significant, and independent, components. Consequently, to further analyze the denaturation trajectories, we use a three-state protein folding model in which a folding intermediate and native and unfolded protein molecules take part. Using a global analysis procedure, we determine the relative concentrations of the species involved and show that the stability of WFF apoflavodoxin against global unfolding is ～4.1 kcal/mol. Analysis of time-resolved anisotropy data of WFF apoflavodoxin unfolding reveals the remarkable observation that W74 is equally well fixed within both the native protein and the molten globule-like folding intermediate. Slight differences between the direct environments of W74 in the folding intermediate and native protein cause different rotameric populations of the indole in both folding species as fluorescence lifetime analysis reveals. Importantly, thermodynamic analyses of the spectral denaturation trajectories of the double-tryptophan-containing protein variants WWF apoflavodoxin and WFW apoflavodoxin show that these variants are significantly more stable (5.9 kcal/mol and 6.8 kcal/mol, respectively) than WFF apoflavodoxin (4.1 kcal/mol) Hence, tryptophan residues contribute considerably to the 10.5 kcal/mol thermodynamic stability of native wild-type apoflavodoxin.     

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.     

Effects of macromolecular crowding agents on protein folding in vitro and in silico

Proteins fold and function inside cells which are environments very different from that of dilute buffer solutions most often used in traditional experiments. The crowded milieu results in excluded-volume effects, increased bulk viscosity and amplified chances for inter-molecular interactions. These environmental factors have not been accounted for in most mechanistic studies of protein folding executed during the last decades. The question thus arises as to how these effects—present when polypeptides normally fold in vivo—modulate protein biophysics. To address excluded volume effects, we use synthetic macromolecular crowding agents, which take up significant volume but do not interact with proteins, in combination with strategically selected proteins and a range of equilibrium and time-resolved biophysical (spectroscopic and computational) methods. In this review, we describe key observations on macromolecular crowding effects on protein stability, folding and structure drawn from combined in vitro and in silico studies. As expected based on Minton’s early predictions, many proteins (apoflavodoxin, VlsE, cytochrome c, and S16) became more thermodynamically stable (magnitude depends inversely on protein stability in buffer) and, unexpectedly, for apoflavodoxin and VlsE, the folded states changed both secondary structure content and, for VlsE, overall shape in the presence of macromolecular crowding. For apoflavodoxin and cytochrome c, which have complex kinetic folding mechanisms, excluded volume effects made the folding energy landscapes smoother (i.e., less misfolding and/or kinetic heterogeneity) than in buffer.     

Design and Structure of an Equilibrium Protein Folding Intermediate: A Hint into Dynamical Regions of Proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.     

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.     

Sequence context modulates the stability of a GxxxG-mediated transmembrane helix-helix dimer

To quantify the relationship between sequence and transmembrane dimer stability, a systematic mutagenesis and thermodynamic study of the protein-protein interaction residues in the glycophorin A transmembrane helix-helix dimer was carried out. The results demonstrate that the glycophorin A transmembrane sequence dimerizes when its GxxxG motif is abolished by mutation to large aliphatic residues, suggesting that the sequence encodes an intrinsic propensity to self-associate independent of a GxxxG motif. In the presence of an intact GxxxG motif, the glycophorin A dimer stability can be modulated over a span of -0.5 kcal mol(-1) to +3.2 kcal mol(-1) by mutating the surrounding sequence context. Thus, these flanking residues play an active role in determining the transmembrane dimer stability. To assess the structural consequences of the thermodynamic effects of mutations, molecular models of mutant transmembrane domains were constructed, and a structure-based parameterization of the free energy change due to mutation was carried out. The changes in association free energy for glycophorin A mutants can be explained primarily by changes in packing interactions at the protein-protein interface. The energy cost of removing favorable van der Waals interactions was found to be 0.039 kcal mol(-1) per A2 of favorable occluded surface area. The value corresponds well with estimates for mutations in bacteriorhodopsin as well as for those mutations in the interiors of soluble proteins that create packing defects.     

Sequence-structure signals of 3D domain swapping in proteins

Three-dimensional domain swapping occurs when two or more identical proteins exchange identical parts of their structure to generate an oligomeric unit. It affects proteins with diverse sequences and structures, and is expected to play important roles in evolution, functional regulation and even conformational diseases. Here, we search for traces of domain swapping in the protein sequence, by means of algorithms that predict the structure and stability of proteins using database-derived potentials. Regions whose sequences are not optimal with regard to the stability of the native structure, or showing marked intrinsic preferences for non-native conformations in absence of tertiary interactions are detected in most domain-swapping proteins. These regions are often located in areas crucial in the swapping process and are likely to influence it on a kinetic or thermodynamic level. In addition, cation-pi interactions are frequently observed to zip up the edges of the interface between intertwined chains or to involve hinge loop residues, thereby modulating stability. We end by proposing a set of mutations altering the swapping propensities, whose experimental characterization would contribute to refine our in silico derived hypotheses.     

Formation of on- and off-pathway intermediates in the folding kinetics of Azotobacter vinelandii apoflavodoxin

The folding kinetics of the 179-residue Azotobacter vinelandii apoflavodoxin, which has an alpha-beta parallel topology, have been followed by stopped-flow experiments monitored by fluorescence intensity and anisotropy. Single-jump and interrupted refolding experiments show that the refolding kinetics involve four processes yielding native molecules. Interrupted unfolding experiments show that the two slowest folding processes are due to Xaa-Pro peptide bond isomerization in unfolded apoflavodoxin. The denaturant dependence of the folding kinetics is complex. Under strongly unfolding conditions (>2.5 M GuHCl), single exponential kinetics are observed. The slope of the chevron plot changes between 3 and 5 M denaturant, and no additional unfolding process is observed. This reveals the presence of two consecutive transition states on a linear pathway that surround a high-energy on-pathway intermediate. Under refolding conditions, two processes are observed for the folding of apoflavodoxin molecules with native Xaa-Pro peptide bond conformations, which implies the population of an intermediate. The slowest of these two processes becomes faster with increasing denaturant concentration, meaning that an unfolding step is rate-limiting for folding of the majority of apoflavodoxin molecules. It is shown that the intermediate that populates during refolding is off-pathway. The experimental data obtained on apoflavodoxin folding are consistent with the linear folding mechanism I(off) <==> U <==> I(on) <== > N, the off-pathway intermediate being the molten globule one that also populates during equilibrium denaturation of apoflavodoxin. The presence of such on-pathway and off-pathway intermediates in the folding kinetics of alpha-beta parallel proteins is apparently governed by protein topology.     

Importance of mutant position in Ramachandran plot for predicting protein stability of surface mutations

Understanding the mechanisms by which mutations affect protein stability is one of the most important problems in molecular biology. In this work, we analyzed the relationship between changes in protein stability caused by surface mutations and changes in 49 physicochemical, energetic, and conformational properties of amino acid residues. We found that the hydration entropy was the major contributor to the stability of surface mutations in helical segments; other properties responsible for size and volume of molecule also correlated significantly with stability. Classification of coil mutations based on their locations in the (phi-psi) map improved the correlation significantly, demonstrating the existence of a relationship between stability and strain energy, which indicates that the role of strain energy is very important for the stability of surface mutations. We observed that the inclusion of sequence and structural information raised the correlation, indicating the influence of surrounding residues on the stability of surface mutations. Further, we examined the previously reported "inverse relationship" between stability and hydrophobicity, and observed that the inverse hydrophobic effect was generally applicable only to coil mutations. The present study leads to a simple method for predicting protein stability changes caused by amino acid substitutions, which will be useful for protein engineering in designing novel proteins with increased stability and altered function.     

Consequences of single-site mutations in the intestinal fatty acid binding protein

The intestinal fatty acid binding protein (IFABP) is a small (15 kDa) protein consisting mostly of 10 antiparallel beta-strands (A-J) and a small helical region that serves as a portal for the ligand. Two beta-sheet structures (strands A-E and F-J) surround a cavity into which the ligand binds. In this work, we investigated how changes in the side chains of specific residues are propagated through the structure. To determine what these changes were and how they relate to changes in stability, (15)N chemical shift perturbations were measured and compared to those of the wild-type protein. Seven mutations, five of which change either valine or leucine to glycine, have been examined. All these mutants were less stable than wild-type IFABP, suggesting some structural changes. For five of the mutants, the data suggest that destabilization of a small region of the protein propagates throughout the structure, resulting in an overall decrease in stability. In two (Leu38Gly and Leu89Gly), the loss of cooperativity in the equilibrium denaturation curves suggests that the destabilization of one region may not be transmitted to other regions in a cooperative manner. It is shown that the effect of mutating hydrophobic residues is much greater than that observed upon mutation of a solvent-exposed polar residue.     

A comparison of the urea-induced unfolding of apoflavodoxin and flavodoxin from Desulfovibrio vulgaris.

The kinetics and thermodynamics of the urea-induced unfolding of flavodoxin and apoflavodoxin from Desulfovibrio vulgaris were investigated by measuring changes in flavin and protein fluorescence. The reaction of urea with flavodoxin is up to 5000 times slower than the reaction with the apoprotein (0.67 s(-1) in 3 m urea in 25 mm sodium phosphate at 25 degrees C), and it results in the dissociation of FMN. The rate of unfolding of apoflavodoxin depends on the urea concentration, while the reaction with the holoprotein is independent of urea. The rates decrease in high salt with the greater effect occurring with apoprotein. The fluorescence changes fit two-state models for unfolding, but they do not exclude the possibility of intermediates. Calculation suggests that 21% and 30% of the amino-acid side chains become exposed to solvent during unfolding of flavodoxin and apoflavodoxin, respectively. The equilibrium unfolding curves move to greater concentrations of urea with increase of ionic strength. This effect is larger with phosphate than with chloride, and with apoflavodoxin than with flavodoxin. In low salt the conformational stability of the holoprotein is greater than that of apoflavodoxin, but in high salt the relative stabilities are reversed. It is calculated that two ions are released during unfolding of the apoprotein. It is concluded that the urea-dependent unfolding of flavodoxin from D. vulgaris occurs because apoprotein in equilibrium with FMN and holoprotein unfolds and shifts the equilibrium so that flavodoxin dissociates. Small changes in flavin fluorescence occur at low concentrations of urea and these may reflect binding of urea to the holoprotein.