On the dynamic nature of the transition state for protein-protein association as determined by double-mutant cycle analysis and simulation

The process of protein-protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein-protein association for three different protein-protein interactions; Barnase-Barstar, TEM1-BLIP and IFNalpha2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant DeltaDeltaG(double dagger)(int) values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant DeltaDeltaG(double dagger)(int) values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental DeltaDeltaG(double dagger)(int) values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental DeltaDeltaG(double dagger)(int) values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNalpha2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.     

Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges

An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-beta-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol(-1 )each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol(-1). Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming net-works of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of DeltaDeltaG(double dagger)K(d) which reflect the strength of short range interactions, while DeltaDeltaG(o)K(D) values which include the effects of long range electrostatic forces that alter specifically DeltaDeltaG(double dagger)k(a) were treated separately.     

Use of protein charge ladders to study electrostatic interactions during protein ultrafiltration

Recent studies have demonstrated the importance of electrostatic interactions in membrane systems, but there is still controversy about the underlying phenomena. Protein charge ladders, consisting of a set of chemical derivatives of a given protein that differ by single charge groups, were used to quantify the electrostatic interactions during protein ultrafiltration. Myoglobin charge ladders were generated by acylation, with the different derivatives analyzed simultaneously by capillary electrophoresis. Filtration experiments were performed using polyethersulfone and composite regenerated cellulose membranes, with the membrane charge determined from the streaming potential. As expected, the rejection increased as the protein became more heavily charged due to the increase in electrostatic repulsion. However, the transmission of the weakly charged myoglobin species increased dramatically at very low ionic strength. This increase in transmission was attributed to a shift in pH within the pore caused by hydrogen ion partitioning into the charged membrane. The sieving data were in good agreement with theoretical calculations accounting for the effects of this pH shift on the electrostatic interactions.     

Fluorous effect in proteins: de novo design and characterization of a four-alpha-helix bundle protein containing hexafluoroleucine

Several studies have demonstrated that proteins incorporating fluorinated analogues of hydrophobic amino acids such as leucine and valine into their hydrophobic cores exhibit increased stability toward thermal denaturation and unfolding by guanidinium chloride. However, estimates for the increase in the thermodynamic stability of a protein (DeltaDeltaG(unfold)) afforded by the substitution of a hydrophobic amino acid with its fluorinated analogue vary quite significantly. To address this, we have designed a peptide that adopts an antiparallel four-helix bundle structure in which the hydrophobic core is packed with leucine, and investigated the effects of substituting the central two layers of the core with L-5,5,5,5',5',5'-hexafluoroleucine (hFLeu). We find that DeltaDeltaG(unfold) is increased by 0.3 kcal/mol per hFLeu residue. This is in good agreement with the predicted increase in DeltaDeltaG(unfold) of 0.4 kcal/mol per residue arising from the increased hydrophobicity of the hFLeu side chain, which we determined experimentally from partitioning measurements on hFLeu and leucine. The increased stability of this fluorinated protein may therefore be ascribed to simple hydrophobic effects, rather than specific "fluorous" interactions between the hFLeu residues.     

The kinetics and thermodynamics of reversible denaturation of crystalline soybean trypsin inhibitor

Crystalline soybean trypsin inhibitor protein undergoes denaturation on heating which is reversed on cooling. In the range of temperature of 35 to 50 degrees C. a solution of the protein consists of a mixture of native and denatured forms in equilibrium with each other. The equilibrium is only slowly established and its final value at any temperature is the same whether a heated, denatured solution of the protein is cooled to the given temperature or whether a fresh solution is raised to that temperature. The kinetics of reversible denaturation of the soybean protein as well as the reversal of denaturation is that of a reversible unimolecular reaction, each process consisting at a given temperature of the same two simultaneous reactions acting in opposite directions. The experimental data on the effect of temperature on the velocity and the equilibrium constants of the opposing reaction were utilized in evaluating the reaction energies and activation energies. The reaction energies for denaturation were found to be as follows:- Change in total heat of reaction DeltaH = 57,000 calories per mole Change in entropy of reaction DeltaS = 180 calories per degree per mole The heat of activation DeltaH(1) (double dagger) for denaturation = 55,000 The heat of activation DeltaH(2) (double dagger) for the reversal of denaturation = -1900 The entropy DeltaS(1) (double dagger) for denaturation = 95 The entropy DeltaS(2) (double dagger) for reversal of denaturation = -84     

Enthalpic denaturation of collagen satisfies the general relation of the change in Delta H(d) from Delta C(p) for proteins at zero hydroxyproline level]

