A native protein in aqueous salt solution usually unfolds reversibly upon heating at a transition temperature (T
m) with unfolding enthalpy (ΔH
m). The method to predict T
m and ΔH
m from the 3D structure assuming simple two-states transition has been developed using experimental values of T
m, ΔH
m and heat capacity change (ΔC
p, u) compiled in the database“ProTherm”. In order to account for dependence of thermostability on pH and salt concentration, the difference in electrostatic interactions between native and denatured states (Δh
u, ele) was computed applying Debye-Huckel type potential between partial charges located on ionizable side chain atoms. Standard unfolding free energy (ΔG
0u), enthalpy (ΔH
0u), and ΔC
p, u at 25°C after subtracting electrostatic energies were obtained using a data set of 63 experimental data on 32 uncharged proteins. Each of these quantities was assumed to be a sum of products of the difference in accessible surface area (ΔASA) of an i-th constituent atomic group between the native and unfolded extended structure times a corresponding adequate constant. Proportional constants for six (aliphatic, aromatic, hydroxyl, amide, carbonyl, and thiol) atomic groups were determined using the data set. Thus, ΔG
0u, ΔH
0u and ΔC
p, u for a protein at uncharged form of known 3D structure can be computed using its ΔASA value of every group, and ΔG
u(T) at T is determined according to thermodynamic relation, adding Δh
u, ele(T). Values of T
m and ΔH
m for the protein may, therefore, be evaluated as values at ΔG
u(T)=0. The agreement of calculated and experimental values for the data set used for the determination of the constants was practically the same as that for another data set of 33 experimental data on T
m and ΔH
m for 19 proteins with standard deviations of about 12°C and 38kcal/mol, respectively. The stability and instability of correctly and misfolded proteins (hemerythrin, and Ig-variable domain fragment) were successfully explained in terms of ΔG
u.
View full abstract