Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.     

Conformational stability and thermodynamic characterization of the lipoic acid bearing domain of human mitochondrial branched chain alpha-ketoacid dehydrogenase

The lipoic acid bearing domain (hbLBD) of human mitochondrial branched chain alpha-ketoacid dehydrogenase (BCKD) plays important role of substrate channeling in oxidative decarboxylation of the branched chain alpha-ketoacids. Recently hbLBD has been found to follow two-step folding mechanism without detectable presence of stable or kinetic intermediates. The present study describes the conformational stability underlying the folding of this small beta-barrel domain. Thermal denaturation in presence of urea and isothermal urea denaturation titrations are used to evaluate various thermodynamic parameters defining the equilibrium unfolding. The linear extrapolation model successfully describes the two-step; native state <-->denatured state unfolding transition of hbLBD. The average temperature of maximum stability of hbLBD is estimated as 295.6 +/- 0.9 K. Cold denaturation of hbLBD is also predicted and discussed.     

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

Osmolyte effects on kinetics of FKBP12 C22A folding coupled with prolyl isomerization

Unfolding and refolding kinetics of human FKBP12 C22A were monitored by fluorescence emission over a wide range of urea concentration in the presence and absence of protecting osmolytes glycerol, proline, sarcosine and trimethylamine-N-oxide (TMAO). Unfolding is well described by a mono-exponential process, while refolding required a minimum of two exponentials for an adequate fit throughout the urea concentration range considered. The bi-exponential behavior resulted from complex coupling between protein folding, and prolyl isomerization in the denatured state in which the urea-dependent rate constant for folding was greater than, equal to, and less than the rate constants for prolyl isomerization within the urea concentration range of zero to five molar. Amplitudes and the observed folding and unfolding rate constants were fitted to a reversible three-state model composed of two sequential steps involving the native state and a folding-competent denatured species thermodynamically linked to a folding-incompetent denatured species. Excellent agreement between thermodynamic parameters for FKBP12 C22A folding calculated from the kinetic parameters and those obtained directly from equilibrium denaturation assays provides strong support for the applicability of the mechanism, and provides evidence that FKBP12 C22A folding/unfolding is two-state, with prolyl isomer heterogeneity in the denatured ensemble. Despite the chemical diversity of the protecting osmolytes, they all exhibit the same kinetic behavior of increasing the rate constant of folding and decreasing the rate constant for unfolding. Osmolyte effects on folding/unfolding kinetics are readily explained in terms of principles established in understanding osmolyte effects on protein stability. These principles involve the osmophobic effect, which raises the Gibbs energy of the denatured state due to exposure of peptide backbone, thereby increasing the folding rate. This effect also plays a key role in decreasing the unfolding rate when, as is often the case, the activated complex exposes more backbone than is exposed in the native state.     

Identification of equilibrium and kinetic intermediates involved in folding of urea-denatured creatine kinase

The unfolding transition and kinetic refolding of dimeric creatine kinase after urea denaturation were monitored by intrinsic fluorescence and far ultraviolet circular dichroism. An equilibrium intermediate and a kinetic folding intermediate were identified and characterized. The fluorescence intensity of the equilibrium intermediate is close to that of the unfolded state, whereas its ellipticity at 222 nm is about 50% of the native state. The transition curves measured by these two methods are therefore non-coincident. The kinetic folding intermediate, formed during the burst phase of refolding under native-like conditions, possesses 75% of the native secondary structure, but is mostly lacking in native tertiary structure. In moderate concentrations of urea, only the initial, rapid change in fluorescence intensity or negative ellipticity is observed, and the final state values do not reach the equivalent unfolding values. The unfolding and refolding transition curves measured under identical conditions are non-coincident within the transition from intermediate to fully unfolded state. It is observed by SDS-PAGE that disulfide bond-linked dimeric or oligomeric intermediates are formed in moderate urea concentrations, especially in the refolding reaction. These rapidly formed, soluble intermediates represent an off-pathway event that leads to the hysteresis in the refolding transition curves.     

Conformational plasticity of cryptolepain: accumulation of partially unfolded states in denaturants induced equilibrium unfolding