One of the noticeable peculiarities of thermodynamics of collagen is an anomalous high magnitude of the enthalpy of denaturation delta Hd at a very low thermostability. Taking into account recent ideas about the role of hydrophobic interactions in determination of the thermodynamic function of protein denaturation, it is shown that the high magnitude of delta Hd of collagen in comparison with those of globular proteins can be explained by two factors: a significant contribution of residues of 4-hydroxyproline and small magnitude of hydrophobic interactions.     

Dual behavior of sodium dodecyl sulfate as enhancer or suppressor of insulin aggregation and chaperone-like activity of camel αS1-casein

Sodium dodecyl sulfate (SDS) at low concentrations considerably enhanced insulin aggregation and reduced the chaperone-like activity of purified camel αS₁-casein (αS₁-CN). These observed changes were the result of repulsive electrostatic interactions between both negative charged head groups of SDS and αS₁-CN, and the net negative charge of insulin molecules, resulting in the greater exposure of hydrophobic patches of insulin and its enhanced aggregation. In contrast, enhanced hydrophobic interactions were primarily responsible for the conformational changes observed in insulin and αS₁-CN at high SDS concentrations, resulting in increased binding of SDS and αS₁-CN to insulin and its reduced aggregation.     

Quantitative membrane electrostatics with the atomic force microscope

The atomic force microscope (AFM) is sensitive to electric double layer interactions in electrolyte solutions, but provides only a qualitative view of interfacial electrostatics. We have fully characterized silicon nitride probe tips and other experimental parameters to allow a quantitative electrostatic analysis by AFM, and we have tested the validity of a simple analytical force expression through numerical simulations. As a test sample, we have measured the effective surface charge density of supported zwitterionic dioleoylphosphatidylcholine membranes with a variable fraction of anionic dioleoylphosphatidylserine. The resulting surface charge density and surface potential values are in quantitative agreement with those predicted by the Gouy-Chapman-Stern model of membrane charge regulation, but only when the numerical analysis is employed. In addition, we demonstrate that the AFM can detect double layer forces at a separation of several screening lengths, and that the probe only perturbs the membrane surface potential by <2%. Finally, we demonstrate 50-nm resolution electrostatic mapping on heterogeneous model membranes with the AFM. This novel combination of capabilities demonstrates that the AFM is a unique and powerful probe of membrane electrostatics.     

Structure of heat-induced beta-lactoglobulin aggregates and their complexes with sodium-dodecyl sulfate

We report on the conformation of heat-induced bovine beta-lactoglobulin (betalg) aggregates prepared at different pH conditions, and their complexes with model anionic surfactants such as sodium dodecyl sulfate (SDS). The investigation was carried out by combining a wide range of techniques such as ultra small angle light scattering, static and dynamic light scattering, small angle neutron scattering, small-angle X-ray scattering, electrophoretic mobility, isothermal titration calorimetry (ITC) and transmission electron microscopy. Three types of aggregates were generated upon heating betalg aqueous dispersions at increasing pH from 2.0 to 5.8 to 7.0: rod-like aggregates, spherical aggregates, and worm-like primary aggregates, respectively. These aggregates were shown not only to differ for their sizes and morphologies, but also for their internal structures and fractal dimensions. The main differences between aggregates are discussed in terms of the ionic charge and conformational changes arising for betalg at different pHs. The formation of complexes between SDS and the various protein aggregates at pH 3.0 was shown to occur by two main mechanisms: at low concentration of SDS, the complex formation occurs essentially by ionic binding between the positive residues of the protein and the negative sulfate heads of the surfactant. At complete neutralization of charges, precipitation of the complexes is observed. Upon further increase in SDS concentration, complex formation of SDS and the protein aggregates occurs primarily by hydrophobic interactions, leading to (i) the formation of an SDS double layer around the protein aggregates, (ii) the inversion of the total ionic charge of each individual protein aggregate, and (iii) the complete redispersion of the protein aggregate-SDS complexes in water. Remarkably, the SDS double layer around the protein aggregates provides an efficient protective shield, preventing precipitation of the aggregates at any possible pH values, including those values corresponding to the isoelectric pH of the aggregates.     

Hydrophobic and ionic effects upon the electrophoretic mobilities of the subunits of coupling factor 1 from mitochondria

