Molecular dynamics simulations reveal structural instability of human trypsin inhibitor upon D50E and Y54H mutations

Serine protease inhibitor Kazal type 1 (SPINK1) plays an important role in protecting the pancreas against premature trypsinogen activation that causes pancreatitis. Various mutations in the SPINK1 gene were shown to be associated with patients with pancreatitis. Recent transfection studies identified intracellular folding defects, probably caused by mutation induced misfolding of D50E and Y54H mutations, as a common mechanism that reduces SPINK1 secretion and as a possible novel mechanism of SPINK1 deficiency associated with chronic pancreatitis. Using molecular dynamics, we investigated the effects of D50E and Y54H mutations on SPINK1 dynamics and conformation at 300 K. We found that the structures of D50E and Y54H mutants were less stable than and were distorted from those of the wild type, as indicated by the RMSD plots, RMSF plots and DSSP series. Specifically, unwinding of the top of helices (the main secondary structures) and the distortion of the loops above the helices were observed. It may be possible that this distorted protein structure may be recognized as “non-native” by members of the chaperone family; it may be further retained and targeted for degradation, leading to SPINK1 secretion reduction and subsequently pancreatitis in patients as Király et al. (Gut 56:1433, 2007) proposed.