Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

Secondary structures and structural fluctuation in a dimeric protein, Streptomyces subtilisin inhibitor.

Based on the nuclear magnetic resonance assignments of a dimeric protein, Streptomyces subtilisin inhibitor (SSI), microscopic details of secondary structures in solution have been elucidated. The chemical shift index of C(alpha) signals, together with information on the hydrogen exchange rates of the backbone amide protons, were used to identify secondary structures. The locations of these secondary structures were found to be different in some critical points from those determined earlier by X-ray crystallography of the crystal. Notably, the beta3 strand is completely missing and the alpha2 helix is extended toward the C-terminus. Furthermore, hydrogen exchange experiments of individual peptide NH protons under strongly folding conditions revealed mechanisms of global and local structural fluctuation within the dimeric structure. It has been suggested that the global fluctuation of the monomeric unit occurs without affecting the accompanying monomer, in contrast to the equilibrium thermal unfolding, which is cooperative. Higher protection against hydrogen exchange for residues in part of the beta4 strand implies that this region might serve as a folding core.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Protein folding kinetics by combined use of rapid mixing techniques and NMR observation of individual amide protons

A method to be used for experimental studies of protein folding introduced by Schmid and Baldwin (J. Mol. Biol. 135: 199-215, 1979), which is based on the competition between amide hydrogen exchange and protein refolding, was extended by using rapid mixing techniques and 1H NMR to provide site-resolved kinetic information on the early phases of protein structure acquisition. In this method, a protonated solution of the unfolded protein is rapidly mixed with a deuterated buffer solution at conditions assuring protein refolding in the mixture. This simultaneously initiates the exchange of unprotected amide protons with solvent deuterium and the refolding of protein segments which can protect amide groups from further exchange. After variable reaction times the amide proton exchange is quenched while folding to the native form continues to completion. By using 1H NMR, the extent of exchange at individual amide sites is then measured in the refolded protein. Competition experiments at variable reaction times or variable pH indicate the time at which each amide group is protected in the refolding process. This technique was applied to the basic pancreatic trypsin inhibitor, for which sequence-specific assignments of the amide proton NMR lines had previously been obtained. For eight individual amide protons located in the beta-sheet and the C-terminal alpha-helix of this protein, apparent refolding rates in the range from 15 s-1 to 60 s-1 were observed. These rates are on the time scale of the fast folding phase observed with optical probes.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

Sequence-specific 1H and 15N resonance assignments and secondary structure of GDP-bound human c-Ha-Ras protein in solution

Summary All the backbone ¹H and ¹⁵N magnetic resonances (except for Pro residues) of the GDP-bound form of a truncated human c-Ha-ras proto-oncogene product (171 amino acid residues, the Ras protein) were assigned by ¹⁵N-edited two-dimensional NMR experiments on selectively ¹⁵N-labeled Ras proteins in combination with three-dimensional NMR experiments on the uniformly ¹⁵N-labeled protein. The sequence-specific assignments were made on the basis of the nuclear Overhauser effect (NOE) connectivities of amide protons with preceding amide and/or Cprotons. In addition to sequential NOEs, vicinal spin coupling constants for amide protons and C protons and deuterium exchange rates of amide protons were used to characterize the secondary structure of the GDP-bound Ras protein; six strands and five helices were identified and the topology of these elements was determined. The secondary structure of the Ras protein in solution was mainly consistent with that in crystal as determined by X-ray analyses. The deuterium exchange rates of amide protons were examined to elucidate the dynamic properties of the secondary structure elements of the Ras protein in solution. In solution, the -sheet structure in the Ras protein is rigid, while the second helix (A66-R73) is much more flexible, and the first and fifth helices (S17-124 and V152-L171) are more rigid than other helices. Secondary structure elements at or near the ends of the effector-region loop were found to be much more flexible in solution than in the crystalline state.     

Multi-state unfolding of the alpha subunit of tryptophan synthase, a TIM barrel protein: insights into the secondary structure of the stable equilibrium intermediates by hydrogen exchange mass spectrometry

