Molecular dynamics simulations identify the topological weak spots of a protease CN2S8A

Thermophilic enzymes are highly desired in industrial applications due to their efficient catalytic activity at high temperature. However, most enzymes exhibit inferior thermostability and it remains challenging to identify the optimal sites for designing mutations to improve protein stability. To tackle this issue, we integrated topological analysis and all-atom molecular dynamics simulations to efficiently pinpoint the thermally-unstable regions in protein structures. Using a protease CN2S8A as the model, we analyzed the intramolecular hydrogen bonding interactions between adjacent secondary structure elements, and then identified the topological weak spots of CN2S8A where weak hydrogen bonding interactions were formed. To examine the role of these sites in protein structural stability, we designed three virtual mutations at different weak spots and characterized the effects of these mutations on the structural properties of CN2S8A. The results showed that all three mutations increased the protein structural stability. In conclusion, these findings provide a novel method to identify the topological weak spots of proteins, with implications in the rational design of biocatalysts with superior thermostability.