Kinetic mechanism and catalysis of a native-state prolyl isomerization reaction

Folding of tendamistat is a rapid two-state process for the majority of the unfolded molecules. In fluorescence-monitored refolding kinetics about 8% of the unfolded molecules fold slowly (lambda=0.083s(-1)), limited by peptidyl-prolyl cis-trans isomerization. This is significantly less than expected from the presence of three trans prolyl-peptide bonds in the native state. In interrupted refolding experiments we detected an additional very slow folding reaction (lambda=0.008s(-1) at pH 2) with an amplitude of about 12%. This reaction is caused by the interconversion of a highly structured intermediate to native tendamistat. The intermediate has essentially native spectroscopic properties and about 2% of it remain populated in equilibrium after folding is complete. Catalysis by human cyclophilin 18 identifies this very slow reaction as a prolyl isomerization reaction. This shows that prolyl-isomerases are able to efficiently catalyze native state isomerization reactions, which allows them to influence biologically important regulatory conformational transitions. Folding kinetics of the proline variants P7A, P9A, P50A and P7A/P9A show that the very slow reaction is due to isomerization of the Glu6-Pro7 and Ala8-Pro9 peptide bonds, which are located in a region that makes strong backbone and side-chain interactions to both beta-sheets. In the P50A variant the very slow isomerization reaction is still present but native state heterogeneity is not observed any more, indicating a long-range destabilizing effect on the alternative native state relative to N. These results enable us to include all prolyl and non-prolyl peptide bond isomerization reactions in the folding mechanism of tendamistat and to characterize the kinetic mechanism and the energetics of a native-state prolyl isomerization reaction.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 β-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 β-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of β-lactamase by almost 20 °C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol^− 1 relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring β-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.     

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.     

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Probing the folding intermediate of Rd-apocyt b562 by protein engineering and infrared T-jump

Small proteins often fold in an apparent two-state manner with the absence of detectable early-folding intermediates. Recently, using native-state hydrogen exchange, intermediates that exist after the rate-limiting transition state have been identified for several proteins. However, little is known about the folding kinetics from these post-transition intermediates to their corresponding native states. Herein, we have used protein engineering and a laser-induced temperature-jump (T-jump) technique to investigate this issue and have applied it to Rd-apocyt b(562) , a four-helix bundle protein. Previously, it has been shown that Rd-apocyt b(562) folds via an on-pathway hidden intermediate, which has only the N-terminal helix unfolded. In the present study, a double mutation (V16G/I17A) in the N-terminal helix of Rd-apocyt b(562) was made to further increase the relative population of this intermediate state at high temperature by selectively destabilizing the native state. In the circular dichroism thermal melting experiment, this mutant showed apparent two-state folding behavior. However, in the T-jump experiment, two kinetic phases were observed. Therefore, these results are in agreement with the idea that a folding intermediate is populated on the folding pathway of Rd-apocyt b(562) . Moreover, it was found that the exponential growth rate of the native state from this intermediate state is roughly (25 microsec)(-1) at 65 degrees C.     

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by F?rster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.     

Kinetic coupling of folding and prolyl isomerization of beta2-microglobulin studied by mutational analysis

β₂-Microglobulin (β2-m), a protein responsible for dialysis-related amyloidosis, adopts a typical immunoglobulin domain fold with the N-terminal peptide bond of Pro32 in a cis isomer. The refolding of β2-m is limited by the slow trans-to-cis isomerization of Pro32, implying that intermediates with a non-native trans-Pro32 isomer are precursors for the formation of amyloid fibrils. To obtain further insight into the Pro-limited folding of β2-m, we studied the Gdn-HCl-dependent unfolding/refolding kinetics using two mutants (W39 and P32V β2-ms) as well as the wild-type β2-m. W39 β2-m is a triple mutant in which both of the authentic Trp residues (Trp60 and Trp95) are replaced by Phe and a buried Trp common to other immunoglobulin domains is introduced at the position of Leu39 (i.e., L39W/W60F/W95F). W39 β2-m exhibits a dramatic quenching of fluorescence upon folding, enabling a detailed analysis of Pro-limited unfolding/refolding. On the other hand, P32V β2-m is a mutant in which Pro32 is replaced by Val, useful for probing the kinetic role of the trans-to-cis isomerization of Pro32. A comparative analysis of the unfolding/refolding kinetics of these mutants including three types of double-jump experiments revealed the prolyl isomerization to be coupled with the conformational transitions, leading to apparently unusual kinetics, particularly for the unfolding. We suggest that careful consideration of the kinetic coupling of unfolding/refolding and prolyl isomerization, which has tended to be neglected in recent studies, is essential for clarifying the mechanism of protein folding and, moreover, its biological significance.     

