Nitrile hydratase (NHase), an excellent bio-catalyst for the synthesis of amide compounds, was composed of two heterologous subunits. A thermoalkaliphilic NHase NHCTA1 (Tm = 71.3°C) obtained by in silico screening in our study exhibited high flexibility of α-subunit but excellent thermostability, as opposed to previous examples. To gain a deeper structural insight into the thermostability of NHCTA1, comparative molecular dynamics simulation of NHCTA1 and reported NHases was carried out. By comparison, we speculated that β-subunit played a key role in adjusting the flexibility of α-subunit and the different conformations of linker in "α5-helix-coil ring" supersecondary structure of β-subunit can affect the interaction between β-subunit and α-subunit. Mutant NHCTA1-α C with a random coil linker and mutant NHCTA1-αβγ with a truncated linker were therefore constructed to understand the impact on NHCTA1 thermostability by varying the supersecondary structure. The varied thermostability of NHCTA1-α C and NHCTA1-αβγ (Tm = 74.4°C, Tm = 65.6°C) verified that the flexibility of α-subunit adjusted by β-subunit was relevant to the stability of NHCTA1. This study gained an insight into the NNHCTA1 thermostability by virtual dynamics comparison and experimental studies without crystallization, and this approach could be applied to other industrial-important enzymes.