Molecular Dynamics Simulations of the Unfolding of the Starch Binding Domain from Aspergillus niger Glucoamylase

Abstract

The 600 ps molecular dynamics simulations to investigate the unfolding of the starch binding domain from Aspergillus niger glucoamylase were conducted in vacuum as well as in an external field with the dielectric constant of 80 with temperature jump technique. Electrostatic interactions play an important role in determining the stability of the β-strands in this domain. The starch binding site 1 is less stable than site 2 since it is more exposed to the surface. The disulfide bond between C509 and C604 is unstable since these two residues are located near the flexible linker domain and in the mobile loop region between β-strands 6 and 7, respectively. The melting temperature, at which the total residual β-strand content is 50% that of the solution structure, is about 544K for the simulations with dielectric constant of 80, leading to the estimated unfolding timescale of 0.48 ms in vitro. In addition, the unfolding of the starch binding domain is proposed to initiate from the interior region by the lost of the integrity of the secondary structure.