Nuclear magnetic resonance characterization of the refolding intermediate of beta2-microglobulin trapped by non-native prolyl peptide bond

beta(2)-Microglobulin (beta2-m), a light chain of the major histocompatibility complex type I, is also found as a major component of amyloid fibrils formed in dialysis-related amyloidosis. Denaturation of beta2-m is considered to initiate the formation of fibrils. To clarify the mechanism of fibril formation, it is important to characterize the intermediate conformational states at the atomic level. Here, we investigated the refolding of beta2-m from the acid-unfolded state by heteronuclear magnetic resonance and circular dichroism spectroscopies. At low temperature, beta2-m refolded slowly, accumulating a rate-limiting intermediate with non-native chemical shift dispersions for several residues, but with compactness and secondary structures similar to those of the native protein. beta2-m has a cis proline residue at Pro32, located on the turn connecting the betaB and betaC strands. The slow refolding phase disappeared upon mutation of Pro32 to Val, indicating that Pro32 is responsible for the accumulation of the intermediate. The distribution of the perturbed residues in the intermediate suggests that the non-native prolyl peptide bond of Pro32 affects large areas of the molecule. A cis proline residue is common to various immunoglobulin domains involved in amyloidosis, implying that a non-native prolyl peptide bond that might occur under physiological conditions is related to the amyloidogenicity of these immunoglobulin domains.     

Downhill versus Barrier-Limited Folding of BBL 1: Energetic and Structural Perturbation Effects upon Protonation of a Histidine of Unusually Low pKa

A dispersion of melting temperatures at pH 5.3 for individual residues of the BBL protein domain has been adduced as evidence for barrier-free downhill folding. Other members of the peripheral subunit domain family fold cooperatively at pH 7. To search for possible causes of anomalies in BBL's denaturation behavior, we measured the pH titration of individual residues by heteronuclear NMR. At 298 K, the pK_a of His142 was close to that of free histidine at 6.47 ± 0.04, while that of the more buried His166 was highly perturbed at 5.39 ± 0.02. Protonation of His166 is thus energetically unfavorable and destabilizes the protein by ∼ 1.5 kcal/mol. Changes in C^α secondary shifts at pH 5.3 showed a decrease in helicity of the C-terminus of helix 2, where His166 is located, which was accompanied by a measured decrease of 1.1 ± 0.2 kcal/mol in stability from pH 7 to 5.3. Protonation of His166 perturbs, therefore, the structure of BBL. Only ∼ 1% of the structurally perturbed state will be present at the biologically relevant pH 7.6. Experiments at pH 5.3 report on a near-equal mixture of the two different native states. Further, at this pH, small changes of pH and pK_a induced by changes in temperature will have near-maximal effects on pH-dependent conformational equilibria and on propagation of experimental error. Accordingly, conventional barrier-limited folding predicts some dispersion of measured thermal unfolding curves of individual residues at pH 5.3.     

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.     

In vitro evolution of a hyperstable Gbeta1 variant

An in-vitro selection strategy was used to obtain strongly stabilized variants of the beta1 domain of protein G (Gbeta1). In a two-step approach, first candidate positions with a high potential for stabilization were identified in Gbeta1 libraries that were created by error-prone PCR, and then, after randomization of these positions by saturation mutagenesis, strongly stabilized variants were selected. For both steps the in-vitro selection method Proside was employed. Proside links the stability of a protein with the infectivity of a filamentous phage. Ultimately, residues from the two best selected variants were combined in a single Gbeta1 molecule. This variant with the four mutations E15V, T16L, T18I, and N37L showed an increase of 35.1 degrees C in the transition midpoint and of 28.5 kJ mol(-1) (at 70 degrees C) in the Gibbs free energy of stabilization. It was considerably more stable than the best variant from a previous Proside selection, in which positions were randomized that had originally been identified by computational design. Only a single substitution (T18I) was found in both selections. The best variants from the present selection showed a higher cooperativity of thermal unfolding, as indicated by an increase in the enthalpy of unfolding by about 60 kJ mol(-1). This increase is apparently correlated with the presence of Leu residues that were selected at the positions 16 and 37.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Over-expression of mitochondrial heat shock protein 70 suppresses programmed cell death in rice