The urea-induced unfolding of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, an eight-stranded (beta/alpha)(8) TIM barrel protein, has been shown to involve two stable equilibrium intermediates, I1 and I2, well populated at approximately 3 M and 5 M urea, respectively. The characterization of the I1 intermediate by circular dichroism (CD) spectroscopy has shown that I1 retains a significant fraction of the native ellipticity; the far-UV CD signal for the I2 species closely resembles that of the fully unfolded form. To obtain detailed insight into the disruption of secondary structure in the urea-induced unfolding process, a hydrogen exchange-mass spectrometry study was performed on alphaTS. The full-length protein was destabilized in increasing concentration of urea, the amide hydrogen atoms were pulse-labeled with deuterium, the labeled samples were quenched in acid and the products were analyzed by electrospray ionization mass spectrometry. Consistent with the CD results, the I1 intermediate protects up to approximately 129 amide hydrogen atoms against exchange while the I2 intermediate offers no protection. Electrospray ionization mass spectrometry analysis of the peptic fragments derived from alphaTS labeled at 3 M urea indicates that most of the region between residues 12-130, which constitutes the first four beta strands and three alpha helices, (beta/alpha)(1-3)beta(4), is structured. The (beta/alpha)(1-3)beta(4) module appears to represent the minimum sub-core of stability of the I1 intermediate. A 4+2+2 folding model is proposed as a likely alternative to the earlier 6+2 folding mechanism for alphaTS.     

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.     

Hydrogen exchange of the tetramerization domain of the human tumour suppressor p53 probed by denaturants and temperature.

We have analysed the hydrogen/deuterium exchange of the tetramerization domain of human tumour suppressor p53 under mild chemical denaturation conditions, and at different temperatures. Exchange behaviour has been measured for 16 amide protons in the chemical-denaturation studies and for seven protons in the temperature-denaturation studies. The exchange rates are in the range observed for other proteins with similar elements of secondary structure. The slowest-exchange core includes contributions from residues in the alpha helix and the beta sheet. However, only some of the slowest-exchanging protons correspond to residues involved in native interactions in the transient intermediate detected during the folding of this domain. The guanidinium-chloride denaturation curves of all residues seem to merge together, although they are well below the main isotherm of global unfolding. Thus, there is no evidence for several subglobal unfolding units. The activation parameters obtained from the temperature-denaturation experiments are similar to those obtained for monomeric proteins, and well below the global unfolding enthalpy obtained by circular dichroism measurements. Thus, the exchange studies at different denaturant concentrations and temperatures indicate that no particular folding intermediate is populated under those conditions.     

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.     

Amide proton hydrogen exchange rates for sperm whale myoglobin obtained from 15N-1H NMR spectra

The hydrogen exchange behavior of exchangeable protons in proteins can provide important information for understanding the principles of protein structure and function. The positions and exchange rates of the slowly-exchanging amide protons in sperm whale myoglobin have been mapped using 15N-1H NMR spectroscopy. The slowest-exchanging amide protons are those that are hydrogen bonded in the longest helices, including members of the B, E, and H helices. Significant protection factors were observed also in the A, C, and G helices, and for a few residues in the D and F helices. Knowledge of the identity of slowly-exchanging amide protons forms the basis for the extensive quench-flow kinetic folding experiments that have been performed for myoglobin, and gives insights into the tertiary interactions and dynamics in the protein.     

Structural and kinetic characterization of early folding events in beta-lactoglobulin

We have defined the structural and dynamic properties of an early folding intermediate of beta-lactoglobulin known to contain non-native alpha-helical structure. The folding of beta-lactoglobulin was monitored over the 100 micros--10 s time range using ultrarapid mixing techniques in conjunction with fluorescence detection and hydrogen exchange labeling probed by heteronuclear NMR. An initial increase in Trp fluorescence with a time constant of 140 micros is attributed to formation of a partially helical compact state. Within 2 ms of refolding, well protected amide protons indicative of stable hydrogen bonded structure were found only in a domain comprising beta-strands F, G and H, and the main alpha-helix, which was thus identified as the folding core of beta-lactoglobulin. At the same time, weak protection (up to approximately 10-fold) of amide protons in a segment spanning residues 12--21 is consistent with formation of marginally stable non-native alpha-helices near the N-terminus. Our results indicate that efficient folding, despite some local non-native structural preferences, is insured by the rapid formation of a native-like alpha/beta core domain.     

Sequence-specific 1H-NMR assignments and folding topology of human CD59.