A sodium dodecyl sulfate (SDS)-urea polyacrylamide gel system was used to investigate certain properties of the subunits of the beef heart mitochondrial ATPase, (native F1, nF1). By examining the affects of urea concentration and acrylamide concentration upon the electrophoretic mobilities of the polypeptides comprising the nF1 enzyme, we have obtained conditions under which all five subunits are simultaneously resolved when the discontinuous buffer system of Laemmli is used (U. K. Laemmli (1970) Nature (London) 277, 680-685). The determination of the apparent molecular weights by analysis of Ferguson plots (K. A. Ferguson (1964) Metabolism 13, 985-1002) revealed that the addition of urea to the SDS gels resulted in a decrease in the apparent molecular weight of the beta subunit. A dramatic increase in the apparent molecular weight of the delta subunit was also brought about by the presence of urea in the SDS gels. In addition, the apparent molecular weight of both the alpha and the beta subunits was dependent upon the acrylamide concentration used, indicating that these subunits contain either areas highly resistant to denaturation by the combined action of urea and SDS, or covalent modifications leading to anomalous electrophoretic mobility. The results of experiments in which urea analogs were used indicate that the interactions of urea with the beta subunit involve the formation of hydrogen bonds between urea and regions of this subunit. On the other hand, the interactions of urea with the delta subunit are primarily of a hydrophobic nature, suggesting that these interactions could involve domains of the delta subunit required for binding of the coupling factor to the mitochondrial membrane.     

Determining the origin of the stabilization of DNA by 5-aminopropynylation of pyrimidines

DNA duplexes are stabilized by aminopropynyl modification of pyrimidines at the 5 position. A combination of thermodynamic analyses as a function of ionic strength, NMR, and molecular modeling has been applied to determine the origin of the stabilization. UV melting studies of a dodecamer bearing one, two, or three nonadjacent modified dU and dC and of a single dU(8) in the Dickerson-Drew dodecamer revealed that the modifications are essentially additive in terms of T(m), DeltaG, and DeltaH, and there is little difference between dU and dC. The free energy change was parsed into electrostatic and nonelectrostatic components, which showed a significant contribution from charge interactions at physiological ionic strength but also a nonelectrostatic contribution that arises in part from hydration. NMR spectroscopy of the modified Dickerson-Drew dodecamer revealed that the conformation of the duplexes is not significantly altered by the modifications, though (31)P NMR shows that the positive charge may affect ionic interactions with the oxygen atoms of the neighboring phosphates. The modified duplex showed significant hydration in both major and minor grooves. The single strands were also analyzed by NMR, which showed evidence of significant stacking interactions in the modified oligonucleotide. Parsing the energy contribution has shown that electrostatics and hydration can produce substantial increases in thermodynamic stability without significant changes in the conformation of the duplex state. These considerations have significance for the design of oligonucleotides used for hybridization.     

Determination of partial unfolding enthalpy for lysozyme upon interaction with dodecyltrimethylammonium bromide using an extended solvation model

The interactions of dodecyltrimethylammonium bromides (DTABs) with hen egg lysozyme have been investigated at pH = 7.0 and 27 degrees C in phosphate buffer by isothermal titration calorimetry. DTAB interacts endothermically and activate lysozyme. The endothermicity of the lysozyme-DTAB interaction is in marked contrast to the exothermic interactions between sodium dodecyl sulphate (SDS) and lysozyme which have been attributed to specific binding between the anionic sulphate head groups and cationic amino acid residues. The enthalpies of interaction between the cationic surfactant (DTAB) and lysozyme are dominated by the endothermic unfolding of the native structure followed by an exothermic solvation of the lysozyme-DTAB complex by the addition of extra DTAB. A new direct calorimetric method to follow protein denaturation, and the effect of surfactants on the stability of proteins was introduced. The extended solvation model was used to reproduce the enthalpies of lysozyme-DTAB interaction over the whole range of DTAB concentrations. The solvation parameters recovered from the new equation, attributed to the structural change of lysozyme and its biological activity. At low concentrations of DTAB, the binding is mainly electrostatic, with some simultaneous interaction of the hydrophobic tail with nearby hydrophobic patches on the lysozyme. These initial interactions presumably cause some protein unfolding and expose additional hydrophobic sites. The DTAB-induced denaturation enthalpy of lysozyme is 86.46 +/- 0.02 kJ mol(-1).     

Comments on selected aspects of nucleic acid electrostatics

Recent experimental, theoretical, and computational developments in the field of nucleic acid electrostatics have brought interesting concepts to the fore. The phosphate charge on the double helix apparently influences its structure. When the charge is neutralized asymmetrically, the resulting force imbalance drives bending toward the neutralized side. When the charge is uniformly neutralized, the force imbalance acts to buckle the helical axis, resulting in a compact tertiary conformation. Sharing of condensed counterions by single strands is a stabilizing factor for formation of the double helix. Sharing of condensed counterions by two double helices causes clustering of DNA and may be a factor in RNA folding. Support for these statements is reviewed.     

Different modes of interaction in cyanobacterial complexes of plastocyanin and cytochrome f

