Kinetic stability is a key parameter to comprehend protein behavior and it plays a central role to understand how evolution has reached the balance between function and stability in cell-relevant timescales. Using an approach that includes simulations, protein engineering, and calorimetry, we show that there is a clear correlation between kinetic stability determined by differential scanning calorimetry and protein thermal flexibility obtained from a novel method based on temperature-induced unfolding molecular dynamics simulations. Thermal flexibility quantitatively measures the increment of the conformational space available to the protein when energy in provided. The (β/α) barrel fold of two closely related by evolution triosephosphate isomerases from two trypanosomes are used as model systems. The kinetic stability-thermal flexibility correlation has predictive power for the studied proteins, suggesting that the strategy and methodology discussed here might be applied to other proteins in biotechnological developments, evolutionary studies, and the design of protein based therapeutics.

No abstract

No abstract

The effect of specific residues on the kinetic stability of two closely related triosephosphate isomerases (from Trypanosoma cruzi, TcTIM and Trypanosoma brucei, TbTIM) has been studied. Based on a comparison of their β-turn occurrence, we engineered two chimerical enzymes where their super secondary β-loop-α motifs 2 ((βα) ) were swapped. Differential scanning calorimetry (DSC) experiments showed that the (βα) motif of TcTIM inserted into TbTIM (2Tc) increases the kinetic stability. On the other hand, the presence of the (βα) motif of TbTIM inserted into TcTIM (2Tb) gave a chimerical protein difficult to purify in soluble form and with a significantly reduced kinetic stability. The comparison of the contact maps of the (βα) of TbTIM and TcTIM showed differences in the contact pattern of residues 43 and 49. In TcTIM these residues are prolines, located at the N-terminal of loop-2 and the C-terminal of α-helix-2. Twelve mutants were engineered involving residues 43 and 49 to study the effect over the unfolding activation energy barrier (E ). A systematic analysis of DSC data showed a large decrease on the E of TcTIM (ΔE ranging from 468 to 678 kJ/mol) when the single and double proline mutations are present. The relevance of Pro43 to the kinetic stability is also revealed by mutation S43P, which increased the free energy of the transition state of TbTIM by 17.7 kJ/mol. Overall, the results indicate that protein kinetic stability can be severely affected by punctual mutations, disturbing the complex network of interactions that, in concerted action, determine protein stability. Proteins 2017; 85:571-579. © 2016 Wiley Periodicals, Inc.