Predicted unfolding order of the 13 alpha-helices in the catalytic domain of glucoamylase from Aspergillus awamori var. X100 by molecular dynamics simulations.

The unfolding mechanism of the 13 alpha-helices in the catalytic domain of Aspergillus awamori var. X100 glucoamylase was investigated by 200 ps molecular dynamics simulations in explicit water with temperature jump technique. Rather than a simultaneous event, the unfolding of these 13 alpha-helices followed a random ordered mechanism as alpha8-->alpha1-->alpha11-->alpha7-->alpha10-->alpha3-->alpha12-->alpha13-->alpha4-->alpha5-->alpha9-->alpha6-->alpha2. No significant relationships were found between the unfolding order and the length and the hydrophobicity of the helix. alpha-Helix 8 located in the inner region of the catalytic domain was predicted to be the first helix to unfold, indicating that the destruction of the secondary structure motif was initiated from the inner region of the catalytic domain. The dynamic behavior of these alpha-helices induced by increased kinetic energy during the unfolding process is considered to be similar to the expansion and compression of a series of springs under the influence of mechanical stress.