pH and chemical denaturant dependent conformational changes of a serine protease cryptolepain from Cryptolepis buchanani are presented in this paper. Activity measurements, near UV, far UV CD, fluorescence emission spectroscopy, and ANS binding studies have been carried out to understand the folding mechanism of the protein in the presence of denaturants. pH and chemical denaturants have a marked effect on the stability, structure, and function of many globular proteins due to their ability to influence the electrostatic interactions. The preliminary biophysical study on cryptolepain shows that major elements of secondary structure are beta-sheets. Under neutral conditions the enzyme was stable in urea while GuHCl-induced equilibrium unfolding was cooperative. Cryptolepain shows little ANS binding even under neutral conditions due to more hydrophobicity of beta-sheets. Multiple intermediates were populated during the pH-induced unfolding of cryptolepain. Temperature-induced denaturation of cryptolepain in the molten globule like state is non-cooperative, contrary to the cooperativity seen with the native protein, suggesting the presence of two parts, possibly domains, in the molecular structure of cryptolepain, with different stability that unfolds in steps. Interestingly, the GuHCl-induced unfolding of A state (molten globule state) of cryptolepain is unique, as lower concentration of denaturant, not only induces structure but also facilitate transition from one molten globule like state (MG(1)) into another (MG(2)). The increase of pH drives the protein into alkaline denatured state characterized by the absence of any ANS binding. GuHCl- and urea-induced unfolding transition curves at pH 12.0 were non-coincidental indicating the presence of an intermediate in the unfolding pathway.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

Effect of Substitutions in Surface Amino Acid on Energy Profile of Apomyoglobin

Studies on the process of spontaneous protein folding into a unique native state are an important issue of molecular biology. Apomyoglobin from the sperm whale is a convenient model for these studies in vitro. Here, we present the results of equilibrium and kinetic experiments carried out in a study on the folding and unfolding of eight mutant apomyoglobin forms of with hydrophobic amino acid substitutions on the protein surface. Calculated values of apparent constants of folding/unfolding rates, as well as the data on equilibrium conformational transitions in the urea concentration range of 0–6 M at 11°C are given. Based on the obtained information on the kinetic properties of the studied proteins, a Φ-value analysis of the transition state has been performed and values of urea concentrations corresponding to the midpoint of the transition from the native to intermediate state have been determined for the given forms of mutant apomyoglobin. It has been found that a significant increase in the stability of the native state can be achieved by a small number of amino acid substitutions on the protein surface. It has been shown that the substitution of only one amino acid residue exclusively affects the height of the energy barrier that separates different states of apomyoglobin.     

Role of a solvent-exposed aromatic cluster in the folding of Escherichia coli CspA

Escherichia coli CspA is a member of the cold shock protein family. All cold shock proteins studied to date fold rapidly by an apparent two-state mechanism. CspA contains an unusual cluster of aromatic amino acids on its surface that is necessary for nucleic acid binding and also provides stability to CspA (Hillier et al., 1998). To elucidate the role this aromatic cluster plays in the determining the folding rate and pathway of CspA, we have studied the folding kinetics of mutants containing either leucine or serine substituted for Phe 18, Phe20, and/or Phe31. The leucine substitutions are found to accelerate folding and the serine substitutions to decelerate folding. Because these residues exert effects on the free energy of the folding transition state, they may be necessary for nucleating folding. They are not responsible, however, for the very compact, native-like transition state ensemble seen in the cold shock proteins, as the refolding rates of the mutants all show a similar, weak dependence of unfolding rate on denaturant concentration. Using mutant cycle analysis, we show that there is energetic coupling among the three residues between the unfolded and transition states, suggesting that the cooperative nature of these interactions helps to determine the unfolding rate. Overall, our results suggest that separate evolutionary pressures can act simultaneously on the same group of residues to maintain function, stability, and folding rate.     

Native state EX2 and EX1 hydrogen exchange of Escherichia coli CspA, a small beta-sheet protein

Escherichia coli CspA is a small all-beta-sheet protein that folds fast (tau = 4 ms) via an apparent two-state mechanism. Our previous studies have shown that a large aromatic cluster on the surface of the protein participates in the rate-limiting step of folding and thus may be part of the folding nucleus of this protein. To obtain a more detailed picture of molecular events at the peptide backbone during unfolding and folding of CspA, we used native state hydrogen exchange and nuclear magnetic resonance spectroscopy (NMR). The experiments with native CspA were performed over a range of pH values from low pH, where exchange is governed by a rapid equilibrium before chemical exchange (EX2 exchange), to high pH, where exchange is dictated by the rate of unfolding (EX1 exchange). Rates of folding and unfolding were determined for 11 residues. The distribution of rates of folding within the structure of CspA suggests that hairpin turns, including one near the aromatic cluster, may nucleate the folding of CspA.     

Measuring the stability of partly folded proteins using TMAO

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.     

Folding studies on a knotted protein

