Insertion of a chaperone domain converts FKBP12 into a powerful catalyst of protein folding

The catalytic activity of human FKBP12 as a prolyl isomerase is high towards short peptides, but very low in proline-limited protein folding reactions. In contrast, the SlyD proteins, which are members of the FKBP family, are highly active as folding enzymes. They contain an extra "insert-in-flap" or IF domain near the prolyl isomerase active site. The excision of this domain did not affect the prolyl isomerase activity of SlyD from Escherichia coli towards short peptide substrates but abolished its catalytic activity in proline-limited protein folding reactions. The reciprocal insertion of the IF domain of SlyD into human FKBP12 increased its folding activity 200-fold and generated a folding catalyst that is more active than SlyD itself. The IF domain binds to refolding protein chains and thus functions as a chaperone module. A prolyl isomerase catalytic site and a separate chaperone site with an adapted affinity for refolding protein chains are the key elements for a productive coupling between the catalysis of prolyl isomerization and conformational folding in the enzymatic mechanisms of SlyD and other prolyl isomerases, such as trigger factor and FkpA.     

Crystal Structure Determination and Functional Characterization of the Metallochaperone SlyD from Thermus thermophilus

SlyD (sensitive to lysis D; product of the slyD gene) is a prolyl isomerase [peptidyl-prolyl cis/trans isomerase (PPIase)] of the FK506 binding protein (FKBP) type with chaperone properties. X-ray structures derived from three different crystal forms reveal that SlyD from Thermus thermophilus consists of two domains representing two functional units. PPIase activity is located in a typical FKBP domain, whereas chaperone function is associated with the autonomously folded insert-in-flap (IF) domain. The two isolated domains are stable and functional in solution, but the presence of the IF domain increases the PPIase catalytic efficiency of the FKBP domain by 2 orders of magnitude, suggesting that the two domains act synergistically to assist the folding of polypeptide chains. The substrate binding surface of SlyD from T. thermophilus was mapped by NMR chemical shift perturbations to hydrophobic residues of the IF domain, which exhibits significantly reduced thermodynamic stability according to NMR hydrogen/deuterium exchange and fluorescence equilibrium transition experiments. Based on structural homologies, we hypothesize that this is due to the absence of a stabilizing β-strand, suggesting in turn a mechanism for chaperone activity by ‘donor-strand complementation.’ Furthermore, we identified a conserved metal (Ni²⁺) binding site at the C-terminal SlyD-specific helical appendix of the FKBP domain, which may play a role in metalloprotein assembly.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

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.     

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.     

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.     

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.     

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.     

How do chemical denaturants affect the mechanical folding and unfolding of proteins?

We present the first single-molecule atomic force microscopy study on the effect of chemical denaturants on the mechanical folding/unfolding kinetics of a small protein GB1 (the B1 immunoglobulin-binding domain of protein G from Streptococcus). Upon increasing the concentration of the chemical denaturant guanidinium chloride (GdmCl), we observed a systematic decrease in the mechanical stability of GB1, indicating the softening effect of the chemical denaturant on the mechanical stability of proteins. This mechanical softening effect originates from the reduced free-energy barrier between the folded state and the unfolding transition state, which decreases linearly as a function of the denaturant concentration. Chemical denaturants, however, do not alter the mechanical unfolding pathway or shift the position of the transition state for mechanical unfolding. We also found that the folding rate constant of GB1 is slowed down by GdmCl in mechanical folding experiments. By combining the mechanical folding/unfolding kinetics of GB1 in GdmCl solution, we developed the “mechanical chevron plot” as a general tool to understand how chemical denaturants influence the mechanical folding/unfolding kinetics and free-energy diagram in a quantitative fashion. This study demonstrates great potential in combining chemical denaturation with single-molecule atomic force microscopy techniques to reveal invaluable information on the energy landscape underlying protein folding/unfolding reactions.     

The equilibrium unfolding intermediate observed at pH 4 and its relationship with the kinetic folding intermediates in green fluorescent protein

Folding mechanisms of a variant of green fluorescent protein (F99S/M153T/V163A) were investigated by a wide variety of spectroscopic techniques. Equilibrium measurements on acid-induced denaturation of the protein monitored by chromophore and tryptophan fluorescence and small-angle X-ray scattering revealed that this protein accumulates at least two equilibrium intermediates, a native-like intermediate and an unfolding intermediate, the latter of which exhibits the characteristics of the molten globule state under moderately denaturing conditions at pH 4. To elucidate the role of the equilibrium unfolding intermediate in folding, a series of kinetic refolding experiments with various combinations of initial and final pH values, including pH 7.5 (the native condition), pH 4.0 (the moderately denaturing condition where the unfolding intermediate is accumulated), and pH 2.0 (the acid-denaturing condition) were carried out by monitoring chromophore and tryptophan fluorescence. Kinetic on-pathway intermediates were accumulated during the folding on the refolding reaction from pH 2.0 to 7.5. However, the signal change corresponding to the conversion from the acid-denatured to the kinetic intermediate states was significantly reduced on the refolding reaction from pH 4.0 to pH 7.5, whereas only the signal change corresponding to the above conversion was observed on the refolding reaction from pH 2.0 to pH 4.0. These results indicate that the equilibrium unfolding intermediate is composed of an ensemble of the folding intermediate species accumulated during the folding reaction, and thus support a hierarchical model of protein folding.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Folding kinetics and thermodynamics of Pseudomonas syringae effector protein AvrPto provide insight into translocation via the type III secretion system