The highly efficient electron-transfer chain in photosynthesis demonstrates a remarkable variation among organisms in the type of interactions between the soluble electron-transfer protein plastocyanin and it partner cytochrome f. The complex from the cyanobacterium Nostoc sp. PCC 7119 was studied using nuclear magnetic resonance spectroscopy and compared to that of the cyanobacterium Phormidium laminosum. In both systems, the main site of interaction on plastocyanin is the hydrophobic patch. However, the interaction in the Nostoc complex is highly dependent on electrostatics, contrary to that of Phormidium, resulting in a binding constant that is an order of magnitude larger at low ionic strength for the Nostoc complex. Studies of the mixed complexes show that these differences in interactions are mainly attributable to the surface properties of the plastocyanins.     

Factors determining the reconstructive denaturation of proteins in sodium dodecyl sulfate solutions. Further circular dichroism studies on structural reorganization of seven proteins

The factors determining the onset and extent of reconstructive denaturation of proteins were considered by comparing circular dichroism (CD) data of seven proteins and previously published findings. The effects of sodium dodecyl sulfate (SDS) on the conformation of the following proteins were tested: lysozyme, the mitogens fromPhytolacca americana (fractions Pa2 and Pa4), lectin fromWistaria floribunda, ovine lutropin, a Bence Jones protein, and histone H2B. While the helix content of lysozyme was raised by SDS slightly, in the Bence Jones protein andW. floribunda lectin it increased from near zero to about 25–30%. In histone H2B the helix content was raised by SDS even to about 48%. However, no clear indication of helix formation could be observed in the mitogens and lutropin, even at low pH or 2.0–2.5. The tertiary structure of the proteins was perturbed by SDS. It was concluded that the reorganization of secondary structure of the proteins was favored by the following factors: (1) presence of helicogenic amino acid sequences in the protein, (2) availability of positively charged sites of the basic amino acids for interactions with the dodecyl ion, (3) absence of a large surplus of negatively charged sites on the surface of protein, and (4) absence of extensive disulfide cross-linking within the macromolecule. Both hydrophobic and electrostatic interactions occur in reconstructive denaturation, and the newly formed helices are stabilized by hydrophobic shielding by the alkyl chains of the alkyl sulfate.     

A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

The Correlation of Cold Denaturation Temperature with Surface Stability Factor of Proteins

Cold denaturation is an intriguing phenomenon in protein denaturation for elucidating protein accessible surface area (ASA). Compared to the impact of protein surface, the importance of protein-water interactions in cold denaturation may be ruled out significantly. Here, based on the ASA, we have defined a new factor, the surface stability factor (SSF). From the SSF, in combination with the cold denaturation temperature (T(g')) or temperature at DeltaS = 0 (T(s)) of a given protein, one can predict the percent of hydrophobic surface area (H), percent of total surface there on positive and negative charge sum (effective charge) be zero (C), percent of patches hydrophobicity (HP) and others critical surface parameters without any need to the crystallographic data.     

The critical role of partially exposed N-terminal valine residue in stabilizing GH10 xylanase from Bacillus sp.NG-27 under poly-extreme conditions

Background

Understanding the mechanisms that govern protein stability under poly-extreme conditions continues to be a major challenge. Xylanase (BSX) from Bacillus sp. NG-27, which has a TIM-barrel structure, shows optimum activity at high temperature and alkaline pH, and is resistant to denaturation by SDS and degradation by proteinase K. A comparative circular dichroism analysis was performed on native BSX and a recombinant BSX (R-BSX) with just one additional methionine resulting from the start codon. The results of this analysis revealed the role of the partially exposed N-terminus in the unfolding of BSX in response to an increase in temperature.

Methodology

We investigated the poly-extremophilicity of BSX to deduce the structural features responsible for its stability under one set of conditions, in order to gain information about its stability in other extreme conditions. To systematically address the role of the partially exposed N-terminus in BSX stability, a series of mutants was generated in which the first hydrophobic residue, valine (Val1), was either deleted or substituted with various amino acids. Each mutant was subsequently analyzed for its thermal, SDS and proteinase K stability in comparison to native BSX.

Conclusions

A single conversion of Val1 to glycine (Gly) changed R-BSX from being thermo- and alkali- stable and proteinase K and SDS resistant, to being thermolabile and proteinase K-, alkali- and SDS- sensitive. This result provided insight into the structure-function relationships of BSX under poly-extreme conditions. Molecular, biochemical and structural data revealed that the poly-extremophilicity of BSX is governed by a partially exposed N-terminus through hydrophobic interactions. Such hitherto unidentified N-terminal hydrophobic interactions may play a similar role in other proteins, especially those with TIM-barrel structures. The results of the present study are therefore of major significance for protein folding and protein engineering.     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.