CD59 is a recently discovered cell-surface glycoprotein that restricts lysis by homologous complement and has limited sequence similarity to snake venom neurotoxins. This paper describes the first results of a two-dimensional NMR study of CD59 prepared from human urine. Nearly complete 1H-NMR assignments were obtained for the 77 amino acid residues and partial assignments for the N-glycan and the glycosylphosphatidylinositol (GPI) anchor. These results together confirm that the C-terminal residue of the mature protein is Asn 77 and that the urine-derived form retains the nonlipid part of the GPI anchor. The data further indicate that the GPI anchor and possibly the N-glycan are structurally inhomogeneous and suggest that the phospholipid present in the intact GPI anchor was removed by phosphatidylinositol-specific phospholipase-D. The folding topology of the protein was determined from NOE enhancements and slowly exchanging backbone amide protons and consists primarily of five extended strands (denoted beta 1-beta 5 in sequence order), arranged into separate two-stranded (beta 1 and beta 2) and three-stranded (beta 3-beta 5) antiparallel beta-sheets. The same folding topology is found in all of the snake venom neurotoxins whose structures have been determined. The region between the beta 4 and beta 5 strands has helical character, a feature that is not present in the neurotoxins but that is seen in the topologically similar wheat germ agglutinin.     

The unfolding of oxidized c-type cytochromes: the instructive case of Bacillus pasteurii

The reversible unfolding of oxidized Bacillus pasteurii cytochrome c(553) by guanidinium chloride under equilibrium conditions has been monitored by NMR and optical spectroscopy. The results obtained indicate that unfolding takes place through a mechanism involving the detachment from heme iron coordination of the sulfur of the Met71 axial ligand and yielding either a high spin (HS) or a low spin (LS(1)) species, depending on the pH value. In the LS(1) form the Met71 is replaced by another protein ligand, possibly Lys. The ligand exchange reaction does not reach completion until the protein backbone reaches a largely unfolded state, as monitored through 1H-15N NMR experiments, thus demonstrating that there is a significant correlation between formation of the Fe-S bond and native structure stability. 1H/2H exchange data, however, show that helix alpha(3), the C-terminal region of helix alpha(4), and helix alpha(5) maintain low exchangeability of the amide protons in the LS(1) form. This finding most likely implies that these regions maintain some ordered non-covalent structure, in which the amide moieties are involved in H-bonds. Finally, a folding mechanism is proposed and discussed in terms of analogies and differences with the larger mitochondrial cytochrome c proteins. It is concluded that the thermodynamic stability of the region around the metal cofactor is determined by the chemical nature of the residues around the axial methionine residue.     

Solution structure of the phosphocarrier protein HPr from Bacillus subtilis by two-dimensional NMR spectroscopy.