In order to infect their hosts, many Gram-negative bacteria translocate agents of infection, called effector proteins, through the type III secretion system (TTSS) into the host cytoplasm. This process is thought to require at least partial unfolding of these agents, raising the question of how an effector protein might unfold to enable its translocation and then refold once it reaches the host cytoplasm. AvrPto is a well-studied effector protein of Pseudomonas syringae pv tomato. The presence of a readily observed unfolded population of AvrPto in aqueous solution and the lack of a known secretion chaperone make it ideal for studying the kinetic and thermodynamic characteristics that facilitate translocation. Application of Nzz exchange spectroscopy revealed a global, two-state folding equilibrium with 16% unfolded population, a folding rate of 1.8 s(-1), and an unfolding rate of 0.33 s(-1) at pH 6.1. TrAvrPto stability increases with increasing pH, with only 2% unfolded population observed at pH 7.0. The R(1) relaxation of TrAvrPto, which is sensitive to both the global anisotropy of folded TrAvrPto and slow exchange between folded and unfolded conformations, provided independent verification of the global kinetic rate constants. Given the acidic apoplast in which the pathogen resides and the more basic host cytoplasm, these results offer an intriguing mechanism by which the pH dependence of stability and slow folding kinetics of AvrPto would allow efficient translocation of the unfolded form through the TTSS and refolding into its functional folded form once inside the host.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Stability of proteins in the presence of polyols estimated from their guanidinium chloride-induced transition curves at different pH values and 25 degrees C

We have recently concluded from the heat-induced denaturation studies that polyols do not affect deltaG(D) degrees (the Gibbs free energy change (deltaG(D)) at 25 degrees C) of ribonuclease-A and lysozyme at physiological pH and temperature, and their stabilizing effect increases with decrease in pH. Since the estimation of deltaG(D) degrees of proteins from heat-induced denaturation curves requires a large extrapolation, the reliability of this procedure for the estimation of deltaG(D) degrees is always questionable, and so are conclusions drawn from such studies. This led us to measure deltaG(D) degrees of ribonuclease-A and lysozyme using a more accurate method, i.e., from their isothermal (25 degrees C) guanidinium chloride (GdmCl)-induced denaturations. We show that our earlier conclusions drawn from heat-induced denaturation studies are correct. Since the extent of unfolding of heat- and GdmCl-induced denatured states of these proteins is not identical, the extent of stabilization of the proteins by polyols against heat and GdmCl denaturations may also differ. We report that in spite of the differences in the structural nature of the heat- and GdmCl-denatured states of each protein, the extent of stabilization by a polyol is same. We also report that the functional dependence of deltaG(D) of proteins in the presence of polyols on denaturant concentration is linear through the full denaturant concentration range. Furthermore, polyols do not affect the secondary and tertiary structures of the native and GdmCl-denatured states.     

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Mutational analysis of the stability of the H2A and H2B histone monomers

The eukaryotic histone heterodimer H2A-H2B folds through an obligatory dimeric intermediate that forms in a nearly diffusion-limited association reaction in the stopped-flow dead time. It is unclear whether there is partial folding of the isolated monomers before association. To address the possible contributions of structure in the monomers to the rapid association, we characterized H2A and H2B monomers in the absence of their heterodimeric partner. By far-UV circular dichroism, the H2A and H2B monomers are 15% and 31% helical, respectively—significantly less than observed in X-ray crystal structures. Acrylamide quenching of the intrinsic Tyr fluorescence was indicative of tertiary structure. The H2A and H2B monomers exhibit free energies of unfolding of 2.5 and 2.9 kcal mol^− 1, respectively; at 10 μM, the sum of the stability of the monomers is ∼ 60% of the stability of the native dimer. The helical content, stability, and m values indicate that H2B has a more stable, compact structure than H2A. The monomer m values are larger than expected for the extended histone fold motif, suggesting that the monomers adopt an overly collapsed structure. Stopped-flow refolding—initiated from urea-denatured monomers or the partially folded monomers populated at low denaturant concentrations—yielded essentially identical rates, indicating that monomer folding is productive in the rapid association and folding of the heterodimer. A series of Ala and Gly mutations were introduced into H2A and H2B to probe the importance of helix propensity on the structure and stability of the monomers. The mutational studies show that the central α-helix of the histone fold, which makes extensive intermonomer contacts, is structured in H2B but only partially folded in H2A.     

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.