In this study, we identified and functionally characterized the mitochondrial heat shock protein 70 (mtHsp70). Over-expression of mtHsp70 suppressed heat- and H₂O₂-induced programmed cell death (PCD) in rice protoplasts, as reflected by higher cell viability, decreased DNA laddering and chromatin condensation. Mitochondrial membrane potential (Δψ_m) after heat shock was destroyed gradually in protoplasts, but mtHsp70 over-expression showed higher Δψ_m relative to the vector control cells, and partially inhibited cytochrome c release from mitochondria to cytosol. Heat treatment also significantly increased reactive oxygen species (ROS) generation, a phenomenon not observed in protoplasts over-expressing mtHsp70. Together, these results suggest that mtHsp70 may suppress PCD in rice protoplasts by maintaining mitochondrial Δψ_m and inhibiting the amplification of ROS.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Purification,crystal structure determination and functional characterization of type III antifreeze proteins from the European eelpout Zoarces viviparus

Antifreeze proteins (AFPs) are essential components of many organisms adaptation to cold temperatures. Fish type III AFPs are divided into two groups, SP isoforms being much less active than QAE1 isoforms. Two type III AFPs from Zoarces viviparus, a QAE1 (ZvAFP13) and an SP (ZvAFP6) isoform, are here characterized and their crystal structures determined. We conclude that the higher activity of the QAE1 isoforms cannot be attributed to single residues, but rather a combination of structural effects. Furthermore both ZvAFP6 and ZvAFP13 crystal structures have water molecules around T18 equivalent to the tetrahedral-like waters previously identified in a neutron crystal structure. Interestingly, ZvAFP6 forms dimers in the crystal, with a significant dimer interface. The presence of ZvAFP6 dimers was confirmed in solution by native electrophoresis and gel filtration. To our knowledge this is the first report of dimerization of AFP type III proteins.     

Reversible control of F1-ATPase rotational motion using a photochromic ATP analog at the single molecule level

Motor enzymes such as F₁-ATPase and kinesin utilize energy from ATP for their motion. Molecular motions of these enzymes are critical to their catalytic mechanisms and were analyzed thoroughly using a single molecule observation technique. As a tool to analyze and control the ATP-driven motor enzyme motion, we recently synthesized a photoresponsive ATP analog with a p-tert-butylazobenzene tethered to the 2′ position of the ribose ring. Using cis/trans isomerization of the azobenzene moiety, we achieved a successful reversible photochromic control over a kinesin-microtubule system in an in vitro motility assay. Here we succeeded to control the hydrolytic activity and rotation of the rotary motor enzyme, F₁-ATPase, using this photosensitive ATP analog. Subsequent single molecule observations indicated a unique pause occurring at the ATP binding angle position in the presence of cis form of the analog.     

Microsecond hydrophobic collapse in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.     

Cooperative unfolding of a metastable serpin to a molten globule suggests a link between functional and folding energy landscapes

Alpha-1 antitrypsin (alpha(1)-AT) is a member of the serpin class of protease inhibitors, and folds to a metastable state rather than its thermodynamically most stable native state. Upon cleavage by a target protease, alpha(1)-AT undergoes a dramatic conformational change to a stable form, translocating the bound protease more than 70 A to form an inhibitory protease-serpin complex. Numerous mutagenesis studies on serpins have demonstrated the trade-off between the stability of the metastable state on the one hand and the inhibitory efficiency on the other. Studies of the equilibrium unfolding of serpins provide insight into this connection between structural plasticity and metastability. We studied equilibrium unfolding of wild-type alpha(1)-AT using hydrogen-deuterium/exchange mass spectrometry to characterize the structure and the stability of an equilibrium intermediate that was observed in low concentrations of denaturant in earlier studies. Our results show that the intermediate observed at low concentrations of denaturant has no protection from hydrogen-deuterium exchange, indicating a lack of stable structure. Further, differential scanning calorimetry of alpha(1)-AT at low concentrations of denaturant shows no heat capacity peak during thermal denaturation, indicating that the transition from the intermediate to the unfolded state is not a cooperative first-order-like phase transition.. Our results show that the unfolding of alpha(1)-AT involves a cooperative transition to a molten globule form, followed by a non-cooperative transition to a random-coil form as more guanidine is added. Thus, the entire alpha(1)-AT molecule consists of one cooperative structural unit rather than multiple structural domains with different stabilities. Furthermore, our results together with previous mutagenesis studies suggest a possible link between an equilibrium molten globule and a functional intermediate that may be populated during the protease inhibition.