The solution structure of the phosphocarrier protein, HPr, from Bacillus subtilis has been determined by analysis of two-dimensional (2D) NMR spectra acquired for the unphosphorylated form of the protein. Inverse-detected 2D (1H-15N) heteronuclear multiple quantum correlation nuclear Overhauser effect (HMQC NOESY) and homonuclear Hartmann-Hahn (HOHAHA) spectra utilizing 15N assignments (reported here) as well as previously published 1H assignments were used to identify cross-peaks that are not resolved in 2D homonuclear 1H spectra. Distance constraints derived from NOESY cross-peaks, hydrogen-bonding patterns derived from 1H-2H exchange experiments, and dihedral angle constraints derived from analysis of coupling constants were used for structure calculations using the variable target function algorithm, DIANA. The calculated models were refined by dynamical simulated annealing using the program X-PLOR. The resulting family of structures has a mean backbone rmsd of 0.63 A (N, C alpha, C', O atoms), excluding the segments containing residues 45-59 and 84-88. The structure is comprised of a four-stranded antiparallel beta-sheet with two antiparallel alpha-helices on one side of the sheet. The active-site His 15 residue serves as the N-cap of alpha-helix A, with its N delta 1 atom pointed toward the solvent to accept the phosphoryl group during the phosphotransfer reaction with enzyme I. The existence of a hydrogen bond between the side-chain oxygen atom of Tyr 37 and the amide proton of Ala 56 is suggested, which may account for the observed stabilization of the region that includes the beta-turn comprised of residues 37-40. If the beta alpha beta beta alpha beta (alpha) folding topology of HPr is considered with the peptide chain polarity reversed, the protein fold is identical to that described for another group of beta alpha beta beta alpha beta proteins that include acylphosphatase and the RNA-binding domains of the U1 snRNP A and hnRNP C proteins.     

Stable intermediate states and high energy barriers in the unfolding of GFP

We present a study of the denaturation of a truncated, cycle3 variant of green fluorescent protein (GFP). Chemical denaturation is used to unfold the protein, with changes in structure being monitored by the green fluorescence, tyrosine fluorescence and far-UV circular dichroism. The results show that the denaturation behaviour of GFP is complex compared to many small proteins: equilibrium is established only very slowly, over the time course of weeks, suggesting that there are high folding/unfolding energy barriers. Unfolding kinetics confirm that the rates of unfolding at low concentrations of denaturant are very low, consistent with the slow establishment of the equilibrium. In addition, we find that GFP significantly populates an intermediate state under equilibrium conditions, which is compact and stable with respect to the unfolded state (m(IU)=4.6 kcal mol(-1) M(-1) and Delta G(IU)=12.5 kcal mol(-1)). The global and local stability of GFP was probed further by measuring the hydrogen/deuterium (H/D) NMR exchange rates of more than 157 assigned amide protons. Analysis at two different values of pH showed that amide protons within the beta-barrel structure exchange at the EX2 limit, consequently, free energies of exchange could be calculated and compared to those obtained from the denaturation-curve studies providing further support for the three-state model and the existence of a stable intermediate state. Analysis reveals that amide protons in beta-strands 7, 8, 9 and 10 have, on average, higher exchange rates than others in the beta-barrel, suggesting that there is greater flexibility in this region of the protein. Forty or so amide protons were found which do not undergo significant exchange even after several months and these are clustered into a core region encompassing most of the beta-strands, at least at one end of the barrel structure. It is likely that these residues play an important role in stabilizing the structure of the intermediate state. The intermediate state observed in the chemical denaturation studies described here, is similar to that observed at pH 4 in other studies.     

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.     

Structure and dynamics of the potato carboxypeptidase inhibitor by 1H and 15N NMR

The solution structure and backbone dynamics of the recombinant potato carboxypeptidase inhibitor (PCI) have been characterized by NMR spectroscopy. The structure, determined on the basis of 497 NOE-derived distance constraints, is much better defined than the one reported in a previous NMR study, with an average pairwise backbone root-mean-square deviation of 0.5 A for the well-defined region of the protein, residues 7-37. Many of the side-chains show now well-defined conformations, both in the hydrophobic core and on the surface of the protein. Overall, the solution structure of free PCI is similar to the one that it shows in the crystal of the complex with carboxypeptidase A. However, some local differences are observed in regions 15-21 and 27-29. In solution, the six N-terminal and the two C-terminal residues are rather flexible, as shown by 15N backbone relaxation measurements. The flexibility of the latter segment may have implications in the binding of the inhibitor by the enzyme. All the remaining residues in the protein are essentially rigid (S2 > 0.8) with the exception of two of them at the end of a short 3/10 helix. Despite the small size of the protein, a number of amide protons are protected from exchange with solvent deuterons. The slowest exchanging protons are those in a small two-strand beta-sheet. The unfolding free energies, as calculated from the exchange rates of these protons, are around 5 kcal/mol. Other protected amide protons are located in the segment 7-12, adjacent to the beta-sheet. Although these residues are not in an extended conformation in PCI, the equivalent residues in structurally homologous proteins form a third strand of the central beta-sheet. The amide protons in the 3/10 helix are only marginally protected, indicating that they exchange by a local unfolding mechanism, which is consistent with the increase in flexibility shown by some of its residues. Backbone alignment-based programs for folding recognition, as opposite to disulfide-bond alignments, reveal new proteins of unrelated sequence and function with a similar structure.     

Hydrogen-deuterium exchange in free and prodomain-complexed subtilisin

Residue-specific exchange rates of 223 amide protons in free and prodomain-complexed subtilisin were determined in order to understand how the prodomain binding affects the energetics of subtilisin folding. In free subtilisin, amide protons can be categorized according to exchange rate: 74 fast exchangers (rates > or = 1 h(-1)); 52 medium exchangers (rates between 1 h(-1) and 1 day(-1)); 31 slow exchangers (rates between 1 day(-1) and 0.001 day(-1)). The remaining 66 amide proteins did not exchange detectibly over 9 months (k(obs) < year(-1)) and were denoted as core protons. Core residues occur throughout the main structural elements of subtilisin. Prodomain binding results in high protection factors (100-1000) in the central beta-sheet, particularly in the vicinity of beta-strands S5, S6, and S7 and the connecting loops between them. These connecting loops provide the ligands to the cation at metal site B. Overall, prodomain binding seems to facilitate the organization of the entire central beta-sheet and alpha-helix C in the left-handed crossover connection between beta-strands two and three. It also appears to facilitate the isomerization of multiple prolines late in folding, allowing the formation of metal site B. The gain of stability region around site B comes at the cost of stability in regions more distal to prodomain binding: the C-terminal alpha-helix H and the N-terminal alpha-helices A and B. The acceleration of exchange in these regions by prodomain binding reveals an antagonism between the folding intermediate and the full native structure. This antagonism helps to explain why the prodomain is needed to stabilize the folding intermediate as well as why the unfolding of free subtilisin seldom occurs via this intermediate.     

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.