Yuzhen Li, Yunqing Wang, Zhenmin Liu, Chunping You
{"title":"关于葡萄糖和半乳糖与赖氨酸的马氏反应的计算研究。","authors":"Yuzhen Li, Yunqing Wang, Zhenmin Liu, Chunping You","doi":"10.1007/s00894-024-06105-w","DOIUrl":null,"url":null,"abstract":"<div><h3>Context</h3><p>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 <i>N</i><sup>α</sup>- and <i>N</i><sup>ε</sup>-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 <i>N</i><sup>ε</sup>-channel is lower than in <i>N</i><sup>α</sup>-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.</p><h3>Methods</h3><p>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 H<sub>2</sub>O and calculated in wB97XD functional with 6–31 + G (d) basis set. All computations were carried out Gaussian 09 suite of quantum chemistry software.</p><h3>Graphical Abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1000,"publicationDate":"2024-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Computational study on the Maillard reactions of glucose and galactose with lysine\",\"authors\":\"Yuzhen Li, Yunqing Wang, Zhenmin Liu, Chunping You\",\"doi\":\"10.1007/s00894-024-06105-w\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><h3>Context</h3><p>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 <i>N</i><sup>α</sup>- and <i>N</i><sup>ε</sup>-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 <i>N</i><sup>ε</sup>-channel is lower than in <i>N</i><sup>α</sup>-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.</p><h3>Methods</h3><p>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 H<sub>2</sub>O and calculated in wB97XD functional with 6–31 + G (d) basis set. All computations were carried out Gaussian 09 suite of quantum chemistry software.</p><h3>Graphical Abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>\",\"PeriodicalId\":651,\"journal\":{\"name\":\"Journal of Molecular Modeling\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":2.1000,\"publicationDate\":\"2024-08-20\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Molecular Modeling\",\"FirstCategoryId\":\"92\",\"ListUrlMain\":\"https://link.springer.com/article/10.1007/s00894-024-06105-w\",\"RegionNum\":4,\"RegionCategory\":\"化学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q4\",\"JCRName\":\"BIOCHEMISTRY & MOLECULAR BIOLOGY\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Molecular Modeling","FirstCategoryId":"92","ListUrlMain":"https://link.springer.com/article/10.1007/s00894-024-06105-w","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"BIOCHEMISTRY & MOLECULAR BIOLOGY","Score":null,"Total":0}
Computational study on the Maillard reactions of glucose and galactose with lysine
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.
期刊介绍:
The Journal of Molecular Modeling focuses on "hardcore" modeling, publishing high-quality research and reports. Founded in 1995 as a purely electronic journal, it has adapted its format to include a full-color print edition, and adjusted its aims and scope fit the fast-changing field of molecular modeling, with a particular focus on three-dimensional modeling.
Today, the journal covers all aspects of molecular modeling including life science modeling; materials modeling; new methods; and computational chemistry.
Topics include computer-aided molecular design; rational drug design, de novo ligand design, receptor modeling and docking; cheminformatics, data analysis, visualization and mining; computational medicinal chemistry; homology modeling; simulation of peptides, DNA and other biopolymers; quantitative structure-activity relationships (QSAR) and ADME-modeling; modeling of biological reaction mechanisms; and combined experimental and computational studies in which calculations play a major role.