YibK is a 160 residue homodimeric protein belonging to the SPOUT class of methyltransferases. Proteins in this group all display a unique topological feature; the backbone polypeptide chain folds to form a deep trefoil knot. Such knotted structures were completely unpredicted, it being thought impossible for a protein to fold efficiently in this way. However, they are becoming more common and there are now a growing number of examples in the Protein Data Bank. These intriguing knotted structures represent a new and significant challenge in the field of protein folding. Here, we present an initial characterisation of the folding of YibK, one of the smallest knotted proteins to be identified. This is the first detailed folding study on a knotted protein to be reported. We have established conditions under which the protein can be denatured reversibly in vitro using urea, thereby showing that molecular chaperones are not required for the efficient folding of this protein. A series of equilibrium unfolding experiments were performed over a 400-fold range of protein concentration. Both secondary and tertiary structural probes show a single, protein concentration-dependent unfolding transition, and data are most consistent with a three-state equilibrium denaturation model involving a monomeric intermediate. Thermodynamic parameters obtained from the fit of the data to this model indicate that the intermediate is a stable species with appreciable secondary and tertiary structure; whether the topological knot remains in the intermediate state is still to be shown. Together, these results demonstrate that, despite its complex knotted structure, YibK is able to fold efficiently and behaves remarkably similarly to other dimeric proteins under equilibrium conditions.     

Investigation of folding/unfolding kinetics of apomyoglobin

Apomyoglobin kinetic and equilibrium unfolding and folding processes were studied at pH 6.2, 11 degrees C by stopped-flow tryptophan fluorescence. There are two distinct consecutive processes in apomyoglobin folding process, namely, the protein fast transition between the unfolded (U) and an intermediate (I) states (U <----> I) and slow transition between the intermediate and the native (N) states (I <----> N). Accumulation of the intermediate state was observed in the wide range of urea concentrations. The presence of the intermediate state was shown even beyond the middle transition on the unfolding limb. The dependence of observed folding/unfolding rates on urea concentration (chevron plot) was obtained. The shape of this dependence was compared with that of two-state proteins, folding from the U to N state.     

Folding and stability of trp aporepressor from Escherichia coli

Equilibrium and kinetic studies of the urea-induced unfolding of trp aporepressor from Escherichia coli were performed to probe the folding mechanism of this intertwined, dimeric protein. The equilibrium unfolding transitions at pH 7.6 and 25 degrees C monitored by difference absorbance, fluorescence, and circular dichroism spectroscopy are coincident within experimental error. All three transitions are well described by a two-state model involving the native dimer and the unfolded monomer; the free energy of folding in the absence of denaturant and under standard-state conditions is estimated to be 23.3 +/- 0.9 kcal/mol of dimer. The midpoint of the equilibrium unfolding transition increases with increasing protein concentration in the manner expected from the law of mass action for the two-state model. We find no evidence for stable folding intermediates. Kinetic studies reveal that unfolding is governed by a single first-order reaction whose relaxation time decreases exponentially with increasing urea concentration and also decreases with increasing protein concentration in the transition zone. Refolding involves at least three phases that depend on both the protein concentration and the final urea concentration in a complex manner. The relaxation time of the slowest of these refolding phases is identical with that for the single phase in unfolding in the transition zone, consistent with the results expected for a reaction that is kinetically reversible. The two faster refolding phases are presumed to arise from slow isomerization reactions in the unfolded form and reflect parallel folding channels.     

The rough energy landscape of superfolder GFP is linked to the chromophore

Many green fluorescent protein (GFP) variants have been developed for use as fluorescent tags, and recently a superfolder GFP (sfGFP) has been developed as a robust folding reporter. This new variant shows increased stability and improved folding kinetics, as well as 100% recovery of native protein after denaturation. Here, we characterize sfGFP, and find that this variant exhibits hysteresis as unfolding and refolding equilibrium titration curves are non-coincident even after equilibration for more than eight half-lives as estimated from kinetic unfolding and refolding studies. This hysteresis is attributed to trapping in a native-like intermediate state. Mutational studies directed towards inhibiting chromophore formation indicate that the novel backbone cyclization is responsible for the hysteresis observed in equilibrium titrations of sfGFP. Slow equilibration and the presence of intermediates imply a rough landscape. However, de novo folding in the absence of the chromophore is dominated by a smoother energy landscape than that sampled during unfolding and refolding of the post-translationally modified polypeptide.     

Unfolding and refolding of bovine beta-lactoglobulin monitored by hydrogen exchange measurements.

