Computational study on the Maillard reactions of glucose and galactose with lysine
Yuzhen Li1, Yunqing Wang1,2, Zhenmin Liu3
1State Key Laboratory of Dairy Biotechnology, Shanghai Engineering Research Center of Dairy Biotechnology, Postdoctoral Workstation of Bright Dairy - Shanghai Jiao Tong University, Dairy Research Institute, Bright Dairy & Food Co., Ltd., Shanghai, 200436, China.
Context:
Milk has nutrient-rich but thermal sensitive matrix that undergoes varying degrees of Maillard reaction (MR) at heating conditions. The MR mainly occurs between lysine residues (Lys) and lactose composed of glucose (Glc) and galactose (Gal), which are abundantly sourced from dairy products. In the present study, the MRs of Glc and Gal with Lys at the initial and intermediate stages have been investigated theoretically using density functional theory (DFT) to simulate the gaseous and aqueous phases. Reaction mechanisms have been proposed, and relative energy changes of different steps were calculated according to the total mass balance. The calculations reveal that both Nα- and Nε-amine groups of Lys can react with the carbonyl functional group of Glc and Gal with the similar potential energy profiles, and Gal is more reactive than Glc. However, the barrier in Nε-channel is lower than in Nα-channel, indicating a faster reaction rate through the former channel compared with the latter. The 5-hydroxymethyl-2-furfural (HMF) and derivative are formed under 3-deoxysone route in the intermediate stage. The calculation results are helpful for proposing a reasonable MR mechanism and suggesting possible control methods of the MRs.
Methods:
In this study, different levels of DFT calculations have been conducted to investigate the mechanisms and favorability of generating MR products in Glc-Lys and Gal-Lys models at initial and intermediate stages in the gaseous and aqueous conditions. In order to elucidate the molecular models from the perspectives of chemistry and geometry, DFT calculations were performed by the mean of B3LYP functional at basis sets of 6-311 + + G (d, p) and 6-311 + + G (2df, 2p) with optional solvation settings. To examine the solvation effect, the study further constructed models with solvent H2O and calculated in wB97XD functional with 6-31 + G (d) basis set. All computations were carried out Gaussian 09 suite of quantum chemistry software.
