QM/MM study on enzymatic mechanism in sinigrin biosynthesis

As the major and abundant type of glucosinolates (GL) in plants, sinigrin has potential functions in promoting health and insect defense. The final step in the biosynthesis of sinigrin core structure is highly representative in GL compounds, which corresponds to the process from 3-methylthiopropyl ds-GL to 3-methylthiopropyl GL catalyzed by sulfotransferase (SOT). However, due to the lack of the crystallographic structure of SOT complexed with the 3-methylthiopropyl GL, little is known about this sulfonation process. Fortunately, the crystal structure of SOT 18 from Arabidopsis thaliana (AtSOT18) containing the substance (sinigrin) similar to 3-methylthiopropyl GL has been determined. To understand the enzymatic mechanism, we employed molecular dynamics (MD) simulation and quantum mechanics combined with molecular mechanics (QM/MM) methods to study the conversion from ds-sinigrin to sinigrin catalyzed by AtSOT18. The calculated results demonstrate that the reaction occurs through a concerted dissociative mechanism. Moreover, Lys93, Thr96, Thr97, Tyr130, His155, and two enzyme peptide chains (Pro92-Lys93 and Gln95-Thr96-Thr97) play a role in positioning the substrates and promoting the catalytic reaction by stabilizing the transition state geometry. Particularly, His155 acts as a catalytic base while Lys93 acts as a catalytic acid in the reaction process. The presently proposed concerted dissociative mechanism explains the role of AtSOT18 in sinigrin biosynthesis, and could be instructive for the study of GL biosynthesis catalyzed by other SOTs.