Bovine beta-lactoglobulin (beta-LG) is a widely studied protein belonging to the lipocalin family, whose structural characterisation has been reported by X-ray crystallography and NMR studies at physiological and acidic pH, respectively. Bovine beta-LG consists of nine antiparallel beta-sheets and a terminal alpha-helix segment. The beta-sheets form a calyx structure with a hydrophobic buried cluster conferring stability to the protein while a hydrophobic surface patch provides stabilising interactions between the barrel and the flanking terminal helix. Here, the stability and the folding properties of bovine beta-LG in the presence of a chemical denaturant is probed. The analysis of the NMR spectra recorded in aqueous solution with increasing amounts of urea revealed that the intensities of the backbone cross-peaks decrease upon increasing urea concentration, while their secondary shifts do not change significantly on going from 0 to 5 M urea, thus suggesting the presence of slow exchange between native and unfolded protein. Hydrogen exchange measurements at different urea concentrations were performed in order to evaluate the exchange rates of individual backbone amide protons. The opening reactions that determine protein exchange can be computed for the most slowly exchanging hydrogen atoms, and the measured exchange rates and the corresponding free energies can be expressed in terms of the equilibrium energetic for the global transition and the local fluctuations. Most of the residues converge to define a common isotherm identifying a unique cooperative folding unit, encompassing all the strands, except strand betaI, and the terminal region of the helix. The amides that do not join the same global unfolding isotherm are characterised by low DeltaGH20op and especially by low m values, indicating that they are already substantially exposed in the native state. A two-state unfolding model N <==> U is therefore proposed for this rather big protein, in agreement with CD data. Renaturation studies show that the unfolding mechanism is reversible up to 6 M urea and suggest a similar unfolding and refolding pathway. Present results are discussed in light of the hypothesis of an alpha-->beta transition proposed for bovine beta-LG refolding.     

Protein Folding Mechanism of the Dimeric AmphiphysinII/Bin1 N-BAR Domain

The human AmphyphisinII/Bin1 N-BAR domain belongs to the BAR domain superfamily, whose members sense and generate membrane curvatures. The N-BAR domain is a 57 kDa homodimeric protein comprising a six helix bundle. Here we report the protein folding mechanism of this protein as a representative of this protein superfamily. The concentration dependent thermodynamic stability was studied by urea equilibrium transition curves followed by fluorescence and far-UV CD spectroscopy. Kinetic unfolding and refolding experiments, including rapid double and triple mixing techniques, allowed to unravel the complex folding behavior of N-BAR. The equilibrium unfolding transition curve can be described by a two-state process, while the folding kinetics show four refolding phases, an additional burst reaction and two unfolding phases. All fast refolding phases show a rollover in the chevron plot but only one of these phases depends on the protein concentration reporting the dimerization step. Secondary structure formation occurs during the three fast refolding phases. The slowest phase can be assigned to a proline isomerization. All kinetic experiments were also followed by fluorescence anisotropy detection to verify the assignment of the dimerization step to the respective folding phase. Based on these experiments we propose for N-BAR two parallel folding pathways towards the homodimeric native state depending on the proline conformation in the unfolded state.     

Revealing a Concealed Intermediate that Forms after the Rate-limiting Step of Refolding of the SH3 Domain of PI3 Kinase

Kinetic and equilibrium studies of the folding and unfolding of the SH3 domain of the PI3 kinase, have been used to identify a folding intermediate that forms after the rate-limiting step on the folding pathway. Folding and unfolding, in urea as well as in guanidine hydrochloride (GdnHCl), were studied by monitoring changes in the intrinsic fluorescence or in the far-UV circular dichroism (CD) of the protein. The two probes yield non-coincident equilibrium transitions for unfolding in urea, indicating that an intermediate, I, exists in equilibrium with native (N) and unfolded (U) protein, during unfolding. Hence, the equilibrium unfolding data were analyzed according to a three-state N ↔ I ↔ U mechanism. An intermediate is observed also in kinetic unfolding studies, and its presence leads to the unfolding reaction in urea as well as in GdnHCl, occurring in two steps. The fast step is complete within the initial 11 ms of unfolding and manifests itself in a burst phase change in fluorescence. At high concentrations of GdnHCl, the entire change in fluorescence during unfolding occurs during the 11 ms burst phase. CD measurements indicate, however, that I retains N-like secondary structure. An analysis of the kinetic and thermodynamic data, according to a minimal three-state N ↔ I ↔ U mechanism, positions I after the rate-limiting transition state, TS1, of folding, on the reaction coordinate of folding in GdnHCl. Hence, I is not revealed when folding is commenced from U, regardless of the nature of the probe used to follow the folding reaction. Interrupted unfolding experiments, in which the protein is unfolded transiently in GdnHCl for various lengths of time before being refolded, showed that I refolds to N much faster than does U, confirms the analysis of the direct folding and unfolding experiments, that I is formed after the rate-limiting step of refolding in GdnHCl.     

Role of solvation barriers in protein kinetic stability

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).