The mechanism of Maillard reaction has been investigated by means of density functional theory calculations in the gaseous phase and aqueous solution. The Maillard reaction is a cascade of consecutive and parallel reaction. In the present model system study, glucose and glycine were taken as the initial reactants. On the basis of previous experimental results, the mechanisms of Maillard reaction have been proposed, and the possibility for the formation of different compounds have been evaluated through calculating the relative energy changes for different steps of reaction under different pH conditions. Our calculations reveal that the TS3 in Amadori rearrangement reaction is the rate-determining step of Maillard reaction with the activation barriers of about 66.7 and 68.8 kcal mol(-1) in the gaseous phase and aqueous solution, respectively. The calculation results are in good agreement with previous studies and could provide insights into the reaction mechanism of Maillard reaction, since experimental evaluation of the role of intermediates in the Maillard reaction is quite complicated.